欢迎来到天天文库
浏览记录
ID:45621560
大小:180.86 KB
页数:3页
时间:2019-11-15
《重大数值分析作业》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验一线性方程组的数值解法1.实验目的:(1)高斯列主元消去法求解线性方程组的过程(2)熟悉用迭代法求解线性方程组的过程(3)设计出相应的算法,编制相应的函数了程序2.实验内容分别用高斯列主元消去法Jacobi迭代&,Gauss-Saidel迭代法,超松弛迭代法求解线性方程组2100-3-4-121234149-3_「10「13x25-4-2-13_^4__7_3、实验原理写出木次实验所用算法的算法步骤叙述或画出算法程序框图4.实验环境及实验文件存档名写出实验环境及实验文件存档名3、实验结果及分析输出计算结果,结果分析和小结等。4.实验结果:高斯列主元消去法■ATLM-JgJX]Fil
2、oEd:tDesktopJfindo*Holp□G»
3、®
4、静旦?CurrentBirectorx)CDocwentsandSettingsA4i.inistrator^®SW8X/^?tShortcutsKxtoAdd[fjVhaCsNewCurrentDirectory-C:DocunentxandSattings...E]er©<©a*Co5、ncO1m)ncO2mjnc036、mincO4mM-fileM-fileM-file1KB2011-10-271KB2011-10-271KB2011-10-271KB2011-10-27<12JConmandHistoryX%A^k+l:n,k:c;=A^k+l:n,k:cj-^A^k+1:n,k)/endendx=xeros(length(b)>1):%回代求解x(n)=A(n,c)/A(run);fork=n-l:-l:!x(k)=(A(k,c)-A(k,k+l:n)*x(k+l:n))/A(k,k)enddispCX=*):disp(x);fornatshort%设矍为默认格式显示,且示6位—11-10-277、下午3:19—%一11-10-27下午3:23—%zlJ^Start105-27詞1.000000000000002.000000000000003.000000000000004.00000000000000»IOVE2100-33-4-1213I23-44149-13Jacobi迭代法:Gauss-Saidel迭代法:超松弛迭代法巧IAT1ABfileEditYie*DebyDesktopWindowHelp二QX电).jk书吋目yCurrantDiractorxShortcutsjdHovtoAdd耳What*sKaw:8、cDocxmgtsandSott】:igvAdiHra9、st"torj£面需斤连文件夹刁10、©CurrentBirect«ry—C:DocsenlsandSettinILastModifiedD«zcr:pti«nKBKBKBKB2011-10-2714:47:282011-10-2715:06:522011-10-2715:10:082011.10-2715:20:30求解例題2.2Jacobi迭代法求解例3.1GaussSeidel.mSOR^f湖例3.32JC^itc11、(k)=(A(k,c)-A(k,k*l:n)sx(k*l:n))/A(k,k);endkdispCX=r);r•disp(x);forna-tshort:%设黃为默认格式显示,显示5位爲“1卜10・27下午3:19—X11-10-27下午3:23—%0・94372745740947-L986?18887395781.17497388025717-L446158637233951.0e+100♦1.08358145418013-2.281137127196451.34909702562603-1.659321327741661.0e+101♦1.24416086299551-2.619112、86149629361.54902403801066-1.905221473794931.0e+102♦1.42853704909609-3.007331740219171.77857887517126-2.187562350656691.0b*103♦1.64023653318175-3.452997869971752.04215217942518-2.511744227026121.Oe*104♦1.88330844235816-3.964
5、ncO1m)ncO2mjnc03
6、mincO4mM-fileM-fileM-file1KB2011-10-271KB2011-10-271KB2011-10-271KB2011-10-27<12JConmandHistoryX%A^k+l:n,k:c;=A^k+l:n,k:cj-^A^k+1:n,k)/endendx=xeros(length(b)>1):%回代求解x(n)=A(n,c)/A(run);fork=n-l:-l:!x(k)=(A(k,c)-A(k,k+l:n)*x(k+l:n))/A(k,k)enddispCX=*):disp(x);fornatshort%设矍为默认格式显示,且示6位—11-10-27
7、下午3:19—%一11-10-27下午3:23—%zlJ^Start105-27詞1.000000000000002.000000000000003.000000000000004.00000000000000»IOVE2100-33-4-1213I23-44149-13Jacobi迭代法:Gauss-Saidel迭代法:超松弛迭代法巧IAT1ABfileEditYie*DebyDesktopWindowHelp二QX电).jk书吋目yCurrantDiractorxShortcutsjdHovtoAdd耳What*sKaw:
8、cDocxmgtsandSott】:igvAdiHra
9、st"torj£面需斤连文件夹刁
10、©CurrentBirect«ry—C:DocsenlsandSettinILastModifiedD«zcr:pti«nKBKBKBKB2011-10-2714:47:282011-10-2715:06:522011-10-2715:10:082011.10-2715:20:30求解例題2.2Jacobi迭代法求解例3.1GaussSeidel.mSOR^f湖例3.32JC^itc11、(k)=(A(k,c)-A(k,k*l:n)sx(k*l:n))/A(k,k);endkdispCX=r);r•disp(x);forna-tshort:%设黃为默认格式显示,显示5位爲“1卜10・27下午3:19—X11-10-27下午3:23—%0・94372745740947-L986?18887395781.17497388025717-L446158637233951.0e+100♦1.08358145418013-2.281137127196451.34909702562603-1.659321327741661.0e+101♦1.24416086299551-2.619112、86149629361.54902403801066-1.905221473794931.0e+102♦1.42853704909609-3.007331740219171.77857887517126-2.187562350656691.0b*103♦1.64023653318175-3.452997869971752.04215217942518-2.511744227026121.Oe*104♦1.88330844235816-3.964
11、(k)=(A(k,c)-A(k,k*l:n)sx(k*l:n))/A(k,k);endkdispCX=r);r•disp(x);forna-tshort:%设黃为默认格式显示,显示5位爲“1卜10・27下午3:19—X11-10-27下午3:23—%0・94372745740947-L986?18887395781.17497388025717-L446158637233951.0e+100♦1.08358145418013-2.281137127196451.34909702562603-1.659321327741661.0e+101♦1.24416086299551-2.6191
12、86149629361.54902403801066-1.905221473794931.0e+102♦1.42853704909609-3.007331740219171.77857887517126-2.187562350656691.0b*103♦1.64023653318175-3.452997869971752.04215217942518-2.511744227026121.Oe*104♦1.88330844235816-3.964
此文档下载收益归作者所有