| 查看: 339 | 回复: 1 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】求助滤波器的幅频响应与实际幅频响应
|
|||
|
设FIR滤波器的脉冲响应为元素均为1/10、长度为10的信号。产生一个长度为1024的带有白噪声的周期信号，并采用上述滤波器进行滤波，采用输入信号和滤波后的输出数据估计该滤波器的幅频响应并与实际幅频响应作比较。 [ Last edited by wenzhenzhong on 2009-12-5 at 01:01 ] |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有49人回复
-
filecode变化情况
已经有22人回复
-
机械工程论文润色/翻译怎么收费?
已经有68人回复
-
filecode与中标关系的预测
已经有5人回复
-
关于filecode
已经有7人回复
-
好奇怪的filecode
已经有5人回复
-
这样的filecode谁见过
已经有22人回复
-
2026国自然放榜时间
已经有9人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有6人回复
-
确定了，国自然21号放榜
已经有7人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
hu99aiai
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 194.4
- 散金: 10
- 帖子: 86
- 在线: 8.2小时
- 虫号: 797085
- 注册: 2009-06-23
- 性别: GG
- 专业: 信号与信息处理
★ ★ ★ ★ ★ ★ ★
1013902040(金币+7,VIP+0):谢谢大侠相助 6-29 22:28
1013902040(金币+7,VIP+0):谢谢大侠相助 6-29 22:28
|
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); 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)); 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 运行一下即可知道结果，可以看到估计的滤波器幅频响应与真实滤波器幅频响应具有较好的一致性 |
2楼2009-06-29 22:11:16
