如何使用 MATLAB Butterworth 濾波

據我所知，MATLAB似乎沒有專用的陷波濾波器功能。 Matalb 中的濾波器功能都是常規的模擬和 FIR 和 IIR 濾波器。 模擬電路中的陷口往往是雙T形的，這是通過相位抵消來實現的。

我不認為MATLAB的設計者會對這些“奇怪的技巧”感興趣。

正式的方法是使用帶阻濾波器函式來實現它，我就不說它的優缺點了。

freqs（） 是將資料從時域傳輸到頻域。

WA 是您的輸入。

b、a]為低通濾波器引數。

h 是輸出。

該程式以 20*log10 的形式繪製。

h = freqs（b， a， w） 根據係數向量計算返回濾波器的復頻域響應。

呼叫此低通濾波器時，請使用以下函式。

sf=filter(a,b,s);s 是要濾波的資料，sf 是通過您設計的低通濾波器後的新資料。

數字還是模擬？ RP是通帶紋波，RS是阻帶衰減，兩者均以dB為單位。

巴特沃斯帶通濾光片設計：

一。 請求。

中心頻率：40MHz

通帶上限：42MHz; 下限：38MHz

頻寬：4MHz

帶內波動AP<

阻帶衰減：>60dB

阻抗：50

二。 設計步驟。

已知 fda=42MHz，fdb=38MHz，BWD=，頻寬=10%

1.首先找到幾何中心頻率：

2.根據公式，我們可以得到kl、kc和頻寬變換常數ka

其中：fc=fb]。

3.選擇過濾器型別：從 P5、A37、A1-12 中選擇 CO5-10-27其中 as=， ap=頻寬： p=10%滿足上述條件。

4.找 bw2， bw4

bw2 = 2 bwd [其中 2 存在於 p37、a1-12、CO5-10-27 中]。

BW4 = 4 BWD [其中 4 存在於 P35、A1-12、CO5-10-27 中]。

5.找到 f2a， f2b， f4a， f4b， [2a]， [2b]， 4a]， [4b]。

6.計算電路的具體引數。

L1 = L1 ka kl（其中 L1 = CO5-10-27 中的 C1）。

L2A = L2 ka KL [2A]（其中 L2 = CO5-10-27 中的 C2）。

l2b=l2ˊ×ka×kl×[2b]

L3 = L3 ka kl（其中 L3 = CO5-10-27 中的 C3）。

L4A = L4 ka Kl [4A]（其中 L4 = CO5-10-27 中的 C4）。

l4b=l4ˊ×ka×kl×[4b]

l5=l5ˊ×ka×kl

自行更改引數。

讓我們用過濾器原型公式做乙個程式模仿。

你有乙個非常高頻的帶通濾波器，所以讓我們用波導來做。

簡介 巴特沃斯濾波器是一種電子濾波器。 巴特沃斯濾波器的特點是通帶中最平滑的頻率響應曲線。 這種濾波器最早是由英國工程師史蒂芬·巴特沃斯（Stephen Butterworth）在1930年發表在英國《無線電工程》雜誌上的一篇文章中提出的。

巴特沃斯濾波器特點巴特沃斯濾波器的特點是在通帶內頻率響應曲線平坦，無波動，而阻帶內的頻率響應逐漸減小到零。 在振幅的對數對角線頻率的波圖上，從一定的邊界角頻率開始，振幅隨著角頻率的增加而逐漸減小，趨於負無窮大。 巴特沃斯低通濾波器用於平滑受雜訊干擾的影象。

例項百分比 by lyqmath

dlut school of mathematical sciences 2008

clc; clear all; close all;

i=imread('');

figure;

subplot(2, 2, 1); imshow(i); title('原始影象');

j1 = imnoise(i, 'salt & pepper');堆疊的鹽和胡椒噪音。

subplot(2, 2, 2); imshow(j1); title('為影象新增噪點');

g = fft2(double(j1));傅利葉變換。

g = fftshift(g);轉換資料矩陣。

m, n]=size(g);

nn = 2;%二階巴特沃茲低通濾波器。

d0 = 20;

m = fix(m/2);

n = fix(n/2);

for i = 1:m

for j = 1:n

d = sqrt((i-m)^2+(j-n)^2);

h = 1/(1+;計算低通濾波器傳遞函式。

result(i,j) = h*g(i,j);

t(i, j) = h;

endend

result = ifftshift(result);

j2 = ifft2(result);

j3 = uint8(real(j2));

subplot(2, 2, 3); mesh(t); title('濾波器示意圖')

box on;

顯示過濾後的影象。

subplot(2, 2, 4); imshow(j3); title('篩選結果'結果。

2 在 WP 30 fs 中，30 是通帶截止頻率，那個改為 2000，fs 略大於 4000。

看到 wp=2*30 fs？

30是截止平流層。

你太棒了，我什至看不懂。