高建军《matlab程序设计》第三章matlab的符号运算2015.10.pptVIP

例3. 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1 syms x y z; g1=x+y+z-1,g2=x-y+z-2,g3=2*x-y-z-1 f=solve(g1,g2,g3,’x’,’y’,’z’) f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1‘x’,’y’,’z’) f.z = 5/6, f.y = -1/2, f.x = 2/3 f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) f = x: [1x1 sym] f.x ans =2/3 y: [1x1 sym] f.y ans =-1/2 z: [1x1 sym] f.z ans =5/6 [x,y,z]=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) x = 2/3 y =-1/2 z =5/6 6. 符号微分方程求解 —用一个函数可以方便地得到微 分方程的符号解 符号微分方程求解指令：dsolve 命令格式：dsolve(f,g) f —— 微分方程，可多至12个微分方程的求 解；g为初始条件 默认自变量为 x,可任意指定自变量t, u等 微分方程的各阶导数项以大写字母D表示 或 或 或 y的一阶导数—— Dy y的二阶导数—— D2y y的 n 阶导数—— Dny [y1,y2…]=dsolve(x1,x2,…xn) —— 返回 微分方程的解 一阶微分方程 dsolve(Dx=y,Dy=x,x(0)=0,y(0)=1) ans = x(t) = (exp(t)-exp(-t))/2, y(t) = (exp(-t)+exp(t))/2 二阶微分方程 dsolve(D2y=-a^2*y,y(0)=1,Dy(pi/a)=0) ans = cos(a*x) 微分方程组求解 注：eq1,eq2,...为微分方程或微分方程组，cond1,cond2,...,是初始条件或边界条件，v是独立变量，默认的独立变量是t。 函数dsolve用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解。 dsolve(eq1,eq2,...,?cond1,cond2,...,?v) [X,Y]=dsolve(Dx+2*x-Dy=10*cos(t),Dx+Dy+2*y=4*exp(-*t),x(0)=2,y(0)=0,t) 例1：求 的通解。 例2. y=dsolve(D2y+2*Dy+2*y=0,y(0)=1,Dy(0)=0) ans = exp(-x)*cos(x)+exp(-x)*sin(x) ezplot(y) ——方程解y(x)的时间曲线图 求该方程的解 注：ezplot()可以在图型窗口绘制出符号函数图形，也可以画隐函数图形，形如f(x,y)=0这种不能写出像y=f(x)这种函数的图形. 小 结 本节介绍了matlab语言的符号运算 功能，通过学习应该掌握： 掌握如何创建、修改符号矩阵 掌握符号运算功能 * 第三章 MATLAB的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 一、符号运算的基本操作 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 特点： ? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解 Symbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算的。 maple软件——主要功能是符号运算，它占据符号软件的主导地位。