我的课题 发表于 2013-4-13 10:55

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=;

ChaChing 发表于 2013-5-5 10:55

求助完整格式:出错代码和出错提示

邢永忠 发表于 2013-9-16 19:34

本帖最后由 邢永忠 于 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');

大喜 发表于 2015-4-19 11:47

同求啊。。。。

liu@yanyong 发表于 2015-6-24 19:47

用支持向量机处理端点效应好吗
页: [1]
查看完整版本: emd端点效应处理程序