<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta charset="utf-8">
<meta name="generator" content="pdf2htmlEX">
<meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1">
<link rel="stylesheet" href="https://static.pudn.com/base/css/base.min.css">
<link rel="stylesheet" href="https://static.pudn.com/base/css/fancy.min.css">
<link rel="stylesheet" href="https://static.pudn.com/prod/directory_preview_static/626bfabd1e41a87e8aaa260f/raw.css">
<script src="https://static.pudn.com/base/js/compatibility.min.js"></script>
<script src="https://static.pudn.com/base/js/pdf2htmlEX.min.js"></script>
<script>
try{
pdf2htmlEX.defaultViewer = new pdf2htmlEX.Viewer({});
}catch(e){}
</script>
<title></title>
</head>
<body>
<div id="sidebar" style="display: none">
<div id="outline">
</div>
</div>
<div id="pf1" class="pf w0 h0" data-page-no="1"><div class="pc pc1 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/626bfabd1e41a87e8aaa260f/bg1.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0"> 2019-2020 <span class="ff2">第<span class="_ _0"> </span></span>1<span class="_ _0"> </span><span class="ff2">学期数值分析上机实习</span></div><div class="t m0 x2 h4 y3 ff2 fs1 fc0 sc1 ls0 ws0">学生姓名： <span class="_ _1"></span> <span class="_ _2"> </span> <span class="sc0">周智<span class="_ _3"></span><span class="sc1"> <span class="_ _4"> </span> <span class="_ _5"></span> <span class="_ _6"> </span> </span></span></div><div class="t m0 x2 h4 y4 ff2 fs1 fc0 sc1 ls0 ws0">学生学号： <span class="_ _7"></span><span class="sc0"> <span class="_ _8"> </span> <span class="ff3">2019201523<span class="_ _9"></span><span class="ff2"> <span class="_ _a"> </span> <span class="ff4"> <span class="_ _b"></span><span class="ff2"> <span class="_ _c"> </span> </span></span></span></span></span></div><div class="t m0 x2 h4 y5 ff2 fs1 fc0 sc1 ls0 ws0">年级专业： <span class="_ _d"></span><span class="sc0"> <span class="_ _d"></span> 牵引动力国家重点实验室<span class="_ _e"></span> <span class="_ _f"> </span> </span></div><div class="t m0 x2 h4 y6 ff2 fs1 fc0 sc1 ls0 ws0">联系电话： <span class="_ _10"></span><span class="sc0"> <span class="_ _11"> </span> <span class="_ _12"></span> <span class="ff3">18402894842<span class="_ _13"></span><span class="ff2"> <span class="_ _14"> </span> <span class="ff4"> <span class="_ _15"></span><span class="ff2"> <span class="_ _16"> </span> </span></span></span></span></span></div><div class="t m0 x3 h5 y7 ff4 fs1 fc0 sc0 ls0 ws0"> </div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>
</body>
</html>
<div id="pf2" class="pf w0 h0" data-page-no="2"><div class="pc pc2 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/626bfabd1e41a87e8aaa260f/bg2.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x3 h3 y8 ff2 fs0 fc0 sc1 ls0 ws0">目录</div><div class="t m0 x3 h6 y9 ff3 fs2 fc0 sc0 ls0 ws0">2019-2020 <span class="ff2">第<span class="_ _17"> </span></span>1<span class="_ _17"> </span><span class="ff2">学期数值分析上机实习</span>.......................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _19"></span>1</div><div class="t m0 x3 h6 ya ff2 fs2 fc0 sc0 ls0 ws0">一、第一题<span class="ff3">......................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>.<span class="_ _1a"> </span>1</span></div><div class="t m0 x4 h6 yb ff3 fs2 fc0 sc0 ls0 ws0">1.1 <span class="ff2">程序代码</span>................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _1a"> </span>1</div><div class="t m0 x5 h6 yc ff3 fs2 fc0 sc0 ls0 ws0">1.1.1Jacobi<span class="_ _17"> </span><span class="ff2">迭代通用程序</span>..........................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>..<span class="_ _1b"></span>1</div><div class="t m0 x5 h6 yd ff3 fs2 fc0 sc0 ls0 ws0">1.1.2 Gauss-Seidel<span class="_ _17"> </span><span class="ff2">迭代通用程序</span>.............................................<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _19"></span>1</div><div class="t m0 x4 h6 ye ff3 fs2 fc0 sc0 ls0 ws0">1.2<span class="_ _17"> </span><span class="ff2">计算结果</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>2</div><div class="t m0 x5 h6 yf ff3 fs2 fc0 sc0 ls0 ws0">1.2.1 Jacobi<span class="_ _17"> </span><span class="ff2">迭代结果</span>.....................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _1b"></span>2</div><div class="t m0 x5 h6 y10 ff3 fs2 fc0 sc0 ls0 ws0">1.2.2 Gauss-Seidel<span class="_ _17"> </span><span class="ff2">迭代结果</span>.........................................<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _1b"></span>3</div><div class="t m0 x3 h6 y11 ff2 fs2 fc0 sc0 ls0 ws0">二、第二题<span class="ff3">......................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>.<span class="_ _1a"> </span>3</span></div><div class="t m0 x4 h6 y12 ff3 fs2 fc0 sc0 ls0 ws0">2.1<span class="_ _17"> </span><span class="ff2">程序代码</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>3</div><div class="t m0 x5 h6 y13 ff3 fs2 fc0 sc0 ls0 ws0">2.1.1 Doolittle<span class="_ _17"> </span><span class="ff2">分解程序</span>..............................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _1b"></span>3</div><div class="t m0 x4 h6 y14 ff3 fs2 fc0 sc0 ls0 ws0">2.2<span class="_ _17"> </span><span class="ff2">计算结果</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>4</div><div class="t m0 x5 h6 y15 ff3 fs2 fc0 sc0 ls0 ws0">2.2.1 <span class="ff2">第一问</span>.............................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _1b"> </span>4</div><div class="t m0 x5 h6 y16 ff3 fs2 fc0 sc0 ls0 ws0">2.2.2 <span class="ff2">第二问</span>.............................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _1b"> </span>4</div><div class="t m0 x5 h6 y17 ff3 fs2 fc0 sc0 ls0 ws0">2.2.3 <span class="ff2">第三问</span>.............................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _1b"> </span>5</div><div class="t m0 x3 h6 y18 ff2 fs2 fc0 sc0 ls0 ws0">三、第三题<span class="ff3">......................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>.<span class="_ _1a"> </span>6</span></div><div class="t m0 x4 h6 y19 ff3 fs2 fc0 sc0 ls0 ws0">3.1 <span class="ff2">程序代码</span>................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _1b"> </span>6</div><div class="t m0 x5 h6 y1a ff3 fs2 fc0 sc0 ls0 ws0">3.1.1<span class="_ _17"> </span><span class="ff2">多项式插值代码</span>.........................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _1b"> </span>6</div><div class="t m0 x5 h6 y1b ff3 fs2 fc0 sc0 ls0 ws0">3.1.2 <span class="ff2">五次多项式拟合代码</span>.........................................................................<span class="_ _1b"></span>7</div><div class="t m0 x4 h6 y1c ff3 fs2 fc0 sc0 ls0 ws0">3.2<span class="_ _17"> </span><span class="ff2">计算结果</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>8</div><div class="t m0 x4 h6 y1d ff3 fs2 fc0 sc0 ls0 ws0">3.3<span class="_ _17"> </span><span class="ff2">结果分析</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>8</div><div class="t m0 x3 h6 y1e ff2 fs2 fc0 sc0 ls0 ws0">四、第四题<span class="ff3">......................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>.<span class="_ _1a"> </span>8</span></div><div class="t m0 x4 h6 y1f ff3 fs2 fc0 sc0 ls0 ws0">4.1<span class="_ _17"> </span><span class="ff2">程序代码</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>8</div><div class="t m0 x5 h6 y20 ff3 fs2 fc0 sc0 ls0 ws0">4.1.1<span class="_ _17"> </span><span class="ff2">非线性方程迭代程序代码</span>.................................................................<span class="_ _18"></span>..<span class="_ _19"></span>8</div><div class="t m0 x4 h6 y21 ff3 fs2 fc0 sc0 ls0 ws0">4.2<span class="_ _17"> </span><span class="ff2">计算结果</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>9</div><div class="t m0 x3 h6 y22 ff2 fs2 fc0 sc0 ls0 ws0">五、第五题<span class="ff3">......................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>.<span class="_ _1a"> </span>9</span></div><div class="t m0 x4 h6 y23 ff3 fs2 fc0 sc0 ls0 ws0">5.1<span class="_ _17"> </span><span class="ff2">程序代码</span>........................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>9</div><div class="t m0 x5 h6 y24 ff3 fs2 fc0 sc0 ls0 ws0">5.1.1 Romberg <span class="ff2">算法代码</span>...........................................................<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _19"></span>9</div><div class="t m0 x5 h6 y25 ff3 fs2 fc0 sc0 ls0 ws0">5.1.2 <span class="ff2">高斯型求积公式程序代码</span>.......................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _1b"></span>10</div></div><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>
<div id="pf3" class="pf w0 h0" data-page-no="3"><div class="pc pc3 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/626bfabd1e41a87e8aaa260f/bg3.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x4 h6 y26 ff3 fs2 fc0 sc0 ls0 ws0">5.2<span class="_ _17"> </span><span class="ff2">计算结果</span>......................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>10</div><div class="t m0 x3 h6 y27 ff2 fs2 fc0 sc0 ls0 ws0">六、第六题<span class="ff3">....................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>.<span class="_ _1b"> </span>11</span></div><div class="t m0 x4 h6 y28 ff3 fs2 fc0 sc0 ls0 ws0">6.1<span class="_ _17"> </span><span class="ff2">程序代码</span>......................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>11</div><div class="t m0 x5 h6 y29 ff3 fs2 fc0 sc0 ls0 ws0">6.1.1<span class="_ _17"> </span><span class="ff2">四阶<span class="_ _17"> </span></span>Runge-Kutta<span class="_ _17"> </span><span class="ff2">法程序代码</span>................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _19"></span>11</div><div class="t m0 x5 h6 y2a ff3 fs2 fc0 sc0 ls0 ws0">6.1.2<span class="_ _17"> </span><span class="ff2">显示<span class="_ _17"> </span></span>4<span class="_ _17"> </span><span class="ff2">阶<span class="_ _1c"> </span></span>Adams<span class="_ _1c"> </span><span class="ff2">程序代码</span>...............................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _19"></span>12</div><div class="t m0 x4 h6 y2b ff3 fs2 fc0 sc0 ls0 ws0">6.2<span class="_ _1c"> </span><span class="ff2">计算结果</span>...................................................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>.<span class="_ _1b"></span>12</div><div class="t m0 x4 h6 y2c ff3 fs2 fc0 sc0 ls0 ws0">6.3 <span class="ff2">结果分析</span>..............................................................<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>....<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>...<span class="_ _18"></span>..<span class="_ _1b"> </span>13</div></div><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a><a class="l" rel='nofollow' onclick='return false;'><div class="d m1"></div></a></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>
<div id="pf4" class="pf w0 h0" data-page-no="4"><div class="pc pc4 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/626bfabd1e41a87e8aaa260f/bg4.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x6 h3 y8 ff4 fs0 fc0 sc0 ls0 ws0">2019-2020 <span class="ff2 sc1">第<span class="_ _0"> </span></span>1<span class="_ _0"> </span><span class="ff2 sc1">学期数值分析上机实习</span></div><div class="t m0 x3 h7 y2d ff2 fs3 fc0 sc0 ls0 ws0">一、第一题</div><div class="t m0 x3 h4 y2e ff3 fs1 fc0 sc0 ls0 ws0">1.1 <span class="ff2">程序代码</span></div><div class="t m0 x3 h8 y2f ff3 fs4 fc0 sc0 ls0 ws0">1.1.1Jacobi<span class="_ _17"> </span><span class="ff2">迭代通用程序</span></div><div class="t m0 x7 h9 y30 ff1 fs5 fc0 sc0 ls0 ws0">function [x] = Jacobi(A,b,x0,epsilon)</div><div class="t m0 x7 ha y31 ff1 fs5 fc0 sc0 ls0 ws0"> %A<span class="_ _1c"> </span><span class="ff2">为系数矩阵；</span>b<span class="_ _1c"> </span><span class="ff2">为等号右端列向量；</span>x0<span class="_ _17"> </span><span class="ff2">为迭代初始向量；</span>epsilon<span class="_ _1c"> </span><span class="ff2">为迭代误差</span></div><div class="t m0 x7 ha y32 ff1 fs5 fc0 sc0 ls0 ws0">Dn=diag(1./diag(A)); %<span class="ff2">原来主对角元素倒数组成的对角矩阵，也就是<span class="_ _1c"> </span></span>D^(-1)</div><div class="t m0 x7 ha y33 ff1 fs5 fc0 sc0 ls0 ws0">B=eye(size(A))-Dn*A; %B=I-D^(-1)*A<span class="_ _1d"></span> ,B<span class="_ _1c"> </span><span class="ff2">是迭代用的矩阵</span></div><div class="t m0 x7 ha y34 ff1 fs5 fc0 sc0 ls0 ws0">g=Dn*b; %g=D^(-1)*b ,<span class="ff2">是迭代用的向量</span></div><div class="t m0 x7 h9 y35 ff1 fs5 fc0 sc0 ls0 ws0">r=max(abs(eig(B)));</div><div class="t m0 x7 h9 y36 ff1 fs5 fc0 sc0 ls0 ws0">while(1)</div><div class="t m0 x7 h9 y37 ff1 fs5 fc0 sc0 ls0 ws0"> if (r>=1)</div><div class="t m0 x7 ha y38 ff1 fs5 fc0 sc0 ls0 ws0"> disp('<span class="ff2">该<span class="_ _1c"> </span></span>Jacobi<span class="_ _17"> </span><span class="ff2">迭代不收敛</span>,<span class="ff2">请尝试<span class="_ _1c"> </span></span>Gauss-Seidel<span class="_ _1c"> </span><span class="ff2">迭代法</span>');</div><div class="t m0 x7 h9 y39 ff1 fs5 fc0 sc0 ls0 ws0"> x=[];</div><div class="t m0 x7 h9 y3a ff1 fs5 fc0 sc0 ls0 ws0"> break</div><div class="t m0 x7 h9 y3b ff1 fs5 fc0 sc0 ls0 ws0"> end</div><div class="t m0 x7 ha y3c ff1 fs5 fc0 sc0 ls0 ws0"> x=B*x0+g; %<span class="ff2">根据上一个值<span class="_ _1c"> </span></span>x<span class="_ _1c"> </span><span class="ff2">计算下一个值<span class="_ _17"> </span></span>x </div><div class="t m0 x7 ha y3d ff1 fs5 fc0 sc0 ls0 ws0"> R=norm(x-x0,2); %<span class="ff2">计算相邻两次迭代结果的二范数</span></div><div class="t m0 x7 ha y3e ff1 fs5 fc0 sc0 ls0 ws0"> if (R<epsilon) %<span class="ff2">判断是否满足误差要求</span></div><div class="t m0 x7 h9 y3f ff1 fs5 fc0 sc0 ls0 ws0"> break</div><div class="t m0 x7 h9 y40 ff1 fs5 fc0 sc0 ls0 ws0"> end</div><div class="t m0 x7 h9 y41 ff1 fs5 fc0 sc0 ls0 ws0"> x0=x;</div><div class="t m0 x7 h9 y42 ff1 fs5 fc0 sc0 ls0 ws0">end</div><div class="t m0 x7 h9 y43 ff1 fs5 fc0 sc0 ls0 ws0">end</div><div class="t m0 x3 h8 y44 ff3 fs4 fc0 sc0 ls0 ws0">1.1.2 Gauss-Seidel<span class="_ _1e"> </span><span class="ff2">迭代通用程序</span></div><div class="t m0 x7 h9 y45 ff1 fs5 fc0 sc0 ls0 ws0">function [x]=GaussSeidel(A,b,x0,epsilon)</div><div class="t m0 x7 ha y46 ff1 fs5 fc0 sc0 ls0 ws0">%A<span class="_ _1c"> </span><span class="ff2">为系数矩阵；</span>b<span class="_ _1c"> </span><span class="ff2">为等号右端列向量；</span>x0<span class="_ _17"> </span><span class="ff2">为迭代初始向量；</span>epsilon<span class="_ _1c"> </span><span class="ff2">为迭代误差</span></div><div class="t m0 x7 ha y47 ff1 fs5 fc0 sc0 ls0 ws0">D=diag(diag(A));U=-triu(A,1);L=-tril(A,-1);%<span class="ff2">分别计算<span class="_ _17"> </span></span>D<span class="_ _1c"> </span><span class="ff2">矩阵，</span>L<span class="_ _1c"> </span><span class="ff2">矩阵和<span class="_ _17"> </span></span>U<span class="_ _1c"> </span><span class="ff2">矩阵</span></div><div class="t m0 x7 ha y48 ff1 fs5 fc0 sc0 ls0 ws0">B=(D-L)\U;%<span class="ff2">计算迭代矩阵（</span>D-L<span class="ff2">）</span>^-1*U</div><div class="t m0 x7 ha y49 ff1 fs5 fc0 sc0 ls0 ws0">g=(D-L)\b;%<span class="ff2">计算（</span>D-L<span class="ff2">）</span>^-1*b</div><div class="t m0 x7 h9 y4a ff1 fs5 fc0 sc0 ls0 ws0">r=max(abs(eig(B)));</div><div class="t m0 x7 h9 y4b ff1 fs5 fc0 sc0 ls0 ws0">while(1)</div><div class="t m0 x7 h9 y4c ff1 fs5 fc0 sc0 ls0 ws0"> if (r>=1)</div><div class="t m0 x7 ha y4d ff1 fs5 fc0 sc0 ls0 ws0"> disp('<span class="ff2">该<span class="_ _1c"> </span></span>Gauss-Seidel<span class="_ _17"> </span><span class="ff2">迭代不收敛</span>,<span class="ff2">请尝试<span class="_ _1c"> </span></span>Jacobi<span class="_ _1c"> </span><span class="ff2">迭代法</span>');</div><div class="t m0 x8 hb y4e ff3 fs6 fc0 sc0 ls0 ws0">1</div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>
