Estoy intentando generar un ajuste de parábola por medio de mínimos cuadrados, entonces estoy tomando 30 números aleatorios y toman valores para x,y y con eso genero mis listas.
Lo que intento hacer es ajustar el polinomio de segundo grado
P_2(x)= a_0 + a_1 X + a_2 x^2
y utilizo 30 datos aleatorios que están dentro de una región y cuando intento graficarlos me aparece la recta porque el x^2 se vuelve 0 pero cuando solo utilizo 9 datos que esten cercanos (en el rango) si se grafica.
El rango para los datos de x es de 153 a 446. El rango para los datos de y es de 117 a 170.
import numpy as np
import matplotlib.pyplot as plt
import scipy.linalg as linalg
dx= [333, 439, 266, 169, 360, 208, 384, 199, 379, 303, 392, 390, 279, 342, 430, 373, 189, 280, 218,
374, 366, 171, 396, 165, 328, 210, 441, 344, 216, 214]
dy= [128, 140, 142, 131, 136, 155, 159, 117, 122, 123, 127, 126, 139, 145, 151, 143, 147, 158, 124,
150, 167, 161, 152, 121, 164, 119, 168, 156, 118, 170]
xi=sorted(dx)
yi=sorted(dy)
x= np.array(xi)
y=np.array(yi)
plt.plot(x,y,'ok',label='data')
n=len(x)
F1=[n,sum(x),sum(x**2)]
F2=[sum(x),sum(x**2),sum(x**3)]
F3=[sum(x**2),sum(x**3),sum(x**4)]
M=np.array([F1,F2,F3])
v=np.array([sum(y),sum(x*y),sum((x**2)*y)])
solution=linalg.solve(M,v)
a0=solution[0]
a1=solution[1]
a2=solution[2]
#print("a0= {:.3f}".format(a0),"a1= {:.3f}".format(a1),"a2= {:.3f}".format(a2))
x_curve=np.linspace(0,5,100)
y_curve=a0+a1*x_curve+a2*(x_curve**2)
plt.plot(x_curve,y_curve,'b--',label='Mejor parabola')
plt.legend()
M
yv
tienes quelinalg.solver()
espera una matriz cuadrada de rango igual al número de betas (en tu caso 3), pero tu le pasasM
que no es una matriz como la anterior que he descrito. Yo lo que haría sería entender bien el problema que quieres resolver matemáticamente, y después buscar las herramientas que se adaptan a tu problema para resolverlo.