mathematica(高等数学实验室word版)全套.doc

准备实验 实验1.1 进入高等数学实验室 实验1.2 变量与函数 实验1.3 分段函数 实验1.4 一元函数的图形 实验1.5 二元函数的图形 实验1.6 表的使用 实验1.7 多项式有理式及方程求解 实验1.8 过程设计  计算实验 实验2.1 极限与导数的计算 实验2.2 泰勒展开式的计算 实验2.3 积分的计算 实验2.4 向量的运算 实验2.5 解常微分方程 实验2.6 偏导数与全微分 实验2.7 多元函数的极值  体验实验 实验3.1 平面上的几何变换 实验3.2 收敛速度与无穷小的阶 实验3.3 机翼轮廓线的加工 实验3.4 数控机床的刀具补偿 实验3.5 迭代式与不动点 实验3.6 牛顿迭代法 实验3.7 卫星轨道的长度 实验3.8 越野赛问题及其他  应用实验 实验4.1 计算机模拟 实验4.2 日常生活中的经济问题 实验4.3 Fibonacci 数列 实验4.4 超音速飞机与马赫锥 实验4.5 分叉与混沌 实验4.6 π的近似计算 实验4.7 万有引力的发现  高等数学实验室教程 李卫国 实验4.7 万有引力的发现 实验目的 本实验将引领同学沿着当年牛顿思考的足迹再进行一次科学探索的实践. 通过实验使学生理解微积分在科学发现中所起过的巨大作用. 实验指导 1. 历史背景 16世纪后期, 丹麦天文学家布拉赫(T. Brache) 对太阳系中的行星运动进行了长达20年之久的观察, 积累了大量丰富的资料. 他的助手, 德国科学家开普勒 (J. Kepler) 参与了一部分观测工作, 并继承了他的全部观测资料. 在此基础上, 开普勒有进行了20年的观测和研究, 总结出关于行星运动的三大规律. 开普勒第一定律 行星绕太阳运行的轨迹是一个椭圆,太阳的位置 正好位于该椭圆的一个焦点上. 开普勒第二定律 将太阳中心指向行星的有向线段作为向径, 当行星运行时, 向径在单位时间内扫过的面积恒为常数. 开普勒第三定律 太阳系中各行星的椭圆轨道半长轴 a 的立方与其公转周期 T 的平方之比 是一个与行星本身无关的常数. 在此基础上, 英国科学家牛顿运用他所创的力学体系和流数法 (微积分), 通过对开普勒三定律进行精致的数理分析, 运用演绎的方法发现了享誉世界的万有引力定律. 本次试验中, 我们要循着牛顿的脚步, 重做一次这样的数理分析. 2. 极坐标下质点的运动分析 (1) 平面质点的运动规律(向量表示) 质点的运动规律可以用向量函数 表示, 其瞬时速度为, 瞬时加速度为. 将视为由坐标原点O指向动点M的向径. 在极坐标下, t 时刻的转角设为, 则质点在极坐标下的运动方程表为 (1) 其中r(t)为向径的模长, 表示向径的方向(径向), 并用 表示与向径垂直的方向 (剪向). (2)速度向量的分解 若运动方程表为 , 则运动速度 速度向量可以分解为: 径向速度, 剪向速度 . (2) (3)加速度向量的分解 由于径向单位向量的导数为 剪向单位向量的导数为 计算可知, 运动的加速度有如下分解: 径向加速度 (3) 剪向加速度 (4) (4) 质点的受力分析 根据牛顿第二定律, t时刻质点受力为 . 由加速度的分解可知, 质点受力同样分为径向力和剪向力. 径向力 (5) 剪向力: (6) 3.万有引力定律的推导 (1) 利用开普勒第一定律----行星轨道的极坐标方程 根据开普勒第一定律知, 行星绕太阳运行的轨迹是一个椭圆, 太阳的位置正好位于该椭圆的一个焦点上. 设椭圆的长轴为2 a, 短轴为2 b , 焦距为2c, 其中. 以椭圆的一个焦点为极点, 以其长轴为起始方向, 建立极坐标系. 设椭圆在极坐标下的方程为 . 由于椭圆上任意动点到两焦点的距离之和等于2 a, 因此 解出 , 这就是行星椭圆轨道的极坐标方程. 绘制椭圆的图形 r[s_]:= 16/(5-3Cos[s]); ParametricPlot[{r[s] Cos[s], r[s] Sin[s]}, {s,0,2Pi}, PlotStyle-RGBC