www.pudn.com > IS95Gausewhitenoise.rar > IS95Gausewhitenoise.m, change:2010-05-07,size:1198b


%高斯分布白噪声经IS-95滤波器 
clear; 
clc; 
detat=0.128; 
t=0:detat:128; 
y0=(randn(1,length(t)));  
y1=xcorr(y0); 
Y0=fft(y1); 
S0=abs(Y0); 
W0=0:1:2000; 
m=[-25288315,-34167931,-35752323,-16733702,21602514,64938487,91002137,81894974,37071157,-21998074,-60716277,-51178658,7874526,84368728,126869306,94528345,-12839661,-143477028,-211829088,-140513128,94601918,441387140,785875640,1000000000,1000000000,785875640,441387140,94601918,-140513128,-211829088,-143477028,-12839661,94528345,126869306,84368728,7874526,-51178658,-60716277,-21998074,37071157,81894974,91002137,64938487,21602514,-16733702,-35752323,-34167931,-25288315]*0.000000001; 
y=filter(m,1,y0); 
y2=xcorr(y); 
Y=abs(fft(y2)); 
S=Y(1:1001); 
W=0:1:1000; 
 
 
 
figure(1); 
subplot(2,1,1),plot(t,y0), 
title('Gause white noise'), 
xlabel('t'),ylabel('y0'); 
grid on; 
subplot(2,1,2),plot(t,y), 
title('noise at the output of IS95 FIR'), 
xlabel('t'),ylabel('y'); 
grid on; 
 
figure(2); 
semilogy(W0/1000,S0/1000), 
title('power spectrum of Gause white noise'), 
xlabel('W'),ylabel('S0'); 
grid on; 
 
figure(3); 
semilogy(W/500,S/1000), 
title('power spectrum at the output of IS95 FIR'), 
xlabel('W'),ylabel('S'); 
grid on;