MATLAB程序设计 课件 第7章.ppt

7.5 符号方程的求解 第7章 符号运算简介 7.5.1 符号代数方程的求解 一般代数方程包括线性（Linear）方程、非线性（Nonlinear）方程和超越（Transcedental）方程等，求解函数是solve。当方程不存在符号解而且无其它自由参数时，solve函数将给出数值解。该函数最常用的调用格式为 S=solve(eq,v) 求方程关于指定变量的解 S=solve(‘eq1’,‘eq2’,…,‘eqn’,‘v1’,‘v2’,…,‘vn’) 求方程组关于指 定变量的解 7.5.1 符号代数方程的求解 【说明】 ? 在第1种格式中，eq可以是含等号的符号表达式的方程，也可以是不含等号的符号表达式，但指的仍然是令eq=0的方程；v是表示求解变量名的字符串，当它省略时，求解对findsym指令确认的变量进行。 ? 第2种格式中输入宗量的含义与第1种格式类似。 在第2种格式中，输出宗量S是结构数组。如果要显示求解结果，必须采用S.v1，S.v2，…，S.vn的提取方式。 ? 第2种格式也可以写成： S=solve(‘eq1’,‘eq2’,…,‘eqn’,‘v1,v2,…,vn’) ? 在得不到“封闭型解析解”时，又不存在其它不确定参数，则给出数值解。 7.5.1 符号代数方程的求解 【例7.26】求方程 ， 以及 的解。 编写M脚本文件exam7_26.m如下： % exam7_26 clear syms a b c x s1=solve(a*x^2+b*x+c) % 无等号，默认求解变量为x s2=solve((x+2)^x=2,x) % 有等号，指定求解变量为x s3=solve(sin(x)+cos(2*x)-1,x) % 无等号，指定求解变量为x 7.5.1 符号代数方程的求解 exam7_26.m的运行结果 【说明】 ? 由于三角函数是周期函数，只能得到在x=0附近的有限几个解。 7.5.1 符号代数方程的求解 % exam7_27 clear S=solve(u*y^2+v*z+w=0,y+z+w=0,y,z); disp(S.y=),disp(S.y),disp(S.z=),disp(S.z) 【例7.27】求方程组 关于y，z的解。 编写M脚本文件exam7_27.m如下： exam7_27.m的运行结果 7.5.1 符号代数方程的求解 【说明】 ? 建议：在使用solve函数求解代数方程组时，最好采用指定变量方式。 7.5 符号方程的求解 7.5.2 符号微分方程的求解 求解符号常微分方程的最常用指令格式为 S=dsolve(eq,cond,v) S=dsolve(eq1,eq2,…,eqn, cond1,cond2,…,condn, v) 【说明】 ? 输入宗量包括三部分：微分方程、初始条件、指定自变量。微分方程是必不可少的输入内容；其余视需要而定，可有可无。输入宗量必须以字符串形式编写。 ? 若不对自变量加以专门的定义，则以小写字母t为自变量。 7.5.2 符号微分方程的求解 微分方程的记述规定：当y是“因变量”时，用“Dny”表示“y的n阶导数”。在t为默认自变量时，Dy表示 ；Dny表示 。 ? 初始条件应写成诸如：‘ y(a)=b,Dy(c)=d’的格式等。a，b，c，d可以是变量使用符以外的其它字符。当初始条件少于微分方程数时，在所得解中将出现常数符C1，C2，…。解中任意常数符的数目等于所缺少的初始条件数。 ? 在第2种格式中，输出宗量S是结构数组。如果y是因变量，则关于它的解在S.y中。 7.5.2 符号微分方程的求解 % exam7_28 clear S_general=dsolve(Dy=a*y) S_special=dsolve(Dy=a*y,y(0)=b) 【例7.28】求微分方程 的通解和在初始条件 下的特解。 编写M脚本文件exam7_28.m如下： exam7_28.m的运行结果 7.5.2 符号微分方程的求解 【说明】 ? 在列写微分方程时，最好遵循“导数在前函数在后，导数阶数降阶”的次序，否则有可能运行出