B øúB^Zã@s<ddlZddlmZddlZddlmZ dd„Z dd„Z dS)éNcCs$d}|||t | |¡}|S)Né)ÚnpÚexp)ÚtÚT_0ÚrÚT_sÚT©r úD/Users/vegardjervell/Documents/6_semester/matmod/matmod4/fit_r_t0.pyr sr c Cs t ¡\}}tdƒd\}}|d|d}}xPt||ƒdkr~t||ƒdkr~tjt||||gd\}}||}}|\}}q0Wt t |ƒt |ƒd¡}tdt |d ƒƒtd t |d ƒƒt|||ƒ} t j ||d d t j|| dd t  d¡t  d¡t  ¡t  ¡dS)NzFitting initial temperature and heat exchange coefficient in analytical model to best fit experimental data using scipy.optimize.curve_fit())éZgš™™™™™¹?égñh㈵øä>g-Cëâ6?)Úp0ièzT0 = ézr = ézExperimental data)ÚlabelzBest fit analytical solutionz Time [min]u Temp [°C])ÚflÚ les_real_cupÚprintÚabsÚoptZ curve_fitr rÚlinspaceÚminÚmaxÚroundÚpltÚscatterÚplotÚxlabelÚylabelÚlegendÚshow) ÚtidÚtempÚT0rZold_T0Zold_rÚcoeffÚcovZtid_akseZTempr r r Ú finn_r_T0s$ "     r') Ú fillesingrZscipy.optimizeÚoptimizerÚnumpyrÚmatplotlib.pyplotÚpyplotrr r'r r r r Ús