fs=1500;??%采樣頻率
wp1=6*2*pi/fs;?ws1=8*2*pi/fs;?
wp2=210*2*pi/fs;?ws2=220*2*pi/fs;?
tr_width=ws1-wp1;?%過度帶寬
M=ceil（6.6*pi/tr_width）+1;?%濾波器階數
wc1=（ws1+wp1）/2;?wc2=（ws2+wp2）/2;?%通帶截止頻率
hd=ideal_lp（wc2M）-ideal_lp（wc1M）;?%帶通濾波器的設計可以看成是兩個低通濾波器想減
w_ham=（hamming（M））‘;
h=hd?.*?w_ham;?
[Hw]=freqz（h1）;
db=20*log10（abs（H））;
n=0:M-1;??
figure（1）;
subplot（221）;stem（nhd）;title（‘理想脈沖響應‘）;
axis（[0?1700?-0.1?0.1]）;xlabel（‘n‘）;ylabel（‘hd（n）‘）;
subplot（222）;stem（nw_ham）;title（‘海明窗‘）;
axis（[0?1700?0?1]）;xlabel（‘n‘）;ylabel（‘w_ham（n）‘）;
subplot（223）;stem（nh）;title（‘實際脈沖響應‘）;
axis（[0?1700?-0.1?0.1]）;xlabel（‘n‘）;ylabel（‘h（n）‘）;
figure（2）;
plot（w*fs/（2*pi）db）;title（‘幅度響應（單位：?dB）‘）;grid??%w*fs/（2*pi）再次將橫軸單位變為HZ，而非弧度，與第二行對應來看
axis（[-10?600?-200?0]）;?xlabel（‘頻率（單位：Hz）‘）;?ylabel（‘分貝‘）
set（gca‘XTickMode‘‘manual‘‘XTick‘[03500600]）
set（gca‘YTickMode‘‘manual‘‘YTick‘[-100]）?
T=2000;%采樣點數
n=0:T-1;
t=n/fs;
y=100*sin（2*pi*1*t）+50*cos（2*pi*2*t）+70*sin（2*pi*100*t）+80*sin（2*pi*150*t）;
figure（3）;
plot（y）;
x=fft（yT）;???%對輸入信號進行快速傅里葉變換，求其頻譜圖
mag=abs（x）;???%求取Fourier變換的振幅
F=n*fs/T;
figure（4）;
plot（Fmag）;?%繪出隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;title（‘N=1024‘）;grid?on;
Y=fftfilt（hy）;%調用fir濾波器
figure（5）;
plot（Y）;
X=fft（YT）;
mag=abs（X）;???%求取Fourier變換的振幅
F=n*fs/T;
figure（6）;
plot（Fmag）;?%繪出隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;title（‘N=1024‘）;grid?on;
??
??
??
??
??
??
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????388??2013-05-22?11:41??matlab設計的四種fir數字濾波器\ideal_lp.m

?????文件???????1881??2013-05-28?16:47??matlab設計的四種fir數字濾波器\gaotong.m

?????文件???????1602??2013-05-30?08:47??matlab設計的四種fir數字濾波器\daitong.m

?????文件???????3510??2013-05-28?15:58??matlab設計的四種fir數字濾波器\ditong.m

?????文件???????1607??2013-05-28?17:38??matlab設計的四種fir數字濾波器\daizu.m

?????目錄??????????0??2013-05-28?17:24??matlab設計的四種fir數字濾波器

-----------?---------??----------?-----??----

?????????????????8988????????????????????6