MATLAB程序设计和教程.pptVIP

第 4 章 MATLAB 符号计算 4.1 符号对象 的创建 4.2 符号 表达式的化简和替换 4.3 符号微积分 4.4 符号方程求解 4.5 简易绘图函数 4.4 符号方程求解 4.4.1 符号代数方程求解 在 MATLAB 中，求解用符号表达式表示的代数方程 可由函数 solve 实现，其调用格式为： solve(s) ：求解符号表达式 s 的代数方程，求解变量 为默认变量。 solve(s,v) ：求解符号表达式 s 的代数方程，求解变量 为 v 。 solve(s1,s2,…,sn,v1,v2,…,vn) ：求解符号表达式 s1,s2,…,sn 组成的代数方程组，求解变量分别 v1,v2,…,vn 。 。 【例 4-21 】分别求解代数方程 ax^2+bx+c 和 cos(2x)+sin(x)=1 。 ? syms a b c x ? s=a*x^2+b*x+c; ? solve(s) ? solve(cos(2*x)+sin(x)=1) ? ans = ? 1/2/a*(-b+(b^2-4*a*c)^(1/2)) ? 1/2/a*(-b-(b^2-4*a*c)^(1/2)) ? ans = ? pi ? 0 ? 1/6*pi ? 5/6*pi 0 50 100 150 200 250 300 350 400 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 0 50 100 150 200 250 300 350 400 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 50 100 150 200 250 300 350 400 -1 -0.5 0 0.5 1 1.5 2 ? 【例 4-22 】求解代数方程组 ? syms x y z ? f=x^2-y^2+z-10; ? g=x+y-5*z; ? h=2*x-4*y+z; ? [x,y,z]=solve(f,g,h) ? S=solve(f,g,h); ? [S.x,S.y,S.z] 4.4.2 符号常微分方程求解 在 MATLAB 中，用大写字母 D 表示导数。例如， Dy 表示 y ， D2y 表示 y ， Dy(0)=5 表示 y(0)=5 。 D3y+D2y+Dy-x+5=0 表 示微分方程 y+y+y-x+5=0 。符号常微分方程求解可以 通过函数 dsolve 来实现，其调用格式为： dsolve(e,c,v) 该函数求解常微分方程 e 在初值条件 c 下的特解。参数 v 描述 方程中的自变量，省略时按缺省原则处理，若没有给出初 值条件 c ，则求方程的通解。 dsolve 在求常微分方程组时的调用格式为： dsolve(e1,e2,…,en,c1,…,cn,v1,…,vn) 该函数求解常微分方程组 e1,…,en 在初值条件 c1,…,cn 下的特 解，若不给出初值条件，则求方程组的通解， v1,…,vn 给 出求解变量。 ? 【例 4-23 】求微分方程的 dy/dt=ay ? 通解和当 y(0)=b 时的特解。 ? dsolve(Dy = a*y) ? dsolve(Dy = a*y, y(0) = b) ? % 【例 4-24 】求微分方程时的特解。 ? dsolve(D2y = -a^2*y) ? dsolve(D2y = -a^2*y, y(0) = 1, Dy(pi/a) = 0) ? ? % 【例 4-25 】绘制函数表达式 x^2-y^4 的二 维图形。 ? syms x y ? ezplot(x^2-y^4) x y x 2 -y 4 = 0 -6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6 ? % 【例 4-26 】绘制误差函数的二维图形。 ? syms x ? ezplot(erf(x)) ? grid -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -1 -0.5 0 0.5 1 x erf(x) ? 【例 4-27 】在极坐标下绘制函数表达式 1+cos(t) 的二维图形。 ? syms t ? ezpolar(1+cos(t)) 0.5 1 1.5 2 30 210 60 240 90 270 120 300 150 330 180 0 r = 1+cos(t) -6 -4 -2 0 2 4