Dwcsoic数值研究实验报告.docVIP

七夕，古今诗人惯咏星月与悲情。吾生虽晚，世态炎凉却已看透矣。情也成空，且作“挥手袖底风”罢。是夜，窗外风雨如晦，吾独坐陋室，听一曲《尘缘》，合成诗韵一首，觉放诸古今，亦独有风韵也。乃书于纸上。毕而卧。凄然入梦。乙酉年七月初七。 -----啸之记。 《数值分析》 实 验 报 告 册 姓名: 学号: 专业: 年级: 计算机科学学院 计算机应用教研室 2008 年 春季 学期 目 录 TOC \o \h \z HYPERLINK \l _Toc191956312 实验一 3 HYPERLINK \l _Toc191956313 实验二 5 HYPERLINK \l _Toc191956314 实验三 7 HYPERLINK \l _Toc191956315 实验四 10 HYPERLINK \l _Toc191956316 实验五 12 HYPERLINK \l _Toc191956317 实验六 15 实验七……………………………………………………………………18  实验一 一、课题名称 非线性方程数值解法 二、目的和意义 学会常用的插值方法，求函数的近似表达式，以解决其它实际问题；明确插值多项式和分段插值多项式各自的优缺点；熟悉插值方法的程序编制；如果绘出插值函数的曲线，观察其光滑性。 三、计算公式 Lagrange插值公式： 牛顿插值公式： 四、结构程序设计 程序设计: #includemath.h float f(float x) { return((x*x*x-1)/3); /*牛顿迭代函数*/ } main() { float x1,x2,eps,d;int k=0; clrscr(); printf(\n input x1=); /*输入迭代初值*/ scanf(%f,x1); printf(\n input eps=); /*输入求解精度eps*/ scanf(%f,eps); do{ k++; x1=x2; x2=f(x1); printf(\n %d %f\n,k,x2); }while(fabs(x2-x1)=eps); printf(the root of f(x)=0 is x=%f,k=%d\n,x2,k); /*输出x和迭代次数k*/ getch(); } 五、结果讨论和分析 计算结果分析： 将六种迭代格式分别代入程序试验： 第一种格式：无论何值都无法求出，即发散 第二种格式：初值为任意的x（x2=1），精度为0.00001 X=-0.347296，k=6 其他值为发散。 第三种格式：初值为任意的x（x0），精度为0.0001 X=1.879372，k=10 其他值为发散。 第四种格式：初值为任意值，精度为0.00001 X=-0.347296,k=5 （5）第五种格式：初值为任意值,精度为0.00001 X=-0.347296,k=4 第六种格式：初值为任意值，精度为0.00001 X=-0.347296,k=4 由此可知不同的初值对公式的计算有影响，当初值不满足函数的收敛条件时，无法计算结果，函数发散。 精度的大小不同也使迭代函数迭代的次数不同，从而影响xn的近似程度。 实验二 一、课题名称 解线性方程组的直接方法 二、目的和意义 掌握线性方程组直接接法的基本思想；了解不同数值方法解线性方程组的原理、实现条件、使用范围、计算公式；培养编程与上机调试能力。 三、计算公式 消去法 设a(k)kk=0,对k=1，2，……，n-1计算 mik=a(k)ik/a(k)kk a(k+1)ij=a(k)ij-mika(k)kj i,j=k+1,k+2,……,n b(k+1)i=b(k)i-mikb(k)k n xn=b(n)n/a(n)nn n j=i+1 xi=(b(i)i-Σa(i)ijxj)/a(i)ii i=n-1,n-2,……,1 j=i+1 平方根法