matlab生猪的出售时机,利用MATLAB进行验证性实验1划艇比赛的成绩2.汽车刹车距离生猪的出售时机模型求解.doc...

河北大学 《数学模型》实验实验报告

班级专业

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=时的