資源簡(jiǎn)介
自己寫的線性相位的FIR低通和帶通的數(shù)字濾波器程序,絕對(duì)沒有錯(cuò)誤,注釋非常詳細(xì)。用的是窗函數(shù)法,畫了單位沖擊響應(yīng)和幅頻響應(yīng)的圖。

代碼片段和文件信息
%?設(shè)計(jì)線性相位的FIR帶通濾波器
%?2014.03.12?
clcclear?all;
wls=0.2*pi;????????????????????????????%?阻帶最小頻率?
wlp=0.35*pi;???????????????????????????%?通帶最小頻率
whp=0.65*pi;???????????????????????????%?通帶最大頻率
whs=0.8*pi;????????????????????????????%?阻帶最大頻率?
wc=[wlp/piwhp/pi];????????????????????%?帶通濾波器的歸一化截止頻率
B=wlp-wls;?????????????????????????????%?過渡帶寬=低端通帶截止頻率-低端阻帶截止頻率?
N=ceil(6.2*pi/B);??????????????????????%?hanning窗階數(shù)N,由過渡帶寬B=8*pi/N算出Nceil為取整函數(shù)
window=hanning(N);?????????????????????%?窗選擇hanning窗,為N*1的矩陣/向量
n=0:1:N-1;?
[h1w]=freqz(window1);????????????????%?window是N*1矩陣,h1是hannning窗的頻率響應(yīng)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(221);????????????????????????%?將多個(gè)圖畫在一個(gè)平面上
stem(window);??????????????????????????%?繪制離散序列圖
axis([0?45?0?1.2]);????????????????????%?定義xy軸范圍?
grid;??????????????????????????????????%?顯示網(wǎng)格線
xlabel(‘n‘);???????????????????????????%?x軸命名
title(‘Hanning窗函數(shù)‘);???????????????
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(222);
plot(w/pi20*log(abs(h1)/abs(h1(1))));?%?w來自freqz函數(shù),為N個(gè)頻率等分點(diǎn)的值(N*1矩陣),后面計(jì)算的是幅度值
axis([0?1?-350?0]);
grid;
xlabel(‘w/pi‘);
ylabel(‘幅度(dB)‘);
title(‘Hanning?窗函數(shù)的頻譜‘);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
hn=fir1(N-1wchanning(N));????????????%?設(shè)計(jì)filter。函數(shù)fir1返回的是filter的系數(shù)h(n)即沖激相應(yīng)
[h2w]=freqz(hn1512);????????????????%?filter的頻率響應(yīng)h2
subplot(223);
stem(nhn);
axis([0?45?-0.25?0.25]);
grid;
xlabel(‘n‘);
ylabel(‘h(n)‘);
title(‘Hanning?窗函數(shù)的單位脈沖響應(yīng)‘);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(224);
plot(w/pi20*log(abs(h2)/abs(h2(1))));??%?w來自freqz(與h2對(duì)應(yīng))
grid;
xlabel(‘w/pi‘);
ylabel(‘幅度(dB)‘);
title(‘Hanning?窗設(shè)計(jì)帶通濾波器的幅頻響應(yīng)‘);
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????2034??2014-03-13?21:33??FIR_BPF.m
?????文件???????1892??2014-03-13?21:34??FIR_LPF.m
-----------?---------??----------?-----??----
?????????????????3926????????????????????2
- 上一篇:移相全橋simuli
nk仿真文件 - 下一篇:汽車二自由度狀態(tài)方程模型
評(píng)論
共有 條評(píng)論