å½
线æ§æ¹ç¨ç»çè§æ¨¡æ¯è¾å¤§æ¶ï¼éç¨
é«æ¯æ¶å
æ³éè¦å¤ªå¤æ¶é´ãè¿æ¶å°±è¦éç¨
è¿ä»£æ³æ±è§£æ¹ç¨ç»äºã
é«æ¯æ¶å
æ³æ¯ä¸ä¸ªO(n^3)ç
æµ®ç¹è¿ç®çæéåºåï¼å¨ç»è¿æéæ¥è®¡ç®ä¹åç论ä¸å¾å°çæ¯ç²¾ç¡®è§£(æ èå
¥è¯¯å·®æ¶)ã
èè¿ä»£æ³å¨ç»è¿æéæ¥è¿ä»£ä¹åä¸è¬ä¸äº§ç精确解ï¼è¿ä»£æ³å¨è®¡ç®è¿ç¨ä¸éæ¸åå°è¯¯å·®ï¼å½è¯¯å·®å°äºå®¹è®¸å¼æ¶åæ¢è¿ä»£è®¡ç®ãæ¹ç¨ç»çç³»æ°ç©éµæ¯ä¸¥æ ¼å¯¹è§å ä¼ç©éµæ¶ï¼è¿ä»£æ»æ¯æ¶æçã
âJacobiè¿ä»£æ³
对äºæ¹ç¨ç»3u+v=5ï¼u+2v=5ï¼å°å
¶æ¹å为å¦ä¸çå½¢å¼
ç±äºæ¹ç¨ç»çç³»æ°ç©éµæ¯ä¸¥æ ¼å¯¹è§å ä¼ç©éµæ¶ï¼è¿ä»£ä¸å®æ¶æã使ç¨åå¼[u0ï¼v0]=[0ï¼0]å¼å§è¿ä»£ï¼ä»¥ä¸æ¯è¿ä»£è¿ç¨ï¼
继ç»è¿ä»£è¿ç¨æç»ä¼æ¶æå°è§£[1ï¼2].è¿ä¸ªè¿ä»£è¿ç¨å°±æ¯Jacobiè¿ä»£ã
对äºæ¹ç¨ç»u+2v=5ï¼3u+v=5ï¼ç±äºæ¹ç¨ç»çç³»æ°ç©éµä¸æ¯ä¸¥æ ¼å¯¹è§å ä¼ç©éµæ¶ï¼å æ¤è¿ä»£ä¸æ¶æãæ¥çè¿ä»£è¿ç¨ï¼
设D表示系æ°ç©éµA ç主对è§é¨åï¼L表示Aç主对è§çº¿ä¸æ¹é¨åï¼U表示Aç主对è§çº¿ä¸æ¹é¨åãåA=D+L+Uï¼AX=bå¯æ¹å为
对äºä¸é¢çæ¹ç¨ç»3u+v=5ï¼u+2v=5ï¼åæç©éµå½¢å¼
è¿ä»£æ ¼å¼ä¸º
è¿ä¸ä¹åçè¿ä»£æ ¼å¼æ¯ä¸è´çã
Fortranæºä»£ç