å·®åæ¹ç¨æ¯å¾®åæ¹ç¨ç离æ£åãä¸ä¸ªå¾®åæ¹ç¨ä¸ä¸å®å¯ä»¥è§£åºç²¾ç¡®ç解ï¼æå®åæå·®åæ¹ç¨ï¼å°±å¯ä»¥æ±åºè¿ä¼¼ç解æ¥ã
æ¯å¦ dy+y*dx=0,y(0)=1 æ¯ä¸ä¸ªå¾®åæ¹ç¨ï¼ xåå¼[0,1]
(注ï¼è§£ä¸ºy(x)=e^(-x));
è¦å®ç°å¾®åæ¹ç¨ç离æ£åï¼å¯ä»¥æxçåºé´åå²ä¸ºè®¸å¤å°åºé´ [0,1/n],[1/n,2/n],...[(n-1)/n,1]
è¿æ ·ä¸è¿°å¾®åæ¹ç¨å¯ä»¥ç¦»æ£å为ï¼
å·®åæ¹ç¨
y((k+1)/n)-y(k/n)+y(k/n)*(1/n)=0ï¼ k=0,1,2,...,n-1 ï¼n 个离æ£æ¹ç¨ç»ï¼
å©ç¨y(0)=1çæ¡ä»¶ï¼ä»¥åä¸é¢çå·®åæ¹ç¨ï¼å°±å¯ä»¥è®¡ç®åº y(k/n) çè¿ä¼¼å¼äºã
è¿æ¯ä¸å¤ªæå~~
å°±æ¯æå¾®åæ¹ç¨æ¹æ两个ä¸è¥¿çå·®çä¸ç§æ¹ç¨
å³dy/dt=f(t)å¨[0,1]
离æ£æn个ç¹ï¼æ¯å¦11个ç¹
y(0),y(1/10),y(2/10),...,y(1)
ç¶å令
[y(1/10)-y(0)]/(1/10-0)=f(0) ï¼çå¼å·¦è¾¹å«å·®åï¼æ¯å¯¼æ°dy/dtçé¼è¿å¼ï¼
[y(2/10)-y(1/10)]/(2/10-1/10)=f(1/10)
追é®ããããæ°å¦ä¸å¥½ãããã
追çåµåµ