《MIMO-OFDM系统原理、应用及仿真》.docx

2.5 仿真实例 实例2-1 瑞利分布与莱斯分布 功能：绘制瑞利分布曲线与莱斯分布曲线 程序名称：Example2_1.m 程序代码： clear, clf N=200000; %产生200000个信道系数供统计使用 level=30; %统计区间被划分的分数。 K_dB=[-40 0 15];%莱斯因子为-40dB、0dB、15dB gss=[k-*; k-o; k-+;k-^];%绘制曲线的颜色、线形与标志符号 % 瑞利模型 Rayleigh_ch=Ray_model(N);%调用Ray_model子程序，产生瑞利分布幅度系数 [temp,x]=hist(abs(Rayleigh_ch(1,:)),level);%统计数据分布 plot(x,temp,gss(1,:)) hold on %莱斯模型 for i=1:length(K_dB);%对不同莱斯因子进行信道模型仿真 Rician_ch(i,:) = Ric_model(K_dB(i),N);%调用Ric_model产生莱斯分布幅度系数 [temp x] = hist(abs(Rician_ch(i,:)),level);%统计数据分布 plot(x,temp,gss(i+1,:)) end xlabel(x), ylabel(Occurrence) legend(Rayleigh,Rician, K=-40dB,Rician, K=0dB,Rician, K=15dB) %瑞利信道模型子程序，子程序程序名称：Ray_model.m function H=Ray_model(L) % 输入参数 L: 仿真信道个数，为N=200000 % 输出参数 H： 返回瑞利信道矩阵 H = (randn(1,L)+j*randn(1,L))/sqrt(2); %产生实部为高斯分布、虚部为高斯分布、包络为瑞利分布的信道系数。实部功率为1/2，虚部功率为1/2，因 %此该行指令返回单位功率的或称归一化功率的瑞利信道幅度系数。 %莱斯信道模型子程序，子程序程序名称：Ric_model.m function H=Ric_model(K_dB,L) % 输入参数 : K_dB 为莱斯因子，L为仿真信道个数 % 输出参数H: 返回莱斯信道矩阵 K = 10^(K_dB/10);%将dB值描述的莱斯因子转换为幅度值 H = sqrt(K/(K+1)) + sqrt(1/(K+1))*Ray_model(L); %产生莱斯信道幅度系数。莱斯信道模型中包含视距通信，收发之间有直通路径。 程序仿真结果见图2-9。 实例2-2 两径信道与指数信道模型 功能：产生一个两径信道和一个指数衰减的多径信道。 程序名称：Example2_2 程序代码： clear, clf scale=1e-9; % 纳秒量级 Ts=10*scale; % 抽样时间间隔为10ns，在这个程序中这个量也为指数信道路径间隔 t_rms=30*scale; % RMS 时延扩展为30ns num_ch=10000; % 仿真信道个数 % 两径信道模型 % 产生并绘制了理想的两径信道模型和瑞利分布两径信道模型。 pow_2=[0.5 0.5]; delay_2=[0 t_rms*2]/scale; %给出理想两径信道功率均为0.5，延时为0和60ns H_2 = [Ray_model(num_ch); Ray_model(num_ch)].*diag(sqrt(pow_2)); %产生瑞利两径信道幅度系数。通过调用子程序Ray_model产生归一化功率的瑞利两径信道幅度系数。 avg_pow_h_2 = mean(H_2.*conj(H_2)); %计算瑞利分布两径信道每一径的平均功率。在这里可以看到上一条语句中diag(sqrt(pow_2))的作用。%当通过对幅度系数进行运算计算功率时，sqrt(pow_2)可以使每一径的功率为pow_2，即每一径的功率为0.5。 subplot(121) stem(delay_2,pow_2,ko), hold on, stem(delay_2,avg_pow_h_2,k.); xlabel(Delay[ns]), ylabel(Channel Power[linear]); title(2-ray Model); legend(Ideal,Simulation); axis([-10 140 0 0.7]); % 指数信道模型 %产生并绘制理想的指数信道模型和瑞利分布的指数信道模型。 pow_e=exp_PDP(t_rms,Ts); %通过调用exp_PDP子程序