来人啊来人啊,很简单的
close all;clc;
x=0:0.02:10;
x=reshape(x,length(x),1);
noise=zeros(size(x));
bgIdx=30;
endIdx=60;
noise(bgIdx:endIdx)=3*randn(endIdx-bgIdx+1,1);
y=0.2*x.^2+2*x+noise;
y1=0.2*x.^2+2*x;
figure,plot(x,y,'r-');
ylim([-5,50]);
%%
alpha=0.5;
lambda=100;
smallNum=0.1;
tic;
y2 = wlsFilter_1d(y,lambda,alpha,y/max(y(:)),smallNum);
toc;
t=toc-tic;
figure,plot(x,y2,'r-');
ylim([-5,50]);
让我知道什么是什么就行,这是信号滤波的代码,定义了一个函数,下面输这个就出图了
close all;clc;
x=0:0.02:10;
x=reshape(x,length(x),1);
noise=zeros(size(x));
bgIdx=30;
endIdx=60;
noise(bgIdx:endIdx)=3*randn(endIdx-bgIdx+1,1);
y=0.2*x.^2+2*x+noise;
y1=0.2*x.^2+2*x;
figure,plot(x,y,'r-');
ylim([-5,50]);
%%
alpha=0.5;
lambda=100;
smallNum=0.1;
tic;
y2 = wlsFilter_1d(y,lambda,alpha,y/max(y(:)),smallNum);
toc;
t=toc-tic;
figure,plot(x,y2,'r-');
ylim([-5,50]);
让我知道什么是什么就行,这是信号滤波的代码,定义了一个函数,下面输这个就出图了