继续接收太多用于插值的数组索引
回答此问题可获得 20 贡献值，回答如果被采纳可获得 50 分。
<h2>interp-使用拉格朗日插值数据的程序</h2>
<p>我无法按照下面的编码顺序完成for循环。我没有看到任何错误，因为我选择NP.WORD（NPUTE）来创建席的一维数组，并且由于某些原因，循环不会填充这些值。p>
<pre><code>def intrpf(xi,x,y):
"""Function to interpolate between data points
using Lagrange polynomial (quadratic)
Inputs
x Vector of x coordinates of data points (3 values)
y Vector of y coordinates of data points (3 values)
xi The x value where interpolation is computed
Output
yi The interpolation polynomial evaluated at xi
"""
#* Calculate yi = p(xi) using Lagrange polynomial
yi = ( (xi-x[1])*(xi-x[2])/((x[0]-x[1])*(x[0]-x[2])) * y[0]
+ (xi-x[0])*(xi-x[2])/((x[1]-x[0])*(x[1]-x[2])) * y[1]
+ (xi-x[0])*(xi-x[1])/((x[2]-x[0])*(x[2]-x[1])) * y[2] )
return yi
#* Initialize the data points to be fit by quadratic
x = np.empty(3)
y = np.empty(3)
print ('Enter data points as x,y pairs (e.g., [1, 2]')
for i in range(3):
temp = np.array(input('Enter data point: '))
x[i] = temp[0]
y[i] = temp[1]
#* Establish the range of interpolation (from x_min to x_max)
xr = np.array(input('Enter range of x values as [x_min, x_max]: '))
</code></pre>
<p><strong>我在这一部分陷入了困境，它似乎设置正确，但for循环中的xi[I]上出现了“数组索引太多”的问题。</strong></p>
<pre><code>#* Find yi for the desired interpolation values xi using
# the function intrpf
nplot = 100 # Number of points for interpolation curve
xi = np.empty(nplot)
yi = np.empty(nplot)
for i in range(nplot) :
xi[i] = xr[0] + (xr[1]-xr[0])* i/float(nplot)
yi[i] = intrpf(xi[i], x, y) # Use intrpf function to interpolate
</code></pre>
请先 登录 后评论
针对此问题: 更多的回答
- 89 关注
- 1 收藏，216 浏览
- 网友 提问于 2天前
相关Python问题
- 如何在乒乓球比赛中预测球的轨迹，对于AI球拍预测？ 7 回答
- 如何在乒乓球游戏中阻止球 6 回答
- 如何在乘法和模中不乘空间？ 5 回答
- 如何在乘法和除以2个不同的数字之间进行交换？ 3 回答
- 如何在也是数据一部分的单个字符上拆分大字符串 1 回答
- 如何在乾草堆中找到針，有更好的解決方案嗎？ 4 回答
- 如何在事件wxWidgets中传递自定义数据 3 回答
- 如何在事件中使用lambda i=i？ 2 回答
- 如何在事件中心只接收最近的数据 2 回答
- 如何在事件发生之前保持云函数运行？ 9 回答
- 如何在事件发生后使页面重定向到同一页面 9 回答
- 如何在事件回调之间保持python生成器的状态 9 回答
- 如何在事件处理程序（pythonsocket、sphinx）中保留docstring 5 回答
- 如何在事件处理程序中更改wxRichTextCtrl的光标位置？ 5 回答
- 如何在事件处理程序中访问外部对象？ 4 回答
- 如何在事件循环中将协程打包为正常函数？ 7 回答
- 如何在事件循环之外运行协同程序？ 5 回答
- 如何在事件循环结束时为并发未来的所有线程调用类方法？ 1 回答
- 如何在事件文件中只保留一份摘要？ 3 回答
- 如何在事件模板中添加事件 6 回答
