梯形积分及辛普生积分被积函数是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.绘出所得
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/26 18:39:56
梯形积分及辛普生积分被积函数是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))
辛普生公式
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))辛普生公式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
大按时代发粉我下 按时打发上;打发上;地方按时地方按时地方;老师傅拉;四六风;啊;色浪费;挖人;‘大师傅山豆根 士大夫敢死队风格