紧急求助：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有些地方有些迷糊
急求高手帮助，导师急着要结果，我已经研究两周了，唉:@( 把程序和数据放上来，让大家看看
哦，好的，都急忘了
我的程序：我的数据太大，上传不了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 我看到论坛里的分析大多针对信号解析式，而对于实际信号的很少，其实在故障诊断中数据是很庞大的，我用了几种方法都算不出来matlab就一直busy没有结果，我的采样数据有8千多行。。。。:@(
用法国人的程序，和nspab和mspc都会出现不同的错误，就是出不来结果，而且对EMD的分解成分也不准确、
高手求助
高手求助：能帮我解释一下我的错误是什么意思吗？大概问题出在哪？是我的实际数据参数过大吗不好意思，程序帖错了
clearclc
%----------------------------------------
%(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))
急求：好奇的小议EDM及附件
在学习论坛知识的同时，在论坛中看到许多帖提到小议EDM及附件，可是找了好长时间都没有找到，谁有这个，能给我提供一份吗 ？万分感谢，再急也要一点点来，啊，真难啊，还是我至今还没找到门道
主要画的HHT谱，三维的，和边际谱都出不来，:@( :'(
辅助：我的邮箱
fgsunny@163.com，如果有的话发我邮箱就行，太感谢了！ 画谱图的时候有个分辨率的问题，一般默认的较高，但是对于数据量这么大的谱图而言，没有那个必要，你调整下你的输出谱图分辨率就可以了。 哦，好的，我试试，谢谢:@) ，我还真没注意这个，因为也是新学有些程序上还有些不清晰，分辨率，好的，我看看，谢谢了:@) :handshake
页:
[1]