




已阅读5页,还剩3页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
河北大学河北大学 数学模型数学模型 实验实验 实验报告实验报告 一 实验目的一 实验目的 1 学会编写程序段 2 能根据 m 文件的结果进行分析 3 根据图像进行比较和分析 二 实验要求二 实验要求 8 18 1 捕鱼业的持续收获捕鱼业的持续收获 运行下面的 m 文件 并把相应结果填空 即填入 clear clc 无捕捞条件下单位时间的增长量 f x rx 1 x N 捕捞条件下单位时间的捕捞量 h x Ex F x f x h x rx 1 x N Ex 捕捞情况下渔场鱼量满足的方程 x t F x 满足 F x 0 的点 x 为方程的平衡点 求方程的平衡点 syms r x N E 定义符号变量 Fx r x 1 x N E x 创建符号表达式 x solve Fx x 求解 F x 0 求根 得到两个平衡点 记为 x0 x1 x0 x 2 x1 x 1 符号变量 x 的结构类型成为 求 F x 的微分 F x syms x 定义符号变量 x 的结构类型为 dF diff Fx x dF simple dF 简化符号表达式 得 F x 求 F x0 并简化 dFx0 subs dF x x0 将 x x0 代入符号表达式 dF dFx0 simple dFx0 得 F x0 求 F x1 dFx1 subs dF x x1 得 F x1 若 E r 有 F x0 0 故 x0 点稳定 x1 点不稳定 根据平衡点稳定 性的准则 若 E r 则结果正好相反 在渔场鱼量稳定在 x0 的前提下 E r 求 E 使持续产量 h x0 达到最大 hm 通过分析 见教材 p216 图 1 只需求 x0 使 f x 达到最大 且 hm f x0 syms r x N fx r x 1 x N df diff fx x x0 solve df x 得 x0 hm subs fx x x0 得 hm 又由 x0 N 1 E r 可得 E 产量模型的结论是 将捕捞率控制在固有增长率的一半 E r 2 时 能够获得最大的持续产量 8 28 2 种群的相互竞争 种群的相互竞争 1 1 补充如下指出的程序段 然后运行该 m 文件 对照教材上的相应结果 clear clc 甲乙两个种群满足的增长方程 x1 t f x1 x2 r1 x1 1 x1 N1 k1 x2 N2 x2 t g x1 x2 r2 x2 1 k2 x1 N1 x2 N2 求方程的平衡点 即解代数方程组 f x1 x2 0 g x1 x2 0 8 38 3 种群的相互竞争 种群的相互竞争 2 2 补充如下指出的程序段 然后运行该 m 文件 对照教材上的相应结果 clear clc 甲乙两个种群满足的增长方程 x1 t f x1 x2 r1 x1 1 x1 N1 k1 x2 N2 x2 t g x1 x2 r2 x2 1 k2 x1 N1 x2 N2 求方程的平衡点 即解代数方程组 f x1 x2 0 g x1 x2 0 编写出该程序段 三 实验内容三 实验内容 8 18 1 捕鱼业的持续收获捕鱼业的持续收获 产量模型产量模型 文件名 p178 m clear clc 无捕捞条件下单位时间的增长量 f x rx 1 x N 捕捞条件下单位时间的捕捞量 h x Ex F x f x h x rx 1 x N Ex 捕捞情况下渔场鱼量满足的方程 x t F x 满足 F x 0 的点 x 为方程的平衡点 求方程的平衡点 syms r x N E 定义符号变量 Fx r x 1 x N E x 创建符号表达式 x solve Fx x 求解 F x 0 求根 得到两个平衡点 记为 x0 N 1 x N x1 0 x0 x 2 x1 x 1 符号变量 x 的结构类型成为 求 F x 的微分 F x syms x 定义符号变量 x 的结构类型为 dF diff Fx x dF simple dF 简化符号表达式 得 F x 求 F x0 并简化 dFx0 subs dF x x0 将 x x0 代入符号表达式 dF dFx0 simple dFx0 得 F x0 E r 求 F x1 r E dFx1 subs dF x x1 得 F x1 r E 若 E r 有 F x0 0 故 x0 点稳定 x1 点不稳定 根据平衡点稳定性的准 则 若 E r 则结果正好相反 在渔场鱼量稳定在 x0 的前提下 E r 求 E 使持续产量 h x0 达到最大 hm 通过分析 见教材 p178 图 1 只需求 x0 使 f x 达到最大 且 hm f x0 syms r x N fx r x 1 x N a df diff fx x x0 solve df x 得 x0 N 2 hm subs fx x x0 得 hm r N 4 又由 x0 N 1 E r 可得 E r 2 产量模型的结论是 将捕捞率控制在固有增长率的一半 E r 2 时 能够获得最大的持续产量 提示 符号简化函数 simple 的格式 simple S 对符号表达式 S 尝试多种不同的算法简化 以显示 S 表达式的长度最短的简化形式 变量替换函数 sub 的格式 Subs S OLD NEW 将符号表达式 S 中的 OLD 变量替换为 NEW 变量 8 28 2 种群的相互竞争种群的相互竞争 文件名 p186 m clear clc 甲乙两个种群满足的增长方程 x1 t f x1 x2 r1 x1 1 x1 N1 k1 x2 N2 x2 t g x1 x2 r2 x2 1 k2 x1 N1 x2 N2 求方程的平衡点 即解代数方程组 f x1 x2 0 g x1 x2 0 得 4 个平衡点 P 1 P1 N1 0 P 2 P2 0 N2 P 3 P3 N1 1 k1 1 k2 k1 N2 1 k2 1 k2 k1 P 4 P4 0 0 平衡点位于第一象限才有意义 故要求 P3 k1 k2 同时小于 1 或同时大于 1 判断平衡点的稳定性 参考教材 p200 fx1 diff f x1 fx2 diff f x2 gx1 diff g x1 gx2 diff g x2 A fx1 fx2 gx1 gx2 syms x1 x2 p subs fx1 gx2 x1 x2 P 1 P 2 p simple p 简化符号表达式 p q subs det A x1 x2 P 1 P 2 q simple q P p q 得到教材 p186 表 2 的前 3 列 经测算可得该表的第 4 列 即稳定条件 8 38 3 种群的相互竞争 种群的相互竞争 2 2 求微分方程组 1 1 2 2 1 1 222 2 2 2 1 1 1 11 1 N x N x xrtx N x N x xrtx 的数值解 分别画出教材 p189 中的图 5 1 图 5 2 图 5 3 有关数据参见教材 p188 中 计算与验证 提示 1 求微分方程组的数值解可参考教材 p139 的程序 2 在 figure 1 中画图 5 1 在 figure 2 中画图 5 2 在 figure 3 中画图 5 3 在 程序中 figure 图形编号 用于定位对应图形 3 使用 text x y 标识文本 坐标点 x y 在 标识文本 的左边 调整 x y 值 使 标识文本 放在图中的适当位置 4 用 axis xmin xmax ymin ymax 控制坐标的刻度范围 5 用 grid on 打开网格 grid off 关闭网格 6 用 hold on 把要画的图形保持在之前在同一 figure 上所画的图形中 同一坐标系 7 图 5 3 中的两 点线 直线 一条的两个端点为 0 1 和 1 0 另一条的两个端点 为 0 2 和 1 6 0 四 实验结果及其分析四 实验结果及其分析 8 18 1 捕鱼业的持续收获捕鱼业的持续收获 产量模型产量模型 8 28 2 种群的相互竞争 种群的相互竞争 1 1 1 实验代码 syms r1 r2 N1 N2 k1 k2 x1 x2 f r1 x1 1 x1 N1 k1 x2 N2 g r2 x2 1 k2 x1 N1 x2 N2 x1 x2 solve f g x1 x2 P x1 x2 fx1 diff f x1 fx2 diff f x2 gx1 diff g x1 gx2 diff g x2 A fx1 fx2 gx1 gx2 syms x1 x2 p subs fx1 gx2 x1 x2 P 1 P 2 p simple p q subs det A x1 x2 P 1 P 2 q simple q P p q 2 实验结果 8 38 3 种群的相互竞争 种群的相互竞争 2 2 1 实验代码 脚本文件 function f fun t x r1 2 5 r2 1 8 N1 1 6 N2 1 k1 0 5 k2 1 6 f r1 x 1 1 x 1 N1 k1 x 2 N2 r2 x 2 1 k2 x 1 N1 x 2 N2 t x ode45 fun 0 8 0 1 0 1 plot t x 1 t x 2 axis 0 8 0
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 展位转合同(标准版)
- 电缆材料知识培训总结
- 2025标准保密合同范本
- 2025短期合同工简易劳动合同模板
- 电线电缆设备知识培训内容课件
- 4.13 罗斯福新政 说课稿 2025-2026学年统编版九年级历史下册
- 三年级英语下册 Module 4 Unit 2 Does Lingling like oranges说课稿2 外研版(三起)
- 二年级下册书法教学计划和教案(2025-2026学年)
- 金融借款合同
- 废料处理合同
- 2025年旅游策划服务框架协议书
- 燃气电厂培训计划
- 2025年人民教育出版社有限公司招聘笔试参考题库含答案解析
- 太空作战指挥控制-洞察分析
- KCA试题库完美版
- 方竹种植和管理技术课件
- 成考语文成人高考(高起专)试题及解答参考(2025年)
- 华为智慧矿山解决方案
- 工业锅炉安装质量证明书
- 铺面装修购销合同模板
- 无人机兴趣班课件:第1节无人机介绍与现代无人机知识
评论
0/150
提交评论