河北大学 《数学模型》实验实验报告
班级专业
15 计科 2 班
姓名
张宇轩
学号 006
实验地点
C1-229
指导老师
司建辉
成绩
实验项目
利用 MATLAB进行验证性实验
1.划艇比赛的成绩
2.汽车刹车距离
3.生猪的出售时机模型求解
一、实验目的
学会利用 MATLAB进行验证性实验,熟练掌握用数据拟合求解模型和参数。了解并使用最小二
乘多项式拟合函数polyfit ,仿照案例今后能够自己解决图形问题。
二、 实验要求
1.划艇比赛的成绩的模型:t=αnβ
其中, t 为比赛成绩(时间),n 为桨手人数,α和β为参数。
为适合数据拟合,将模型改为:log t=log α+βlog
桨手人数n比赛平均成绩t
1
2
4
8
1>. 参数 α 和β估计
程序如下:
clear;clc;
n=[1 2 4 8];
t=[];
logt=log(t);logn=log(n);
p=polyfit(logn,logt,1);
beta=p(1)
alfa=exp(p(2))
2>. 实际值与计算值比较(数据比较和和拟合图形)
参考数据结果:
ans =
1
2
4
8
参考图形结果:
图 1:题给拟合图形结果
要求:
运行以上程序。
编程:实际值与计算值比较(数据比较和和拟合图形)。
用 help 查询函数 polyfit 的用法。
2.汽车刹车距离的模型: d=t1v+kv2
其中, d 为刹车距离,变量v 为车速,参数t1 为反应时间,参数k 为比例系数。取经验值
t 1=秒。
实际数据表
车速实际刹车距离
(英里 / 小时)(英尺 / 秒)(英尺)
20
44
30
78
40
124
50
186
60
268
70
372
80
506
①用数据拟合求参数k
为适合数据拟合,将模型改为:y=k
其中 y=/v2
程序如下:
clear;clc;
v=[]; %英尺 / 秒
d=[44 78 124 186 268 372 506]; %最大实际刹车距离(英尺)
y=*v)./v.^2;
k=polyfit(v,y,0)
②用所得模型计算刹车距离和刹车时间(数据比较)
程序如下:
clear;clc;
k= ; %输入上题所求得的结果
v=[]; %英尺 / 秒
d=[44 78 124 186 268 372 506]; %最大实际刹车距离(英尺)
dd=*v+k*v.^2; %计算刹车距离
t=d./v; %计算刹车时间
format short g;
[v',d',round(10*[dd',t'])/10]
③实际和计算刹车距离的比较(拟合图形)
程序如下:
clear;clc;
k= ; %输入题1 所求得的结果
vh=[20 30 40 50 60 70 80]; %英里 / 小时
v=[];% 英尺 / 秒
d=[44 78 124 186 268 372 506];%最大实际刹车距离(英尺)
dd=*v+k*v.^2;%计算刹车距离
plot(vh,d,'r+',vh,dd,'b-');
title(' 实际和计算刹车距离的比较');
axis([20,80,0,510]);
xlabel('v 英里 / 小时 ');
ylabel('d 英尺 ');
要求:
运行以上程序,结果与教材相应内容比较。
题 2 和题 3 中要求输入题 1 所求得的 k 值。
理解程序。
3.生猪的出售时机模型求解
目标函数(生猪出售纯利润,元):
Q(t)=(8-gt)(80+rt)-4t-640
其中, t≥0 为第几天出售, g 为每天价格降低值(常数,元/ 公斤),r 为每天生猪体
重增加值(常数,公斤)。
t 使 Q(t)最大。
图解法
绘制目标函数
Q(t)=(8-gt)(80+rt)-4t-640
的图形( 0≤t ≤20)。其中,g=, r=2。
程序如下:
clear;clc;
g=;r=2;
fplot(@(t)(8-g*t)*(80+r*t)-4*t-640,[0,20]);
grid;
xlabel('t'); ylabel('Q');
②代数法
对目标函数
Q(t)=(8-gt)(80+rt)-4t-640
MATLAB 求 t 使 Q(t)最大。其中, r,g 是待定参数。程序如下:
clear;clc;
syms t;%定义符号变量t
Q=sym('(8-g*t)*(80+r*t)-4*t-640')%建立符号表达式
dQ=diff(Q,'t') %求微分dQ/dt
t=solve(dQ,t)%求 dQ=0 的解 t
r=2;g=;
t=eval(t) %求 r=2,g=时的