




已阅读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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 8.4.1 澳大利亚 说课稿2025-2026学年人教版地理七年级下册
- 小螺号说课稿-2025-2026学年小学音乐人音版五线谱四年级上册-人音版(五线谱)
- Unit 2 Let's play sports教学设计初中英语七年级上册牛津译林版
- 2025-2030基因编辑技术临床应用伦理审查与产业化监管趋势
- 2025-2030基因编辑作物商业化种植监管框架与国际市场准入策略
- 2025-2030基因测序仪技术迭代与临床应用价值评估报告
- 2025-2030基因治疗行业深度分析及技术突破与投资策略报告
- 2025-2030基因治疗临床试验进展与产业化障碍研究
- 2025-2030基因检测技术在儿童天赋潜能开发领域的应用现状与政策风险研判
- 2025-2030基于图像识别的投币洗衣机筒自洁检测系统研发报告
- 化工仪表检修与维护课件
- 2025年华为软件开发工程师招聘面试题库及答案解析
- 程序化广告课件
- 电工基础课件
- 四川省石渠县2025年上半年公开招聘辅警试题含答案分析
- 真菌生物膜毒力因子-洞察及研究
- 基孔肯雅热危害及预防课件
- 副校长在任职宣布会上的表态发言材料
- 同学互助基金管理办法
- 农行对导盲犬管理办法
- 2024年了解脑退化症及各项可用服务和支持的指南-澳大利亚脑退化症支持协会
评论
0/150
提交评论