设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 980|回复: 0

[讨论] 数字滤波器程序问题,麻烦大家帮看看

[复制链接]

该用户从未签到

尚未签到

发表于 2010-3-30 09:25:41 | 显示全部楼层 |阅读模式

马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

您需要 登录 才可以下载或查看,没有账号?立即加入

×
clear all' E" K+ |' w  p$ B$ C
%-------------------------------------------------------------* u. G7 b4 e. X. s) f# P7 c
N=500;
8 `! K$ ^3 \  Y3 F  w' l9 gp=0.05;
6 B$ g% q6 u) ?' D" ~* Tf=1/50;
  u9 ~  D' R3 tn=1:N;4 M3 V4 D$ ^; c! Y' a( B0 E" K0 A
u=rand(1,N);
6 Q, F  n2 J" _& Ks=sin(2*pi*f*0:N-1);# Z) ]" w! L# u' p$ t
x=u(1:N)+s;
8 n* ~6 V& Y- O9 }) d/ t4 L& |subplot(311)7 d' r% v* r+ a4 x! a3 N- R1 K' q
plot(x(1:500));grid on;
$ I! x4 Z" O6 y' ?% D7 `3 ]8 ~ylabel('x(n)'). N( ^, Q* j. X6 I
%---------------------------------------------------------------
& I; w6 _2 P0 l- ifp=200;8 ]0 y, R- u0 [
fs=400;Fs=1200;" c3 e( K# H, R' g: p5 X3 b
rp=3;  %通带衰减不大于3dB( Z6 a, ^& @* [0 o
rs=35;  %阻带衰减不大于35dB. M8 ?2 _) ^  }/ j  z( `; c
wp=2*fp/Fs;
  x  l- P4 p6 ^0 x& n0 u1 Gws=2*fs/Fs;& Z0 _4 `8 i+ O. m  u
syms a w real
0 |! J/ U1 N8 T) s( o  V3 O[nz,wn]=buttord(wp,ws,rp,rs);
. d2 k- A# f, H2 T* a[bz,az]=butter(nz,wp,'low');: Y$ C2 e5 E9 I; F
[h,w]=freqz(bz,az,1024,Fs);
. }/ g9 n+ H/ Z  ~  B1 Jih=ifourier(h);' t2 z5 _! [& S/ `
subplot(312)
2 k" w' ]  i2 O8 c6 _plot(x,h);grid;
. B+ X+ q  F5 r: J$ M+ Vylabel('Low-pass DF')$ J2 Y  J( z: ~1 C
%--------------------------------------------------------" _1 A7 k5 I9 g  E' q2 S
y=fftfilt(h,x);* I1 G% p! o7 f* ?2 q  j
subplot(313)# h7 l. x% c3 U
plot(y);* z$ r' |# `( Q2 _! }; h7 a
4 y2 \2 P1 G4 n+ O/ U) U
画出的图不对,好像是逆傅里叶变换用的不太对似的。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
您需要登录后才可以回帖 登录 | 立即加入

本版积分规则

招聘斑竹

小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

GMT+8, 2026-3-19 08:22

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表