第六讲定积分与定积分的近似计算..docVIP

第六讲 定积分与定积分的近似计算 实验目的 １．通过本实验加深理解积分理论中分割、近似、求和、取极限的思想方法. ２．学习并掌握用matlab求不定积分、定积分、二重积分、曲线积分的方法. ３．学习matlab命令sum、symsum与int. 4. 了解定积分近似计算的矩形法、梯形法。（***） 实验内容 学习matlab命令 （１）求和命令sum调用格式. sum(x),给出向量x的各个元素的累加和,如果x是矩阵,则sum(x)是一个元素为x的每列列和的行向量. 例4.1.x=[1,2,3,4,5,6,7,8,9,10]; sum(x) ans=55 例4.2.x=[1,2,3;4,5,6;7,8,9] x= 1 2 3 4 5 6 7 8 9 sum(x) ans=12 15 18 （２）求和命令symsum调用格式. symsum(s,n), 求 symsum(s,k,m,n),求 当x的元素很有规律,比如为表达式是的数列时,可用symsum求得x的各项和,即 symsum= symsum 例4.3.syms k n symsum(k,1,10) ans=55 symsum(k^2,k,1,n) ans=1/3*(n+1)^3-1/2*(n+1)^2+1/6*n+1/6 （３）matlab积分命令int调用格式 int(函数) 计算不定积分 int(函数,变量名) 计算不定积分 int(函数) 计算定积分 int(函数变量名) 计算定积分 计算不定积分 例4.4.计算 解:输入命令: int(x^2*log(x)) 可得结果: ans=1/3*x^3*log(x)-1/9*x^3 注意设置符号变量. 例4.5.计算下列不定积分: 1. 2. 3. 解:首先建立函数向量. syms x syms a real y=[sqrt(a^2-x^2),(x-1)/(3*x-1)^(1/3),x^2*asin(x)]; 然后对y积分可得对y的每个分量积分的结果. int(y,x) ans = [1/2*x*(a^2-x^2)^(1/2)+1/2*a^2*asin((1/a^2)^(1/2)*x), -1/3*(3*x-1)^(2/3)+1/15*(3*x-1)^(5/3), 1/3*x^3*asin(x)+1/9*x^2*(1-x^2)^(1/2)+2/9*(1-x^2)^(1/2)] ３．定积分的概念. 定积分是一个和的极限.取,积分区间为,等距划分为20个子区间. x=linspace(0,1,21); 选取每个子区间的端点,并计算端点处的函数值. y=exp(x); 取区间的左端点乘以区间长度全部加起来. y1=y(1:20);s1=sum(y1)/20 s1=1.6757 s1可作为的近似值. 若选取右端点: y2=y(2:21);s2=sum(y2)/20 s2=1.7616 s2也可以作为的近似值.下面我们画出图象. plot(x,y);hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],’b’) end 如果选取右端点,则可画出图象. for i=1:20; fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i+1),y(i+1),0],’b’) hold on end plot(x,y,’r’) 在上边的语句中,for … end是循环语句,执行语句体内的命令20次,fill命令可以填充多边形,在本例中,用的是兰色(blue)填充.从图上看,当分点逐渐增多时,的值越来越小,读者可试取50个子区间看一看结果怎样.下面按等分区间计算