emd端点效应处理程序
我是一emd新手,大家帮帮忙,如何用支持向量机程序来抑制信号的分解的端点效应,程序总出错,实在没有头绪,不知道怎么改,程序如下:clc;clear;close all;
fs=1024;n=1024;
%t=linspace(0,1,n);
t=0:1/fs:(n-1)/fs;
x=sin(60*pi*t)+sin(50*pi*t);
imf=emd(x);
L_IMF1=length(imf(:,1));
figure;
fori=1:L_IMF1
subplot(L_IMF1,1,i);
plot(t,imf(i,:));
set(gca,'xlim',);
set(gca,'fontname','Times New Roman','Fontsize',9');
end
a=imf(1,:);
b=imf(2,:);
num=100;%设置要延拓的点数
l=100; %设置训练样本数
fprintf('向右延拓中……\n');
for i=1:l
X(i,1:N-l)=x(i:N-l+i-1);
Y(i,1)=x(N-l+i);
end
= svr(X,Y,'linear',Inf,'einsensitive',0);
for j=0:num-1
m=0;
XX=x(l+1+j:N+j);
for i=1:l
m=m+beta(i)*(X(i,1:N-l)*XX');
end
s(N+j+1)=m+bias;
end
extend_right=x(N+1:N+num);
extend=;
求助完整格式:出错代码和出错提示 本帖最后由 邢永忠 于 2013-9-16 19:47 编辑
ChaChing 发表于 2013-5-5 10:55 http://forum.chinavib.com/static/image/common/back.gif
求助完整格式:出错代码和出错提示
clc;clear;close all;
fs=1024;n=1024;
t=linspace(0,1,n);
t=0:1/fs:(n-1)/fs;
x=sin(60*pi*t)+sin(50*pi*t);
num=100;%设置要延拓的点数
l=100; %设置训练样本数
N=100;
fprintf('向右延拓中……\n');
for i=1:l
X(i,1:N-l)=x(i:N-l+i-1);
Y(i,1)=x(N-l+i);
end
= svr(X,Y,'linear',Inf,'einsensitive',0);
for j=0:num-1
m=0;
XX=x(l+1+j:N+j);
for i=1:l
m=m+beta(i)*(X(i,1:N-l)*XX');
end
s(N+j+1)=m+bias;
end
extend_right=x(N+1:N+num);
extend=;
plot(x,'-b');
hold on;
plot(extend,'*r');
同求啊。。。。 用支持向量机处理端点效应好吗
页:
[1]