2009-数学类作业(综合实验).pptVIP

计算 * 某厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60台、80台．每季度的生产费用为 （元），其中x是该季生产的台数．若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c元．已知工厂每季度最大生产能力为100台，第一季度开始时无存货，设a=50、b=0.2、c=4，问工厂应如何安排生产计划，才能既满足合同又使总费用最低．讨论a、b、c变化对计划的影响，并作出合理的解释． 综合题1 综合题2 绘制山区地貌图 要在某山区方圆大约27平方公里范围内修建一条公路，从山脚出发经过一个居民区，再到达一个矿区。横向纵向分别每隔400米测量一次，得到一些地点的高程：(平面区域0=x=5600,0=y=4800)，需作出该山区的地貌图和等高线图。 Matlab shanqu 返回 综合题3 用 MATLAB 求解问题时，一般要经历建模和编程两个过程，只有在建模正确的前提下，方能得出正确的结果。 一、单自由度系统有阻尼自由振动 1.建立计算模型 由动力学可知，单自由度有阻尼自由振动的振动方程为： 无量刚化后有： 其中 , 上述方程的解为： 其中 x0 表示初始位置， v0 表示初始速度。 综合题4 参数ωn ＝10， x0 ＝1， v0 ＝0，计算的终止时间 t=2。试求ξ 从 0.1 到 1运动方程的解，并画出波形。 2. MATLAB 编程 编写 M 文件 ex1.m 新建一个图形窗口，绘制三维网格图 如果改变初始条件令x0＝0，v0＝1，其运动曲线实际上就是系统的脉冲过渡函数。 二、气体分子运动的麦克斯韦分布曲线 通过本例说明如何用复杂的数学公式绘制曲线。 利用气体分子运动的麦克斯韦速度分布律，求氯分子运动的速度分布曲线，并讨论温度T及分子量mu对速度分布曲线的影响。 1. 建立计算模型 麦克斯韦速度分布律为: 其中，m---分子质量, m=mu/NA, mu---分子量, NA---阿伏加德罗数 k---波尔茨曼常数 T----气体的绝对温度 v----分子速度 综合题5 为研究单个参数的影响，先把麦克斯韦分布律编为一个函数子程序，以便重复调用，同时将常数项也放在子程序中。 需要强调的是：子程序不得与主程序放在同一个 M 文件中，只能将子程序单独做成 M 文件，并放在与主程序同一个工作路径中。 2. MATLAB 编程 首先建立计算麦克斯韦分布律的子程序 mxw.m 给出T和mu的 绘制分布曲线 为了看出不同的T和mu对曲线形状的影响，再次给定T和mu，在同一幅图中绘制分布律曲线的图形 三、方波的分解 在连续信号系统中，方波可以用相应频率的基波及其奇次谐波合成，这也是将方波展开为正弦级数的出发点。本节将演示这一现象。 1. 建立计算模型 一个以原点为奇对称中心的方波 y(t)可以用奇次正弦波的叠加来逼近： 方波的宽度为π，周期为 2π。 综合题6 2. MATLAB编程 建立 M 文件 ex3.m % 演示基波和奇次谐波合成方波 t=0:0.1:10; %首先设定一个有101个点的时间数组 %绘制频率w＝1（ f=1/2π）的正弦基波，并设置暂停 %叠加3次谐波，绘图并设置暂停 %叠加1、3、5、7、9次谐波，绘图并设置暂停 %为了绘制三维曲面，需要将各次波形数据存储为一个三维数组，因此需要重新定义y，重新编程，本例将求至 19 次谐波 %将各个波形叠合绘出，并设置暂停 %将各个波形绘制成三维网格图 考虑人工捕获 设表示捕获能力的系数为e，相当于食饵的自然增长率由r1 降为r1-e，捕食者的死亡率由r2 增为 r2+e 设战前捕获能力系数e=0.3, 战争中降为e=0.1, 则战前与战争中的模型分别为: 综合题7 模型求解: 1、分别用m-文件shier1.m和shier2.m定义上述两个方程 2、建立主程序shark1.m, 求解两个方程，并画出两种情况下鲨鱼数在鱼类总数中所占比例 x2(t)/[x1(t)+x2(t)] 实线为战前的鲨鱼比例，“*”线为战争中的鲨鱼比例 结论：战争中鲨鱼的比例比战前高！ 返 回 地中海鲨鱼问题 意大利生物学家Ancona曾致力于鱼