xindao.m

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

clc;clear;         
Nm=100;
t1=5e-6;t2=2.5e-9;
x=zeros(1,1000);
y=zeros(1,1000);
h=zeros(1,1000);
H=zeros(1,1000);
rand('seed',1);%把均匀分布伪随机发生器设为1。 
tn=unifrnd(t1-t2,t1+t2,1,Nm);%取Nm个服从均匀分布随机的数。
for f=3:3:3000
    i=f/3;
     x(i)=sum(cos(2*pi*1e+006*f*tn));
     y(i)=sum(sin(2*pi*1e+006*f*tn));
    h(i)=sqrt(x(i).^2+y(i).^2);
    H(i)=10*log10(h(i));
end; 
f=3:3:3000;
figure,plot(f,H)
title('Nm=100的频率响应曲线')
grid on
xlabel('frequency(MHz)')
ylabel('H(f) in dB')