3MATLAB的符号运算.ppt

例3. 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1 g1=x+y+z=1,g2=x-y+z=2,g3=2*x-y-z=1 f=solve(g1,g2,g3) f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) f = z = 5/6, y = -1/2, 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 5. 符号微分方程求解 —— 用一个函数可以方便地得到微 分方程的符号解 符号微分方程求解指令：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) = sin(t), y(t) = cos(t) 二阶微分方程 dsolve(D2y=-a^2*y,y(0)=1,Dy(pi/a)=0) ans = cos(a*x) 例3. 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(t)的时间曲线图 求该方程的解 三、maple函数——符号运算的扩展 maple——是专门进行数学运算的软件工具， 具有超强的符号运算能力，提供了 几乎包括所有数学领域的专用函数 matlab——依赖于maple的内核与函数库，扩 展了自己的符号运算功能。 matlab还设计了对maple库函数的调用功能 使得已有的maple数学功能，可以扩充matlab 中,作为自身符号运算能力的扩展。 1. maple内核访问函数 可以访问maple内核的matlab函数: maple ——— 访问maple内核函数 mapleinit —— maple函数初始化 mpa ———— maple函数定义 mhelp ——— maple函数帮助命令 procread —— maple函数程序安装 ?. maple 的调用格式 maple(表达式） —— 将表达式送至maple内核， 返回符号表达式结果。 maple （函数，变量1，变量2） ——调用maple函数，传递给定 变量。 例1. 展开5阶 bernoulli 多项式，计算 x=3 时bernoulli 数。 a=maple(bernoulli(5,x)) a = -1/6*x+5/3*x^3+x^5-5/2*x^4 a=maple(bernoulli(5,3)) a = 85 例2. 化简三角函数式sin2x+cos2x a=maple(simplify(sin(x)^2+cos(x)^2);) a = 1 例4. 求f(t)=e-3tsint的拉式变换 f=maple(laplace(exp(-3*t)*sin(t),t,s);) f = 1/((s+3)^2+1) 例4. 寻找二次多项式的完全平方 f (x) = x2+2x+2 a=maple(completesquare(x^2+2*x+2)) a = completesquare(x^2+2*x+2) 将工具包装入内存 maple(with(s