声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1733|回复: 9

[HHT] 紧急求助:HHT谱分析问题

[复制链接]
发表于 2008-8-27 19:13 | 显示全部楼层 |阅读模式

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

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

x
高手帮帮忙,我是新手,
   通过在论坛上的学习,应用论坛中的HHT程序,其分析的是发动机转子实验测得的一个数据序列,能够得出EDM结构,但是好像不对,运行到HHT谱分析是就运出现如下错误:
??? Error using ==> accumarray
Out of memory. Type HELP MEMORY for your options.
Error in ==> toimage at 98
im = accumarray([indf(:),indt(:)],A(:),[sply,splx]);
Error in ==> Untitled at 45
[im,tt]=toimage(A,f);
不知怎么改,对toimage有些地方有些迷糊

急求高手帮助,导师急着要结果,我已经研究两周了,唉:@(
回复
分享到:

使用道具 举报

发表于 2008-8-27 19:27 | 显示全部楼层
把程序和数据放上来,让大家看看
 楼主| 发表于 2008-8-27 19:47 | 显示全部楼层

哦,好的,都急忘了

我的程序:我的数据太大,上传不了
clear
clc
%----------------------------------------
%(1) read data
%
%----------------------------------------
% case 1: x1(2,:), non-stationary data
% case 2: x2(3,:), non-stationary data
%  data format:   
        fid1=fopen ('E:\Study Sources\important\文章\time_del1_run01.dat','r');
           x=fscanf(fid1,'%e  ',[3 inf]);
        fclose(fid1);
        Ns=length(x);
        x1=x(2,:);
        x2=x(3,:);
        fs=1000;  % Sampling frequency
        t=0:(1/fs):((Ns-1)/fs);
         freq=[1:Ns/2]*fs/Ns;      % frequency axis
     
   figure (1)
        subplot(2,1,1);plot(t,x1);
               ylabel ('x1 /um ','fontsize',14,'fontweight','bold')
        subplot(2,1,2);plot(t,x2);
               ylabel ('x2 /um ','fontsize',14,'fontweight','bold')
               xlabel ('Time/s ','fontsize',14,'fontweight','bold')
   pause
%----------------------------------------  


   
   
   
   
   
   
   
   imf=emd(x1);
emd_visu(x1,1:Ns,imf)
for i=1:nIMF
[A,f,tt]=hhspectrum(c(:,i));
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))
    end
 楼主| 发表于 2008-8-27 19:52 | 显示全部楼层
我看到论坛里的分析大多针对信号解析式,而对于实际信号的很少,其实在故障诊断中数据是很庞大的,我用了几种方法都算不出来matlab就一直busy没有结果,我的采样数据有8千多行。。。。:@(
用法国人的程序,和nspab和mspc都会出现不同的错误,就是出不来结果,而且对EMD的分解成分也不准确、
 楼主| 发表于 2008-8-27 19:54 | 显示全部楼层

高手求助

高手求助:能帮我解释一下我的错误是什么意思吗?大概问题出在哪?是我的实际数据参数过大吗
 楼主| 发表于 2008-8-27 20:05 | 显示全部楼层

不好意思,程序帖错了

clear
clc
%----------------------------------------
%(1) read data
%
%----------------------------------------
% case 1: x1(2,:), non-stationary data
% case 2: x2(3,:), non-stationary data
%  data format:   
        fid1=fopen ('E:\Study Sources\important\文章\time_del1_run01.dat','r');
           x=fscanf(fid1,'%e  ',[3 inf]);
        fclose(fid1);
        Ns=length(x);
        x1=x(2,:);
        x2=x(3,:);
        fs=1000;  % Sampling frequency
        t=0:(1/fs):((Ns-1)/fs);
         freq=[1:Ns/2]*fs/Ns;      % frequency axis
     
   figure (1)
        subplot(2,1,1);plot(t,x1);
               ylabel ('x1 /um ','fontsize',14,'fontweight','bold')
        subplot(2,1,2);plot(t,x2);
               ylabel ('x2 /um ','fontsize',14,'fontweight','bold')
               xlabel ('Time/s ','fontsize',14,'fontweight','bold')
   pause
  %----------------------------------------

imf=emd(x1);
emd_visu(x1,1:Ns,imf)

[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))
 楼主| 发表于 2008-8-27 20:16 | 显示全部楼层

急求:好奇的小议EDM及附件

在学习论坛知识的同时,在论坛中看到许多帖提到小议EDM及附件,可是找了好长时间都没有找到,谁有这个,能给我提供一份吗 ?
万分感谢,再急也要一点点来,啊,真难啊,还是我至今还没找到门道
主要画的HHT谱,三维的,和边际谱都出不来,:@( :'(
 楼主| 发表于 2008-8-27 20:17 | 显示全部楼层

辅助:我的邮箱

fgsunny@163.com,如果有的话发我邮箱就行,太感谢了!
发表于 2008-8-28 13:43 | 显示全部楼层
画谱图的时候有个分辨率的问题,一般默认的较高,但是对于数据量这么大的谱图而言,没有那个必要,你调整下你的输出谱图分辨率就可以了。
 楼主| 发表于 2008-8-28 16:33 | 显示全部楼层
哦,好的,我试试,谢谢:@) ,我还真没注意这个,因为也是新学有些程序上还有些不清晰,分辨率,好的,我看看,

谢谢了:@) :handshake
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-26 17:00 , Processed in 0.063670 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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