计算方法第一次上机报告.docVIP

《计算方法》实验报告一 题目： 班级： 计1402 学号：1403140126 姓名： 应成龙 2016.3.17 程序名： 一、问题和要求： 1. 土木工程和环境工程师在设计一条排水渠道时必须考虑渠道的各种参数（如宽度，深度，渠道内壁光滑度）及水流速度、流量、水深等物理量之间的关系。 假设修一条横断面为矩形的水渠，其宽度为B，假定水流是定常的，也就是说水流速度不随时间而变化。为了不同的工业目的（比如说要把污染物稀释到一定的浓度以下，或者为某工厂输入一定量的水），需要指定流量Q和B，求出水的深度。这样，就需要求解 其中Q 是水的流量（），U是流速（），H是水的深度（）。 一个具体的案例是求出渠道中水的深度H。 的收敛性，并求出使 || x(k + 1) – x(k) ||1 ≤0.0001 的近似解及相应的迭代次数。考虑用 雅可比迭代法；（2）赛德尔迭代法；（3）超松驰迭代法（( 取1.2，1.3，1.9，0.9）。 3. 对下面的9阶对称矩阵分别做LU分解和LDLT分解 二、设计总体方案及分析： 2.1 问题分析 第 第2题考查的是线性代数方程组的迭代法，分别是（1）雅可比迭代法，（2）赛德尔迭代法， （3）超松驰迭代法 第3题考查的是线性代数方程组的矩阵三角分解法，分别是LU分解和LDLT分解 2.2 概要设计 2.3 详细设计（二合一） 第’(x)是否收敛，从x=0开始+1-Xk的值小于0.000001，才停止。并输出所有的迭代值 1题用雅可比迭代法时，由于 公式中，迭代的次数是k次，所以需要用一个额外的数组来保存之前第k次迭代结果，我用了x[6][N]={ 0 }，来存放那个数组，通过控制每个次迭代与 k次迭代差=0.001控制循环次数，当每个差都同时=0.001，循环停止，迭代结束 2题用的是由于已经进行过的迭代，直接覆盖原值，所以 第3题用的是超松弛迭代法， 超松弛法在就行了。 第1题用的是 这题我用了三个嵌套循环，在最内层循环有两个并列循环，先给,u[][]的遍历，最终实现全部赋值。 第2题用的是 实际用到的是这3个公式 这一题我通过修改我自己写的 2.4 调试分析 第1题结果我就算出来两个206213，估计是 第2-1这 第2,2-3，都是在1的基础上稍作修改就很快算出了结果。 第1也是从公式入手，通过观察公式，发现规律，有了前面的经验之后很快就能写出来。 第2是由于公式选择的问题，导致一直编译出错，最后在借鉴了同学的代码之后完成。 第 第1 第2-2题 第2-3题 第3-1 第2 三、源程序及注释 第1题 #includeiostream #includeiomanip #includemath.h #define n 0.03 #define B 20 #define S 0.0002 #define Q 5 using namespace std; bool shoulian(double H) { if (0.8 / B*pow(n*Q / pow(S, 0.5) / (B + 2 * H), 0.6) 1) return true; else return false; } double diedai(double H) { return (pow(n*Q / pow(S, 0.5)*pow(B + 2 * H, 2 / 3), 0.6) / B); } int main() { double H; for (double h = 0; h = 100; h = h + 0.5) { if (shoulian(h)) { H = h; break; } } while (fabs(H - diedai(H))=0 { cout H \t diedai(H) \t endl; H = diedai(H); } cout H \t diedai(H) \t endl; } 第2-1题 #includeiostream #includeiomanip #includemath.h using namespace std; #define N 1000 int main() { int counter = 0; int a[][6] = { 4, -1, 0, -1, 0, 0, -1, 4, -1, 0, -1, 0, 0, -1, 4, -1, 0, -1, -1, 0, -1, 4, -1, 0