%程式設計
%任務書中給出的要求為中心頻率200hz,頻寬150hz。
%故設上通帶截止頻率為110hz,下通帶截止頻率290hz,阻帶上限頻率140hz,阻帶下限頻率260hz。
%此處僅以boxcar窗為示例,其他窗函式的程式**基本相同,
%只是在window=
boxcar
(n)、n=
ceil
(1.8
*pi/delta_w)兩處作出各個窗函式相應的修改即可。
flp=
490;
fhp=
510;
fls=
499;
fhs=
501;
fs=2048
;wlp=
2*pi*flp/fs;
whp=
2*pi*fhp/fs;
wls=
2*pi*fls/fs;
whs=
2*pi*fhs/fs;
wc=[
(wlp+wls)/(
2*pi)
,(whp+whs)/(
2*pi)];
delta1=wls-wlp;
delta2=whp-whs;
delta_w=
min(delta1,delta2);n=
ceil
(1.8
*pi/delta_w);%
//不同的窗要選擇係數不同//
n=n+
rem(n,2)
;n=0:n-1;
window=
boxcar
(n+1);
%//選擇窗函式//
[h1,w]
=freqz
(window,1)
;subplot(2
,2,1
)stem
(window,
'.')
;xlabel
('n');
title
('boxcar窗函式');
subplot(2
,2,2
)plot
(w*fs/(2
*pi),20
*log
(abs
(h1)
/abs(h1
(1))
));grid;
xlabel
('f/hz');
ylabel
('幅度(db)');
title
('boxcar窗函式的頻譜');
hn=fir1
(n,wc,
'stop'
,window)
;[h2,w]
=freqz
(hn,1,
512)
;subplot(2
,2,3
)stem
(hn,
'.')
;xlabel
('n');
ylabel
('h(n)');
title
('boxcar窗函式的單位脈衝響應');
subplot(2
,2,4
)plot
(w*fs/(2
*pi),20
*log
(abs
(h2)
/abs(h2
(1))
));grid;
xlabel
('f/hz');
ylabel
('幅度(db)');
title
('boxcar帶阻濾波器的幅頻特性'
);
matlab設計模擬帶阻濾波器
簡單記錄下在matlab上如何設計出模擬的帶阻濾波器,包括 巴特沃斯濾波器 切比雪夫i型濾波器 切比雪夫ii型濾波器 橢圓型濾波器。設計帶阻濾波器 巴特沃斯 切比雪夫i型 切比雪夫ii型 橢圓型濾波器 clear all wp和ws分別是通帶和阻帶的頻率 截止頻率 當wp ws時,為高通濾波器 當w...
論高通濾波器,帶阻濾波器,陷波濾波器
首先,對一副影象進行如下二維傅利葉變換。我們將u 0和v 0帶上式,我們可以得到如下式子。根據上式,可以到f 0,0 的值是非常大的。這裡,我們將 f 0,0 稱為直流分量,直流分量比其他的成分要大好幾個數量級。所以,這也就是傅利葉譜為什麼需要使用對數變換才能看清楚的原因。這裡,對於高通濾波器而言,...
matlab設計濾波器
陳懷琛的 數字訊號處理教程 matlab釋義與實現 從效能上來說,iir濾波器傳遞函式包括零點和極點兩組可調因素,對極點的惟一限制是在單位圓內。因此可用較低的階數獲得高的選擇性,所用的儲存單元少,計算量小,效率高。但是這個高效率是以相位的非線性為代價的。選擇性越好,則相位非線性越嚴重。fir濾波器傳...