华北科技学院。
上机报告 系(部基础部。
专业、班级计算b121
姓名房超学号 201209014119
课程名称数学建模算法与应用。
上机题目数学模型。
任课教师谭立云。
指导教师谭立云。
成绩(优、良、中、及格、不及格)
华北科技学院基础部。
一.上机目的。
第15章典型例题。
二.上机环境。
等。三.上机内容。
第十五章**方法习题。
15.2 已知兰彻斯特的游击战争模型如下:
其中:参数c,d,α,的值未知,现在有连续20天的观测数据如表15.2所列,拟合参数c,d,α,
解:微分方程对应的差分方程为。
可以改写成。
cx(k)y(k) -x(k) 00 d
00 -x(k)y(k) -y(k) α
上述所有的差分方程可以改写成矩阵形式,即。
-x(1)y(1) -x(1) 00c
d-x(19)y(19) -x(19) 00
00 -x(1)y(19) -y(1
00 -x(19)y(19) -y(19
利用最小二乘法,求得c=0,α=0.0393,d=0,β=0.0931。
拟合的matlab程序如下:
clc,clear
a=[1500 1400 1320 1100 1000 950 880 800 700 680
x=a([1,4],:x=x';x=x(:)
y=a([2,5],:y=y';y=y(:)
dx=diff(x);%求一阶向前差分。
dy=diff(y);
temp=x(1:end-1).*y(1:end-1);%构造分块矩阵的子矩阵。
a=[-temp -x(1:end-1) zeros(19,2);zeros(19,2),-temp -y(1:end-1)];
b=[dx;dy];
solution=a\b
运行结果:solution =
数学建模上机实验报告
西华大学数学建模上机实验报告。一 实验目的。1.掌握建立和执行m文件的方法 2.掌握实现选择结构的方法 3.掌握实现循环结构的方法。二 实验内容。1 几种常用循环语句的练习。1 对n 1,2,10,求的值,并记录在x中。2 设银行年利率为11.25 将10000元钱存入银行,问多长时间会连本带利翻一...
数学建模上机实验报告
西华大学数学建模上机实验报告。一 实验目的。熟悉matlab软件的一些功能以及编程的过程,同时能解决部分的基本问题,掌握相关语法的应用等等。二 实验内容。对于此问题,编写函数m文件 1 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。2 有一个矩阵,编程求出其最大值及其所处的位置...
《数学建模与数学实验》上机报告
数学建模与数学实验 上机报告 第 1 次 一 上机训练目的 题目或内容 简述综述 等。题目一 数学软件 mathtype5.2 matlab maple mathematica4.0 lingo8.0 安装调试 基本命令使用 变量赋值 定义函数 过程控制 绘图命令 拟合 线性规划 非线性规划 高等数...