arryf.m

内燃机转子仿真

function [m,c,k]=arryf(M,C,kt,n);
m=zeros(n);
c=zeros(n);
k=zeros(n);
for i=1:n
    m(i,i)=M(i);
end
for i=1:n
   k(i,i)=kt(i)+kt(i+1);
    c(i,i)=C(i)+C(i+1);
end
for i=1:n-1
           k(i,i+1)=-kt(i+1);
            c(i,i+1)=-C(i+1);
            k(i+1,i)=-kt(i+1);
            c(i+1,i)=-C(i+1);
end