第三章(符号计算与微积分实验).pptVIP

* 符号计算与微积分实验 符号表达式及其应用 微分和积分符号计算 台劳展开式符号计算 旋转曲面绘制方法 ? ? ? ? 在符号计算中, 符号表达式是主要操作对象. 符号表达式——符号变量、运算符、函数、数字组成 在定义符号表达式之前,首先要创建符号变量. 创建符号变量的常方法是在 syms 空格 跟符号变量。 符号计算又称为计算机代数，以符号形式处理数学表达式，关注准确的计算和公式推导。符号计算不仅用于数学研究，还可以用于工程计算 。1993年MathWorks公司购买Maple的使用权,开发了符号计算工具箱——Symbolic Math Toolbox. syms 符号变量1 符号变量2 … syms 符号变量1 符号变量2 … 例. 用ezplot命令绘制 函数 y = e – x2的图形 syms x y y=exp(-x^2); ezplot(y) 符号变量的定义 例. 转换数值变量为符号变量 A=[1/3,1/4;1/5,1/7] B=sym(A) B = [ 1/3, 1/4] [ 1/5, 1/7] syms x f=2*x^2+3*x-5; g=x^2+x-7; h1 = f+g,h2 = expand(f*g ) factor(h2) h1 = 3*x^2+4*x-12 h2 =2*x^4+5*x^3-16*x^2-26*x+35 ans =(x-1)*(2*x+5)*(x^2+x-7) expand() —— 多项式展开 factor() —— 多项式因式分解 例.多项式的符号表达式创建与运算 syms x ; f = x*sin(1/x);figure(1),ezplot(f,[-0.1, 0.1]) figure(2) ezplot(exp(-0.5*x)*sin(5*x),[0,10,-.8,.9]) 例 .分别对表达式 x sin(1/x) 和e-0.5xsin5x 绘图. S1=subs(S, ‘old’, ‘new’) 用‘new’置换符号表达式S中的’old’ 例 . 输入不同的参数a,b绘制函数的图形 f(x)=exp(a x)sin(b x) syms a b x f=exp(a*x)*sin(b*x); f1=subs(f,a,-0.6); f2=subs(f1,b,0.5); ezplot(f2,[0,2*pi]) 符号表达式中变量替换 调用格式：double(A) numeric(A) 例5. A=sym([1/3,2/5;10/7,2/5]) numeric(A) A = [ 1/3, 2/5] [ 10/7, 2/5] ans = 0.3333 2.5000 1.4286 0.4000 将符号矩阵转化为数值矩阵 复杂表达式的化简 syms x y z a b c f=(x+y)*(a+b^c)^z/(x+a)^2 pretty(f) 常用化简命令: 降幂排列:collect(P,x); 展开:expand(P); 重叠: horner(P); 因式分解: factor(P); 5. 化简: simplify(P) 微积分符号计算 diff(f) — 对缺省变量求微分 diff(f,v) — 对指定变量v求微分 diff(f,v,n) —对指定变量v求n阶微分 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b) 区间求定积分 例4.计算 f = 1/(5+4cos(x)) 关于x的导数 syms x f=1/(5+4*cos(x)) ezplot(f) f1=diff(f,x,1) ezplot(f1) f1 = 4/(5+4*cos(x))^2*sin(x) 例5 计算不定积分 syms x int(exp(a*x)*sin(b*x)) g=simplify(ans) g = exp(a*x)*(-b*cos(b*x)+a*sin(b*x))/(a^2+b^2) diff(g) f=simplify(ans) f = exp(a*x)*sin(b*x) 例. 绘函数 (a = 1, b = 3 )在 [0,3.2]上的图形. 并计算 syms a b x f=exp(a*x)*sin(b*x) f1=subs(f,a,1),f1=subs