MateWiki - Contribuciones del usuario [es]

Trabajo sobre algoritmo de la biseccion grupo 23

2019-12-13T12:37:49Z

C.novero:

{{ TrabajoED | Aproximación de raíces por el método de bisección. Grupo 23 | [[:Categoría:Matemáticas I|Matemáticas I]]|[[:Categoría:MatI/19|Curso 2019-20]] | Covadonga Corral Martín, Carlos Novero Garrido, Bárbara Herraez Arribas }}

Vamos a proceder a encontrar un valor de x tal que tan(x)=4x

== Planteamiento ==

Queremos averiguar el punto de corte de las funciones tan (x) y 4x

== Método ==

El método de la bisección consiste en aplicar el teorema de Bolzano, para ellos partimos de un intervalo (a,b) que verifica f(a)*f(b)<0 y se procede a dividir el intervalo en dos mitades y ver cual de los dos subintervalos cumple dicha condición
== Aplicación ==

Para averiguarlo vamos a aplicar el método de la bisección a la función f(x)=tan(x)-4x, de manera que encontramos una raiz que se corresponde con el punto de corte de ambas funciones, lo haremos en el intervalo (-0.3,0.5) y el error máximo que vamos a admitir es de 10^-3 por lo que el programa se detendrá cuando el intervalo sea inferior, en longitud, a 10^-3

EL valor de la aproximación es -3.906249999999930e-04, valor muy cercano a cero por lo que podemos considerarlo el punto que buscamos, el error es menor de 10^-3
== Programa ==

{{matlab|codigo=m
f=@(x) tan(x)-(4*x); %definimos la función en función de X
ei=-0.3; %establecemos el extremo inferior
ed=0.5; %establecemos el extremo superior
while (ed-ei)>1.e-3 %mientras el intervalo es mayor que 10^-3 el programa sigue realizando el bucle
if f(ei)*f((ei+ed)/2)<0 %si se cumple el teorema de bolzano entre el extremo inicial y la mitad del intervalo, establecemos la mitad del intervalo como el extremo final
ed=(ei+ed)/2;
else %de lo contrario establecemos la mitad del intervalo como extremo inicial
ei=(ei+ed)/2;
end
end %se para el bucle
sol=(ei+ed)/2 %la solución es la mitad del último intervalo
}}

[[Categoría:Matemáticas I]]
[[Categoría:MatI/19]]

Trabajo sobre algoritmo de la biseccion grupo 23

2019-12-13T12:36:41Z

C.novero: Página creada con «{{ TrabajoED | Aproximación de raíces por el método de bisección. Grupo 23 | Matemáticas I|Curso 2019-20 | Nuest...»

{{ TrabajoED | Aproximación de raíces por el método de bisección. Grupo 23 | [[:Categoría:Matemáticas I|Matemáticas I]]|[[:Categoría:MatI/19|Curso 2019-20]] | Nuestros nombres }}

Vamos a proceder a encontrar un valor de x tal que tan(x)=4x

== Planteamiento ==

Queremos averiguar el punto de corte de las funciones tan (x) y 4x

== Método ==

El método de la bisección consiste en aplicar el teorema de Bolzano, para ellos partimos de un intervalo (a,b) que verifica f(a)*f(b)<0 y se procede a dividir el intervalo en dos mitades y ver cual de los dos subintervalos cumple dicha condición
== Aplicación ==

Para averiguarlo vamos a aplicar el método de la bisección a la función f(x)=tan(x)-4x, de manera que encontramos una raiz que se corresponde con el punto de corte de ambas funciones, lo haremos en el intervalo (-0.3,0.5) y el error máximo que vamos a admitir es de 10^-3 por lo que el programa se detendrá cuando el intervalo sea inferior, en longitud, a 10^-3

EL valor de la aproximación es -3.906249999999930e-04, valor muy cercano a cero por lo que podemos considerarlo el punto que buscamos, el error es menor de 10^-3
== Programa ==

{{matlab|codigo=m
f=@(x) tan(x)-(4*x); %definimos la función en función de X
ei=-0.3; %establecemos el extremo inferior
ed=0.5; %establecemos el extremo superior
while (ed-ei)>1.e-3 %mientras el intervalo es mayor que 10^-3 el programa sigue realizando el bucle
if f(ei)*f((ei+ed)/2)<0 %si se cumple el teorema de bolzano entre el extremo inicial y la mitad del intervalo, establecemos la mitad del intervalo como el extremo final
ed=(ei+ed)/2;
else %de lo contrario establecemos la mitad del intervalo como extremo inicial
ei=(ei+ed)/2;
end
end %se para el bucle
sol=(ei+ed)/2 %la solución es la mitad del último intervalo
}}

[[Categoría:Matemáticas I]]
[[Categoría:MatI/19]]

Aproximación por minimos cuadrados grupo 23

2019-11-12T11:18:31Z

C.novero:

{{ Trabajo | Aproximación por mínimos cuadrados grupo 23 | [[:Categoría:Matemáticas I|Matemáticas I]]|[[:Categoría:MatI/19|Curso 2019-20]] | Covadonga Corral Martín, Mohamed Ayoub Ould el Makak Berjalí, Carlos Novero Garrido }}
Se está testeando un nuevo metamaterial sintético en el laboratorio frente a la fatiga. Para ello se golpea periódicamente el material con un dispositivo automático, manteniendo la fuerza pero cambiando la frecuencia. Se registran los datos de tiempo de rotura tr en horas para cada frecuencia de golpeo ω (número de golpes por minuto).
Los datos obtenidos son:
{| class="wikitable"
|+'''Frecuencia/tiempo de rotura'''
|-
|1.33
|1040
|-
|2.67
|1036
|-
|4
|999
|-
|5.33
|1017
|-
|6.67
|1049
|-
|8.67
|1000
|-
|10.67
|1007
|-
|12
|NaN
|-
|14.67
|NaN
|-
|16
|NaN
|-
|17.3
|NaN
|-
|18.67
|943
|-
|20
|916
|-
|21.33
|944
|-
|22
|929
|-
|22.67
|955
|-
|24
|940
|-
|24,67
|911
|-
|25.33
|785
|-
|25.67
|703
|-
|26
|634
|-
|26.33
|326
|-
|26.67
|7
|}
Algunos ensayos no tienen datos por un error en el laboratorio.
== Dibujar los puntos y generar la gráfica que los representa ==
Comenzamos por introducir en Matlab los datos y a continuación generar la gráfica que los representa para su fácil visualización:
{{matlab|codigo=
w=[1.33 2.67 4 5.33 6.67 8.67 10.67 18.67 20 21.33 22 22.67 23.33 24 24.67 25.33 25.67 26 26.33 26.67] %vector de la frecuencia
tr=[1040 1036 999 1017 1049 1000 1007 943 916 944 929 955 931 940 911 785 703 624 326 7] %vector del tiempo de rotura
plot(w.',tr.','k--o')
xlabel('frecuencia')
ylabel('tiempo de rotura')
title('ensayo fatiga')
grid %dibujamos la gráfica tiempo de rotura/frecuencia y le damos formato
}}
== Cálculo y representación de una recta por mínimos cuadrados ==
=== Cálculo de la recta ===
Procedemos a ajustar los datos a una recta, para ello introducimos en Matlab:
{{matlab|codigo=
% añadir al programa anterior
m1r=[ones(20,1),w.'] %matriz de los coeficientes, tamaño 20x2, siendo la primera fila una columna de unos y la segunda el vector de las frecuencias traspuesto
m2r=[m1r,tr.'] %matriz ampliada de la matriz anterior junto con el vector de los tiempos de rotura traspuesto
rank(m1r)
rank(m2r) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento obtenemos que el rango de m1 es 2 y el de m3 es 3, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la recta y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
% añadir al programa anterior
xr=m1r\tr.' %definimos xr como el vector con las soluciones del sistema de la recta
yr=xr(1)+xr(2)*w
hold on %le pedimos que mantenga la gráfica anterior
plot(w,yr,'b-') %añade la recta calculada por mínimos cuadrados
legend('datos','recta de minimos') %añade una leyenda
}}

=== Cálculo del error cuadrático medio de la recta ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
% añadir al programa anterior
errorr=(sum((tr-yr).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
% añadir al programa anterior
zr=xr(1)+xr(2)*16
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xr=[1.146876510867884e+03;-16.980810431368106], por lo que la ecuación de la recta es: tr=1146.9-17*ω

errorr=4.310542091992706e+04, el error cuadrático medio es de 43105

zr=8.751835439659939e+02, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 875 horas 11 minutos
y además representa la siguiente gráfica:
[[Archivo:Grafica_recta_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación lineal de mínimos cuadrados]]

== Cálculo y representación de una parábola por mínimos cuadrados ==
=== Cálculo de la parábola ===
Procedemos a ajustar los datos a una parábola, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1p=[ones(20,1),w.',w.'.*w.'] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de las frecuencias al cuadrado
m2p=[m1p,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1p)
rank(m2p) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices, siendo 3 para m1 y 4 para m2, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
xp=m1p\tr.'
y1p=xp(1)+xp(2)*w+xp(3)*w.^2 %ecuación de la parábola calculada por mínimos cuadrados
h=1:2:27
y2p=xp(1)+xp(2)*h+xp(3)*h.^2 %representamos una parábola equivalente usando más datos
hold on
plot(h,y2p,'b-')
legend('datos','parabola de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
errorp=(sum((tr-y1p).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
Zp=xp(1)+xp(2)*16+xp(3)*16.^2
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xp=[8.301635504889812e+02;57.398707088415634;-2.553441460170276] por lo que la ecuación de la parábola es: tr=830.16-57.39*ω-2.55*ω^2

errorp=2.947828440943336e+04 por lo que el error cuadrático medio es de 29478

zp=1.094861850100041e+03, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 1094 horas, y además representa la siguiente gráfica:
[[Archivo:Gráifca_parabola_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación parabólica de mínimos cuadrados]]

== Cálculo y representación de una ecuación exponencial por mínimos cuadrados ==
=== Cálculo de la ecuación exponencial ===
Procedemos a ajustar los datos a una ecuación exponencial, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1e=[ones(20,1),w.',exp(w.')] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de los exponentes de las frecuencias
m2e=[m1e,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1e)
rank(m2e) %calculamos los rangos de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices y obtenemos que es 3 para m1 y 4 para m2 por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
xe=m1e\tr.' %definimos x como el vector con las soluciones del sistema
ye=xe(1)+xe(2)*w+xe(3)*exp(w) %ecuación exponencial calculada por mínimos cuadrados
hold on
plot(w,ye,'b-')
legend('datos','exponencial de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
errore=sum((tr-ye).^2)/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
ze=xe(1)+xe(2)*16+xe(3)*exp(16)
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xe=[1.025976427703857e+03;-1.974268082492748;-2.319344866513919e-09] por lo que la ecuación exponencial es tr=1025.97-1.97*w-2.319e-9*exp(w)

errore=1.908149369400242e+03 el error cuadrático medio vale 1908.1

ze=9.943675284291536e+02 para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 994.36 horas y además genera la siguiente gráfica:
[[Archivo:Gráfica_exponencial_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación exponencial de mínimos cuadrados]]

== Conclusiones ==
La distribución con el menor error cuadrático medio es la función exponencial, además es la que más coincide, visualmente, con los datos de partida por lo que nos parece que 1098.1 es la predicción más fiable de las 3

[[Categoría:Matemáticas I]]
[[Categoría:MatI/19]]

Archivo:Grafica recta minimos grupo 23.png

2019-11-12T11:15:55Z

C.novero:

Aproximación por minimos cuadrados grupo 23

2019-11-12T11:00:34Z

C.novero:

{{ Trabajo | Aproximación por mínimos cuadrados grupo 23 | [[:Categoría:Matemáticas I|Matemáticas I]]|[[:Categoría:MatI/19|Curso 2019-20]] | Covadonga Corral Martín, Mohamed Ayoub Ould el Makak Berjalí, Carlos Novero Garrido }}
Se está testeando un nuevo metamaterial sintético en el laboratorio frente a la fatiga. Para ello se golpea periódicamente el material con un dispositivo automático, manteniendo la fuerza pero cambiando la frecuencia. Se registran los datos de tiempo de rotura tr en horas para cada frecuencia de golpeo ω (número de golpes por minuto).
Los datos obtenidos son:
{| class="wikitable"
|+'''Frecuencia/tiempo de rotura'''
|-
|1.33
|1040
|-
|2.67
|1036
|-
|4
|999
|-
|5.33
|1017
|-
|6.67
|1049
|-
|8.67
|1000
|-
|10.67
|1007
|-
|12
|NaN
|-
|14.67
|NaN
|-
|16
|NaN
|-
|17.3
|NaN
|-
|18.67
|943
|-
|20
|916
|-
|21.33
|944
|-
|22
|929
|-
|22.67
|955
|-
|24
|940
|-
|24,67
|911
|-
|25.33
|785
|-
|25.67
|703
|-
|26
|634
|-
|26.33
|326
|-
|26.67
|7
|}
Algunos ensayos no tienen datos por un error en el laboratorio.
== Dibujar los puntos y generar la gráfica que los representa ==
Comenzamos por introducir en Matlab los datos y a continuación generar la gráfica que los representa para su fácil visualización:
{{matlab|codigo=
w=[1.33 2.67 4 5.33 6.67 8.67 10.67 18.67 20 21.33 22 22.67 23.33 24 24.67 25.33 25.67 26 26.33 26.67] %vector de la frecuencia
tr=[1040 1036 999 1017 1049 1000 1007 943 916 944 929 955 931 940 911 785 703 624 326 7] %vector del tiempo de rotura
plot(w.',tr.','k--o')
xlabel('frecuencia')
ylabel('tiempo de rotura')
title('ensayo fatiga')
grid %dibujamos la gráfica tiempo de rotura/frecuencia y le damos formato
}}
== Cálculo y representación de una recta por mínimos cuadrados ==
=== Cálculo de la recta ===
Procedemos a ajustar los datos a una recta, para ello introducimos en Matlab:
{{matlab|codigo=
% añadir al programa anterior
m1r=[ones(20,1),w.'] %matriz de los coeficientes, tamaño 20x2, siendo la primera fila una columna de unos y la segunda el vector de las frecuencias traspuesto
m2r=[m1r,tr.'] %matriz ampliada de la matriz anterior junto con el vector de los tiempos de rotura traspuesto
rank(m1r)
rank(m2r) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento obtenemos que el rango de m1 es 2 y el de m3 es 3, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la recta y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
% añadir al programa anterior
xr=m1r\tr.' %definimos xr como el vector con las soluciones del sistema de la recta
yr=xr(1)+xr(2)*w
hold on %le pedimos que mantenga la gráfica anterior
plot(w,yr,'b-') %añade la recta calculada por mínimos cuadrados
legend('datos','recta de minimos') %añade una leyenda
}}

=== Cálculo del error cuadrático medio de la recta ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
% añadir al programa anterior
errorr=(sum((tr-yr).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
% añadir al programa anterior
zr=xr(1)+xr(2)*16
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xr=[1.146876510867884e+03;-16.980810431368106], por lo que la ecuación de la recta es: tr=1146.9-17*ω

errorr=4.310542091992706e+04, el error cuadrático medio es de 43105

zr=8.751835439659939e+02, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 875 horas 11 minutos
y además representa la siguiente gráfica:
[[Archivo:Gráfica_recta_mínimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación lineal de mínimos cuadrados]]

== Cálculo y representación de una parábola por mínimos cuadrados ==
=== Cálculo de la parábola ===
Procedemos a ajustar los datos a una parábola, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1p=[ones(20,1),w.',w.'.*w.'] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de las frecuencias al cuadrado
m2p=[m1p,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1p)
rank(m2p) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices, siendo 3 para m1 y 4 para m2, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
xp=m1p\tr.'
y1p=xp(1)+xp(2)*w+xp(3)*w.^2 %ecuación de la parábola calculada por mínimos cuadrados
h=1:2:27
y2p=xp(1)+xp(2)*h+xp(3)*h.^2 %representamos una parábola equivalente usando más datos
hold on
plot(h,y2p,'b-')
legend('datos','parabola de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
errorp=(sum((tr-y1p).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
Zp=xp(1)+xp(2)*16+xp(3)*16.^2
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xp=[8.301635504889812e+02;57.398707088415634;-2.553441460170276] por lo que la ecuación de la parábola es: tr=830.16-57.39*ω-2.55*ω^2

errorp=2.947828440943336e+04 por lo que el error cuadrático medio es de 29478

zp=1.094861850100041e+03, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 1094 horas, y además representa la siguiente gráfica:
[[Archivo:Gráifca_parabola_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación parabólica de mínimos cuadrados]]

== Cálculo y representación de una ecuación exponencial por mínimos cuadrados ==
=== Cálculo de la ecuación exponencial ===
Procedemos a ajustar los datos a una ecuación exponencial, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1e=[ones(20,1),w.',exp(w.')] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de los exponentes de las frecuencias
m2e=[m1e,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1e)
rank(m2e) %calculamos los rangos de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices y obtenemos que es 3 para m1 y 4 para m2 por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
xe=m1e\tr.' %definimos x como el vector con las soluciones del sistema
ye=xe(1)+xe(2)*w+xe(3)*exp(w) %ecuación exponencial calculada por mínimos cuadrados
hold on
plot(w,ye,'b-')
legend('datos','exponencial de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
errore=sum((tr-ye).^2)/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
ze=xe(1)+xe(2)*16+xe(3)*exp(16)
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xe=[1.025976427703857e+03;-1.974268082492748;-2.319344866513919e-09] por lo que la ecuación exponencial es tr=1025.97-1.97*w-2.319e-9*exp(w)

errore=1.908149369400242e+03 el error cuadrático medio vale 1908.1

ze=9.943675284291536e+02 para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 994.36 horas y además genera la siguiente gráfica:
[[Archivo:Gráfica_exponencial_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación exponencial de mínimos cuadrados]]

== Conclusiones ==
La distribución con el menor error cuadrático medio es la función exponencial, además es la que más coincide, visualmente, con los datos de partida por lo que nos parece que 1098.1 es la predicción más fiable de las 3

[[Categoría:Matemáticas I]]
[[Categoría:MatI/19]]

Aproximación por minimos cuadrados grupo 23

2019-11-11T12:20:19Z

C.novero:

{{ Trabajo | Aproximación por mínimos cuadrados grupo 23 | [[:Categoría:Matemáticas I|Matemáticas I]]|[[:Categoría:MatI/19|Curso 2019-20]] | Covadonga Corral Martín, Mohamed Ayoub Ould el Makak Berjalí, Carlos Novero Garrido }}
Se está testeando un nuevo metamaterial sintético en el laboratorio frente a la fatiga. Para ello se golpea periódicamente el material con un dispositivo automático, manteniendo la fuerza pero cambiando la frecuencia. Se registran los datos de tiempo de rotura tr en horas para cada frecuencia de golpeo ω (número de golpes por minuto).
Los datos obtenidos son:
{| class="wikitable"
|+'''Frecuencia/tiempo de rotura'''
|-
|1.33
|1040
|-
|2.67
|1036
|-
|4
|999
|-
|5.33
|1017
|-
|6.67
|1049
|-
|8.67
|1000
|-
|10.67
|1007
|-
|12
|NaN
|-
|14.67
|NaN
|-
|16
|NaN
|-
|17.3
|NaN
|-
|18.67
|943
|-
|20
|916
|-
|21.33
|944
|-
|22
|929
|-
|22.67
|955
|-
|24
|940
|-
|24,67
|911
|-
|25.33
|785
|-
|25.67
|703
|-
|26
|634
|-
|26.33
|326
|-
|26.67
|7
|}
Algunos ensayos no tienen datos por un error en el laboratorio.
== Dibujar los puntos y generar la gráfica que los representa ==
Comenzamos por introducir en Matlab los datos y a continuación generar la gráfica que los representa para su fácil visualización:
{{matlab|codigo=
w=[1.33 2.67 4 5.33 6.67 8.67 10.67 18.67 20 21.33 22 22.67 23.33 24 24.67 25.33 25.67 26 26.33 26.67] %vector de la frecuencia
tr=[1040 1036 999 1017 1049 1000 1007 943 916 944 929 955 931 940 911 785 703 624 326 7] %vector del tiempo de rotura
plot(w.',tr.','k--o')
xlabel('frecuencia')
ylabel('tiempo de rotura')
title('ensayo fatiga')
grid %dibujamos la gráfica tiempo de rotura/frecuencia y le damos formato
}}
== Cálculo y representación de una recta por mínimos cuadrados ==
=== Cálculo de la recta ===
Procedemos a ajustar los datos a una recta, para ello introducimos en Matlab:
{{matlab|codigo=
% añadir al programa anterior
m1r=[ones(20,1),w.'] %matriz de los coeficientes, tamaño 20x2, siendo la primera fila una columna de unos y la segunda el vector de las frecuencias traspuesto
m2r=[m1r,tr.'] %matriz ampliada de la matriz anterior junto con el vector de los tiempos de rotura traspuesto
rank(m1r)
rank(m2r) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento obtenemos que el rango de m1 es 2 y el de m3 es 3, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la recta y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
% añadir al programa anterior
xr=m1r\tr.' %definimos xr como el vector con las soluciones del sistema de la recta
yr=xr(1)+xr(2)*w
hold on %le pedimos que mantenga la gráfica anterior
plot(w,yr,'b-') %añade la recta calculada por mínimos cuadrados
legend('recta de mínimos','datos') %añade una leyenda
}}

=== Cálculo del error cuadrático medio de la recta ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
% añadir al programa anterior
errorr=(sum((tr-yr).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
% añadir al programa anterior
zr=x(1)+x(2)*16
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xr=[1.146876510867884e+03;-16.980810431368106], por lo que la ecuación de la recta es: tr=1146.9-17*ω

errorr=4.310542091992706e+04, el error cuadrático medio es de 43105

zr=8.751835439659939e+02, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 875 horas 11 minutos
y además representa la siguiente gráfica:
[[Archivo:Gráfica_recta_mínimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación lineal de mínimos cuadrados]]

== Cálculo y representación de una parábola por mínimos cuadrados ==
=== Cálculo de la parábola ===
Procedemos a ajustar los datos a una parábola, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1p=[ones(20,1),w.',w.'.*w.'] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de las frecuencias al cuadrado
m2p=[m1p,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1p)
rank(m2p) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices, siendo 3 para m1 y 4 para m2, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
y1p=x(1)+x(2)*w+x(3)*w.^2 %ecuación de la parábola calculada por mínimos cuadrados
h=1:2:27
y2p=x(1)+x(2)*h+x(3)*h.^2 %representamos una parábola equivalente usando más datos
hold on
plot(h,y2p,'b-')
legend('datos','parabola de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
errorp=(sum((tr-y1p).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
Zp=x(1)+x(2)*16+x(3)*16.^2
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xp=[8.301635504889812e+02;57.398707088415634;-2.553441460170276] por lo que la ecuación de la parábola es: tr=830.16-57.39*ω-2.55*ω^2

errorp=2.947828440943336e+04 por lo que el error cuadrático medio es de 29478

zp=1.094861850100041e+03, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 1094 horas, y además representa la siguiente gráfica:
[[Archivo:Gráifca_parabola_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación parabólica de mínimos cuadrados]]

== Cálculo y representación de una ecuación exponencial por mínimos cuadrados ==
=== Cálculo de la ecuación exponencial ===
Procedemos a ajustar los datos a una ecuación exponencial, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1e=[ones(20,1),w.',exp(w.')] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de los exponentes de las frecuencias
m2e=[m1e,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1e)
rank(m2e) %calculamos los rangos de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices y obtenemos que es 3 para m1 y 4 para m2 por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
xe=m1e\tr.' %definimos x como el vector con las soluciones del sistema
ye=x(1)+x(2)*w+x(3)*exp(w) %ecuación exponencial calculada por mínimos cuadrados
hold on
plot(w,ye,'b-')
legend('datos','exponencial de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
error=sum((tr-ye).^2)/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
ze=x(1)+x(2)*16+x(3)*exp(16)
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xe=[1.025976427703857e+03;-1.974268082492748;-2.319344866513919e-09] por lo que la ecuación exponencial es tr=1025.97-1.97*w-2.319e-9*exp(w)

errore=1.908149369400242e+03 el error cuadrático medio vale 1908.1

ze=9.943675284291536e+02 para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 994.36 horas y además genera la siguiente gráfica:
[[Archivo:Gráfica_exponencial_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación exponencial de mínimos cuadrados]]

== Conclusiones ==
La distribución con el menor error cuadrático medio es la función exponencial, además es la que más coincide, visualmente, con los datos de partida por lo que nos parece que 1098.1 es la predicción más fiable de las 3

[[Categoría:Matemáticas I]]
[[Categoría:MatI/19]]

Aproximación por minimos cuadrados grupo 23

2019-11-11T12:17:15Z

C.novero: Página creada con «{{ Trabajo | Aproximación por mínimos cuadrados grupo 23 | Matemáticas I|Curso 2019-20 | Covadonga Corral Martín,...»

{{ Trabajo | Aproximación por mínimos cuadrados grupo 23 | [[:Categoría:Matemáticas I|Matemáticas I]]|[[:Categoría:MatI/19|Curso 2019-20]] | Covadonga Corral Martín, Mohamed Ayoub Ould el Makak Berjalí, Carlos Novero Garrido }}
Se está testeando un nuevo metamaterial sintético en el laboratorio frente a la fatiga. Para ello se golpea periódicamente el material con un dispositivo automático, manteniendo la fuerza pero cambiando la frecuencia. Se registran los datos de tiempo de rotura tr en horas para cada frecuencia de golpeo ω (número de golpes por minuto).
Los datos obtenidos son:
{| class="wikitable"
|+'''Frecuencia/tiempo de rotura'''
|-
|1.33
|1040
|-
|2.67
|1036
|-
|4
|999
|-
|5.33
|1017
|-
|6.67
|1049
|-
|8.67
|1000
|-
|10.67
|1007
|-
|12
|NaN
|-
|14.67
|NaN
|-
|16
|NaN
|-
|17.3
|NaN
|-
|18.67
|943
|-
|20
|916
|-
|21.33
|944
|-
|22
|929
|-
|22.67
|955
|-
|24
|940
|-
|24,67
|911
|-
|25.33
|785
|-
|25.67
|703
|-
|26
|634
|-
|26.33
|326
|-
|26.67
|7
|}
Algunos ensayos no tienen datos por un error en el laboratorio.
== Dibujar los puntos y generar la gráfica que los representa ==
Comenzamos por introducir en Matlab los datos y a continuación generar la gráfica que los representa para su fácil visualización:
{{matlab|codigo=
w=[1.33 2.67 4 5.33 6.67 8.67 10.67 18.67 20 21.33 22 22.67 23.33 24 24.67 25.33 25.67 26 26.33 26.67] %vector de la frecuencia
tr=[1040 1036 999 1017 1049 1000 1007 943 916 944 929 955 931 940 911 785 703 624 326 7] %vector del tiempo de rotura
plot(w.',tr.','k--o')
xlabel('frecuencia')
ylabel('tiempo de rotura')
title('ensayo fatiga')
grid %dibujamos la gráfica tiempo de rotura/frecuencia y le damos formato
}}
== Cálculo y representación de una recta por mínimos cuadrados ==
=== Cálculo de la recta ===
Procedemos a ajustar los datos a una recta, para ello introducimos en Matlab:
{{matlab|codigo=
% añadir al programa anterior
m1r=[ones(20,1),w.'] %matriz de los coeficientes, tamaño 20x2, siendo la primera fila una columna de unos y la segunda el vector de las frecuencias traspuesto
m2r=[m1r,tr.'] %matriz ampliada de la matriz anterior junto con el vector de los tiempos de rotura traspuesto
rank(m1r)
rank(m2r) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento obtenemos que el rango de m1 es 2 y el de m3 es 3, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la recta y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
% añadir al programa anterior
xr=m1r\tr.' %definimos xr como el vector con las soluciones del sistema de la recta
yr=xr(1)+xr(2)*w
hold on %le pedimos que mantenga la gráfica anterior
plot(w,yr,'b-') %añade la recta calculada por mínimos cuadrados
legend('recta de mínimos','datos') %añade una leyenda
}}

=== Cálculo del error cuadrático medio de la recta ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
% añadir al programa anterior
errorr=(sum((tr-yr).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
% añadir al programa anterior
zr=x(1)+x(2)*16
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xr=[1.146876510867884e+03;-16.980810431368106], por lo que la ecuación de la recta es: tr=1146.9-17*ω

errorr=4.310542091992706e+04, el error cuadrático medio es de 43105

zr=8.751835439659939e+02, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 875 horas 11 minutos
y además representa la siguiente gráfica:
[[Archivo:Gráfica_recta_mínimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación lineal de mínimos cuadrados]]

== Cálculo y representación de una parábola por mínimos cuadrados ==
=== Cálculo de la parábola ===
Procedemos a ajustar los datos a una parábola, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1p=[ones(20,1),w.',w.'.*w.'] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de las frecuencias al cuadrado
m2p=[m1p,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1p)
rank(m2p) %usamos la función que nos calcula el rango de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices, siendo 3 para m1 y 4 para m2, por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
y1p=x(1)+x(2)*w+x(3)*w.^2 %ecuación de la parábola calculada por mínimos cuadrados
h=1:2:27
y2p=x(1)+x(2)*h+x(3)*h.^2 %representamos una parábola equivalente usando más datos
hold on
plot(h,y2p,'b-')
legend('datos','parabola de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
errorp=(sum((tr-y1p).^2))/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
Zp=x(1)+x(2)*16+x(3)*16.^2
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xp=[8.301635504889812e+02;57.398707088415634;-2.553441460170276] por lo que la ecuación de la parábola es: tr=830.16-57.39*ω-2.55*ω^2

errorp=2.947828440943336e+04 por lo que el error cuadrático medio es de 29478

zp=1.094861850100041e+03, es decir, para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 1094 horas, y además representa la siguiente gráfica:
[[Archivo:Gráifca_parabola_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación parabólica de mínimos cuadrados]]

== Cálculo y representación de una ecuación exponencial por mínimos cuadrados ==
=== Cálculo de la ecuación exponencial ===
Procedemos a ajustar los datos a una ecuación exponencial, para ello introducimos en Matlab:
{{matlab|codigo=
%añadir al programa anterior
m1e=[ones(20,1),w.',exp(w.')] %matriz de los coeficientes, tamaño 20x3, siendo la primera fila una columna de unos, la segunda el vector de las frecuencias traspuestas y la tercera es el vector de los exponentes de las frecuencias
m2e=[m1e,tr.'] %matriz ampliada con los datos del tiempo de rotura
rank(m1e)
rank(m2e) %calculamos los rangos de las matrices
}}
Si ejecutásemos el programa en este momento nos calcularía el rango de cada una de las matrices y obtenemos que es 3 para m1 y 4 para m2 por lo que sabemos que el sistema es incompatible, pasamos a calcular la resolución del sistema, que Matlab calcula automáticamente por mínimos cuadrados y a continuación definimos la parábola y la añadimos en el gráfico creado anteriormente escribiendo:
{{matlab|codigo=
%añadir al programa anterior
xe=m1e\tr.' %definimos x como el vector con las soluciones del sistema
ye=x(1)+x(2)*w+x(3)*exp(w) %ecuación exponencial calculada por mínimos cuadrados
hold on
plot(w,ye,'b-')
legend('datos','exponencial de minimos') %añadimos a la gráfica anterior la parábola junto con una leyenda
}}

=== Cálculo del error cuadrático medio de la parábola ===

Queremos tener una referencia del error cometido por la aproximación, por lo que calculamos el error cuadrático medio de la siguiente forma
{{matlab|codigo=
%añadir al programa anterior
error=sum((tr-ye).^2)/20
}}

=== Aproximación del valor 16 golpes/min ===

Para saber qué valor predice la recta para una frecuencia de 16 golpes por minuto ejecutamos:
{{matlab|codigo=
%añadir al programa anterior
ze=x(1)+x(2)*16+x(3)*exp(16)
}}

=== Resultado ===

una vez escrito el código lo ejecutamos y Matlab nos da los siguientes resultados:

xe=[1.025976427703857e+03;-1.974268082492748;-2.319344866513919e-09] por lo que la ecuación exponencial es tr=1025.97-1.97*w-2.319e-9*exp(w)

errore=1.908149369400242e+03 el error cuadrático medio vale 1908.1

ze=9.943675284291536e+02 para una frecuencia de 16 golpes por minuto se estima un tiempo de rotura de 994.36 horas y además genera la siguiente gráfica:
[[Archivo:Gráfica_exponencial_minimos_grupo_23.png|600px|thumb|left|frame|gráfica de los datos de partida junto con la ecuación exponencial de mínimos cuadrados]]

== Conclusiones ==
La distribución con el menor error cuadrático medio es la función exponencial, además es la que más coincide, visualmente, con los datos de partida por lo que nos parece que 1098.1 es la predicción más fiable de las 3

Archivo:Gráifca parabola minimos grupo 23.png

2019-11-11T11:49:42Z

C.novero:

Archivo:Gráfica recta mínimos grupo 23.png

2019-11-11T11:49:21Z

C.novero:

Archivo:Gráfica exponencial minimos grupo 23.png

2019-11-11T11:48:53Z

C.novero: