第3章 MALB符号计算.pptVIP

【例3.23】 求方程ax2+bx+c=0和sinx=0的解。 f1=sym(a*x^2+b*x+c) %无等号 f1 = a*x^2+b*x+c solve(f1) %求方程的解x ans = -(b + (b^2 - 4*a*c)^(1/2))/(2*a) -(b - (b^2 - 4*a*c)^(1/2))/(2*a) f2=sym(sin(x)) f2 = sin(x) solve(f2,x) ans = 0 程序分析：当sinx=0有多个解时，只能得出0附近的有限几个解。 【例3.24】 求三元非线性方程组 的解。 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、y、z x = ?1 y = ?3/5 z = 5/3 程序分析：输出结果为“结构对象”，如果最后一句为“S=solve(eq1,eq2,eq3) ”，则结果如下。 S = x: [1x1 sym] y: [1x1 sym] z: [1x1 sym] MATLAB提供了dsolve命令，可以用于对符号常微分方程进行求解。 语法： dsolve(eq, con, v)