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;?%濾波器階數(shù)
wc1=（ws1+wp1）/2;?wc2=（ws2+wp2）/2;?%通帶截止頻率
hd=ideal_lp（wc2M）-ideal_lp（wc1M）;?%帶通濾波器的設(shè)計(jì)可以看成是兩個低通濾波器想減
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（‘理想脈沖響應(yīng)‘）;
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（‘實(shí)際脈沖響應(yīng)‘）;
axis（[0?1700?-0.1?0.1]）;xlabel（‘n‘）;ylabel（‘h（n）‘）;
figure（2）;
plot（w*fs/（2*pi）db）;title（‘幅度響應(yīng)（單位：?dB）‘）;grid??%w*fs/（2*pi）再次將橫軸單位變?yōu)镠Z，而非弧度，與第二行對應(yīng)來看
axis（[-10?600?-200?0]）;?xlabel（‘頻率（單位：Hz）‘）;?ylabel（‘分貝‘）
set（gca‘XTickMode‘‘manual‘‘XTick‘[03500600]）
set（gca‘YTickMode‘‘manual‘‘YTick‘[-100]）?
T=2000;%采樣點(diǎn)數(shù)
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）;???%對輸入信號進(jìn)行快速傅里葉變換，求其頻譜圖
mag=abs（x）;???%求取Fourier變換的振幅
F=n*fs/T;
figure（4）;
plot（Fmag）;?%繪出隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;title（‘N=1024‘）;grid?on;
Y=fftfilt（hy）;%調(diào)用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;
??
??
??
??
??
??
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1602??2013-05-30?08:47??matlab四種fir數(shù)字濾波器\daitong.m
?????文件????????1607??2013-05-28?17:38??matlab四種fir數(shù)字濾波器\daizu.m
?????文件????????3510??2013-05-28?15:58??matlab四種fir數(shù)字濾波器\ditong.m
?????文件????????1881??2013-05-28?16:47??matlab四種fir數(shù)字濾波器\gaotong.m
?????文件?????????388??2013-05-22?11:41??matlab四種fir數(shù)字濾波器\ideal_lp.m
?????文件???????????0??2019-06-14?20:38??matlab四種fir數(shù)字濾波器\四種濾波器.docx
?????目錄???????????0??2019-06-14?20:38??matlab四種fir數(shù)字濾波器\