数值方法实验报告模板..docVIP

2010-2011学年 《数值分析》实验报告 学 院：信息与数学学院 专业班级：信计 10802 班 学生姓名：施耀凯、宋龙飞、李聪、 陈慧、沈江炎 班级序号：10、8、16、33、21 题目序号（组）：N0.1 组长：陈慧 实验一 1. 题目 用追赶法解方程组： 2. 理论分析及算法描述 2.1 理论分析 利用矩阵的直接三角分解法来推导解对角占优的三对角方程组的计算公式．由系数阵的特点,可以将分解为两个三角阵的乘积,即 　　 (1) 其中为下三角阵, 为单位上三角阵．设 (2) 其中，，为待定系数．比较式(2)两端即得   　　 (3) 由式(3)得到 ． 这就是说,由的假设条件完全确定了,实现了的分解． 求解等价于解两个三角方程组与,先后求与,从而得到以下解三对角方程组的追赶公式： 步1 计算的递推公式 ； 步2 解：； 步3 解：． 2.2 算法描述 追赶法解三对角方程组,其步骤如下： 步1 计算 （1）冲掉 （2）冲掉 　 ． 步2 解 （1）解出放到里 （2）解出放到里 ． 步3 解 （1） ， （2） 放到里 　 ． 3. 源程序 #include stdio.h #include math.h /*追赶法解三对角方程组*/ void main() { int i; float a[3],b[4],c[3]; float B[4]={5,-12,11,-1}; for(i=0;i3;i++) a[i]=-1; for(i=0;i4;i++) b[i]=2; for(i=0;i3;i++) c[i]=-1; c[0]/=b[0]; for(i=1;i3;i++) c[i]/=(b[i]-a[i-1]*c[i-1]); B[0]/=b[0]; for(i=1;i4;i++) B[i]=(B[i]-a[i-1]*B[i-1])/(b[i]-a[i-1]*c[i-1]); for(i=2;i=0;i--) B[i]=B[i]-c[i]*B[i+1]; printf(输出结果为:); for(i=0;i4;i++) printf(%5.0f,B[i]); } 4. 运行结果及其分析 结果分析：对计算的结果，冲掉原有的数据存储单元，大大节省了存储空间． 实验二 1. 题目 试求三次样条差值，并满足： （1）， （2） 已知： 0．25 0．30 0．39 0．45 0．53 0.5000 0.5477 0.6245 0.6708 0.7280 2. 理论分析及算法描述 2.1 理论分析 由于s（x）在区间[]上是三次多项式，故在[]上是线性函数，可表示为： 对积分两次并利用及定出积分常数得： （1） 说明：将上式分解展开后各项系数分别为： 三次项：； （2） 二次项：； （3） 一次项： （4） 常数项： （5） 对S(x)求导得： 所以 由(i=1,2,...,n-1) 得： 其中 边界条件为 由公式 得 即 从中解出Mi (i=0,1,...,n) 带入方程（1）得三次样条S(x).也可通过式（2）（3）（4）（5）求出其各项系数从而得到S(x). 2、边界条件为已知 从中解出Mi(i=1,2,...,n-1)带入方程（1）得 三次样条S(x)，也可通过式（2）（3）（4）（5）求出其各项系数从而得到S(x). 2.2 算法描述 三弯矩方程解三次样条函数算法： 步一：输入初始数据，（j=0,1,..n）及，和n 步二：j从0到n-1计算= 步三：j从1到n-1由式