紧急求助：HHT谱分析问题

fgsunny

高手帮帮忙，我是新手，
   通过在论坛上的学习，应用论坛中的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有些地方有些迷糊

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

hering

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

fgsunny

我的程序：我的数据太大，上传不了
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

fgsunny

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

fgsunny

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

fgsunny

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))

fgsunny

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

fgsunny

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

wuqiong

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

fgsunny

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

谢谢了:@) :handshake