fgsunny 发表于 2008-8-27 19:13

紧急求助:HHT谱分析问题

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

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

hering 发表于 2008-8-27 19:27

把程序和数据放上来,让大家看看

fgsunny 发表于 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',);
      fclose(fid1);
      Ns=length(x);
      x1=x(2,:);
      x2=x(3,:);
      fs=1000;% Sampling frequency
      t=0:(1/fs):((Ns-1)/fs);
         freq=*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
=hhspectrum(c(:,i));
=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))
    end

fgsunny 发表于 2008-8-27 19:52

我看到论坛里的分析大多针对信号解析式,而对于实际信号的很少,其实在故障诊断中数据是很庞大的,我用了几种方法都算不出来matlab就一直busy没有结果,我的采样数据有8千多行。。。。:@(
用法国人的程序,和nspab和mspc都会出现不同的错误,就是出不来结果,而且对EMD的分解成分也不准确、

fgsunny 发表于 2008-8-27 19:54

高手求助

高手求助:能帮我解释一下我的错误是什么意思吗?大概问题出在哪?是我的实际数据参数过大吗

fgsunny 发表于 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',);
      fclose(fid1);
      Ns=length(x);
      x1=x(2,:);
      x2=x(3,:);
      fs=1000;% Sampling frequency
      t=0:(1/fs):((Ns-1)/fs);
         freq=*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)

=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))

fgsunny 发表于 2008-8-27 20:16

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

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

fgsunny 发表于 2008-8-27 20:17

辅助:我的邮箱

fgsunny@163.com,如果有的话发我邮箱就行,太感谢了!

wuqiong 发表于 2008-8-28 13:43

画谱图的时候有个分辨率的问题,一般默认的较高,但是对于数据量这么大的谱图而言,没有那个必要,你调整下你的输出谱图分辨率就可以了。

fgsunny 发表于 2008-8-28 16:33

哦,好的,我试试,谢谢:@) ,我还真没注意这个,因为也是新学有些程序上还有些不清晰,分辨率,好的,我看看,

谢谢了:@) :handshake
页: [1]
查看完整版本: 紧急求助:HHT谱分析问题