定稿第4章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 f3= matlab4.2的解 [ 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)) f2 = matlab4.2的解 [ 0] [acos(1/2+1/2*3^(1/2))] [acos(1/2 -1/2*3^(1/2))] numeric(f3) ans = 0 3.1416 0 + 0.8314i 0 - 0.8314i 1.9455 -1.9455 numeric(f2) ans = 0 0 + 0.8314i 1.9455 matlab4.2与6.1的对比 例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——依