数值分析教教案11.docVIP

2.1 计算积分的MATLAB符号法 使用MATLAB的符号计算功能，可以计算出许多积分的解析解和精确解，只是有些精确解显示冗长烦杂，这时可以用vpa或eval函数把它转换成位数有限的数字，有效数字的长度可按需选取。符号法计算积分非常方便，常常用它得到的结果跟近似计算结果进行比较，所以先予以介绍。 求积分的符号运算命令int（取自integrate前三个字)，调用格式： s=int(fun,v,a,b) ① 输入参量fun是被积函数的符号表达式，可以是函数向量或函数矩阵。 ② 输入参量v是积分变量，必须被界定成符号变量。如果被积函数中只有一个变量时可以缺省。 ③ 输入参数a、b为定积分的积分限，缺省时输出被积函数fun的一个原函数。 ④ 输入参量s为积分结果。s为有理表达式而过于冗长时，可在fun两端加引号，使它自动转换成默认的32位有效数字，或者用vpa或eval命令把它转换成有限长度的小数。 【例2-1】计算。 解：该函数的图像如图2-1所示，生成方法如下： y=1:10; ff=exp(-y.^2)+log(y); plot(ff,linewidth,2),title(exp(-y.^2)+log(y)),grid 图2-1 函数的图像 在编辑窗口输入： int(exp(-y^2)+log(y),1,10) 运行得到： ans = 1/2*pi^(1/2)*erf(10)+