計算方法A上机报告.docxVIP

计算方法A：上机报告姓名：学号：电气工程学院电力系统及其自动化共轭梯度法：算法流程：程序使用说明：运行gongetidu.m文件首先出现输入提示‘1-example3.2,2-other examples:’如果输入为1则计算书上习题3.2实例，接下来的提示输入为‘N:’，代表要计算的矩阵阶数，按习题3.2计算可以输入100或者200或者400.输入后出线提示符’Eps:’，用来设置精度，例如可以输入’1E-3’.如果第一次选择输入为2则可以计算任意实例，接下来输入提示为’A:’,输入方程组的系数矩阵，输入格式例如：[1,2,3;4,5,6;7,8,9]接下来的输入提示’b:’,为方程的右端常数项，输入格式例如:[1;2;3]接下来的输入提示’x0:’ ，为方程解迭代的初始值，可以输入:[0;0;0]接下来提示’Eps:’,输入精度，例如：1E-3最后的解存储与’x0’中在命令行输入: disp(x0) 就可以查看算法原理:其中：搜索步长的计算公式：残向量：更新解向量的公式：更新搜索方向的公式：算例：例题3.2其中：，分别计算N=100，200，400时的解均为当取初始向量为时，循环分别执行了50，100，200次。三次样条插值：算法流程：程序使用说明：打开文件sanciyangtiao.m会有输入提示 ‘Enter the number of points:’输入数据点的个数，例