% *****************************************************%將白噪聲過程經過低通濾波器來仿真***濾波器是IIR濾波器% *****************************************************echo onp=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.9plot(1:N,new_b1,'-',1:N,new_b2,'--',1:N-d,new_c,':') ;