xiangguan11.m

来自「多模光纤的信道特性研究」· M 代码 · 共 21 行

clc;clear;
Nm=100;
r=zeros(400,Nm);
t1=5e-6;t2=5e-9;
rand('seed',1);%把均匀分布伪随机发生器置为1状态
tn=unifrnd(t1-t2,t1+t2,1,Nm);%生成均匀分布随机数
for v=1:1:400,
    i=v;  
    for j=1:Nm,
    h(i,j)=cos(2*pi*v*1e+6*tn(j));
    r(i,j)=xcorr(h(i,j));
    end;
    H=sum(r,2)-50;
    R=H/50;
end;
v=1:1:400;
figure,plot(v,R,'k-','LineWidth',3) 
grid on
xlabel('v(MHz)'),ylabel('R(v)')