关于MATLAB中求固有频率以及振型的问题

springen

在论坛里看到前辈给出的求固有频率以及振型的程序，但是运行不成功，不知道为什么。我才开始学习matlab，所以请各位前辈指点一下，多谢。程序如下：
输入M，D，K

1. function [v,w,zeta]=vbr_sf(m,d,k)
   
2. %vbr_sf  vbr_sf(m,d,k)
   
3. %        [v,w,zeta]=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.     zeta=(v'*m*v)\(v'*d*v)/2/w;
    
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.     zeta=-real(w1)/w;
    
25.   end
    
26. end
    
27. w=diag(w);zeta=diag(zeta);
    

复制代码

ChaChing

求助完整格式：出错代码和出错提示

springen

ChaChing 发表于 2014-11-30 20:17
求助完整格式：出错代码和出错提示

多谢，我已经解决了

springen

ChaChing 发表于 2014-11-30 20:17
求助完整格式：出错代码和出错提示

前辈您好，我还有一个问题想请教您。我在论坛里问了一个问题，您能不能帮我看看如何解决呢http://forum.vibunion.com/thread-134903-1-1.html       十分感谢！