齐次线性微分方程组的特解怎么求
matlab求齐次线性微分方程组的基解矩阵:某过程求解释 A=[2,0,0,0,0;4,2,4,0,-4;4,-1,6,1,-5;12,-2,12,4,-14;8,-1,8,1,-7]下面这句求矩阵A的特征值和特征向量,以特征值构成对角阵lambda[v,lambda]=eig(sym(A));n本质上是系统的阶次n=length(lambda);这个判断有点莫名其妙,我想不出什么情况下if的条件会不满足if length(v(1,:))=nsyms t%定义符号变量tPhi=[];置初始的Φ(t)矩阵为空,后面依次添加各列for j=1:n通过循环依次添加Φ(t)矩阵各列,函数eval纯属多余Phi=[Phi,eval(v(:,j)*exp(lambda(j,j)*t))];end下面这一句又是莫名其妙,如果想显示变量Phi,直接写变量名即可,没必要给自身赋值Phi=Phi通过处以Φ(0),对Φ(t)进行类似归一化的操作。eval同样是多余的Psi=eval(Phi/subs(Phi,t,'0'))else这个分支根本似乎不可能执行到disp('False');end其实求状态转移矩阵更简单的做法如下:[P,D]=eig(sym(A));Psi1=P*expm(D*t)*P^-1可以验证,这样求出的基解矩阵和上面相同Psi-Psi1
已知齐次线性微分方程的基本解组x1=t,x2=e∧t求非齐次线性方程组x''+t/(1-t)x 下面是解题方法,并附加该方微分程齐次的解法.以上使解题过程,该方微分程齐次的解法如下