原文: http://math.mit.edu/~djk/calculus_beginners/chapter11/section01.html
如果我们有一个线性方程,例如
因此,在这个例子中,您可以向两侧添加
然而,假设我们有一个更复杂的等式,例如
我们的任务是找到这样一个等式的解决方案或所有解决方案。 我们假设我们方程中的函数在我们感兴趣的领域中是连续的和可微的。
首先要注意的是,在这里绘制左侧是一个好主意,并粗略地观察它改变符号的位置或非常接近
在过去,这是一项非常繁琐的任务,总的来说,人们试图在没有绘图的情况下解决方程,这有点像飞行盲目。如果你能做到这没关系,但为什么不尝试呢?
有一种解决这些方程式的标准技术显然可以追溯到牛顿。在这里。
你开始猜测你寻求的解决方案,选择一个参数,称之为
现在你做同样的事情,从
在过去,对于任何函数而言,这是一件非常繁琐的事情。从
现在有了电子表格,您可以在一分钟内完成设置,并通过练习找到解决方案。您只需执行一次每个步骤,然后复制。
怎么样?
具有通过适当地从两侧分开和减去而获得的解决方案
那我该怎么做电子表格呢?
假设我们在框 A1 中进行了第一次猜测。我们将把它和随后的猜测放在 A 栏开头说,用 A3(只是留下标签的空间)。
为此,我们需要进行以下输入:
在 A3 中,输入= A1(这将在 A3 中开始猜测
如果现在复制 A4(而不是 A3!)以及 A,B 和 C 列中的 B3 和 C3,则已实现该算法。
您可以通过更改 A1 来更改起始猜测,并通过适当更改 B3 和 C3 来更改您的函数,并将结果复制下来。
这真的有用吗?
这种方法在大多数时间内收敛得非常快。如果你从
什么是“好的一面”?
假设您从解决方案开始,调用解决方案
为什么?
在
所有这些意味着
假设我们想要求解方程
填写完成后,我们的电子表格说明应如下所示:
在 A1 中,输入 0.3 在 A2 中,输入 xj。在 B2 中,f(xj)。在 C2 中,f'(xj)。 在 A3 中,= A1。在 B3 中,= sin(A3)-exp(A3)+2。在 C3 中,= cos(A3)-exp(A3) 在 A4 中,= A3-B3 / C3。在 B4 中,= sin(A4)-exp(A4)+2。在 C4 中,= cos(A4)-exp(A4) 向下复制列 A,B 和 C.
<button aria-controls="newtons-method-spreadsheet" aria-expanded="false" class="btn bg-light border-secondary" data-target="#newtons-method-spreadsheet" data-toggle="collapse" id="toggle-spreadsheet-table" type="button">显示表</button>
Number of steps255 10 25 50Number of digits after decimal point105 10 15
练习:
11.1 假设
11. 2 如果
尽管如此,这种方法可以做出奇怪的事情。如果
以下小程序允许您只需输入函数即可绘制和查看方法。 (这比使用电子表格从头开始稍微简单一些)。
<iframe frameborder="0" height="620" src="../mathlets/newtons-method.html" width="100%"></iframe>
练习:
11.3 如果您寻找
不!您可以选择相对于函数比例非常小的 d 值,并将=(f(A3 + d)-f(A3-d))/(2 * d)放在 C3 而不是= f'( A3)。
这几乎与常规的牛顿方法一样。
练习 11.5 在 C3 中重做练习 11.5,=(f(A3 + B $ 1)-f(A3-B $ 1)/(2 * B $ 1)。你的答案如何受到影响?
可能出什么问题?
首先,我们的等式可能没有真正的解决方案。在那种情况下,这两种方法都找不到。针对
当你的方程有多个解时,会出现另一个问题。那么你得到的取决于你从哪里开始。 Applet 说明了这一点。
一般来说,如果你到达
但是这种方法无论如何都很有趣,你可以很容易地判断它什么时候不起作用。
分而治之
还有另一种求解方程的方法,在每次迭代时通过
练习:弄清楚如何在电子表格上实现此方法。 (提示:您可以输入= if(D5 * F5> 0,C5,A5),如果 D5 和 F5 具有相同的符号则给出 C5,否则输入 A5。)