Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
 
 
主题工具 显示模式
旧 2010-06-27, 16:03   #1
穿球鞋的猫
初级会员
 
注册日期: 2010-06-27
帖子: 1
声望力: 0
穿球鞋的猫 正向着好的方向发展
默认 求助 路过帮忙看一下

已经有100万个点,要求FFT,

A1=load('D:\。。。\sc1004.dat');
A1=A1-mean(A1);

A2=load('D:\。。。\sc2004.dat');
A2=A2-mean(A2);

A3=load('D:\。。。\SC1005.dat');
A3=A3-mean(A3);

A4=load('D:\。。。\SC2005.dat');
A4=A4-mean(A4);



B1=abs(fft(A1(99900:200000),1000000));
这里计算A1的幅值,但是没有明白fft((99900:2000000),1000000)啥意思,下面的又不要,而且说是取8192个点 ———


B2=abs(fft(A2,1000000));
B3=abs(fft(A3,1000000));
B4=abs(fft(A4,1000000));




%t=0:0.01:10000.03;
%z=(t');

%plot(z,A1)




%B1=B1(1:500000);
%B2=B2(1:500000);
%B3=B3(1:125000);
%B4=B4(1:125000);

U1=B2./B1;
U2=B4./B3;

F=0100000000)/999999100000000);
f=(F');

U1=U1/max(U1);
U2=U2/max(U2);

%Plot(f,U1,'r',f,U2,'k')
%semilogx(f,U1,'r',f,U2,'k')
%loglog(F,UF1,'r',F,U2,'k')

B1=B1/max(B1);
B2=B2/max(B2);
B3=B3/max(B3);
B4=B4/max(B4);

%plot (f, B1,'b',f,B2,'r'),grid
figure (3)
subplot (2,1,1);
loglog(f,B1,'b',f,B2,'r'),grid
title('')
xlabel ('')
ylabel ('')

subplot (2,1,2);
loglog(f,U1,'g'),grid
title ('')
xlabel ('')
ylabel ('')

figure (4)
subplot (2,1,1);
loglog(f,B3,'b',f,B4,'r'),grid
title('')
xlabel ('')
ylabel ('')

subplot (2,1,2);
loglog(f,U2,'g'),grid
xlabel ('')
ylabel ('')
title ('')
穿球鞋的猫 当前离线   回复时引用此帖
 


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 10:58


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.