雅可比法实验报告matlab,1200710222李卫康数值分析试验一.doc

1200710222李卫康数值分析试验一.doc

附录Ⅰ实验报告的一般格式桂林电子科技大学数学与计算科学学院实验报告实验室实验日期年月日院(系)数学与计算科学学院年级、专业、班1200710222姓名李卫康成绩课程名称数值分析实验项目名称指导教师一，实验目的通过实验加深理解线性方程组的雅可比迭代法，高斯赛德尔迭代法和SOR迭代法的算法思想，掌握三个迭代法MATLAB程序的编写和使用。二，实验原理MATLAB中有三个内置函数可以提取矩阵中的部分元素，DDIAGA当A为方阵时，输出D为向量，它由A的对角线元素构成；当A为向量时，输出D为对角阵，它以A为对角线元素。LTRILA,K输入A为方阵，K为对角线序号，主对角线K0，向右上移动依次K递增1，向左下移动K依次递减1，输出L为第K条对角线及其以下元素构成的三角阵。UTRIUA,K输入A为方阵，K为对角线序号，输出U为第K条对角线及其以上元素构成的三角阵。根据线性方程组的雅可比迭代法的迭代公式111KKKXDBLXDLUXB编写MATLAB程序FUNCTIONX,KJACOBIA,B,X0,E,N输入A为系数矩阵,B为常数的列向量,X0为迭代初始向量,N为最大迭代次数输出为解X和迭代次数KATLONGDDIAGDIAGAL1TRILA,1U1TRIUA,1K0XD\LUX0D\BWHILENORMXX0EX0XXD\LUX0D\BKK1DISPK,X1,X2,X3IFKNBREAKENDEND三，使用仪器，材料MATLAB编译器四，实验内容与步骤1编写高斯赛德尔迭代法和SOR迭代法的MATLAB程序。A)FUNCTIONX,KGAOSISAIDEA,B,X0,E,N输入A为系数矩阵,B为常数的列向量,X0为迭代初始向量,N为最大迭代次数输出为解X和迭代次数KATLONGDDIAGDIAGAL1TRILA,1U1TRIUA,1K0XDL\UX0DL\BWHILENORMXX0EX0XXDL\UX0DL\BKK1DISPK,X1,X2,X3IFKNBREAKENDENDA20231812315B241230 X0000 X,KGAOSISAIDEA,B,X0,000001,100100000000000000074850000000000114268750000000212873750000000200000000000000076642062500000113810523437500212543163020833300000000000000076737473203125113839920472005212536321000651400000000000000076735559802702113841014899581212536795006223500000000000000076735379259109113840978216834212536811742152X076735379259109113840978216834212536811742152K5BFUNCTIONX,KSORA,B,X0,E,N输入A为系数矩阵,B为常数的列向量,X0为迭代初始向量,N为最大迭代次数输出为解X和迭代次数KATLONGW15DDIAGDIAGAL1TRILA,1U1TRIUA,1K0XDWL\1WDWUX0WBWHILENORMXX0EX0XXDWL\1WDWUX0WBKK1DISPK,X1,X2,X3IFKNBREAKENDENDA20231812315B241230 X0000 X,KSORA,B,X0,000001,10010000000000000001099687500000007117910156250016286560546875020000000000000013817681103515613296499612426823082133389587430000000000000003903204493721010791999340666320915892208662040000000000000009723522105090611359110145838221005082518402450000000000000006708228858938411624199043829521643072682160060000000000000008032564360470411123718532592120909066346603670000000000000007610620111890211590699522736421500552661141180000000000000007618460666887711244835239626421099482107939890000000000000007756660906325711467055565012021339043434268510000000000000000760032643937491134034036618592121251510484651100000000000000077259698267942114038638922244212697076498852120000000000000007640751281545111377357004244621250203020021813000000000000000769172512908591138470996992002125196645514801400000000000000076647384875608113857628382821212567979263985150000000000000007676986797037611382033945214721250813860957616000000000000000767276839098351138581135515372125578279787061700000000000000076731929717143113829113656260212523434164096180000000000000007674189540606811384819387746421254436199997819000000000000000767293417653511138370836052692125330757285212000000000000000076739824537617113842789397465212538334047455210000000000000007673254416089511384031563720221253641883525422000000000000000767369863360401138410787117812125367169286992300000000000000076734583716255113841091773180212537052324353240000000000000007673570760291611384081162704721253657580535425000000000000000767352948982