Пример использования np.linalg.solve в решении уравнения

python > Пример использования np.linalg.solve в решении уравнения
05.04.2020 15:18:23


Наиболее часто встречающиеся слова в статье:

[системы] [matplotlib] [optimize] [curve_fit] [Матрица] ["Матрица] ["Вектор] ["Решение]


Статья:

from matplotlib import pylab as plt
from scipy.optimize import curve_fit
import math  as m
import numpy as np
 
def f(x):
    return np.sin(x / 5.) * np.exp(x / 10.) + 5. * np.exp(-x / 2.)
 
x = np.arange(0, 16, 1)
y = f(x)
plt.plot(x, y)
plt.show()
 
a = np.array([[1, 1,1,1], [1, 4,16,64], [1, 10,100,1000], [1, 15,225,3375]]) # Матрица (левая часть системы)
b = np.array([f(1), f(4), f(10),f(15)]) # Вектор (правая часть системы)
x = np.linalg.solve(a, b)
print(a.dot(x))
print("Матрица A:\n", a)
print("Вектор b:\n", b)
print("Решение системы:\n", x)