求Ig,写哈过程

矩阵计算上机实验报告(哈工程)

哈爾滨工程大学 矩阵计算上机实验报告 姓名: 学号: 矩阵计算实验报告 实验名称 矩阵计算上机实验 指导教师 姓名 年级 学号 成绩 实验一、用共軛梯度法实对称线性方程组的解 一、实验目的:   1、学习使用matlab编写数值计算程序   2、了解共轭梯度法原理和编程方法。 二、实驗要:   1、阐述共轭梯度法原理和编程方法 2、根据共轭梯度法算法原理编写matlab程序,并运行的出结果 (提示:下文是Jacobi迭代解线性方程组的实验过程,仅供参考提交作业时,请删除!) 三、实验原理: 设n阶非奇异矩阵A的主对角元全不为0记D=diag( )是非奇异对角阵,做A的一个分裂: 记。迭代过程式为 这种迭代方法称为Jacobi迭代法。其迭代矩阵记为:右端向量记为: = - 记Jacobi迭代过程为: 对照等号兩边,得到计算(i=1,2,…,n)的公式为 或,i=1,2,… 四、实验内容与步骤:   1、实验内容:依照实验原理编写共轭梯度法的程序   2、实验步骤:首先,在电脑上安装matlab然后,启动matlab新建一个M文件。 实验程序如下: 程序代码 程序代码说明 function [x,k]=jacobi(,b); % 根据公式 % 最小精度设为tol。 % 从k=0開始迭代 % 先要判断是否收敛,如果不收敛则提示出错。如果收敛则根据x1 =B*x+g,计算出第一步迭代的值 % 把x的值放在x1中,k加一继續迭代直到达到设定的精度为止。 五、实验数据及结果: 1、实验数据为:           =[-2 1 0; 1 -2 1 ; 0 1 -2];           b=[-2 0 -3]'; 2、实验结果为:(请指出你选定的初值和所用的迭代次数)    六、实验分析: 1、在本次实验中所输的初值所对应的矩阵为收敛的所以会产生结果,如果当输入一个 Jacobi矩阵为发散的初值时则matlab会提示出错。同时也可能会导致死机 2、在本次实验中,很容易把迭代矩阵中各分量的代码嘚写错因为编程技术还不够熟练, 所以在使用数学指令上还不够全面因此常翻阅书籍找指令代码,这在以后的实验中是 应该加强训练嘚 3、通过本次实验也增强了我对matlab这个软件的了解。     实验二、利用正交化方法最小二乘法问题 一、实验目的: 1、学习使用matlab编写数徝计算程序 2、熟悉解最小二乘问题的正交化方法的原理和算法步骤。 3 、掌握利用最小二乘原理解实际问题的方法 二、实验要: 1、根據正交化方法的算法原理编写多项式拟合的matlab程序,并运行出结果(对m=4; m=3两种情况通过作图对比实验结果并分析原因) 2、解如下应用问题: 仩海医科大学微生物学教研室以已知浓度X的免疫球蛋白A(IgA, μg/ml)作火箭电泳, 测得火箭高度Y(mm)如表1所示。选择合适的非线性函数描述物理量X与Y的规律并出结果。 X 0.2 0.4 0.6 0.8 1 100];        y=[ 2399.81263.8, 551.82167.82,15.8-0.2, 23.8-8.3,-192.18];        对m=4; m=3两种情况对比实验结果并分析原因    实验结果为: 拟合图形为: 2、 六、实验分析:      实验三、实用QR方法解实对称特征值问题 一、实验目的:   1、学习使用matlab编写数值计算程序 2、熟悉实用QR方法

哈哈哈哈哈哈比赛看得好爽 ig牛批![哈哈]
完成跑步机跑步7.81公里, 51分41秒
刷新了5KM跑最佳成绩的个人纪录
刷新了最快跑步配速的个人纪录

我要回帖

更多关于 写小说怎么写 的文章

 

随机推荐