matlabä¸å次æ°çææ ¼ææ¥æå¼ï¼y=1/(1+25x^2)ï¼å¾åæä¹æ¾å¨ä¸ä¸ªåæ ç³»ä¸ï¼è¿ä¸ªé®é¢å¯ä»¥è¿æ ·å¤çï¼plotï¼ï¼å½æ°åé¢ï¼å ä¸ä¸ªhold on å½ä»¤ã
é£ææ ·æ¥å®ç°ï¼å¯ä»¥ç¨ä¸å代ç æ¥å®æï¼
1ãé¦å
ç¡®å®xçåºé´ï¼å¦ xå¨ã-1ï¼1ãï¼åx=[-1:0.1:1]åéæ°æ®
2ãå©ç¨ææ ¼ææ¥æå¼å
¬å¼ï¼è¿è¡æå¼è®¡ç®ï¼å³
w_x_diff(j)=subs(w_x_diff(n+1),x_para,x(j));
ln_x(j)=(1/(1+x(j)^2))*w_x(n+1) / ((x_para-x(j))*(w_x_diff(j)));
3ãå©ç¨å¾ªç¯è¯å¥ï¼è®¡ç®n=5,10,20次çæå¼
for n=[5,10,20]
ããã
w_x(n+1)=prod(part_x);
w_x_diff(n+1)=diff(w_x(n+1));
ããã
Ln_x=sum(ln_x);
x_vector=-0.99:0.037:0.99;
Ln_x=subs(Ln_x,x_para,x_vector);
ããã
end
4ãå©ç¨plotå½æ°ï¼ç»å¶å次æ°çæå¼å¤é¡¹å¼çæ²çº¿å¾
plot(x0,y),hold on
plot(x_vector,Ln_x),hold on
5ãå®å代ç ï¼å¹¶è¿è¡å¾å°n=5,10,20次æå¼å¤é¡¹å¼çæ²çº¿å¾
6ãä»å¾åä¸ï¼å¯ä»¥è§å¯å°ï¼éç¨ææ ¼ææ¥å¤é¡¹å¼æå¼æ¶ï¼éçæå¼æ¬¡æ°è¶æ¥è¶å¤§ï¼æå¼å¤é¡¹å¼ä¸è§å¾æ¬¡æ°è¶é«è¶åç¡®ï¼Rungeç°è±¡ï¼ï¼Ln(x)并ä¸ä¸å®æ¶æäºf(x)ï¼èåºç°å¤é¶éè¡æ²çº¿å¨æå®åºé´ä¸å¾å¾ã
追ç