Báo cáo thí nghiệm thông tin vô tuyến 2011

Báo cáo thí nghiệm thông tin vô tuyến Bách Khoa khóa K53
Nội dung trích xuất từ tài liệu:
Báo cáo thí nghiệm thông tin vô tuyến 2011 Phụ lục Mã nguồn Matlab và các kết quả mô phỏng Bài số 1: Mô phỏng kênh Rayleigh theo phương pháp RiceBài 1.1: Tính các thông số f, C, θ:%================================%Tinh cac thong so f, C, th%================================clearf_m=91; %Tan so Doppler lon nhatb=1; %Phuong sai cua qua trinh g(t)N1=9; %So phan tu g(t) phan thuc mo phongN2=10; % So phan tu g(t) phan ao mo phongfor n=1:1:N1; c1(n)=sqrt(2*b/N1); f1(n)=f_m*sin(pi*(n-0.5)/(2*N1)); th1(n)=2*pi*n/(N1+1);endfor n=1:1:N2; c2(n)=sqrt(2*b/N2); f2(n)=f_m*sin(pi*(n-0.5)/(2*N2)); th2(n)=2*pi*n/(N2+1);endsave ex4p1_Res f1 f2 c1 c2 th1 th2Bài 1.2: Lập hàm g.m để tính toán%=======================================%Ham tao qua trinh g(t)%=======================================function y=g(c,f,th,t) y=zeros(size(t)); for n=1:length(f); y=y+c(n)*cos(2*pi*f(n).*t+th(n)); end;Bài 1.3: Viết hàm để đưa ra biên độ kênh và vẽ hàm :%=====================================%Tinh bien do kenh va ham alpha%=====================================clear;load ex4p1_Res f1 f2 c1 c2 th1 th2f_s=270800; %Tan so lay mauT_sim=0.4; %Thoi gian mo phongt=0:1/f_s:T_sim; % Thoi gian mo phong tg1=g(c1,f1,th1,t);g2=g(c2,f2,th2,t);g=g1+j*g2;alpha=abs(g);alpha_dB=20*log10(alpha);Hoàng Sỹ Tân – MSSV: 20083473 – Lớp ĐT3-K53plot(t,alpha_dB);title(Bien do kenh (dB));xlabel(alpha(t));legend(alpha(t)in dB,0);Đồ thị:Bài 1.4: xác định các thông số giá trị trung bình, độ biến đổi, hàm mật độ xác suất (PDF):%======================================%Ham tinh gia tri trung binh va% do bien doi cua g1(t), g(t), alpha(t)%Ham mat do xac suat cua cac qua trinh% bang bieu do%======================================clear;load ex4p1_Res f1 f2 c1 c2 th1 th2f_s=50000; %Tan so lay mauT_sim=20; %Thoi gian mo phongt=0:1/f_s:T_sim;g1=g(c1,f1,th1,t);g2=g(c2,f2,th2,t);g=g1+j*g2;Hoàng Sỹ Tân – MSSV: 20083473 – Lớp ĐT3-K53alpha=abs(g);mean_g=mean(g); %Gia tri trung binh cua gmean_g1=mean(g1); % Gia tri trung binh cua g1mean_alpha=mean(alpha); % Gia tri trung binh cua alphavar_g=var(g); %Do bien doi cua gvar_g1=var(g1); % Do bien doi cua g1var_alpha=var(alpha); % Do bien doi cua alphan=length(alpha);x=0:0.1:3; %Thoi gian tinh bang giayb=hist(alpha,x); %Bieu do mat do xac suat cua alphafigure(1);stem(x,b/n/(x(2)-x(1)));%Ve tin hieu roi rachold on;k=0; %He so k=s^2/(2*b_0)total_p=2; %Cong suat nhan duocalpha_p=(2.*x.*(k+1)/total_p).*exp(-k-((k+1).*x.^2/total_p)).*besseli(0,(2.*x.*sqrt(k*(k+1)/total_p)));plot(x,alpha_p,r);title(Mat do xac suat cua alpha(x));xlabel(x);ylabel(P_{alpha}(x));legend(p_{alpha}(x),Phan bo Rayleigh);hold off;figure(2);n1=length(g1);x1=-4:0.1:4;c=hist(g1,x1); %Bieu do mat do xac suat cua g1stem(x1,c/n1/(x1(2)-x1(1)));hold on;p=(1/sqrt(2*pi))*exp(-x1.^2/2);plot(x1,p,r);title(Mat do xac suat PDF cua g1);xlabel(x);ylabel(P_{g1}(x));legend(p_{g1}(x),Phan bo Gaussian);hold off;Đồ thị:Hoàng Sỹ Tân – MSSV: 20083473 – Lớp ĐT3-K53Hoàng Sỹ Tân – MSSV: 20083473 – Lớp ĐT3-K53Hoàng Sỹ Tân – MSSV: 20083473 – Lớp ĐT3-K53 Bài số 2: Mô hình kênh Fadding phân tập đa đường theo phương pháp Monte CarloBài 2.1: Mô hình kênh phân tập đa đường :Hàm mô hình Monte Carlo%==========================================================================%Ham mo hinh kenh phan tap da duong theo phuong phap Monte Carlo%==========================================================================function [h,t_next]=MCM_channel_model(u,initial_time,number_of_summations,symbol_duration,f_dmax,channel_coefficients);%u la bien ngau nhien duoc dung de roi rac tan so, pha cua hieu uong Doppler%number_of_summation: so luong hoac pha roi rac duoc su dung trong%phuong phap Monte Carlo%symbol+duration la khoang thoi gian lay mau cua tin hieu%f_fmax La tan so Doppler lon nhatt=initial_time; %thoi gian bat dau quan sat kenhChannel_Length=length(channel_coefficients); % Chieu dai dap ung kenhh_vector=[]; %Vector mo hinh CIRfor k=1:Channel_Length; u_k=u(k,:); %Mot bien ngau nhien cua k phi=2*pi*u_k; %tan so Doppler roi rac su dung trong Monte Carlo f_d=f_dmax*sin(2*pi*u_k); h_tem=channel_coefficients(k)*1/(sqrt(number_of_summations))*sum(exp(j*phi).*exp(j*pi*f_d*t)); h_vector=[h_vector,h_tem];end;h=h_vector;t_next=initial_time+symbol_duration;endChương trình chính:%==========================================================================%Chuong trinh chinh mo phong phuong phap Monte Carlo%==========================================================================clear all;%Dat cac thong soNFFT=64; %Chieu dai FFTt_a=8.0000e-007; %Th ...