华北水院第2章MATLAB符号运算.pptVIP

例2. 符号方程cos(x)=sin(x) tan(2*x)=sin(x)求解 f1=solve(cos(x)=sin(x)), f1 = 1/4*pi solve(f , b ) —— 对指定变量b求解 ans = -(a*x^2+c)/x f2= [ 0] [ pi] [ atan(1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))] [ atan(-1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))] [ atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))+pi] [ -atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))-pi] f2=solve(tan(2*x)=sin(x)) 例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 中,作为自身符号运算能力的扩展。 maple内核访问函数 可以访问maple内核的matlab函数: maple ——— 访问maple内核函数 mapleinit —— maple函数初始化 mpa ———— maple函数定义 mhelp ——— maple函数帮助命令 procread —— maple函数程序安装 maple 的调用格式 maple(表达式） —— 将表达式送至maple内核， 返回符号表达式结果。 maple （函数，变量1，变量2） ——调用maple函数，传递给定变量。 例1. 展开5阶 bernoulli 多项式