UC知道梯形积分及辛普生积分
来源：百度知道 编辑：UC知道 时间：2024/09/21 18:45:44
被积函数是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.绘出所得图形及误差图形
高手们帮帮忙吧，小弟在此谢过了
梯形公式
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.绘出所得图形及误差图形
高手们帮帮忙吧，小弟在此谢过了
看看是不是这样的：
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
大按时代发粉我下 按时打发上；打发上；地方按时地方按时地方；老师傅拉；四六风；啊；色浪费；挖人；‘大师傅山豆根 士大夫敢死队风格