数学建模第五节符号计算资料.ppt

【例】：求定积分 x = sym(x); int(abs(1-x),1,2) ans = 1/2 级数求和 symsum 函数用于求无穷级数的和，格式为： symsum(S,v,a,b) S 表示级数通项，是一个符号表达式 v 是求和变量 a 和 b 是求和的首项和末项。 【例】：求下列级数的和 (1) (2) syms n s1 = symsum(1/n^2,n,1,inf) s1 = 1/6*pi^2 s2 = symsum((-1)^(n+1)/n,1,inf) s2 = log(2) Taylor 级数 taylor 函数将函数展开为幂级数，格式为： taylor(f,v,n,a) 该函数将 f 按变量 v 展开为 Taylor 级数，展开到第 n 项，n 默认为 6。参数 a 指定在变量 v = a 处展 开，a 默认为 0。 【例】：求函数的 Taylor 级数展开式 (1) 求 的 5 阶 Taylor 级数展开 (2) 将 在 x = 1 处按 5 次多项式展开 syms x f1 = sqrt(1-2*x+x^3)-(1-3*x+x^2)^(1/3); f2 = (1+x+x^2)/(1-x+x^2); taylor(f1,x,5) ans = 1/6*x^2+x^3+119/72*x^4 taylor(f2,6,1) ans = 3-2*(x-1)^2+2*(x-1)^3-2*(x-1)^5 四、符号方程求解 代数方程是指未涉及微积分运算的方程， 可用 solve 函数求解。调用格式为： solve(eq1,eq2,...,eqN,‘v1,‘v2,...‘vN) 【例】：解下列方程 (1) (2) x = solve(x-(x^3-4*x-7)^(1/3)=1) x = 3 x = solve(2*sin(3*x-pi/4)-1) ％仅标出方程左端 x = 5/36*pi 【例】：求方程组的解。 [x y] = solve(1/x^3+1/y^3=28,1/x+1/y=4) x = 1 1/3 y = 1/3 1 符号微分方程求解 MATLAB 中，用大写字母 D 表示导数。例如Dy 表示 y‘，D2y 表示 y’‘，Dy(0) = 5 表示y(0)=5。 dsolve 函数可以求解符号常微分方程，格式为： dsolve(eqn,c,v) 求解常微分方程 eqn 在初值条件 c 下的特解。 【例】：求微分方程 的通解。 y =dsolve(Dy*x^2 + 2*x*y = exp(x),x) y = (exp(x)+C1)/x^2 【例】：求微分方程 的特解， 。 y = dsolve((Dy-x^2)/(1+y^2),y(2)=1,x) y = 1/3*x^3-5/3 【例】：求微分方程 的通解。 [x,y] = dsolve(Dx = 4*x - 2*y,Dy = 2*x - y,t) x = 2*C2*exp(3*t)+1/2*C1 y = C1+C2*exp(3*t) 练习 1 1. 已知 x = 6, y = 5, 利用符号表达式求 2. 分解因式: (a) (b) 2200 3. 化简表达式: (a) (b) 4. 已知 完成下列运算: (1) (2) 计算 B 的逆矩阵并验证结果 (3) 取出包括 B 的主对角线元素的下三角阵作为 L 练习 2 1. 用符号方法求下列极限或导数 (1)