뉴튼-랩슨 법
수학 혹은 물리학에 대해 다루고 있는 문서입니다. 오사위키는 사용자 여러분에게 유머가 아니더라도 불알을 탁 치게 만드는 일목요연하고 명쾌한 설명을 해준다면, 언제든지 환영입니다. 비전문가도 쉽게 이해할 수 있는 서술을 부탁드립니다 |
이 문서는 진지병에 걸렸습니다. 노잼을 감수하실 분만 스크롤을 내리십시오. 물론 진지병에 걸렸으니 보다가 토하셔도 오사위키는 책임 안 집니다. |
길이가 짧은 문서 문서가 너무 짧습니다. 마치 튀김기에서 갓 나온 치즈스틱처럼요! 식어버리기 전에 얼른 깨물고 잡아당겨 쭈욱 길이를 늘려주세요 |
개요[편집]
뉴튼-랩슨 법은 어떤 연속함수의 방정식을 구할 때 사용하는 방법으로 다음의 점화식을 말한다.
점화식[편집]
<math>x_{n}=x_{n-1}-\frac{f\left ( x_{n-1} \right )}{f'\left ( x_{n-1} \right )}</math>
초기값[편집]
초기값을 적당히 잘 잡아야 몇단계 내로 끝낼 수 있지만 이 초기값을 잘 잡는 게 오사스럽다.
간단한 방정식의 경우에는 굳이 초기값을 쫌 크게 잡아도 대부분 몇번 깔짝거리면 해당 방정식의 해로 수렴한다.
예제[편집]
ln 2[편집]
우리는 방정식 <math>f\left ( x \right )=e^{x}-2</math>의 근을 구할 것이다.
그리고 <math>f\left ( x \right )=e^{x}-2,f'\left ( x \right )=e^{x}</math>이다.
이제 점화식을 구해보면,
<math>x_{n}=x_{n-1}-\frac{e^{x_{n-1}}-2}{e^{x_{n-1}}}</math>
초기값을 1로 잡고 그 값을 표로 정리하면 다음과 같다.
비교후 맞는 자리까지 입력 후 그 값을 아래 점화식에 그대로 대입하였다.
예를 들어 0.693까지 맞으면 그 0.693이란 값을 그 점화식에 그대로 대입하였다.
횟수 | 값 | 소수점 이하 몇째 자리까지 맞는가. |
---|---|---|
1 | 1 | 0 |
2 | <math>x_{2}=1-\frac{e^{1}-2}{e^{1}}\approx0.73</math> | 0 |
3 | <math>x_{3}=0.73-\frac{e^{0.73}-2}{e^{0.73}}\approx0.693</math> | 3 |
4 | <math>x_{4}=0.691-\frac{e^{0.691}-2}{e^{0.691}}\approx0.6931471</math> | 7 |
5 | <math>x_{5}=0.6931471-\frac{e^{0.6931471}-2}{e^{0.6931471}}\approx0.69314718055994</math> | 14 |
6 | <math>x_{6}=0.69314718055994-\frac{e^{0.69314718055994}-2}{e^{0.69314718055994}}\approx0.6931471805599453094172321214</math> | 28 |
대충 이렇게 실제 계산기로 일고리즘을 적용하면 상당히 빨리 수렴한다고 한다..
물론 초기값을 잘 잡아야 몇번안에 수렴한다.
\(\sqrt{2}\)[편집]
이 수는 방정식 <math>x^{2}-2=0</math>의 근을 구하는 것과 같으므로 다음과 같다.
<math>\begin{cases}
& f\left ( x \right )=x^{2}-2 \\ & f'\left ( x \right )=2x
\end{cases}</math>
점화식은 다음과 같다.
<math>x_{n}=x_{n-1}-\frac{{x_{n-1}}^{2}-2}{2x_{n-1}}</math>
이제 구해보면...
초기값은 1이다.
횟수 | 근사값 | 소수점 이하 몇번째 자리까지 맞는가. |
---|---|---|
1 | 1 | 0 |
2 | <math>x_{2}=1-\frac{1^{2}-2}{2*1}=1.5</math> | 0 |
3 | <math>x_{3}=1.5-\frac{1.5^{2}-2}{2*1.5}\approx1.41</math> | 2 |
4 | <math>x_{4}=1.41-\frac{1.41^{2}-2}{2*1.41}\approx1.41421</math> | 5 |
5 | <math>x_{5}=1.41421-\frac{1.41421^{2}-2}{2*1.41421}\approx1.41421356237</math> | 11 |
6 | <math>x_{6}=1.41421356237-\frac{1.41421356237^{2}-2}{2*1.41421356237}\approx1.4142135623730950488016</math> | 22 |
7 | <math>x_{7}=1.4142135623730950488016-\frac{1.4142135623730950488016^{2}-2}{2*1.4142135623730950488016}\approx1.41421356237309504880168872420969807856967187</math> | 44 |
방정식 \(10^{x}+5^{x}-25=0\)의 근 구하기[편집]
일단 점화식을 세워 놓자.
<math>x_{n}=x_{n-1}-\frac{10^{x_{n-1}}+5^{x_{n-1}}-25}{10^{x_{n-1}}\ln 10+5^{x_{n-1}}\ln 5}</math>
역시 표를 만들자.
참값은 모르므로 대충 해당 수를 위 방정식에 대입하고 그것과 25의 오차를 계산하여야 한다.
초기값은 이번엔 2로 해보자.
횟수 | 값 | 방정식에 대입시 25와의 오차. |
---|---|---|
1 | 2 | 100 |
2 | 1.6303066574143 | 34.224894714379 |
3 | 1.3690541667307 | 9.7539643088519 |
4 | 1.2602343153801 | 7.4471099668641 |
5 | 1.2453191381825 | 0.012787933776494 |
6 | 1.2450753285419 | 3.342886010671E-6 |
7 | 1.2450752647744 | -1.5560885913146E-12 |
아무튼 오차의 절댓값이 빠르게 0으로 수렴한다.
방정식 \(\cos x+x^3=0\)의 근 구하기[편집]
<math>\begin{cases}
& f\left ( x \right )=\cos x+{x}^{3} \\ & f'\left ( x \right )=-\sin x+3{x}^{2}
\end{cases}</math>
점화식을 세워보면.
<math>x_{n}=x_{n-1}-\frac{\cos x_{n-1}-{{x_{n-1}}}^{3}}{-\sin x_{n-1}+3{{x_{n-1}}}^{2}}</math>
초기값은 2로 하겠습니다.
횟수 | 값 | 원래 방정식 대입 시 값 |
1 | 2 | 7.5838531635 |
2 | 1.316197227955930 | 2.5320048119 |
3 | 0.717524104322385 | 1.1228467036 |
4 | -0.548366659918160 | 0.6884799691 |
5 | -1.032049489939260 | -0.5862021602 |
6 | -0.887441424951214 | -0.0675083594 |
7 | -0.865929037655441 | -0.0013695035 |
8 | -0.865474233953938 | -6.04272077819878E-07 |
9 | -0.865474033101654 | -1.17794662912729E-13 |
10 | -0.865474033101614 | 대략 <math>{10}^{-15}</math>정도. |
방정식 \({x}^{x}=2\)의 근 구하기[편집]
<math>\begin{cases}
& f\left ( x \right )={x}^{x}-2 \\ & f'\left ( x \right )={x}^{x}\left ( \ln x+1 \right )
\end{cases}</math>
<math>x_{n}=x_{n-1}-\frac{{x_{n-1}}^{x_{n-1}}-2}{{x_{n-1}}^{x_{n-1}}\left ( \ln{x_{n-1}}+1 \right )}</math>
횟수 | 값 | 원래 방정식에 대입한 값-2 |
---|---|---|
1 | 2 | 2 |
2 | 1.70469194542517937512809654533837422144167347556305996627986 | 약 0.48 |
3 | 1.577944557476270445692208233738232151068642360866752760706135853 | 0.05 |
4 | 1.5599245375170789924282304443524568239218251485248606005822021465993 | 0.0009 |
5 | 1.55961056257717667811411825090521591061446411465339074078530052555797972 | 2.6E-7 |
6 | 1.559610469462377536255432533823718726734151756806798834775151244831110372965 | 2.3E-14 |
7 | 1.5596104694623693499703887688282766760754367277887488683140719268229542994126414 | 1.8E-28 |
8 | 1.55961046946236934997038876876500299328488351184309142472337460260886493677807203430 | 1E-56 |
9 | 1.559610469462369349970388768765002993284883511843091424719594569413973034549590587105413 | -1.2E87 |