Increasing accuracy of eigen/eig function

Question

Osama Anwar 2021-1-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/729533-increasing-accuracy-of-eigen-eig-function

评论： Osama Anwar 2021-1-28

采纳的回答： Bjorn Gustavsson

phi.png

在 MATLAB Online 中打开

How can I increase decimal places of the result obtained by using eig function

Here is my code I need it for phi

clear
clc
syms t
m=45000;
k=57*10^6;
h=3;
n=5;
M=zeros(n,n);
K=zeros(n,n);
Phi=zeros(n,n);
H=zeros(n,1);
w=zeros(n,1);
Q=zeros(n,1);
% applicable if all stories' m and k are same
for i=1:n
    
    for j=1:n
        
        if i==j
            
            M(i,j)=m;
            K(i,j)=2*k;
            
            if and(i==n,j==n)
                
                K(i,j)=k;
                 
            end
            
        end
        
        if or(i==j+1,j==i+1)
                
                K(i,j)=-k;
                 
        end
            
    end
    
    H(i,1)=i*h;
    
end
M
K
[phi,d]=eig(K,M);
phi