MATLAB求解复系数特征方程的实数根的方法

        已知一个复系数的特征方程，我们规定其特征根为实数，那么求根的过程其实就是将该方程展开，分别列出其实部方程和虚部方程，然后依次求解。这里需要用到的MATLAB函数有：

expand（扩展表达式并简化函数输入）、 sym2poly（提取多项式系数）、 solve（求解方程）

        代表性问题：针对二元机翼的动气弹问题。

        运用Theodorsen方法求解二元机翼的颤振速度和颤振频率，在求解过程中我们会遇到一个复杂的行列式求解的问题，行列式如下：

        其中M和L是气动力矩和气动力，注意是复数，可以查表得到，而X是我们要求解的根，且前提条件是它是实数，这就需要将这个特征方程展开，并且拆分为实部方程和虚部方程分别求解，手算太过于复杂，我们想到了用MATLAB代替铅笔进行数值计算。

        查表得Lh, La, Mh, Ma的值分别为 0.5407-1.9293i，-3.1749-2.8305i，0.5，0.375-1.6667i，求解代码如下：

clc;
clear;
Lh=0.5407-1.9293i; La=-3.1749-2.8305i; Mh=0.5; Ma=0.375-1.6667i;  %气动力参数值
syms x xr xi %设置符号变量x及其实部和虚部
p=(5*(1-0.25*x)+Lh)*(5*0.5*(1-x)+(Ma-0.4*(La+Mh)+0.16*Lh))-(5*0.25+(Mh-0.4*Lh))*(5*0.25+(La-0.4*Lh)); %用f表示复系数行列式
ep=expand(p); %将行列式展开为标准多项式形式
xishu=sym2poly(ep); %提取行列式系数
shi=real(xishu); %提取实部
xu=imag(xishu); %提取虚部
eqn1= shi(1)*xr^2+shi(2)*xr+shi(3)==0; %列出实部方程
eqn2= xu(1)*xi+xu(2)*xi+xu(3)==0; %列出虚部方程
%我们知道列出的方程最高次数是2
s1=double(solve(eqn1,xr)); %求解实部方程对应的根
s2=double(solve(eqn2,xi)); %求解虚部方程对应的根
%solve求解得到的不是最终的数字，需要再用double简化一下
s1
s2

 运行后得到的结果为：

 由于虚部方程二次项系数为0，所以只得到一个根。

注意在用solve求解后需要再使用double命令化简一下结果，不然得到的结果是这样的：

        对于更一般的复数方程求根问题，若方程的根是一个复数，那么求解过程就更简单了，列出方程后我们直接使用solve命令求解，而不需要再提取实部虚部系数分别列方程求解。