08实验3插值与数值积分2教材.ppt

本次课的主要内容;;trapz(y) 按梯形公式计算定积分（单位步长）。 trapz(x,y) x , y同长度，输出 y 对 x 的按梯形公式计算的积分 （变步长）。;function s=simp(x,y) if mod(length(x),2)==0 error(数据点必须为奇数) end n=length(x);m=(n-1)/2; h=(x(n)-x(1))/2/m; s1=0; s2=0; for i=1:m s1=s1+y(2*i); end for j=1:m-1 s2=s2+y(2*j+1); end s=(y(1)+y(n)+4*s1+2*s2)*h/3;将[a,b]区间等分成2m个小区间：;令f(x)= xk ，用(11)式计算;将 f (x) =1, x , x 2 , x 3 依次代入，得：;将区间(a, b)作m等分，记h=(b-a)/m，xk=a+kh, k=0,1,...,m, 作变换;function s=gaussinteg(f,a,b,m) s=0; h=(b-a)/m;k=1:m+1; x(k)=a+(k-1)*h;j=1:m; z1(j)=(x(j)+x(j+1))/2+h/2/sqrt(3); z2(j)=(x(j)+x(j+1))/2-h/2/sqrt(3); for j=1:m s=s+feval(f,z1(j))+feval(f,z2(j)); end s=s*h/2;;三 用MATLAB实现定积分计算：;⑵ 用辛甫森公式求定积分;一 求卫星轨道长度;c; 我们用梯形公式和辛甫生公式来求卫星轨道的长度。;广义积分、二维数值积分;p65 3 , 4 2) ,10;四 蒙特卡罗方法(用随机模拟方法计算数值积分);若随机变量X 的概率分布密度是 p(x), a? x? b，则;通常还要先设法使被积函数在[0, 1] 范围取值再应用随机投点法。;4、运用蒙特卡罗方法计算重积分;计算二重积分：;§3 数值微分; 将区间[a, b] n等分，步长;二 应用实例——人口增长率; 1 记时刻t 的人口为x(t)，人口相对增长率为 ;表3 某地区20世纪70年代人口增长率数据;p93 5.2 ２) b，d 5) 10)