一个通用数值积分的matlab实现(带附录)(5)

2019-08-03 14:14

end

if nargout==2 out1=djf; out2=n; end case 'mjx'

x=linspace(a,b,n+1)'; %n等分点记入数组中 h=(b-a)/n;

djf=0; %累加器清零

for i=1:1:n

djf=djf+h*f((x(i)+x(i+1))/2); end

if nargout==0 out1=djf; end

if nargout==1 out1=djf; end

if nargout==2 out1=djf; out2=n; end case 'tx'

x=linspace(a,b,n+1)'; %n等分点记入数组中

h=(b-a)/n;

djf=0; %累加器清零

for i=1:1:n

djf=djf+(f(x(i))+f(x(i+1)))*h/2; end

if nargout==0

17

out1=djf; end

if nargout==1 out1=djf; end

if nargout==2 out1=djf; out2=n; end case 'sps'

x=linspace(a,b,n+1)'; %n等分点记入数组中

h=(b-a)/n;

djf=0; %累加器清零

for i=1:1:n

djf=djf+(h/6)*(f(x(i))+4*f((x(i)+x(i+1))/2)+f(x(i+1))); end

if nargout==0 out1=djf; end

if nargout==1 out1=djf; end

if nargout==2 out1=djf; out2=n; end case 'nc'

x=linspace(a,b,n+1)'; %n等分点记入数组中

h=(b-a)/n;

18

djf=0; %累加器清零

for i=1:1:n

djf=djf+(h/90)*(7*f(x(i))+32*f(x(i)+h/4)+12*f(x(i)+2*h/4)+32*f(x(i)+3*h/4)+7*f(x(i)+h));

end

if nargout==0 out1=djf; end

if nargout==1 out1=djf; end

if nargout==2 out1=djf; out2=n; end case 'zsy'

n=1;s1=(f(a)+4*f((a+b)/2)+f(b))*(b-a)/n/6; n=2; s2=s1+1;

while abs(s2-s1)>wc s1=s2; n=n*2; s=0; h=(b-a)/n; for i=1:n

s=s+4*f(a+(i-1/2)*(h)); end; for i=1:n-1 s=s+2*f(a+i*h);

19

end;

s=(f(a)+f(b)+s)*h/6; s2=s; end; if nargout==0 out1=s2; end

if nargout==1 out1=s2; end

if nargout==2 out1=s2; out2=n; end end

20

致 谢

这篇论文是在段敏老师细心的指导下完成的。导师渊博的专业知识,严谨的治学态度,精益求精的工作作风,诲人不倦的高尚师德,严以律己、宽以待人的崇高风范,朴实无华、平易近人的人格魅力对我影响深远。不仅使我树了远大的学术目标、掌握了基本的研究方法,还使我明白了许多待人接物与为人处世的道理。本论文从选题到完成,每一步都是在导师的指导下完成的,倾注了导师的心血。至此,我要向我的导师表示衷心的感谢!四年学业即将完成,人生学习的道路才刚刚开始,以后的道路上我将努力的学习,好好地工作,给老师争光,给学校争气!

黄 凯

2015年5月20日于合肥学院

21


一个通用数值积分的matlab实现(带附录)(5).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:PLC期末复习题

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: