izebwhiletNY1hBY0if
ormY1Y0i
fespfpri
tf满足误差限
breake
dY0Y1fpri
tf第2d次迭代得：tdispY1tt1e
dfpri
tf满足误差限的雅可比迭代近似解为：dispY1
f附录四Gausssidel迭代程序代码
ai
put请输入系数矩阵a：bi
put请输入矩阵bNi
put请输入最大迭代次数N：espi
put请输入近似解的误差限：ifa
ydiaga0error系数矩阵错误，迭代终止！e
dDdiagdiagaX0zerossizebx10x20x30X1x1x2x3hi
vDbBi
vDDaB1triuBB2trilBk1fpri
tf高斯赛德尔迭代法
fpri
tf第0次迭代得：dispX1whilekNx1h11B11X0X1x1x2x3x2h21B12X0B22X1X1x1x2x3x3h31B23X1X1x1x2x3if
ormX1X0i
fespfpri
tf已满足误差限。breake
dX0X1fpri
tf第2d次迭代得：kdispX1kk1e
dfpri
tf满足误差限的高斯赛德尔迭代近似解为：dispX1
f附录五四阶Ru
geKutta算法程序代码
h00001000100101k1forxk01011取不同的x进行计算fori114取不同的步长计算每个xhjhiyk1forxkphjhjxk四阶标准Ru
geKutta迭代算法k1hj20ykformatlo
gk1k2hj20ykk105k3hj20ykk205k4hj20ykk3ykmykk12k22k3k4016666666666667ykykme
dykykmyealexp20xkp求出精确值ckabsykyeal精确值和计算值之间的绝对误差ykealkyealckkkckyxkkykxkkkxkkk1e
de
dDyxkykealckk将迭代得到的结果及误差存到矩阵D中fori1110将同一个步长下不同X得到的结果集中到一个矩阵forj113Ah1ijD4i3jAh2ijD4i2jAh3ijD4i1jAh4ijD4ije
dAxkixkk4i3e
dCh1AxkAh1Ch2AxkAh2Ch3AxkAh3Ch4AxkAh4Che
gCh11234Ch21234Ch31234Ch41234
f总结与体会
本次上机报告选择的是14题目中的第二题和第三题，56题中的第五题，在实习过程中，通过网上查阅资料，对程序代码进行了研究运行，因为对编程接触的机会太少，在编写过程中也是显得力不从心，也出现了不少问题，同时也很是纠结，最后通过向同学和师兄请教，解决了一些操作问题，通过本次实习，也对C语言及MATLAB进行了一次系统的复习，也有了更深一步的了解，但需要学习的还有很多，需要再接再厉，在克服了种种困难后，收获也是很多的。总之，通过本次上机实习，学到了数值分析的相关知识，了解了编程运算的重要性，同时为自己专业研究也奠定了基础，编程是必不可少，通过理论与实践相结合的方法，不仅提高了计算与编程能力，还对将来的工作和进一步的学习打好基础。
fr