第3次MATLAB符号计算3第三次课.pptVIP

例(2) 求三元非线性方程组 的解。 eq1=sym(x^2+2*x+1); eq2=sym(x+3*z=4); eq3=sym(y*z=-1); [x,y,z]=solve(eq1,eq2,eq3) x = -1 y = -3/5 z = 5/3 3.6.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(eq,c,v) 该函数求解常微分方程eq在初值条件c下的特解。参数v描述方 程中的自变量，省略时按缺省原则处理，若没有给出初值条件 c，则求方程的通解。 dsolve在求常微分方程组时的调用格式为： dsolve(eq1,eq2,…,eqn,c1,…,cn,v1,…,vn) 该函数求解常微分方程组eq1,…,eqn在初值条件c1,…,cn下的特 解，若不给出初值条件，则求方程组的通解，v1,…,vn给出求 解变量。 例(1) 求微分方程 , y(1)=0，y(5)=0的解。 y=dsolve(x*D2y-3*Dy=x^2,x) %求微分方程的通解 y = -1/3*x^3+C1+C2*x^4 y=dsolve(x*D2y-3*Dy=x^2,y(1)=0,y(5)=0,x) %求微分方程的特解 y = -1/3*x^3+125/468+31/468*x^4 例(2) 求微分方程组 的解。 [x,y]=dsolve(Dx=y’,’Dy=-x) x = cos(t)*C1+sin(t)*C2 y = -sin(t)*C1+cos(t)*C2 练习： 1、求解方程组 3、求解微分方程 2、求函数 的极值点。 3.7符号函数的可视化 3.7.1符号函数的绘图命令 1. ezplot和 ezplot3命令 ezplot(F,[xmin,xmax],fig)%画符号表达式的图形 ezplot3(x,y,z,[tmin,tmax],’animate’) 【例3.25】画出上一小节【例3.23】中y(x)特解的图形，如图3.2所示。 y =sym(-1/3*x^3+1/3*x^4) y = -1/3*x^3+1/3*x^4 ezplot(y) ezplot(y,[0,100]) %绘制符号函数y在[0,100]中的图形 图3.2 (a) y(x)在[-2л，2л]的图形 (b) y(x)在[0，100]的图形 【例3.26】用ezplot3绘制三维符号表达式曲线，如图所示。 x=sym(sin(t)); z=sym(t); y=sym(cos(t)); ezplot3(x,y,z,[0,10*pi],animate) %绘制t在[0,10*pi]范围的三维曲线 图3.3 符号表达式绘制三维曲线 2. 其它绘图命令 命令名 含义 举例 ezcontour 画等高线 ezcontour(x*sin(t),[-4,4]) ezcontourf 画带填充颜色等高线 ezcontourf(x*sin(t),[-4,4]) ezmesh 画三维网线图 ezmesh(sin(x)*exp(-t),cos(x)*exp(-t),x,[0,2*pi]) ezmeshc 画带等高线的三维网线图 ezmeshc(sin(x)*t,[-pi,pi]) ezpolar 画极坐标图 ezpolar(sin(t),[0,pi/2]) ezsurf 画三维曲面图 ezsurf(x*sin(t),x*cos(t),t,[0,10*pi]) ezsurfc 画带等高线的三维曲面图 ezsurfc(x*sin(t),x*cos(t),t,[0,pi,0,2*pi]) 3.7.2图形化的符号函数计算器由funtool.m文件生成，在MATLAB命令窗口输入命令“funtool”，就会出现该图形化函数计算器。 图3.4 (a) Figure No. 1窗口 (b) Figure No. 2窗口 图3.4 (c) Figure No. 3窗口 *3.8 MuPAD的使用 __自学P110-P112 21 * g=