3-matlab符号运算与数值积分‘程序.ppt

大连海事大学数学系 第三讲 Matlab符号运算与数值微积分 数学软件与实验 大连海事大学数学系 本讲主要内容 符号运算 一、Matlab 符号运算 符号运算举例 符号对象 符号对象的建立 符号对象的建立 符号表达式 基本符号运算 符号对象的基本运算 findsym findsym 举例 subs subs 举例 符号矩阵 符号矩阵与数值矩阵的转化 将数值矩阵转化成符号矩阵 任意精度的数学运算 在symbolic中有三种不同的算术运算： 浮点算术运算 1/2+1/3 －－(定义输出格式format long) ans = 0.83333333333333 符号运算 sym(1/2)+(1/3) ans = 5/6 －－精确解 任意精度算术运算 digits(n) —— 设置可变精度，缺省16位 vpa(x,n) —— 显示可变精度计算 digits(25) vpa(1/2+1/3) ans = .8333333333333333333333333 vpa(5/6,40) ans = .8333333333333333333333333333333333333333 a=sym([1/4,exp(1);log(3),3/7]) a = [ 1/4,exp(1)] [log(3), 3/7] vpa(a,10) ans = [.2500000000, 2.718281828] [1.098612289, .4285714286] 二、 符号极限 limit(F,x,a) 返回函数F在符号变量x趋于a的极限 limit(F,a) 返回函数F在符号变量findsym(F)趋于a 的极限 limit(F) 默认a=0 limit(F,x,a,’left’) 返回函数F在符号变量x趋于a的左极限 limit(F,x,a,’right’) 返回函数F在符号变量x趋于a的右极限 例 求 实验过程： syms x t; limit((1+2*t/x)^(3*x),x,inf) 执行后得 ans = exp(6*t) 验证重要极限 实验过程： syms x limit(sin(x)/x) 执行后得： ans = 1 例：求exp(1/x)在x=0点的极限以及左右极限。 三、 符号微分 diff(f) — 对缺省变量求微分 diff(f,v) — 对指定变量v求微分 diff(f,v,n) —对指定变量v求n阶微分 例 已知 求y的一阶和三阶导数。 实验过程： syms x y=exp(2*x)*log(x^2+1)*tan(-x) diff(y,x) diff(y,x,3) 例 已知 求z 关于变量x的一至二阶偏导，以及关于x,y的二阶混合偏导。 实验过程：syms x y; z=(x^2+y^2)*exp((x^2+y^2)/(x*y)) diff(z,x) diff(z,x,2) diff(diff(z,x),y) 四、 符号积分 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b) 区间求定积分 例1.计算二重不定积分 解法二： syms x y F=int(int(x*exp(-x*y),x),y) F = 1/y*exp(-x*y) 例：计算 实验过程 syms x int(exp(-x^2),x,0,Inf) 输出结果： ans = 1/2*pi^(1/2) 五、其他符号计算 taylor(f,n) —— 泰勒级数展开 ztrans(f) —— Z变换 Invztrans(f) —— 反Z变换 Laplace(f) —— 拉氏变换 Invlaplace(f) —— 反拉氏变换 fourier(f) —— 付氏变换 Invfourier(f) —— 反付氏变换 六、 数值积分 1. 数值积分基本原理 简单的梯形法、辛普生(Simpson)法。 它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。这样求定积分问题就分解为求和问题。 2. 矩形数值积分公式及其实现 （左矩形公式） （右矩形公式）