Problema con mínimos cuadrados y librerías para cálculo

[huevofritopamojarpan]

Hola! Este es mi primer mensaje.
Estoy aprendiendo python. Ya sé algo de C y bueno, no me está costando, la verdad.
El problema es que he hecho un programa en C al que se le pasan una serie de puntos del plano bidimensional y calcula la recta de aproximación mediante mínimos cuadrados y he intentado hacer lo mismo en Python, pero no funcona correctamente.

************************************************** **************
i=0
n=0
while n<3:
n = int(raw_input("\n\nCuantos puntos tiene la funcion a lineaizar: "))

horizontal=vertical=range(n)
while (i<n):
print "\n"
horizontal[i]=float(raw_input("x["+str(i)+"]="))
print horizontal[i]
vertical[i]=float(raw_input("y["+str(i)+"]="))
print vertical[i]

i=i+1

print "\n"
i=0
while (i<n):
print horizontal[i], vertical[i], "\n"
i = i + 1


i=0
SUMx=SUMy=SUMxy=SUMxx=a=b=0.0

while (i<n):
print SUMx, SUMy, SUMxy, SUMxy, "\n"
SUMx = SUMx + horizontal[i]
SUMy = SUMy + vertical[i]
SUMxy = SUMxy + (horizontal[i]*vertical[i])
SUMxx = SUMxx + (horizontal[i]*horizontal[i])
i = i + 1

print SUMx, SUMy, SUMxy, SUMxy, "\n"
a = ( (SUMx*SUMy) - (n*SUMxy) ) / ( (SUMx*SUMx) - (n*SUMxx) )
b = ( SUMy - (a*SUMx) ) / n

print a,"x",
if b>=0:
print "+",
print b
************************************************** **************


El problema es que el resultado que muestra es erróneo porque le asigna a x[i] el contenido de y[i].
Alguien sabe cuál donde está el problema?

Otra cosa es que me gustaría saber como funcionan las "librerías" numpy, scipy y demás, cómo se instalan y eso. Porque me interesaba sobre todo que los scripts que hago se puedan utilizar directamente en otros sistemas y ordenadores, pero creo que si utilizo una de estas librerías, el ordenador que la ejecuta tendría que tener estas librerías instaladas.

Muchas gracias.
Un saludo.

[venkman]

Cuando publiques código en el foro ponlo entre etiquetas de código como:

[highlight=python]
# aquí tu código
[/highlight]

Así saldrá como:

Código python:

# aquí tu código
 

En el caso de Python esto es particularmente importante porque los espacios importan. Y si no lo metes entre esas etiquetas, los espacios no se van a ver y el código no es igual.

[huevofritopamojarpan]

Ouch! Cómo dije antes, era mi primer mensaje, jeje.

Bueno, ahí va:

Código python:

i=0
n=0
while n<3:
    n = int(raw_input("\n\nCuantos puntos tiene la funcion a lineaizar: "))
 
horizontal=vertical=range(n)
while (i<n):
    print "\n"
    horizontal[i]=float(raw_input("x["+str(i)+"]="))
    print horizontal[i]
    vertical[i]=float(raw_input("y["+str(i)+"]="))
    print vertical[i]
    i=i+1
 
print "\n"
i=0
while (i<n):
    print horizontal[i], vertical[i], "\n"
    i = i + 1
 
i=0
SUMx=SUMy=SUMxy=SUMxx=a=b=0.0
 
while (i<n):
    print SUMx, SUMy, SUMxy, SUMxy, "\n"
    SUMx = SUMx + horizontal[i]
    SUMy = SUMy + vertical[i]
    SUMxy = SUMxy + (horizontal[i]*vertical[i])
    SUMxx = SUMxx + (horizontal[i]*horizontal[i])
    i = i + 1
 
print SUMx, SUMy, SUMxy, SUMxy, "\n"
a = ( (SUMx*SUMy) - (n*SUMxy) ) / ( (SUMx*SUMx) - (n*SUMxx) )
b = ( SUMy - (a*SUMx) ) / n 
 
print a,"x",
if b>=0:
    print "+",
print b
 

[yuo2]

Iniciado por huevofritopamojarpan

Otra cosa es que me gustaría saber como funcionan las "librerías" numpy, scipy y demás, cómo se instalan y eso. Porque me interesaba sobre todo que los scripts que hago se puedan utilizar directamente en otros sistemas y ordenadores, pero creo que si utilizo una de estas librerías, el ordenador que la ejecuta tendría que tener estas librerías instaladas.

Muchas gracias.
Un saludo.

Primero debes descargar la libreria y cuando lo instales selecciona la misma carpeta donde esta tu python y listo.
Una vez q estes en el GUI (consola) scribes
import tulibreria

Antes de practicar con la libreria busca documentacion sobre ella para poder usarla.

Te recomiendo la libreria matplotlib
http://matplotlib.sourceforge.net/

te la puedes descargar desde aqui:
http://sourceforge.net/project/showfiles.php?group_id=80706

y Si, tienes razon si quieres ejecutar tu programita esa pc tmb debe tener esa libreria, pero es rapido de instalar..


Aver si te pones un ejemplo correcto de la operacion de tu codigo para intentar hacerlo, aunq sea hecho a mano.. ;)
suerte...

[huevofritopamojarpan]

Hola! Muchas gracias por contestar. He visto (muy por encima) la de la librería matplotlib y me parece interesante, pero cuando intente hacer las gráficas utilizaré algo parecido al popen de c para interactuar con gnuplot, que es lo que uso normalmente.

Lo que pretendo hacer con este miniprogramilla es obtener la recta que mejor interpola a una serie de puntos.
Ejermplo:

x y
1 2
2 3
3 4

Si estos son los puntos la recta y=x+1 pasa por todos ellos. El objetivo del programa es obtener esta recta. Si no existiese una recta que pasa por todos ellos se obtendría la que está más cerca de todos ellos. Ese es el fin.
Pero el problema no viene de un mal algoritmo, ya que lo he implementado en C y funciona perfectamente. El problema tiene que estar en la toma de datos, ya que si los imprimo obtengo que x[i] es igual a y[i], cuando al introducir los valores por teclado yo no he puesto eso.

Muchas gracias por el interés.
Un saludo!

[venkman]

El problema de tu programa está claro. Es la línea 6. Cuando haces esto:

Código python:

horizontal=vertical=range(n)
 

Más concretamente, cuando haces:

Código python:

horizontal=vertical
 

No estás haciendo 2 listas distintas sino que lo que haces es que tanto vertical como horizontal apunten a la misma lista; por eso cuando cambias una cambias la otra (porque son la misma) y te salen los valores iguales. La solución es fácil: crear dos listas diferentes:

Código python:

horizontal=range(n)
vertical=range(n)
 

[huevofritopamojarpan]

Ajá! Joer, qué estupidez! No lo había pensado. Con esto de los tipos dinámicos me olvidé completamente de fijarme en eso.
Muchísimas gracias!!
Agur!