马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clc;
clear;
f=[1:1536];
t=[1:50000]/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=[0:1536];
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
% [y,t0]=wdf2(x,t,f,3);
% toc
% f1=[0:1024];
% tic
% [y1,t0]=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个大小的数据
再大的话就出现上述情况 |