xiangguan21.m

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

clc;clear;
Nm=150;
t1=5e-6;t2=5e-9;
tm=unifrnd(t1-t2,t1+t2,1,Nm);%生成均匀分布随机数
tr=unifrnd(t1-t2,t1+t2,1,Nm);%生成均匀分布随机数
for v=1:1:400,i=v;
    z=2*pi*v*1e+006;
    for j=1:Nm
        for k=1:Nm
            a(j,k)=cos(z*(tm(j)+tr(k)));
            b(j,k)=cos(z*(tm(j)-tr(k)));
            c(j,k)=a(j,k)+b(j,k);
        end;
     h(j)=sum(xcorr(c(j,:)))/2;   
    end;
    R(i)=sum(h(j))+2e+004;
    H(i)=R(i)/4e+004;
end;
v=1:1:400;
figure,plot(v,H,'k-','LineWidth',3)
axis([0,400,0.4,1])
grid on
xlabel('v(MHz)'),ylabel('<|H(f)H(f+v)|^2>')