第6章节 Matlab符号计算.pptVIP

第6章 MATLAB符号计算 例 f = ax2+bx+c = 0 求解 f=a*x^2+b*x+c; solve(f) —— 对缺省变量x求解 ans = [1/2/a*(-b+(b^2-4*a*c)^(1/2))] [1/2/a*(-b-(b^2-4*a*c)^(1/2))] solve(f , b ) —— 对指定变量b求解 ans =(a*x^2+c)/x 例 符号方程cos(x)=sin(x)， tan(2*x)=sin(x)求解 f1=solve(cos(x)=sin(x)), f1 =1/4*pi f2=solve(tan(2*x)=sin(x)) f2 = [ 0] %f2(1) [ pi] %f2(2) [ atan(1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))] %f2(3) [ atan(-1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))] %f2(4) [ atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))+pi] %f2(5) [ -atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))-pi] %f2(6) 例 解方程组 g1=x+y+z=1; g2=xy+z=2; g3=2*x-y-z=1; f=solve(g1,g2,g3) f.z = 5/6, f.y = -1/2, f.x = 2/3 y=dsolve(D2y+2*Dy+2*y=0,y(0)=1,Dy(0)=0) y =exp(-t)*cos(t)+exp(-t)*sin(t) ezplot(y) —— 方程解y(t)的时间曲线图 例 求微分方程 ， y(1)=0，y(5)=0的解。 y=dsolve(x*D2y-3*Dy=x^2,x) %求微分方程的通解 y =-1/3*x^3+C1+C2*x^4 y=dsolve(x*D2y-3*Dy=x^2, y(1)=0,y(5)=0, x) %求微分方程的特解 y =-1/3*x^3+125/468+31/468*x^4 例 求微分方程组 的解。 [x,y]=dsolve(Dx=y,Dy=-x) x =cos(t)*C1+sin(t)*C2 y =-sin(t)*C1+cos(t)*C2 程序分析：默认的自由变量是t，C1、C2为任意常数，程序也可指定自由变量，结果相同： [x,y]=dsolve(Dx=y,Dy=-x,t) 例 求线性方程组AX=b的解。 解方程组(1)的命令如下： A=[34, 8, 4; 3, 34, 3; 3, 6, 8]; b=[4; 6; 2]; X=linsolve(A,b) %调用linsolve函数求(1)的解 A\b %用另一种方法求(1)的解 解方程组(2)的命令如下： syms a11 a12 a13 a21 a22 a23 a31 a32 a33 b1 b2 b3; A=[a11,a12,a13;a21,a22,a23;a31,a32,a33]; b=[b1;b2;b3]; X=linsolve(A,b) %调用linsolve函数求(2)的解 XX=A\b %用左除运算求(2)的解 6.5 代数方程的符号求解 6.5 代数方程的符号求解 6.5 代数方程的符号求解 6.5 代数方程的符号求解 常微分方程符号求解函数dsolve，其调用格式为： dsolve(eqn1,condition,var) 求解微分方程eqn1在初值条件condition下的特解，若没有给出初值条件，求方程的通解 dsolve(eqn1,eqn2,…,eqnN,condition1,…,conditionN,var1,…,varN) 求解微分方程组eqn1、…、eqnN在初值条件conditoion1、…、conditionN下的解，若没有给出初值条件，求方程组的通解 6.6 常微分方