数学建模实验报告解析.docVIP

南通大学实验报告 系 （院）： 理学院 专 业： 信息与计算科学 班 级： 131 学 号： 1302022023 姓 名： 李泽奥 课程名称： 数学实验与数学建模 指导教师： 刘晓惠 实验名称 数列极限与函数极限 实验目的 从刘徽的割圆术、裴波那奇数列研究数列的收敛性并抽象出极限的定义；理解数列收敛的准则；理解函数极限与数列极限的关系。 实验内容 割圆术 中国古代数学家刘徽在《九章算术注》方田章圆田术中创造了割圆术计算圆周率。刘徽先注意到圆内接正多边形的面积小于圆面积；其次，当将边数屡次加倍时，正多边形的面积增大，边数愈大则正多边形面积愈近于圆的面积。 “割之弥细，所失弥少。割之又割以至不可割，则与圆合体而无所失矣。”这几句话明确地表明了刘徽的极限思想。 以表示单位圆的圆内接正多边形面积，则其极限为圆周率。用下列Mathematica程序可以从量和形两个角度考察数列{}的收敛情况： m=2;n=15;k=10; For[i=2,i=n,i++, l[i_]:=N[2*Sin[Pi/(3*2^i)],k]; (圆内接正多边形边长) s[i_]:=N[3*2^(i-1)*l[i]*Sqrt[1-(l[i])^2/4],k]; (圆内接正多边形面积) r[i_]:=Pi-s[i]; d[i_]:=s[i]-s[i-1]; Print[i, ,r[i], ,l[i], ,s[i], ,d[i]] t=Table[{i,s[i]},{i,m,n}] (数组) ListPlot[t] (散点图) —————————————————————————————————————— 裴波那奇数列和黄金分割 由有著名的裴波那奇数列。 如果令，由递推公式可得出 ，； 。用下列Mathematica程序可以从量和形两个角度考察数列{}的收敛情况：n=14,k=10; For[i=3,i=n,i++, t1=(Sqrt[5]+1)/2; t2=(1-Sqrt[5])/2; f[i_]:=N[(t1^(i+1)-t2^(i+1))/Sqrt[5],k]; (定义裴波那奇数列通项) rn=(5^(1/2)-1)/2-f[i-1]/f[i];Rn=f[i-1]/f[i];dn=f[i-1]/f[i]-f[i-2]/f[i-1]; Print[i, ,rn, ,Rn, ,dn]; t=Table[{i,f[i-1]/f[i]},{i,3,n}] ListPlot[t] 3、函数极限与数列极限的关系 用Mathematica程序 m=0;r=10^m;x0=0; f[x_]=x*Sin[1/x] Plot[f[x],{x,-r,r}] Limit[f[x],x-x0] 结论 1、 数据：（m=2;n=15;k=10;）{{2,3.000000000},{3,3.105828541},{4,3.132628613},{5,3.139350203},{6,3.141031951},{7,3.141452472}, {8,3.141557608},{9,3.141583892},{10,3.141590463},{11,3.141592106},{12,3.141592517},{13,3.141592619}, {14,3.141592645},{15,3.141592651}} 做散点图为： —————————————————————————————————————————————— 数据：（m=2;n=20;k=10;） {10,3.141590463},{11,3.141592106},{12,3.141592517},{13,3.141592619},{14,3.141592645},{15,3.141592651}, {16,3.141592653},{17,3.141592653},{18,3.141592654},{19,3.141592654},{20,3.141592654}} 做散点图为： 数据：（m=2;n=50;k=30;） {40,3.14159265358979323846264148344},{41,3.14159265358979323846264290832}, {42,3.14159265358979323846264326454},{43,3.14159265358979323846264335359}, {44,3.141