22/09/2018, 05:02:59 am *
Bienvenido(a), Visitante. Por favor, ingresa o regístrate.
¿Perdiste tu email de activación?

Ingresar con nombre de usuario, contraseña y duración de la sesión
Noticias: ¡Atención! Hay que poner la matemática con LaTeX, y se hace así (clic aquí):
 
 
Páginas: [1]   Ir Abajo
  Imprimir  
Autor Tema: Cuadrados mínimos - MATLAB  (Leído 364 veces)
0 Usuarios y 1 Visitante están viendo este tema.
Émilie.du
Junior
**

Karma: +0/-0
Desconectado Desconectado

Sexo: Femenino
France France

Mensajes: 19


Ver Perfil
« : 17/10/2017, 10:24:03 pm »

Hola, ¿como se podría dar una justificación válida para la siguiente afirmación?
Dados dos vectores [texx]x[/texx] e [texx]y[/texx] de la misma longitud [texx]m[/texx], el comando Matlab polyfit determina los coeficientes [texx]c_i[/texx] del polinomio de grado [texx]n[/texx]: [texx]p(x) = c_1x^n + · · ·+ c_nx + c_{n+1}[/texx] cuya gráfica ajusta por cuadrados mínimos los puntos [texx](x_1, y_1), . . . ,(x_m, y_m)[/texx]. En particular, cuando [texx]n=m−1[/texx], [texx]p(x)[/texx] es el polinomio de interpolación determinado por esos puntos.
En línea
delmar
Moderador Global
Pleno*
*

Karma: +0/-0
Desconectado Desconectado

Sexo: Masculino
Perú Perú

Mensajes: 1.399


Ver Perfil
« Respuesta #1 : 18/10/2017, 12:25:41 am »

Hola

La afirmación es cierta. ¿Por que?   [texx]n=m-1\Rightarrow{m=n+1}[/texx], es decir se tienen n+1 puntos de una función, desconocida en el resto de puntos. polyfit(x,y,n), retorna un vector constituido por los coeficientes de un polinomio de grado n,  [texx]W=(C_1,C_2,C_3,...C_{n+1})[/texx], que minimizan el error cuadrático. El polinomio sería : [texx]p(t)=C_1 t^n +C_2 t^{n-1}+C_3 t^{n-2}+...C_nt+C_{n+1}[/texx]

Es cierto que, si existiese un polinomio, respecto al cual el error cuadrático es cero (polinomio pasa por todos los puntos de la función), polyfit(x,y,n), retornaría sus coeficientes como vector. Para el caso n=m-1, existe, por la razón, que  se puede establecer un sistema de n+1 ecuaciones , con n+1 incógnitas (los coeficientes del polinomio), aquí el detalle :

[texx]p(x_1)=y_1[/texx]

[texx]p(x_2)=y_2[/texx]

....

[texx]p(x_n)=y_n[/texx]

[texx]p(x_{n+1})=y_{n+1}[/texx]

Esto se traduce en :

[texx]C_1 x_1^n +C_2 x_1^{n-1}+C_3 x_1^{n-2}+...C_nx_1+C_{n+1}=y_1[/texx]

[texx]C_1 x_2^n +C_2 x_2^{n-1}+C_3 x_2^{n-2}+...C_nx_2+C_{n+1}=y_2[/texx]

...

[texx]C_1 x_{n+1}^n +C_2 x_{n+1}^{n-1}+C_3 x_{n+1}^{n-2}+...C_nx_{n+1}+C_{n+1}=y_{n+1}[/texx]


Sistema de n+1 ecuaciones y n+1 incógnitas (coeficientes del polinomio), que se puede resolver, cuando las ecuaciones son linealmente independientes. Evidentemente, el polinomio pasa por los puntos de la función, y la aproxima en el resto de puntos, en consecuencia es un polinomio de interpolación,  con error cuadrático cero (pasa por todos los puntos), al tener error cuadrático cero, será también el que tiene el mínimo error cuadrático y en consecuencia, sus coeficientes serán los retornados, como un vector,  por la función polyfit

Saludos
En línea
Ignacio Larrosa
Moderador Global
Pleno*
*

Karma: +0/-0
Desconectado Desconectado

Sexo: Masculino
España España

Mensajes: 2.240


Ver Perfil WWW
« Respuesta #2 : 18/10/2017, 03:52:03 am »

Hola

La afirmación es cierta. ¿Por que?   [texx]n=m-1\Rightarrow{m=n+1}[/texx], es decir se tienen n+1 puntos de una función, desconocida en el resto de puntos. polyfit(x,y,n), retorna un vector constituido por los coeficientes de un polinomio de grado n,  [texx]W=(C_1,C_2,C_3,...C_{n+1})[/texx], que minimizan el error cuadrático. El polinomio sería : [texx]p(t)=C_1 t^n +C_2 t^{n-1}+C_3 t^{n-2}+...C_nt+C_{n+1}[/texx]

Es cierto que, si existiese un polinomio, respecto al cual el error cuadrático es cero (polinomio pasa por todos los puntos de la función), polyfit(x,y,n), retornaría sus coeficientes como vector. Para el caso n=m-1, existe, por la razón, que  se puede establecer un sistema de n+1 ecuaciones , con n+1 incógnitas (los coeficientes del polinomio), aquí el detalle :

[texx]p(x_1)=y_1[/texx]

[texx]p(x_2)=y_2[/texx]

....

[texx]p(x_n)=y_n[/texx]

[texx]p(x_{n+1})=y_{n+1}[/texx]

Esto se traduce en :

[texx]C_1 x_1^n +C_2 x_1^{n-1}+C_3 x_1^{n-2}+...C_nx_1+C_{n+1}=y_1[/texx]

[texx]C_1 x_2^n +C_2 x_2^{n-1}+C_3 x_2^{n-2}+...C_nx_2+C_{n+1}=y_2[/texx]

...

[texx]C_1 x_{n+1}^n +C_2 x_{n+1}^{n-1}+C_3 x_{n+1}^{n-2}+...C_nx_{n+1}+C_{n+1}=y_{n+1}[/texx]


Sistema de n+1 ecuaciones y n+1 incógnitas (coeficientes del polinomio), que se puede resolver, cuando las ecuaciones son linealmente independientes.

Esto último ocurre si y solo si los [texx]x_i[/texx] son todos distintos. En ese caso, el determinante de la matriz de coeficientes, determinante de Vandermonde, es distinto de cero.

Saludos,
En línea

Daría todo lo que se por la mitad de lo que ignoro (R. Descartes)
O incluso por muchísimo menos ...  (yo)
Páginas: [1]   Ir Arriba
  Imprimir  
 
Ir a:  

Impulsado por MySQL Impulsado por PHP Powered by SMF 1.1.4 | SMF © 2006, Simple Machines LLC XHTML 1.0 válido! CSS válido!