[求助]lms图像滤波程序调试……

[edward2000]

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('预测信号,实际信号');




我调试不出来了，咋弄啊，费劲了，我就是想把一个加入了高斯噪声的图像初步滤波，然后我好进行我的仿真，咋就出不来呢！我都疯掉了