45827《MIMO-OFDM系统原理、应用及仿真》李莉(实例代码).pdfVIP

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