声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2313|回复: 5

[HHT] 那位师兄有SVR回归支持向量机预测EMD处理端点效应的程序

[复制链接]
发表于 2008-4-7 02:04 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
那位师兄有SVR回归支持向量机预测EMD处理端点效应的程序?

谢谢,大家!
回复
分享到:

使用道具 举报

发表于 2008-4-14 21:24 | 显示全部楼层
以前编的一个程序,今天拿出来给大家看看,请大家指教。不过挺不喜欢这种求程序的做法,应该自己去想想才对,其实挺简单的,
%支持向量回归机在处理HHT端点效应中的应用(只给出SVR向右预测部分)
%分以下几个步骤进行:
%1、将原信号向两端延拓;2、将延拓后的信号EMD分解
%3、将分解后的各IMF分量舍去两端超出原信号长度部分
%4、再将3中的各IMF延拓,并求其Hilbert谱,结果中再舍去延拓出的那部分就OK了
clear;
clc;
N=1024;         %采样点数
fs=5120;        %采样频率
dt=1/fs;        %采样时间间隔
t=(0:N-1)*dt; %产生时间序列
s=(1+sin(2*pi*10*t)).*sin(2*pi*80*t)+sin(2*pi*150*t);%产生信号
num=100;%设置要延拓的点数
l=100; %设置训练样本数

fprintf('向右延拓中……\n');
for i=1:l
    X(i,1:N-l)=s(i:N-l+i-1);
    Y(i,1)=s(N-l+i);
end
[nsv beta bias] = svr(X,Y,'linear',Inf,'einsensitive',0);
   
for j=0:num-1
    m=0;
    XX=s(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=s(N+1:N+num);

extend=[s,extend_right];

[ 本帖最后由 wangkun 于 2008-4-14 21:43 编辑 ]
 楼主| 发表于 2008-4-14 21:29 | 显示全部楼层
谢谢您,我知道了以后要多动脑,真的很感谢你!:@)
发表于 2010-5-17 21:58 | 显示全部楼层
谢谢,2楼,我试过,但是效果不太好
发表于 2010-7-26 17:21 | 显示全部楼层

回复

还有人做这个吗
发表于 2013-8-18 17:25 | 显示全部楼层

我运行有错啊,什么情况?向右延拓中……
??? Attempt to execute SCRIPT svr as a function:
E:\结课作业\机械故障诊断学\duandianchuli\svr.m

Error in ==> svr at 16
[nsv beta bias] =
svr(X,Y,'linear',Inf,'einsensitive',0);
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-16 01:23 , Processed in 0.065943 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表