| 查看: 346 | 回复: 1 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】求助滤波器的幅频响应与实际幅频响应
|
|||
|
设FIR滤波器的脉冲响应为元素均为1/10、长度为10的信号。产生一个长度为1024的带有白噪声的周期信号，并采用上述滤波器进行滤波，采用输入信号和滤波后的输出数据估计该滤波器的幅频响应并与实际幅频响应作比较。 [ Last edited by wenzhenzhong on 2009-12-5 at 01:01 ] |
回复此楼» 猜你喜欢
各位大佬，面上项目预算多少合适？预算有用吗？
已经有21人回复
-
工学085800能动专硕，初试307分求调剂，满足B区国家线
已经有1人回复
-
金属材料论文润色/翻译怎么收费?
已经有163人回复
-
0856材料306分调剂
已经有4人回复
-
0858专硕.339分求调剂
已经有0人回复
-
非晶合金
已经有2人回复
-
面上参与人员重要吗？找不到参与人
已经有7人回复
-
【博士招生】英国拉夫堡大学博士招生 （LU-CSC scholarship, 暂免ATAS）
已经有0人回复
-
北京理工大学 空天学院-招收发动机相关专业（力学、燃烧传热等方向）博士生
已经有2人回复
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事(长期有效)
已经有10人回复
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事
已经有10人回复
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
5