Mathematica软件及其应用.ppt

* * 动画演示3： For[i=1,i=16,i=i+1, ParametricPlot3D[{Sin[v]*Cos[u],Sin[v]*Sin[u],v}, {u,0,2Pi*i/16},{v,0,Pi}, AspectRatio-1,PlotRange-{{-1,1},{-1,1},{0,Pi}}] ] * * 动画演示4： For[i=0,i2Pi,i=i+Pi/4, Plot3D[Sin[Sqrt[x^2+y^2]+i],{x,-10,10},{y,-10,10}, Axes-False,Boxed-False, PlotRange-{-2,2},PlotPoints-30] ] 或用Table产生图形序列： Table[Plot3D[Sin[Sqrt[x^2+y^2]+i],{x,-10,10},{y,-10,10}, Axes-False,Boxed-False, PlotRange-{-2,2},PlotPoints-30],{i,0,2Pi,Pi/4} ] * * (2) 动画的保存 可用导出命令Export将若干个图形保存为GIF格式文件。 Export[“文件名.gif”,图形表,ConversionOptions-{选项}] 常用选项： “Loop”-True：是否循环播放 “AnimationDisplayTime”-t：每帧显示时间(t秒) * * 如： For[i=1,i2Pi,i=i+Pi/4, t[i]=Plot3D[Sin[Sqrt[x^2+y^2]+i], {x,-15,15},{y,-15,15}, Axes-False,Boxed-False, PlotRange-{-4,4},PlotPoints-40] ] Export[D:\dh1.gif, Table[t[i],{i,1,2Pi-0.1,Pi/4}], ConversionOptions-{Loop-True, AnimationDisplayTime-0.02}] * * 或： p=Table[Plot3D[Sin[Sqrt[x^2+y^2]+i], {x,-15,15},{y,-15,15}, Axes-False,Boxed-False, PlotRange-{-4,4},PlotPoints-40, ViewPoint-{2,1,8}], {i,0,2Pi,Pi/4}] Export[D:\dh2.gif,p, ConversionOptions-{Loop-True, AnimationDisplayTime-0.02}] * * 第四章 数学建模中的应用 一、梯子的最大长度 如图，室内一条a=1米宽的通道与b=2米宽的通道相交成直角，一个梯子需水平绕过拐角，试问梯子的最大长度是多少？ * * 设梯子长度为L，则 ，问题 转化为求L(x)在区间 上的最小值。 先画出函数L(x)的图形，以观察极小值点的大概位置，编写程序如下： a=1;b=2;L=a/Cos[x]+b/Sin[x]; Plot[L,{x,0,Pi/2}] 由图可知，极小值点在 x =1 附近。 * * 由程序的计算结果可知能通过拐角的最大梯子长度为4.16米。 下面求L(x)在 x =1 附近的极小值，编写程序如下： m=FindMinimum[L,{x,1}] Print[梯子的最大长度是：,m[[1]],米] 思考：1、若拐角不是直角，而是120°，则能通过拐角的最大梯子长度为多少米？ 2、若拐角为直角，将梯子换为宽度为 k = 0.4 米的手推车，则车的最大长度为多少米？ * * 二、空中电缆的长度 在相距100米的两个塔上高度相等的两个点之间悬挂一根电缆，在电缆的中间下垂10米，计算两塔之间电缆的长度。 以电缆的最低点为坐标原点建立如图坐标系。 由条件可知，点(-50,0)、 (0,0)、(50,0)在电缆形成的曲线上，用抛物线模拟电缆曲线，然后计算这段抛物线的长度。 * * 通过数据拟合求出抛物线的方程，再利用定积分求曲线的长度。编写程序如下： d={{-50,10},{0,0},{50,10}}; y=Fit[d,{1,x,x^2},x]; L=Integrate[Sqrt[1+(D[y,x])^2],{x,-50,50}]; Print[L