prohl.m

来自「计算关于转子动力学中」· M 代码 · 共 32 行

function prohl(m,l,K,J,v,ei,N)
%求传递矩阵T
Tn=ones(4);
syms w 
    for n=1:N-1
        T(1,1)=1+l(n)^3*(1-v(n))*(m(n)*w*w-K(n))/(6*ei(n));
        T(1,2)=1+l(n)^2*J(n)*w*w/(2*ei(n));
        T(1,3)=l(n)^2/(2*ei(n));
        T(1,4)=l(n)^3*(1-v(n))/(6*ei(n));
        T(2,1)=l(n)^2*(m(n)*w*w-K(n))/(2*ei(n));
        T(2,2)=1+l(n)*J(n)*w*w/ei(n);
        T(2,3)=l(n)/ei(n);
        T(2,4)=l(n)^2/(2*ei(n));
        T(3,1)=l(n)*(m(n)*w*w-K(n));
        T(3,2)=J(n)*w*w;
        T(3,3)=1;
        T(3,4)=l(n);
        T(4,1)=m(n)*w*w-K(n);
        T(4,2)=0;
        T(4,3)=0;
        T(4,4)=1;
        Tn=T.*Tn;  
    end
    Tn=[1 0 0 0 ;0 1 0 0;0 J(N)*w*w 1 0;m(N)*w*w-K(N) 0 0 1].*Tn;
    r=Tn(3,1)*Tn(4,2)-Tn(3,2)*Tn(4,1);
    w=solve(r)
    wn=size(w)