【发布时间】:2017-08-20 19:17:55
【问题描述】:
我有一个功能
f(x) = sin(x/5.0)*exp(x/10.0) + 5*exp(-x/2.0)
我需要求解线性方程组
w0 + w1x1 + w2(x1)**2 + ... + wn(x1)**n = f(x1)
我解决了这个问题,但我在绘制它时遇到了问题
from math import sin, exp
from scipy import linalg
import numpy as np
b = []
def f(x):
return sin(x/5.0)*exp(x/10.0) + 5*exp(-x/2.0)
for i in [1, 15]:
b.append(f(i))
A = []
for i in [1, 15]:
ij = []
x0 = i ** 0
x1 = i ** 1
ij.append(x0)
ij.append(x1)
A.append(ij)
matrix = np.array(A)
b = np.array(b).T
x = linalg.solve(matrix, b)
from matplotlib import pyplot as plt
plt.plot(x, f(x))
但它会返回
TypeError: only length-1 arrays can be converted to Python scalars
我该如何解决这个问题?
【问题讨论】:
标签: python matplotlib scipy