频谱图与响应问题
本帖最后由 wdhd 于 2016-9-13 13:31 编辑=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(,sf(1:Nw,i))
title('傅氏频谱图')
%xlim()
xlabel('Frequency (rad/sec)')
ylabel('H(w) (m/N)')
%ylabel('|Fx(w)|')
%xlim()
%legend('1-DOf','2-DOF','3-DOF')
xlim()
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：
要得到一组正确的响应图与可以对应的频谱图该改哪些东西呢? 我
知道我问题很多很杂
请版友帮帮忙
谢谢
页:
[1]