四种低通滤波器

107 阅读 0 评论 71 点赞

我是靠谱客的博主精明牛排，这篇文章主要介绍四种低通滤波器，现在分享给大家，希望可以做个参考。

今天介绍四种低通滤波器对信号频率成分进行滤波，首先进行正傅里叶变换，设计一个滤波器用点操作的方法加工频谱数据，然后进行反傅里叶变换。关键是求出传递函数H（u,v）。

一.理想低通滤波器

matlab代码：

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
clc;
close all;
clear all;
I = rgb2gray(imread('plane.jfif'));
subplot(131);imshow(I);title('原图');
J = fftshift(fft2(I));
subplot(132);imshow(log(abs(J)),[8,10]);title('傅里叶变换后');
[M,N] = size(I);
H = zeros(M,N);
for u =1:M
    for v = 1:N
        if sqrt((u-M/2)^2+(v-N/2)^2)<=30
            H(u,v) = 1;
        else
            H(u,v) = 0;
       
        end;
    end;
end;
I1 = H.*J;
I1 = ifftshift(I1);
I2 = real(ifft2(I1));
subplot(133);imshow(I2,[]);title('理想低通滤波器');

二.巴特沃斯低通滤波器

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
clc;
close all;
clear all;
I = rgb2gray(imread('plane.jfif'));
subplot(131);imshow(I);title('原图');
J = fftshift(fft2(I));
subplot(132);imshow(log(abs(J)),[8,10]);title('傅里叶变换后');
[M,N] = size(I);
H = zeros(M,N);
n1 = floor(M/2);n2 = floor(N/2);
for i =1:M
    for j = 1:N
        d = sqrt((i-n1)^2+(j-n2)^2);
        H = 1/(1+(d/30)^(2*2));
        I1(i,j)=H*J(i,j);
    end;
end;
I1 = ifftshift(I1);
I2 = real(ifft2(I1));
subplot(133);imshow(I2,[]);title('巴特沃斯低通滤波器');

三.指数低通滤波器

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
clc;
close all;
clear all;
I = rgb2gray(imread('plane.jfif'));
subplot(131);imshow(I);title('原图');
J = fftshift(fft2(I));
subplot(132);imshow(log(abs(J)),[8,10]);title('傅里叶变换后');
[M,N] = size(I);
H = zeros(M,N);
n1 = floor(M/2);n2 = floor(N/2);
for i =1:M
    for j = 1:N
        d = sqrt((i-n1)^2+(j-n2)^2);
        o = -(d/30)^2;
        H = exp(1)^o;
        I1(i,j)=H*J(i,j);
    end;
end;
I1 = ifftshift(I1);
I2 = real(ifft2(I1));
subplot(133);imshow(I2,[]);title('指数低通滤波器');

四.梯形低通滤波器

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
clc;
close all;
clear all;
I = rgb2gray(imread('plane.jfif'));
subplot(131);imshow(I);title('原图');
J = fftshift(fft2(I));
subplot(132);imshow(log(abs(J)),[8,10]);title('傅里叶变换后');
[M,N] = size(I);
H = zeros(M,N);
for u =1:M
    for v = 1:N
        if sqrt((u-M/2)^2+(v-N/2)^2)<30
            H(u,v) = 1;
        elseif 30<= sqrt((u-M/2)^2+(v-N/2)^2)<=180
            H(u,v) = (sqrt((u-M/2)^2+(v-N/2)^2)-180)/(30-180);
        else
            H(u,v)=0;
        end;
    end; 
 
end;
I1 =H.*J;
I1 = ifftshift(I1);
I2 = real(ifft2(I1));
subplot(133);imshow(I2,[]);title('梯形低通滤波器');