5符号计算程序.ppt

* * * * * * * * 由线性代数我们知道A非奇异时,A的行列式不为0,此时方程的解是唯一的。 在实际应用中,除法解方程的速度要比求逆法快2.5倍精确度更高,明显优于求逆法,所以推荐尽量使用除运算,少用逆运算. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 1. 一般符号代数方程的求解 求解函数是solve，调用格式为： solve(eqn1,eqn2,…,eqnN,var1,var2,…,varN) 其中eqn1,eqn2,eqnN等为符号表达式或不带符号的字符串；var1,var2,…,varN 为自变量 * 例 解方程。 命令如下： x1=solve(1/(x+2)+4*x/(x^2-4)=1+2/(x-2),x) %解方程(1) f=sym(x-(x^3-4*x-7)^(1/3)=1); x2=solve(f) %解方程(2) x3=solve(2*sin(3*x-pi/4)=1) %解方程(3) x4=solve(x+x*exp(x)-10,x) %解方程(4)。仅标出方程的左端 x1 =1 x2 =3 x3 =5/36*pi x4 =1.633506170155846384193165178978923428636 * 例 求方程组 的解。 命令如下： [x y]=solve(1/x^3+1/y^3=28,1/x+1/y=4,x,y) x = 1 1/3 y = 1/3 1 * 2. 符号线性方程组的求解 求解线性方程组的函数是linsolve，调用格式为： linsolve(A,B) 它的结果与sym(A)/sym(B)相同 * 3. 符号非线性方程组的求解 求解函数是fsolve，调用格式为： fsolve(f,X0) 以X0为初始矩阵来求解方程f 例：求解下列方程组 3x2-y2=0 3xy2-x2-1=0 * 解微分方程的基本操作命令见下表 在求解微分方程中，用Dy表示 ，D2y表示 ． 输 入 格 式 含　　　义 dsolve(Dy=f(x,y),x) 求一阶微分方程 的通解 dsolve(Dy=f(x,y),y(0)=a,x) 求一阶微分方程 的特解 dsolve(D2y=f(x,y,Dy),y(0)=a,Dy(0)=b,x) 求二阶微分方程 的特解 4. 常微分方程的符号求解 * 例 求以下微分方程和微分方程组： (1)求 的通解； (2)求 的通解； (3)求 满足初始条件y(0)=1，Dy(0)=0 的特解。 * 解 输入及结果如下： syms x y y=dsolve(Dy=x,x) y =1/2*x^2+C1 dsolve (‘D2y= Dy+1’,’x’) ans =-x+C1+C2*exp(x) dsolve(D2y= Dy+1,y(0)=1,Dy(0)=0,x) ans =-x+exp(x) (1) (2) (3) * 1. 级数的符号求和 级数符号求和函数symsum，调用格式为： symsum(s,v,a,b) 符号表达式S中的变量v从a到b的有限和 四、符号级数运算 * 例 求级数 之和。 命令如下： n=sym(n); s1=symsum(1/n^2,n,1,inf) %求s1 s1 = 1/6*pi^2 * 函数的泰勒级数 MATLAB提供了taylor函数将函数展开为幂级数，其调用格式为： taylor(f,v,n,a) 该函数将函数f按变量v展开为泰勒级数，展开到第n项(即变量v的n-1次幂)为止，n的缺省值为6。v的缺省值与diff函数相同。参数a指定将函数f在自变量v=a处展开，a的缺省值是0。 * 例 求函数在指定点的泰勒展开式。 命令如下： x=sym(x); f2=sqrt(1-2*x+x^3)-(1-3*x+x^2)^(1/3); taylor(f2,6) %求泰勒展开式 ans = 1/6*x^2+x^3+119/72*x^4+239/72*x^5 * 1. 函数曲线图? ezplot(F, [a,b]) 作函数F