第5章 数积分.pptVIP

数值计算与MATLAB 第五章 数值积分 我们知道,若函数f(x)在区间[a,b]上连续且其原函数为F(x),则可用Newton-Leibnitz公式 求定积分的值 , Newton- Leibnitz公式 无论在理论上还是在解决实际问题上都起了很大作用，但它并不能完全解决定积分的计算问题，因为积分学涉及的实际问题极为广泛，而且极其复杂，在实际计算中经常遇到以下三种情况： (3) 被积函数f(x)没有具体的解析表达式, 其函数 关系由表格或图形表示。 对于这些情况, 要计算积分的准确值都是十分困难的。由此可见, 通过原函数来计算积分有它的局限性, 因而研究一种新的积分方法来解决Newton-Leibniz公式所不能或很难解决的积分问题, 这时需要用数值解法来建立积分的近似计算方法。 将积分区间细分,在每一个小区间内用简单函数代替复杂函数进行积分，这就是数值积分的思想。 5.1计算积分的matlab符号法 使用matlab的符号计算功能，可以计算出许多积分的解析解和精确解，只是有些精确解显示冗长繁杂，这时可以用vpa或eval函数把它转换成位数有限的数字，有效数字的长度可按需选取。 求积分的符号运算命令为int（取自integerate的前3个字母），调用格式为 s=int(fun,v,a,b) 1)输入参量fun是被积函数的符号表达式，可以是函数向量或函数矩阵。 2）输入参量v是积分变量，必须被界定成符号变量。如果被积函数中只有一个变量时可以缺省。 3）输入参量a、b为定积分的积分限，缺省时输出被积函数fun的一个原函数。 4）输入参量s为积分结果。s为有理表达式有时过于冗长时，可用vpa或eval命令把它转换成有限长度的小数。 例5.1 计算 解：该函数图线的生成方法如下： y=1:10; ff=exp(-y.^2)+log(y); plot(ff,linewidth,2),title(exp(-y.^2)+log(y)),grid 在matlab命令窗口中输入 int(exp(-y^2)+log(y),1,10) 回车得到 ans = 1/2*pi^(1/2)*erf(10)-9+10*log(2)+10*log(5)-1/2*pi^(1/2)*erf(1) 输出的结果比较复杂，下面是用两种方式进行转换的输出结果 是比较他们的差别 输入 eval(1/2*pi^(1/2)*erf(10)+10*log(2)+10*log(5)-9-1/2*pi^(1/2)*erf(1)) ans = 14.1653 输入 vpa(1/2*pi^(1/2)*erf(10)+10*log(2)+10*log(5)-9-1/2*pi^(1/2)*erf(1)) ans = 14.165253722580787828429530852382 例5.2 计算 解：该函数的图像生成方法如下 x=-10:1:10; ff=1./(x.^2+1); plot(ff,linewidth,2),title(1./(x.^2+1)),grid 在matlab命令窗口中输入 syms x; f=1/(x^2+1); a=int(f,-inf,inf) 回车得到 a = pi 得到精确值 例5.3计算 解：在matlab命令窗口中输入 syms x t; f=2*t*x; a=int(f,x,1,sin(t)) 回车得到 a = t*(sin(t)^2-1) 例5.4求 ，准确结果是 解：在matlab命令窗口中输入 syms x; Is=int(sqrt(log(1/x)),x,0,1) 回车的到 Is = int(log(1/x)^(1/2),x = 0 .. 1) 输入 vpa(Is) ans = .88622692545275801364908374167057 再输入 vpa(Is,5) ans = .88623 例5.5计算单位域上的积分 先把二重积分 转换为二次积分的形式 解在matlab命令窗口中输入 syms x y; Q=int(int(exp(-x^2/2)*sin(x^2+y),x,-sqrt(1-y^2),sqrt(1-y^2)),y,-1,1) 回车得到 Q = int(-i*pi^(1/2)*(exp(i*y)^2*erf(1/2*(1-y^2)^(