大学生科技创新论文.doc

大学生科技创新论文题目：MATLAB与计算方法关于数值计算问题姓名：张欣（1008300062） 杜昕阳（1008300061）班级：10级信息班专业：信息与计算科学内蒙古包头师范学院二零一二 年 十一 月 一 日一、 题目：MATLAB与计算方法关于数值计算问题二、 摘要：MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。1 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。计算方法比较全面地介绍了科学与工程计算中常用的计算方法，具体介绍了这些计算方法的基本理论与实际应用，同时对这些数值计算方法的计算效果、稳定性、收敛效果、适用范围以及优劣性与特点也作了简要的分析。全书共9章，内容包括引论、线性代数方程组求解方法、非线性方程求根、函数插值、函数逼近、矩阵特征值与特征向量的数值算法、数值积分与数值微分、常微分方程初值问题的数值解法、自治微分方程稳定区域的计算等。计算方法概念清晰，语言叙述通俗易懂，理论分析严谨，结构编排由浅入深，在分析问题时注重启发性，例题选择具有针对性且注重实际应用。前8章附有一定数量的习题，供读者学习时进行练习。计算方法可作为高等院校数学与应用数学、信息与计算科学、应用物理学、计算机科学等专业的高年级本科生和工科硕士研究生使用，也可供从事科学与工程计算的科技工作者参考。三、 引言：计算方法又称“数值分析”。是为各种数学问题的数值解答研究提供最有效的算法。主要内容为函数逼近论，数值微分，数值积分，误差分析等。常用方法有迭代法、差分法、插值法、有限元素法等。现代的计算方法还要求适应电子计算机的特点。数值分析即“计算方法”.四、 正文：插值与拟合 已知a) 求函数在处的函数值；（先编函数，再求值保存到向量中）b) 对上述数据进行多项式插值，作出多项式的图像，与原函数图象比较；（先列出差分表，再用牛顿插值公式编写出多项式函数）c) 对上述数据做线性拟合，作出多项式的图像；（先定义内积函数，再列出法方程，然后求解，最后编出多项式函数）d) 构造在区间内的5次切比雪夫多项式，并作出图像。（先生成切比雪夫点，再列出差分表，再插值）（1）函数代码为 结果为 （2）差分表代码为 结果为 牛顿插值公式代码为 结果为 （3）对上述数据做线性拟合，关键在于找出法方程。定义如下的内积函数：，这里取权重为1，然后再构造如下的法方程，（k=0,1）,由法方程得到线性方程组，从而通过解线性方程组得到线性拟合曲线的系数。线性拟合程序代码为 结果为 （4）利用的零点作插值节点，构造在区间内的5次多项式。由于切比雪夫零点是在区间之间，对于一般区间只要利用如下的变换即可。 所以首先要生成切比雪夫点。，利用上述区间变换的公式可得节点： ，进一步得到差分表，找出插值函数，并作出相关图像。代码为 结果为4、线性方程组的迭代解法（1）构造4阶Hilbert矩阵H，求； （2）分裂，求，； （3）设，找使得迭代法收敛速度最快；（4）分别用（3）中的方法，雅可比迭代法，高斯塞德尔迭代法，使 时终止，比较迭代次数与（1）（2）题中条件数的关系。（5）用共轭梯度法分别求解和，终止，比较迭代次数与系数矩阵条件数的关系。（1）代码为 结果为 （2）代码为结果为 （3）要使迭代速度最快，只要有w满足迭代矩阵的谱半径最小，则代码为： 结果为（4）雅阁比迭代法和高斯赛德尔迭代法程序为代码为结果为可以看出，条件数小于上述迭代法，结论成立。（5）共轭梯度法是构造一组正交向量，使得搜索方向快速下降，对于最速下降法的一种改进方法。程序代码为 结果为 条件数为 五、 结论：由于该题所给数据局限，两种迭代方法无法迭代出满意的结果，尤其是雅阁比迭代法，竟然是不收敛的，但是可以知道一个结论，就是条件数越小，迭代速度越快。所以，如果需要开发或者衡量一种迭代方法的话，可以考虑相应的条件数，令其足够小。而第三种迭代法，条件数小于上述迭代法，且迭代结果较为合理。六、 参考文献：1. 李庆阳、王能超、易大义.数值分析，第三版。武汉：华中理工大学出版社，19862. 冯康等编，数值计算方法，北京：国防工业出版社，19783. 黄友谦，李岳生，数值逼近，第二版，北京：高等教育出版社，19874. J.Stoer,R.Bulirsch.Introduction to Numerical Analysis ,Spring-Verlay, New York,19805. De Boor carl.A Practical Guide to Spline .Spring-Verlag,New York ,19786. J.M奥特加，数值分析，张丽君等译，北京：高等教育出版社，19737. 菜大用，数值代数，北京：清华大学出版社，19878. 路易斯.汉格曼,大卫.杨，实用迭代法.菜大用等译，北京：清华大学出版社，19819. 易大义、