多径信道仿真.m

来自「多径信道仿真以及产生GOLD序列程序,需要的朋友下来看看,资源互换最总要.」· M 代码 · 共 32 行

% *****************************************************
%将白噪声过程经过低通滤波器来仿真***滤波器是IIR滤波器
% *****************************************************
echo on
p=0.99;
N=1000;
d=5;
A=[1 -2*p p^2];
B=(1-p)^2;
for i=1:N,
   [white_noise_seq1(i) white_noise_seq2(i)] = gngauss;
   echo off;
end;
echo on;
b1=filter(B,A,white_noise_seq1);
b2=filter(B,A,white_noise_seq2);
c=b1(d+1:N)+b2(1:N-d);
%p的值控制着滤波器的带宽当p靠近1时窄0宽
p=0.9;
A=[1 -2*p p^2];
B=(1-p)^2;
for i=1:N,
   [white_noise_seq1(i) white_noise_seq2(i)] = gngauss;
   echo off;
end;
echo on;
new_b1=filter(B,A,white_noise_seq1);
new_b2=filter(B,A,white_noise_seq2);
new_c=new_b1(d+1:N)+new_b2(1:N-d);
plot(1:N,b1,'-',1:N,b2,'--',1:N-d,c,':') ; 
pause % p=0.99 ; Press a key to see results with pole at p=0.9
plot(1:N,new_b1,'-',1:N,new_b2,'--',1:N-d,new_c,':') ;