1、二维卷积矩阵

T=convmtx2(H,m,n)

H=magic(3);

>> T=convmtx2(H,1,2)

T =

(1,1)        8

(2,1)        3

(3,1)        4

(4,1)        1

(5,1)        5

(6,1)        9

(7,1)        6

(8,1)        7

(9,1)        2

(4,2)        8

(5,2)        3

(6,2)        4

(7,2)        1

(8,2)        5

(9,2)        9

(10,2)        6

(11,2)        7

(12,2)        2

2、filter2函数

filter2函数主要进行二维线性滤波操作

Y=filter2(h,x)

函数表示使用二维FIR滤波器h对矩阵X进行滤波操作，X为矩阵时，h通常有fspecial函数输出得到。

Y=filter2(h,X,shape)

shape表示Y的形式，full表示Y的维数大于X的维数，same表Y的维数等于x的维数，valid表示Y的维数小于X的维数，

X=magic(3);

>> h=fspecial('motion',20,45);

>> Y=filter2(h,X)

Y =

0.4094 0.4793 0.9311

0.4793 0.9311 0.9365

0.9311 0.9365 0.2985

3、fspecial函数

h=fspecial(type)函数，产生呢type类型的二维滤波器H

h=fspecial(type,parameters)函数,设定滤波器的参数

I=imread('lena.jpg');

subplot(2,2,1);

imshow(I);

title('originnal image');

H=fspecial('motion',20,45);

M=imfilter(I,H,'replicate');

subplot(2,2,2);

imshow(M);

H=fspecial('disk',10)%圆形均值滤波器

B=imfilter(I，H,'replicate');

subplot(2,2,3);

imshow(B);

H=fspecial('unsharp');%生成锐化滤波器

S=imfilter(I,H,'replicate');

subplot(2,2,4);

imshow(S)

4、imfilter函数

B=imfilter(A,H),使用多维滤波器H对图像A进行滤波   H通常通过fspecial函数产生

代码上面的

5、freqz2  二维频率响应、

功能计算滤波器的二维平率响应、

【 H f1 f2]=freqz2(h,n1,n2);

计算二维FIR滤波器h的n2*n1频率响应H,以及频率向量f2 f1

H=zeros(16,16);

H(5:12,5:12)=1;

H(7:10,7:10)=0;

h=fwind1(H,bartlett(16))；%用一维bartlett床涉及二维FIR滤波器

freqz2(h,[32 32]);

axis([-1 1 -1 1 0 1]);

6\dct2函数  二维离散余弦变换

B=dct2(A,m,n)

计算图像A的离散余弦变换系数矩阵B,m n指将A增0补到m*n维

I=imread('lena.jpg');

I=rgb2gray(I);

J=dct2(I);

imshow(log(abs(J)),[]);

colormap(jet(64))

colorbar;

7\idct2(J)  表示反变换，用于图像的重构

I=imread('lena.jpg');

I=rgb2gray(I);

J=dct2(I);

J(abs(J)<10)=0;

K=idct2(J);

imshow(I);

figure;

imshow(K,[]);