实验指导书(2012.4)new.doc

《数学实验》实验指导书 龚 劬 温罗生 肖 剑 重庆大学数学实验教学示范中心 2010年8月目 录 预备实验——桥梁分析 3 实验1 MATLAB软件入门 8 实验2 方程模型及其求解算法 24 实验3 收敛与混沌——迭代 28 实验4 微分方程模型、求解及稳定性分析 31 实验5 插值方法 34 实验6 数据拟合及参数辨识方法 37 实验7 回归分析模型、求解及检验 40 实验8 连续系统与离散系统的计算机模拟 43 实验9 线性规划模型、求解及灵敏度分析 45 实验10 非线性规划与多目标规划模型及其求解 49 实验11 如何表示二元关系—图的模型及矩阵表示 52 实验12 改进技术的最佳实施问题——综合实验 55 实验13 人口增长模型及其数量预测——综合实验 57 实验14 River-bay系统水污染问题_____综合实验 59 实验15 炮弹发射角的确定———综合实验 61 实验16 探究实验 62 实验17 开采沙子——综合实验 63 实验18 海水中提取淡水——综合实验 66 实验19 氯仿污染——综合实验 70 实验20 机动车尾气排放——综合实验 80 实验21 计算机断层扫描图像——综合实验 88 实验22 一个有年龄结构的人口增长问题——综合实验 96 预备实验——桥梁分析 教学目的和要求： 通过桥梁分析问题，使学生： 了解线性代数在土木工程中的应用； 了解如何通过做一些使问题简化的假设，建立实际问题的数学模型； 体会学好线性代数知识的重要性； 激发学习线性代数的兴趣。 知识点：线性方程组 向量分解 必备技能： 1. 力的平衡分析； 2. 向量分解； 3. 求解线性方程组。 主要内容 应用场景 问题分析 建立数学模型 实验任务 1．应用场景 解方程组在许多领域都有应用。下面给出一个在土木工程中的应用例子，虽然加入了一些幽默元素，但类似的情形土木工程师会经常遇到。 一位货运司机正驾着卡车为一个数学家聚会运送物资，但他的卡车超载了。途中，他遇到了一个古老的桁架桥  3.建立数学模型 为了判断该桥是否会断裂，可以用节点法。假设，当桥中桁架的任何组件断裂时，整座桥就会断裂。节点法是分析桁架每个节点的平衡条件。节点平衡要求该点所受到的所有内力和外力相互平衡，否则，桁架将会移动。 如图3所示，有三个未知的外力,和，作用在桁架与地面接触的节点上。这三个外力与4280磅和5270磅的负载共同作用导致桁架组件的内力。桁架每个组件的内力可以通过对每个节点运用节点法来确定。图3是节点法在每个节点应用的示意图。注意各组件已经被“割断”，假设内力方向如箭头所示。被“割断”的组件的所有力都按张力来标示，所谓张力就是作用在桁架组件的拉力。 为了把卡车的重量从桥的中间转化为桥两端的支撑力，一些组件一定会被压缩，另一些一定会被拉伸。所以如果最终得出的力为负，则说明假设为张力是错误的，该力实际上应该是压力。在图3中具有负值的力的箭头应向内。 为了使桁架满足平衡条件，在每一节点处，x 轴方向和y轴方向的合力应分别为0。对应的数学表达式为 和。大小为4280 lb的力作用在B点和大小为5720 lb的力作用在C点说明，卡车重量直接作用在节点处。让我们来看一下如何根据节点法分析任一节点处的受力情况 节点A 每个桁架用桁架两端节点处的字母表示，如图3，对支点A有四个力的作用, , , 和: ·表示作用在A点作用在x轴方向上的力，表示从A点指向B点的力。 · 表示力的方向为从A指向E。 注意，当分析A点时，将力分解在X和Y方向上。 图4是对力的分析，演示了如何对力进行分解。在在Y方向上的分力： 在X轴方向的分力大小： 为了使A点保持平衡状态，该点处X轴方向和Y轴方向的合力分别为0。A点处Y轴方向的合力： 接着你需要求出A点处X轴方向的合力。然后，通过在每一个支点处重复这样的过程，你可以建立一个方程组以求得各个桁架组件间未知的力。这个方程组有唯一解（因为未知数的个数和方程的个数相同）。 4. 任务 1）确定作用桁架为了评估桁架的安全，需要确定土木工程专业实验1 MATLAB软件入门 实验目的及意义 [1] 熟悉MATLAB软件的用户环境； [2] 了解MATLAB软件的一般目的命令； [3] 掌握MATLAB数组操作与运算函数； [4] 掌握MATLAB软件的基本绘图命令； [5] 掌握MATLAB语言的几种循环、条件和开关选择结构。 通过该实验的学习，使学生能灵活应用MATLAB软件解决一