clear?all;close?all;
xt=xtg;N=1000;
fp=120;fs=150;Rp=0.1;As=60;Fs=1000;%指標
%======================窗函數========================
B=2*pi*（fs-fp）/Fs;
M=ceil（11*pi/B）;%blackman窗過渡帶寬度Bt精確值為11pi/N
wc=（fp+fs）/Fs;?????????????????????%書7.2.20
hn=fir1（Mwcblackman（M+1））;
HW=fft（hn1024）;
Hw=abs（HW）;
ywt=fftfilt（hnxtN）;
%畫圖
figure（2）
subplot（211）
n=0:M;stem（nhn‘.‘）;grid?on;
axis（[0M+11.2*min（hn）1.2*max（hn）]）;
xlabel（‘n‘）;ylabel（‘h（n）‘）;title（‘（a）單位脈沖響應‘）
subplot（212）;
f=[0:1023]*Fs/1024;
plot（f20*log10（Hw））;grid?on;
axis（[0Fs/2-12010]）;
xlabel（‘f/Hz‘）;ylabel（‘幅度（dB）‘）;title（‘（b）低通濾波器幅頻特性‘）
figure（3）
subplot（211）
t=[0:N-1]/Fs;Tp=N/Fs;
plot（tywt）;grid;
axis（[0Tp/2-11]）;xlabel（‘t/s‘）;ylabel（‘y_w（t）‘）;
title（‘（c）濾除噪聲后的信號波形‘）
subplot（212）
YWn=abs（fft（ywt1024））;
f=[0:1023]*Fs/1024;
stem（fYWn‘.‘）;grid
axis（[0Fs/201.2*max（YWn）]）;
xlabel（‘f‘）;ylabel（‘|Y_w（n）|‘）;title（‘（d）輸出信號頻譜‘）
%======================等波紋最佳逼近========================
fb=[fpfs];m=[10];?		%?確定remezord函數所需參數fmdev
dev=[（10^（Rp/20）-1）/（10^（Rp/20）+1）10^（-As/20）];
[MefomoW]=remezord（fbmdevFs）;	%?確定remez函數所需參數
hn=remez（MefomoW）;		%?調用remez函數進行設計
Hw=abs（fft（hn1024））;		%?求設計的濾波器頻率特性
yet=fftfilt（hnxtN）;???????%?調用函數fftfilt對xt濾波
YEn=fft（yet1024）;
%畫圖
figure（4）;
subplot（211）
n=0:Me;stem（nhn‘.‘）;grid?on;
axis（[0M+11.2*min（hn）1.2*max（hn）]）;
xlabel（‘n‘）;ylabel（‘h（n）‘）;title（‘（e）單位脈沖響應‘）
subplot（212）
f=[0:1023]*Fs/1024;
plot（f20*log10（Hw））;grid;title（‘（f）低通濾波器幅頻特性‘）
axis（[0Fs/2-8010]）;
xlabel（‘f/Hz‘）;ylabel（‘幅度‘）
figure（5）
subplot（211）;
plot（tyet）;grid;
axis（[0Tp/2-11]）;xlabel（‘t/s‘）;ylabel（‘y_e（t）‘）;
title（‘（g）濾除噪聲后的信號波形‘）
subplot（212）
stem（fYEn‘.‘）;grid;
xlabel（‘f‘）;ylabel（‘|Y_e（n）|‘）;title（‘（h）輸出信號頻譜‘）
axis（[0Fs/201.2*max（YEn）]）;

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

?????文件???????1948??2014-12-12?21:01??實驗5\exp5.m

?????文件???????1269??2014-12-08?09:45??實驗5\xtg.m

?????文件?????222087??2014-12-22?20:28??實驗5\實驗五：FIR熟悉濾波器設計與軟件實現.docx

?????目錄??????????0??2014-12-22?20:29??實驗5

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

???????????????225304????????????????????4