西安交通大学计算方法实验.docxVIP

目录实验一：共轭梯度法求解线性方程1一、实验原理1二、程序框图3三、程序代码及说明4四、使用案例分析5实验二： 三次样条插值6一、实验原理6二、程序框图7三、程序代码及说明8四、使用案例分析11实验三： Romberg 积分14一、实验原理14二、程序框图15三、程序代码及说明16四、使用案例分析17实验四： 四阶 Runge-Kutta 法求解常微分方程的初值问题19一、实验原理19二、程序框图19三、程序代码及说明20四、使用案例分析20实验一：共轭梯度法求解线性方程实验原理本次实验研究的是共轭梯度法解线性方程组的使用方法，其待求解的方程组的系数矩阵必须是对称正定的矩阵A。同时对于良态大型稀疏性方程组，该法只需经过比方程阶数n小很多的迭代次数就能获得满足精度要求的近似解。同所有反演问题都是寻找目标函数的极值来获取对应的极小值作为反演结果一样，共轭梯度法的基本原理就是把求解线性方程组的问题转化成求解一个与之等价的二次函数极小化的问题。然后从任意给定的初始点出发，沿一组关于矩阵A的共轭方向进行线性搜索，在无舍入误差的假定下，最多迭代n次，就可以求得二次函数的极小点，也即求得了线性方程组的解。定理：设A是n阶对称正定矩阵，则x*是方程的解的充分必要条件是x*是二次函数的极小值点，即 （1.1）共轭梯度法在形式上具有迭代法的特征，即给定初始向量，由迭代格式 （1.2）产生的迭代序列,