tianyw 发表于 2009-5-12 12:02

急需帮助 时频分析源程序帮忙分析一下

clc;
clear;
f=;
t=/6144;
x=load('20071031C_单参数导出_VY16s.txt');%exp(j*pi*t.^2/2);%
x=x';
%x=sin(2*pi*500*t);
figure(1);
%plot(t,x);
plot(t,x(1:length(t)));
xlabel('时间(s)');
ylabel('振动码值');
%title('时域信号波形');
tic
y=wdf(x,t,f);
toc
figure(2);
image(t,f,abs(y)*200);
xlabel('时间(s)');
ylabel('频率(Hz)');
%title('WVD频谱分布—WVD谱1');
colormap(gray(256));
% figure(3);
% mesh(t,f,abs(y)*200);
% xlabel('时间(s)');
% ylabel('频率(Hz)');
% title('WVD频谱分布—WVD谱1');

f1=;
tic
y1=wdf1(x,t,f1);
toc;
figure(4);
image(t,f,abs(y1)*200);
xlabel('时间(s)');
ylabel('频率(Hz)');
%title('WVD频谱分布—WVD谱2');
colormap(gray(256));
% tic
% =wdf2(x,t,f,3);
% toc
% f1=;
% tic
% =wdf3(x,t,f1,3);
% toc
% figure(4)
% subplot(1,2,1)
% image(t0,f,abs(y)*100)
% colormap(gray(256))
% subplot(1,2,2)
% image(t0,f,abs(y1)*100)
% colormap(gray(256))
L = length(x);
Fs = 6144;
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(x,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
figure(5);
plot(f,2*abs(Y(1:NFFT/2+1))) ;
%title('飞机机翼振动单边幅度谱');
xlabel('频率 (Hz)');
ylabel('幅值|X(f)|');

我在处理 的时候,总是出现内存溢出“out of memory”提示
有没有人帮忙解决一下

还有这个程序只能处理4000个大小的数据
再大的话就出现上述情况

quit 发表于 2009-5-14 21:06

将其做1024点的fft

xiaotian99 发表于 2009-6-1 08:47

你的数据文件??

你的数据文件??
页: [1]
查看完整版本: 急需帮助 时频分析源程序帮忙分析一下