积分和级数的MATLAB实现.docVIP

积分与级数 1 积分问题 求解方法： Q=quad(F,0,2); %自适应法Simpson法 F=int(f,x); F=int(f,x,a,b); 1.1 定积分的数值解 例：绘制积分区域 f=inline(x.^3-x.^2+2*x+1); figure; hold on fplot(f,[-1000,1000]) patch([-600 -600:800 800],[0 f(-600:800) 0],r); 例：计算x.^3-x.^2+2*x+1的定积分 内联函数 M函数 f=inline(x.^3-x.^2+2*x+1); Q = quad(f,0,2); function y = myfun(x) y= x.^3-x.^2+2*x+1; Q = quad(@myfun,0,2); 针对函数值序列的数值积分 trapz(sin(0:0.01:pi))*0.01 % 梯形数值积分 例：编程实现上例功能 自定义数值积分函数 function area=getarea(f,x1,x2) delta=0.01; % 子区间的宽度 x=x1:delta:x2; y=feval(f,x); area=cumsum(y)*delta; area=getarea(f,0,2); area(end) area=get