资源描述：
《MATLAB编程经验》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、MATLAB编程经验1.由imread读进来的彩色图像一般是RGB格式，是一个三维矩阵，其每个元素都uint8类型，而matlab的基本矩阵运算和傅立叶变换都不支持uint8类型。因此必须将矩阵转化为double类型才能参与运算。用到的函数：rgb2gray：RGB转256阶灰度double()：强制类型转换2.Matlab的imshow只能显示像素值为uint8的彩色或黑白图像，如果用double类型会显示出一片白色，如果用uint16类型会显示出一片黑色用到的函数：uint8()：强制类型转换3.利用傅立叶变换通过快速傅立叶变换进行矩阵运算（1）计算
2、各子矩阵与模板矩阵的元素乘积和（即互相关）：abw=abs(ifft2(fft2(bw).*fft2(rot90(a,2),m,n)));所得矩阵abw的各元素abw(i,j)的值为以bw(i,j)为右下角元素的子矩阵与模板矩阵的元素乘积和（2）计算各子矩阵的元素平方和：bw2=abs(ifft2(fft2(bw.^2).*fft2(ones(size(a)),m,n)));bwsqrt=sqrt(bw2);所得矩阵bw2的各元素bw2(i,j)的值为以bw(i,j)为右下角元素的子矩阵的各元素的平方和的平方根（3）计算模板矩阵的元素平方和的平方根：as
3、qrt=sqrt(sumsqr(a));（1）计算各子矩阵与模板矩阵的相关系数：T=abw./bwsqrt./asqrt;所得矩阵T的各元素T(i,j)的值为以bw(i,j)为右下角元素的子矩阵与模板矩阵的相关系数（2）找出T中元素的最大值所在的位置(i,j)即为模板在大图中的右下角的位置Tmax=max(max(T));[r,c]=find(T>=Tmax);（3）计算矩阵中以(i,j)为右下角，大小为N*N的子矩阵中各元素的和T=real(ifft2(fft2(result(:,:,dimension)).*fft2(ones(N,N),m,n)))
4、;（4）计算矩阵中以(i,j)为“中心”，大小为N*N的子矩阵中各元素的和T=real(ifft2(fft2(result(:,:,dimension)).*fft2(ones(N,N),m,n)));fori=1:1:mforj=1:1:n%result(i,j)=floor((T(i,j)+T(mod(i+N,m),j)+T(i,mod(j+N,n))+T(mod(i+N,m),mod(j+N,n)+1))/N^2/4);irange=i+N;jrange=j+N;ifirange>mirange=mod(irange,m);endifjrange>
5、njrange=mod(jrange,n);endresult(i,j,dimension)=floor((T(i,j)+T(irange,j)+T(i,jrange)+T(irange,jrange))/N^2/4);endend
