Técnicas Experimentales de Termodinámica. Curso 2009-2010
Práctica 19. Temperatura de Curie de la aleación monel    


Para ejecutar un comando, posiciona sobre él el cursor  y aprieta mayúsculas - ENTER.

Primero introduce todos los datos experimentales en forma de pares de valores {V/mV,t/ºC}

datosVt={{65.2,21.0},{53.66,23.0},{44.47,24.0},{37.60,25.0},{33.80,26.0},{29.16,27.0},{26.94,28.0},{23.02,30.0},{20.89,32.0},{19.33,34.0},{18.51,36.0}, {17.99,38.0},{17.62,40.0}};Needs["NonlinearRegression`"]

Después elige qué datos vas a emplear en el ajuste (en el mismo formato anterior) . Conviene eliminar algunos de los extremos y una vez que este programa consiga hacer bien el ajuste puedes probar a incluir alguno de los datos eliminados y repetir el ajuste. ¡ENSAYO Y ERROR!

datos={{53.66,23.0},{44.47,24.0},{37.60,25.0},{33.80,26.0},{29.16,27.0},{26.94,28.0},{23.02,30.0},{20.89,32.0},{19.33,34.0},{18.51,36.0}, {17.99,38.0},{17.62,40.0}};Clear[tc,hr,k,Vinf];FindFit[datos,-273.15+(tc+273.15)(hr+k (V-Vinf))/ArcTanh[k (V-Vinf)],{{tc,25},{hr,.001},{k,.001},{Vinf,15}},V]

{tc,hr,k,Vinf}={tc,hr,k,Vinf}/.%;

Show[Plot[-273.15+(tc+273.15) (hr+k (V-Vinf))/ArcTanh[k (V-Vinf)],{V,Vinf+.01,87},PlotRange->{{0,80},{10,50}},AxesLabel->{"V/mV","t /°C"}],ListPlot[datosVt]]

datostV=Transpose[RotateLeft[Transpose[datosVt],1]];Show[ParametricPlot[{-273.15+(tc+273.15) (hr+k (V-Vinf))/ArcTanh[k (V-Vinf)],V},{V,Vinf+.01,80},PlotRange->{{10,50},{0,80}},AxesLabel->{"t /°C","V/mV"}],ListPlot[datostV],AspectRatio->1/GoldenRatio]

Clear[tc,hr,k,Vinf];NonlinearRegress[datos,-273.15+(tc+273.15)(hr+k (V-Vinf))/ArcTanh[k (V-Vinf)],{{tc,100},{hr,.001},{k,.001},{Vinf,15}},V,RegressionReport->{ParameterTable}]

Para solicitar ayuda sobre el uso de este notebook Mathematica puedes contactarme en : Jose.A.Manzanares@uv.es