matlab-(含教程)基于Jakes信道的MATLAB性能仿真

function [chann,fm,doppler_rate]=modified_jakes(v,fc,fs,U,M) % v: 物体运动速度（单位：km/h） % fc: 载波频率（单位：Mhz） % fs: 采样频率（单位：Khz） % U: 路径的个数 % M: Number of channel coefficients to generate % chann: channels coefficients matrix of size U by M % fm: 最大多普勒频率（单位：Hz） % doppler_rate: 多普勒率 fm=(1000/3600)*v*fc/300; %最大多普勒频率（单位：Hz） doppler_rate=fm/(fs*1000); No=24; % 不同振荡器的数目 NN=4*No; %振荡器的数目 wm=2*pi*fm; %最大多普勒频率 H=hadamard(No); % 用Walsh_Hadamard码来产生不相关的信道 for n=1:No theta(n)=rand*2*pi; %附加相位，相互独立，且在(0,2 ]均匀分布 end for u=1:U sum=0; t=1/(fs*1000):1/(fs*1000):(1/(fs*1000))*M; % 采样时间 for n=1:No wn=wm*cos(2*pi*(n-0.5)/NN); sum=sum+H(n)*(cos(pi*n/No)+j*sin(pi*n/No))*cos(wn*t+theta(n)); % sum of No distinct oscillators end R=sqrt(2/No)*sum; chann(u,:)=R; %第u条路径 end