Matlab解决弗兰克—赫兹实验的实验数据处理.doc

《数学实验》报告 实验名称?? ?MATLAB与数学实验_ 学 院 专业班级??? 姓??名??? ?? 学??号??? ? ?? ?? 年 月 工科物理实验 一、【作业题目】 弗兰克—赫兹实验的实验数据处理； 二、【作业程序】 x=0:0.5:53; y=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.001,0.012,0.027,0.051,0.068,0.085,0.095,0.105,0.113,0.121,0.126,0.127,0.128,0.125,0.12,0.111,0.101,0.088,0.075,0.062,0.048,0.037,0.027,0.021,0.019,0.023,0.041,0.074,0.114,0.154,0.183,0.211,0.231,0.242,0.247,0.243,0.234,0.214,0.191,0.158,0.128,0.095,0.078,0.065,0.054,0.047,0.04,0.039,0.056,0.084,0.148,0.213,0.263,0.313,0.347,0.372,0.388,0.39,0.381,0.364,0.332,0.289,0.252,0.198,0.146,0.102,0.076,0.054,0.045,0.044,0.045,0.066,0.142,0.221,0.295,0.348,0.405,0.437,0.47,0.492,0.5,0.498,0.483,0.445,0.417,0.356,0.285]; p2=polyfit(x,y,2) p7=polyfit(x,y,7) p10=polyfit(x,y,10) p20=polyfit(x,y,20) disp(二次拟合函数),f2=poly2str(p2,x) disp(七次拟合函数),f7=poly2str(p7,x) disp(十次拟合函数),f10=poly2str(p10,x) disp(二十次拟合函数),f20=poly2str(p20,x) x1=0:0.5:53; y2=polyval(p2,x1); y7=polyval(p7,x1); y10=polyval(p10,x1); y20=polyval(p20,x1); plot(x,y,rp,x1,y2,--,x1,y7,bd,x1,y10,y-.,x1,y20,k-) 三、【作业结果】 p2 = 0.0001 0.0023 0.0025 p7 = -0.0000 0.0000 -0.0000 0.0001 -0.0024 0.0230 -0.0854 0.0648 p10 = -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0001 -0.0014 0.0112 -0.0458 0.0727 -0.0202 p20 = Columns 1 through 13 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0001 Columns 14 through 21 -0.0012 0.0084 -0.0398 0.1268 -0.2473 0.2547 -0.1007 0.0027 二次拟合函数 f2 = 8.8472e-05 x^2 + 0.0023245 x + 0.0024688 七次拟合函数 f7 = -2.4449e-10 x^7 + 4.5621e-08 x^6 - 3.3585e-06 x^5 + 0x^4- 0.0023934 x^3 + 0.022986 x^2 - 0.085418 x + 0.06481 十次拟合函数 f10 = -8.8668e-14 x^10 + 2.1806e-11 x^9 - 2.2792e-09 x^8 + 1.3223e-07 x^7- 4.6647e-06 x^6 + 0x^5 - 0.0014065 x^4 + 0.011236 x^3- 0.045832 x^2 + 0.072654 x - 0