llgd.net
当前位置:首页 >> mAtlAB维纳滤波 >>

mAtlAB维纳滤波

%实例一程序_维纳滤波的计算机实现 %初步处理,并接受输入数据,包括信号样本个数L和滤波器阶数N clear all close all L=input('L='); N=input('N='); a=0.95; %定义w,v,u w=sqrt(3*(1-a^2))*(2*rand(1,L)-1); v=sqrt(3)*(2*rand(1,L)-1); u=o...

这是matlab图像处理工具箱中的二维维纳滤波函数 其中J1表示待滤波图像矩阵 [5,5]表示滤波器窗口大小 k表示滤波后图像矩阵 noise表示噪声功率估计值 进一步可以查看matlab自带帮助

clear; I=imread('img_orignal.tif'); figure; subplot(2,2,1);imshow(I);title(' 原图像 '); [m,n]=size(I);F=fftshift(fft2(I)); k=0.005; for u=1:m for v=1:n H(u,v)=exp((-k)*(((u-m/2)^2+(v-n/2)^2)^(5/6))); end end G=F.*H; I0=real(iff...

%实例一程序_维纳滤波的计算机实现 %初步处理,并接受输入数据,包括信号样本个数L和滤波器阶数N clear all close all L=input('L='); N=input('N='); a=0.95; %定义w,v,u w=sqrt(3*(1-a^2))*(2*rand(1,L)-1); v=sqrt(3)*(2*rand(1,L)-1); u=o...

clc;clear all; %读原始图像% format long Blurred=imread('fig525(b).bmp'); subplot(1,2,1);imshow( Blurred);title('原图像'); %自编函数进行维纳滤波% k=0.0025; [m,n]=size(Blurred); spectrum=zeros(m,n); H=zeros(m,n); for u=1:m for v=...

PSF = fspecial('motion',len,ang); %建立扩散子,其中len是模糊长度,ang是模糊角度 img2=deconvlucy(img,PSF,n); %用lucy-richardson方法复原图像,其中img是运动模糊图像,PSF是扩散子,n是迭代次数,img2是复原图像

初步看了下,是LMS算法

可能应该改一下21和11吧 一个是模糊尺度 一个是模糊方向 。。但是如果你不知道PSF 应该用盲复原啊 用deconvblind函数

已经是维纳滤波了

clc;clear all; %读原始图像% format long Blurred=imread('fig525(b).bmp'); subplot(1,2,1);imshow( Blurred);title('原图像'); %自编函数进行维纳滤波% k=0.0025; [m,n]=size(Blurred); spectrum=zeros(m,n); H=zeros(m,n); for u=1:m for v=...

网站首页 | 网站地图
All rights reserved Powered by www.llgd.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com