数值分析课程设计 沈阳航空航天大学 课程设计用纸 正 文
C(n)knnhnnt?j(?1)n?k?(t?j)dt?dt?nk!(n?k)!?0?....................(2) b?a?0j?0k?jj?0j?kj?k由于是多项式的积分，柯特斯系数的计算不会遇到实质性的困难。 当n=1时，C(1)0?C(1)1?1，这时的求积公式就是我们所熟悉的梯形公式。 2当n=2时，按（2）式，这时的柯特斯系数为
(2)C0?C1(2)(2)C2121(t?1)(t?2)dt?4?06124???t(t?2)dt?206 121??t(t?1)dt?4062b?aa?b[f(a)?4f()?f(b)]就为辛普森求积公式。?f(x)k
162 将区间[a,b]划分为n等份，分点xk=a+kh,h=(b-a)/n,k=0,1,...,n,在每个子区间
相应的求积公式S?[xk,xk+1](k=0,1,?,n-1)上采用辛普森公式，若记
则有
记
=为复合辛普森公式，
余项为,。
于是当f（x）
，与复合辛普森公式余项=,
第 6 页
数值分析课程设计 沈阳航空航天大学 课程设计用纸 正 文
2.程序流程图
开始 输入n,a,b,g
定义未知变量x，将g定义为函数f(x)，计算h的值 调用函数利用for循环计算和号的值 利用复合辛普森公式，带入数值，计算S的值
输出S值
结束
3.计算?sin?x?dx
0?> Untitled5(2,0,pi,@(x)sin(x)) f =
@(x)g(x) ans =
2.0046
>> Untitled5(3,0,pi,@(x)sin(x)) f =
@(x)g(x)
第 7 页
数值分析课程设计 沈阳航空航天大学 课程设计用纸 正 文
ans =
2.0009
>> Untitled5(5,0,pi,@(x)sin(x)) f =
@(x)g(x) ans =
2.0001
>> Untitled5(7,0,pi,@(x)sin(x)) f =
@(x)g(x) ans =
2.0000
>> Untitled5(10,0,pi,@(x)sin(x)) f =
@(x)g(x) ans =
2.0000 >>
第 8 页
数值分析课程设计 沈阳航空航天大学 课程设计用纸 正 文
二、总结
通过这两周的MATLAB课程设计，我学到了很多知识，也了解到MATLAB是
一个运用于很多领域的软件。初次接触MATLAB，感到自己还有很多不足，对一些程序的设计还没有足够精通，但是通过两周的学习，在MATLAB设计方面，觉得自己还是有所进步。在课设过程中，也遇到过很多问题，有时候程序运行不出来或者结果错误，仅仅只是因为忘记加分号，或者是多加了括号，或者是字母的大小写不一致等，这让我意识到MATLAB的严谨性，也在此方面更加注意，避免此类错误发生。有时是因为函数编写思路不正确，通过上网查找，去图书馆查找，以及请教同学等方法来解决这些问题。MATLAB以矩阵计算为基础，在科学研究以及工程设计邻域有着十分广泛的应用，我以后会继
续学习，掌握更多的知识，更好的运用MATLAB。
完稿日期:2016 年1 月 4日
第 9 页
数值分析课程设计 沈阳航空航天大学 课程设计用纸 参考文献
参考文献
[1] 李庆阳,王能超,易大义.数值分析.5版.北京:清华大学出版社,2008.12
[2] 张霞萍,MATLAB 8.X程序设计及典型应用.西安:西安电子科技大学出版社, 2014.7
第 10 页