´óÁ¬Ãñ×åѧԺ
Êý ѧ ʵ Ñé ±¨ ¸æ
¿Î³Ì£º ÊýÖµ·ÖÎöÓëÈí¼þ ʵÑéÌâÄ¿: ·Ö±ðÓÃJacobi¡¢Gauss-Seidel¡¢SORµü´ú·¨Çó½âÏßÐÔ·½³Ì×é ϵ±ð£º ÀíѧԺ רҵ£º ÊýѧÓëÓ¦ÓÃÊýѧ ÐÕÃû£º ÕÔÓ¢´ä °à¼¶£º Êýѧ101 Ö¸µ¼½Ìʦ£ºÅ£´óÌï Íê³ÉѧÆÚ£º 2012 Äê 4 Ô 14 ÈÕ
1 / 10
ʵÑéÄ¿µÄ£º 1. ÊìϤmatlabÉÏ»ú»·¾³£¬ÕÆÎÕmatlabÏà¹ØÖªÊ¶£» 2. ÕÆÎÕÇó½âÏßÐÔ·½³Ì×éµÄÈýÖÖµü´ú·¨£ºJacobi¡¢Gauss-Seidel¡¢SORµü´ú·¨£» 3. ±È½ÏÈýÖÖµü´ú·¨µÄ²îÒ죻 ʵÑéÄÚÈÝ£º£¨ÎÊÌâ¡¢ÊýѧģÐÍ¡¢ÒªÇ󡢹ؼü´Ê£© 1.·Ö±ðÔËÓÃJacobi¡¢Gauss-Seidel¡¢SORµü´ú·¨Çó½âÈçÏÂÏßÐÔ·½³Ì×é ?4?100??x1??5.84????x????14?10?5.5???2?????0?14?1??x3??7.33? ???????00?14??x4??3.35? È¡³õʼµãx(0)?(0,0,0,0)T,¾«¶ÈÒªÇó??10?6£»SORµü´ú·¨ÖÐËɳÚÒò×Ó??1.05¡£ 2.Çó½âÉÏÊö·½³Ì×é׼ȷ½â£¬ÓëÉÏÊöÈýÖÖ·½·¨µÄ½á¹û½øÐбȽϡ£ ¹Ø¼ü´Ê£ºJacobi¡¢Gauss-Seidel¡¢SORµü´ú·¨ ÏßÐÔ·½³Ì×é 2 / 10
ʵÑé·½·¨ºÍ²½Ö裨°üÀ¨ÊýÖµ¹«Ê½¡¢Ëã·¨²½Öè¡¢³ÌÐò£©£º 1. ÀûÓÃJacobiµü´ú·¨Çó½â£º ±àÖÆÃûΪmajacobi.mµÄÎļþ£¬ÄÚÈÝÈçÏ£º function x=majacobi(A,b,x0,ep,N) n=length(b); if nargin>5,Warning('´«µÝµÄ²ÎÊý¸öÊý´íÎó');end if nargin<5,N=500;end if nargin<4,ep=1e-6;end if nargin<3,x0=zeros(n,1);end x=zeros(n,1);k=0; while k function x=maseidel(A,b,x0,ep,N) n=length(b); if nargin<5,N=500;end if nargin<4,ep=1e-6;end if nargin<3,x0=zeros(n,1);end x=zeros(n,1);k=0; while k ±àÖÆÃûΪmasor.mµÄÎļþ£¬ÄÚÈÝÈçÏ£º function x=masor(A,b,omega,x0,ep,N) n=length(b); if nargin<6,N=500;end if nargin<5,ep=1e-6;end if nargin<4,x0=zeros(n,1);end if nargin<3,omega=1.5;end x=zeros(n,1);k=0; while k