samp9_14.m

MATLAB7.x数字信号处理 光盘内容

%Samp9_14
clf;Fs=1000;  %信号采样频率
N=1024;Nfft=256;n=0:N-1;t=n/Fs; %数据长度、分段长度和时间序列
window=hanning(256);   %采用的窗口
noverlap=128;    %重叠数据个数
dflag='none';
randn('state',0);  %设置产生随机数的状态
xn=sin(2*pi*50*t)+randn(1,N);   %输入信号
h=ones(1,10)/10;   %设计一个滤波器的脉冲响应
yn=filter(h,1,xn);  %滤波后信号
%采用自功率谱和互功率谱估计的方法
[Pxx,f1]=psd(xn,Nfft,Fs,window,noverlap,dflag);  %估计自功率谱
[Pxy,f1]=csd(xn,yn,Nfft,Fs,window,noverlap,dflag);  %估计互功率谱
He=Pxy./Pxx;   %求得系统传递函数
[HEST,f]=tfe(xn,yn,Nfft,Fs,window,noverlap,dflag); %采用tfe函数估计传递函数
%比较结果
H=freqz(h,1,f,Fs);   %求得系统的频率响应
subplot(3,1,1),plot(f,abs(H));  %绘制原滤波器幅频特性
ylabel('振幅');
title('实际模型的幅频响应');
axis([0 500 0 1]);grid on;
subplot(3,1,2),plot(f,abs(HEST)); %绘制用tfe函数估计的传递函数幅频特性
ylabel('振幅');
title('TFE函数估计的幅频响应');
axis([0 500 0 1]);grid on
subplot(3,1,3),plot(f1,abs(He)); %绘制用功率谱估计的传递函数幅频特性
xlabel('频率/Hz');ylabel('振幅');
title('功率谱估计的幅频响应');
axis([0 500 0 1]);grid on