作业帮 > 综合 > 作业

梯形积分及辛普生积分被积函数是exp(-x) 上限a=0,b=1梯形公式f(x)=(b-a)/2 * (f(a)+f(b

来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/11/11 01:21:22
梯形积分及辛普生积分
被积函数是exp(-x) 上限a=0,b=1
梯形公式
f(x)=(b-a)/2 * (f(a)+f(b))
辛普生公式
f(x)=(b-a)/6 * (f(a)+4*f((a+b)/2)+f(b))
要求1.用matlab编程并给出说明
2.输出结果的误差
3.绘出所得图形及误差图形
高手们帮帮忙吧,小弟在此谢过了
梯形积分及辛普生积分被积函数是exp(-x) 上限a=0,b=1梯形公式f(x)=(b-a)/2 * (f(a)+f(b
看看是不是这样的:
function zint
clear all;clc;
a=0;b=1;d=0.05;
x=a:d:b;y=func(x);
format long
y_trapz=trapz(y)*d %梯形积分
y_simps=quad(@func,a,b) %辛普生积分
y_accurt=eval(int('exp(-t)',a,b))% 精确值
%误差
y_trapz_error=y_trapz-y_accurt
y_simps_error=y_simps-y_accurt
%图形
plot(x,exp(-x),'o-')
function y=func(x)
y=exp(-x);
结果:
y_trapz = 0.632252245124816
y_simps = 0.632120560205377
y_accrt = 0.632120558828558
y_trapz_error = 1.316862962582599e-004
y_simps_error = 1.376819436238463e-009