Matlab如何求频响函数

springen

小弟初学Matlab，现在做一个关于模态分析的程序。现在程序基本都写好了，但是不知道如何编写频响函数的程序，希望各位前辈指点一下。

程序如下：

1. %二自由度可解耦系统的振动模态
   
2. clc;
   
3. clear;
   
4. 
   
5. m1=1;
   
6. m2=2;
   
7. k1=2;
   
8. k2=4;
   
9. c1=1;
   
10. c2=2;
    
11. 
    
12. x0=[0;1];    %初始位移
    
13. xd0=[0;-1];  %初始速度
    
14. 
    
15. tf=50;       %时间
    
16. dt=0.1;      %时间间隔
    
17. 
    
18. M=[m1,0;0,m2];
    
19. K=[k1+k2,-k2;-k2,k2];
    
20. C=[c1+c2,-c2;-c2,c2];
    
21. A=[zeros(2,2),eye(2);-M\K,-M\C];
    
22. 
    
23. y0=[x0;xd0];
    
24. 
    
25. for i=1:round(tf/dt)+1       %设定计算点，作循环计算
    
26.     t(i)=dt*(i-1);
    
27.     y(:,i)=expm(A*t(i))*y0;    %循环计算矩阵指数
    
28. end
    
29. 
    
30. figure(2);
    
31. subplot(2,1,1);plot(t,y(1,:));xlabel('Zeit');ylabel('Verschiebung1');grid on
    
32. subplot(2,1,2);plot(t,y(2,:));xlabel('Zeit');ylabel('Verschiebung2');grid on
    
33. 
    
34. [v,w]=vbr_sf(M,K,C)

复制代码

函数程序如下：

1. function [v,w]=vbr_sf(m,d,k)
   
2. %vbr_sf  vbr_sf(m,d,k)
   
3. %        [v,w]=vbr4(m,d,k)
   
4. %        function vbr_sf finds the mode shapes and natural frequencies of
   
5. %        a linear second order matrix equation.  
   
6. %        [v,w]=vbr_sf(m,k) finds the mode shapes and natural frequencies
   
7. %        for the undamped case.
   
8. if nargin==2
   
9.   k=d;
   
10.   [v,w]=eig(m\k);
    
11.   w=sqrt(w);
    
12. end
    
13. if nargin==3
    
14.   if norm(d/m*k-k/m*d) < 1e-8*norm(k/m*d)
    
15.     %disp('Damping is proportional, eigenvectors are real.')
    
16.     [v,w]=eig(m\k);
    
17.     w=sqrt(w);
    
18. 
    
19.    else
    
20.     %disp('Damping is non-proportional, eigenvectors are complex.')
    
21.     a=[0*k eye(length(k));-m\k -m\d];
    
22.     [v,w1]=eig(a);
    
23.     w=abs(w1);
    
24. 
    
25.   end
    
26. end
    
27. w=diag(w);
    

复制代码

wenwuyang

chybeyond

自功率谱和互功率谱的比值

zswseu

记得有个专门的函数，直接把输入和输出相比。。。实测的数据。

danzhilan

顶  看看大神回答

HHT

继续， 讲明白一点.......

sealight

help TFESTIMATE

memo1221

zswseu 发表于 2014-12-21 20:07
记得有个专门的函数，直接把输入和输出相比。。。实测的数据。

tfestimate

mzbys

chybeyond 发表于 2014-12-17 10:01
自功率谱和互功率谱的比值

是互谱比自谱吧

刘长

你的程序是响应函数吧，频响函数如果你要从模态分析的方法做，建议参考曹树谦编著的振动结构模态分析，很详细。希望对你有帮助

jiangwj

刘长 发表于 2016-6-7 15:51
你的程序是响应函数吧，频响函数如果你要从模态分析的方法做，建议参考曹树谦编著的振动结构模态分析，很详 ...

有代码吗  可以分享一下