![]() |
[求助]lms图像滤波程序调试……
function LMS_adaptive_filter;%定义功能函数
% samples_per_period=40;%周期采样频率 % N=10*samples_per_period-1;%总采样次数 % signal=sin(2*pi*0.02*[0:N-1]);%初始输入信号 x=imread('e:\22.jpg'); signal= rgb2gray(x) figure;%作图1的第一子图 subplot(2,1,1); imshow(signal); grid; title('自适应滤波器的理想输入'); % nvar=0.5;%噪声方差 % noise=nvar*randn(1,N);%噪声信号 % X=signal+noise;%带加性噪声的输入信号 % delayX=[0 X];%输入信号延时 % SIGNAL=[signal 0]; X=imnoise(x,'gaussian'); M=32;%滤波器长度 mu=0.035;%迭代步长 subplot(2,1,2);%作图1的第二子图 imshow(X); grid; title('自适应滤波器的噪声输入'); initial_status=initlms(zeros(1,M),mu);%滤波器设初值 [Y,e,initial_status]=adaptlms(X,SIGNAL,initial_status);%LMS自适应滤波器输出 figure;%作图2的第一子图 subplot(2,1,1); plot(0:N-1,Y,'.',0:N,SIGNAL,'k'); grid; title('预测信号与实际信号的比较'); legend('预测信号,实际信号'); 我调试不出来了,咋弄啊,费劲了,我就是想把一个加入了高斯噪声的图像初步滤波,然后我好进行我的仿真,咋就出不来呢!我都疯掉了 |
回复: [求助]lms图像滤波程序调试……
initial_status=initlms(zeros(1,M),mu);%滤波器设初值
[Y,e,initial_status]=adaptlms(X,SIGNAL,initial_status); 尤其是这两句,不知道什么意思,也不知道用法 |
所有时间均为北京时间。现在的时间是 12:40。 |
Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.