数值分析C语言超长编程,数值分析编程题c语言.doc

数值分析编程题c语言.doc

数值分析实习报告1数值分析实习报告2上机实习题一一、题目已知A与B1238412,2115237,1061074,1112336,0113584,0718719,1742382,3067813,20317432115237,19141823,3125432,1012345,2189736,1563849,0784165,1112348,31231241061074,3125432,15567914,3123848,2031454,1836742,1056781,0336993,10101031112336,1012345,3123848,27108437,4101011,3741856,2101023,071828,0037585A0113584,2189736,2031454,4101011,19897918,0431637,3111223,2121314,17841370718719,1563849,1836742,3741856,0431637,9789365,0103458,1103456,02384171742382,0784165,1056781,2101023,3111223,0103458,147138465,3123789,22134743067813,1112348,0336993,071828,2121314,1103456,3123789,30719334,44467822031743,3123124,1010103,0037585,1784317,0238417,2213474,4446782,4000001B{21874369,33992318,25173417,084671695,1784317,86612343,11101230,4719345,56784392}1．用HOUSEHOLD变换，把A化为三对角阵B(并打印B)。2．用超松弛法求解BXB(取松弛因子Ω14,X00,迭代9次)。3．用列主元素消去法求解BXB。二、解题方法的理论依据1、用HOUSEHOLDER变换的理论依据﹝1﹞令A0A,AIJ1AIJ,已知AR_1即AR_1AIJR﹝2﹞SRSQRTPOWA,2﹝3﹞ARSRSRABSAR1,RSR﹝4﹞YRAR_1U®/A®﹝5﹞KR/2UR的转置YR/A®﹝6﹞QRYRKRUR﹝7﹞ARAR1QRUR的转置URQR的转置R1,2,,,N22、用超松弛法求解其基本思想在高斯方法已求出XM,XM1的基础上，组合新的序列，从而加快收敛速度。其算式0101IXWIIBBIBIBIIX其中Ω是超松弛因子，当Ω1时，可以加快收敛速度3、用消去法求解用追赶消去法求BXB的方法,,1IBID11IBIA1IIXIIXIII数值分析实习报告3,,1IBIB11IBICQ100,U100,8,2,IQIABICIQ9,21,11}1IQABUDUX9U19,78,IIXQI三、1．计算程序INCLUDE“MATHH“INCLUDE“STDIOH“DEFINEGE8VOIDMAIN{INTSIGNDOUBLEXDOUBLEA9{{1238412,2115237,1061074,1112336,0113584,0718719,1742382,3067813,2031743},{2115237,19141823,3125432,1012345,2189736,1563849,0784165,1112348,3123124},{1061074,3125432,15567914,3123848,2031454,1836742,1056781,0336993,1010103},{1112336,1012345,3123848,27108437,4101011,3741856,2101023,071828,0037585},{0113584,2189736,2031454,4101011,19897918,0431637,3111223,2121314,1784317},{0718719,1563849,1836742,3741856,0431637,9789365,0103458,1103456,0238417},{1742382,0784165,1056781,2101023,3111223,0103458,14713846,3123789,2213474},{3067813,1112348,0336993,071828,2121314,1103456,3123789,30719334,4446782},{2031743,3123124,1010103,0037585,1784317,0238417,2213474,4446782,4000001}}DOUBLEK,H,S,WINTI,J,N,M,GDOUBLEU9,X19,Y9,Q9,B1910,X9DOUBLEB9{21874369,33992318,25173417,084671695,1784317,86612343,11101230,4719345,56784392}FORJ0J0SSSAJ1JSSSAJ1JFORG0G0IXIYIAII1XI1/UIFORI0I1E611RETURNZ}2．打印结果HOUSEHOLD1238412489308000000000000000000000000000000000000000000数值分析实习报告64893082539842649410000000000000000000000000000000000000000000649410206115082439300000000000000000000000000000000000000000082439323422841388007000000000000000000000000000000000000000000138800729698284534500000000000000000000000000000000000000000004534501600612488144000000000000000000000000000000000000000000488144260133245036300000000000000000000000000000000000000000045036321254064504500000000000000000000000000000000000000000004504501453412X01073409X12272579X22856601X32292514X42112165X56422586X61357802X70634259X80587042X01075799X12275744X22855515X32293099X42112634X56423838X61357923X70634244X80587266四、问题讨论此程序具有很好的通用性。在GS方法的基础上，已经求出X的第M解，第M1基础上，经过重新组合得新的序列，而在此新序列收敛速度加快。上机实习题二一、题目已知函数值如下表X12345FX00693147810