数值分析程序代码（MATLAB）+实验报告.rar

本文件针对于数值分析课程，主要内容是数值分析课程实验，包括：牛顿法求函数零点、牛顿插值法、求三次样条插值多项式、通用多项式拟合、插值型求积公式、Runge-Kutta 4阶算法等。本文件仅为个人课程实验程序代码及报告，仅供参考！

function [Snn]=ComSimpson（fabe）
% Made by 付林
% Input
% f:需要计算的积分表达式，如y=4/（1+x^2）;
% a:为待求积分区间的下界
% b：为待求积分区间的上界
% e:估计的精度要求 |En（f）|% Output
% Sn:估计值
% n:节点个数

syms x;
% 计算复合辛普森公式的误差：En（f）=-（b-a）^5/2880*n^4*f（4）（ζ）
% 计算函数表达式f的四阶导数f4
f4=diff（f4）;
% 计算f（4）（ζ）时，x在区间[ab]上的间距为s
s=（b-a）/10000;
% 计算f（4）（ζ）在区间[ab]上的最大值
Max=abs（max（subs（f4xa:s:b）））;
% |En（f）|% 将为sym类型的Max转化为double型
Max=double（Max）;
% 计算划分节点数
n=ceil（（（b-a）^5*Max/（2880*e））^（1/4））;
% 每个划分节点的x值为X
X=linspace（abn）;
% 计算每个间距
h=（b-a）/（n-1）;
% 将函数表达式f转化为matlab函数Y进行直接调用
Y=matlabFunction（f）;
% 计算每个X节点的y值
y=Y（X）;
% 复合辛普森计算公式
Sn=（y（1）+4*sum（y（2:2:n-1））+2*sum（y（3:2:n-1））+y（n））*h/3;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        861  2017-12-21 20:38  ComSimpson\ComSimpson.m

     文件        845  2017-12-17 09:43  ComSimpson\ComSimpson_1.m

     文件        786  2017-12-21 20:35  ComTrapz\ComTrapz.m

     文件        810  2017-12-16 23:44  ComTrapz\ComTrapz_1.m

     文件        618  2017-12-16 23:29  ComTrapz\ComTrapz_2.m

     文件        584  2017-12-16 22:46  ComTrapz\ComTrapz_3.m

     文件       1062  2017-11-29 00:44  NewtonInt\NewtonInt.m

     文件       1064  2017-11-29 00:40  NewtonInt\NewtonInt_1.m

     文件       1039  2017-12-21 18:55  NewtonRoots\NewtonRoots.m

     文件      14337  2017-12-02 11:13  PolyFit\PolyFit.jpg

     文件       1234  2017-12-21 20:34  PolyFit\PolyFit.m

     文件       1220  2017-12-21 20:30  PolyFit\PolyFit_1.m

     文件       1288  2017-12-21 20:39  Romberg\Romberg.m

     文件       1292  2017-12-17 15:07  Romberg\Romberg_1.m

     文件      22320  2017-12-26 17:43  Runge-Kutta4\Runge-Kutta4.xlsx

     文件        780  2017-12-22 21:00  Runge-Kutta4\Runge_Kutta4.m

     文件      15193  2017-12-22 20:36  Runge-Kutta4\Step_0.01.jpg

     文件      15816  2017-12-22 21:00  Runge-Kutta4\Step_0.05（1）.jpg

     文件      15309  2017-12-22 20:31  Runge-Kutta4\Step_0.05.jpg

     文件      16160  2017-12-22 20:32  Runge-Kutta4\Step_0.1.jpg

     文件      16307  2017-12-22 20:59  Runge-Kutta4\Step_0.15（1）.jpg

     文件      15480  2017-12-22 20:35  Runge-Kutta4\Step_0.15.jpg

     文件      15790  2017-12-22 21:00  Runge-Kutta4\Step_0.2（1）.jpg

     文件      15434  2017-12-22 20:37  Runge-Kutta4\Step_0.2.jpg

     文件       3451  2017-12-21 20:19  Spline3\Spline3.m

     文件      19317  2017-12-20 21:50  Spline3\Spline3_1.jpg

     文件       3367  2017-12-26 15:10  Spline3\Spline3_1.m

     文件      19398  2017-12-26 15:01  Spline3\Spline3_11.jpg

     文件      21110  2017-12-20 21:51  Spline3\Spline3_2.jpg

     文件      21127  2017-12-26 14:59  Spline3\Spline3_21.jpg

............此处省略13个文件信息