概率论与随机信号分析3

号经过系统 ys=filter(b,a,ss); subplot(2,3,2) plot(t,ys) title(' 输出信号'); grid on subplot(2,3,5); fys=abs(fft(ys)); plot(f1,fys(1:n1)); title(' 输出信号频谱'); grid on; % 信号直方图统计 subplot(2,3,3); hs=ss; hist(hs,20); Title(' 滤波前信号大小直方图'); subplot(2,3,6); hs=ys; hist(hs,20); title(' 滤波后信号大小直方图') % 噪声及其功率谱 figure; s n=rand(1,n); sn=2*sn-1; subplot(2,3,1) plot(t,sn); title(' 输出噪声 N(0,1) '); grid on; subplot(2,3,4); rn=xcorr(sn,sn,'biased'); fn=abs(fft(rn)); f2=(0:n-1)/(2*n+1)*fs; plot(f2,fn(1:n)); title(' 输入噪声功率谱'); grid on; yn=filter(b,a,sn); subplot(2,3,2) plot(t,yn); title(' 输出噪声'); grid on subplot(2,3,5) rn=xcorr(yn,yn,'biased'); fn=abs(fft(rn)); plot(f2,fn(1:n)); title(' 输出噪声功率谱'); grid on; subplot(2,3,3); hs=sn; hist(hs,20); title(' 滤波前噪声大小直方图') subplot(2,3,6); hs=yn; hist(hs,20); title(' 滤波后噪声大小直方图')