|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 wdhd 于 2016-9-13 13:31 编辑
[num dof]=size(disp);
for i=1:dof;
f(:,i)=disp(:,i);
dt=0.01;
w=pi/dt;
Nt=num;
Nw=Nt/2;
dw=w/Nw;
%ff=fft(f)*dt;
T=Nt*dt;
sf(:,i)=1/(2*pi*T)*(abs(fft(f(:,i)))*dt).^2;
%sf1=sum(sf,2);
%cc=ifft(sqrt(sf*2*pi*T)/dt);
figure(i)
plot([0:dw:(Nw-1)*dw],sf(1:Nw,i))
title('傅氏频谱图')
%xlim([dw Nw*dw])
xlabel('Frequency (rad/sec)')
ylabel('H(w) (m/N)')
%ylabel('|Fx(w)|')
%xlim([0*dw 10000*dw])
%legend('1-DOf','2-DOF','3-DOF')
xlim([0 400])
end
请问以上这些要怎么解释?? 例如 sf(:,i)=1/(2*pi*T)*(abs(fft(f(:,i)))*dt).^2
这是甚么意思?为什么有2*pi*T? 还有 dt 跟 dw 有甚么不同跟影响 ? 拜托请帮我! 还有我用这两组档案下去跑响应与频谱图
用newmark法的dt设定有限制吗? 我发现Q1:newmark的dt只有带0.01才是比较像正确的响应,但频谱图只能对应到w1。 Q2:两边 dt 同时带0.0001时,响应图是错的,但可以对应到完整的频谱图。 Q3:
要得到一组正确的响应图与可以对应的频谱图该改哪些东西呢? 我
知道我问题很多很杂
请版友帮帮忙
谢谢
|
|