bzhjjjzh.m

来自「这是一些matlab函数：欧拉公式一、格斯比函数等等 它是一些已知的m文件 可以」· M 代码 · 共 33 行

function Phi=bzhjjjzh(A,n,k,lambda,m)
syms t xi; f=[ ];
for j=1:k 
    for i=1:m(j)
        if imag(lambda(j))==0
        v=t^(i-1)*exp(lambda(j)*t)/prod(1:(i-1));
        elseif imag(lambda(j))>0
        v=t^(i-1)*exp(real(lambda(j))*t)*cos(imag(lambda(j))*t)/prod(1:(i-1));
        else v=t^(i-1)*exp(real(lambda(j))*t)*sin(-imag(lambda(j))*t)/prod(1:(i-1));
        end; 
        f=[f;v];
    end; 
end
Ft=f;f1=f;
for i=1:n-1 
    df=diff(f1);Ft=[Ft df];f1=df; 
end
F0=subs(Ft,t,0)
E=sym(eye(n));
Phi=sym(zeros(n));
for l=1:n;
    G=sym(zeros(n)); 
    G(:,1)=E(:,l);
    for i=2:n; 
        Apower=A^(i-1); 
        G(:,i)=Apower(:,l);
    end
    Ge=G; 
    Re=Ge/F0; 
    Phi(:,l)=(Ge/F0)*f;
end
Phi=simplify(Phi);