第二次数学模型的作业

第二次数学模型的作业

第一题 最小二乘法拟合二次曲线

首先，用最小二乘的原理¹简单推导一下。得到要计算的公式。

下面两张图片是我针对这一题的推导（可能有些乱，如果看懂了这个博客应该可以看懂）


由于接下来的矩阵和幂次函数的计算比较麻烦，并且要求画图，下面matlab来做

v=(20:20:140);%
d=[6.5,17.8,33.6,57.1,83.4,118.0,153.5]*(0.001);%单位转化为km
A(1,1)=sum(v.^2);
A(1,2)=sum(v.^3);
A(2,1)=A(1,2);
A(2,2)=sum(v.^4);
C(1,1)=sum(d.*v);
C(2,1)=sum(d.*(v.^2));
B=A\C;
c1=B(1)*3600
c2=B(2)*10^3

X=(20:0.01:140);
D=polyval([B(2),B(1),0],X);
plot(v,d,'r*',X,D,'b');
xlabel('时间(Km/h)','FontSize',16);
ylabel('刹车距离d(km)','FontSize',16);

代码效果是这样的：

所以，由四舍五入可以得到c1的结果大约是0.7s,c2大约是0.01$m{h}^2/k{m}^2$.

第二题 五种模型估算出租车的数量

通过ppt上的结论，我们得到通过五种模型来估算的公式。公式都很简单，这里就不介绍了。看代码，应该也能看懂出来公式是什么。

x=sort(randi([1 8900],3,200),2);
Ans=[];
for j=1:3
Ans(j,1)=round(sum(x(j,:))/200)*2;
Ans(j,2)=2*(round(median(x(j,:))))-1;
Ans(j,3)=x(j,200)+x(j,1)-1;
Ans(j,4)=round((1+1/200)*x(j,200));
Ans(j,5)=round((1+1/(2*200-1))*(x(j,200)-1/(2*200)));
Ans(j,6)=max(Ans(j,1:5))-min(Ans(j,1:5));%%最大相差
end

Ans(j+1,1:5)=max(Ans(:,1:5))-min(Ans(:,1:5));%最大相差
Ans

运算结果

最后，老师的要求是手抄写下来三个不同样本的估算结果，直接运行上面的代码，记录下来表格的结果就可以了。（注意！！随机函数每次的结果都是不同，最好不要抄别人的结果，很容易被发现是抄别人的）

关于作业内容的音频（老师上课的音频）
链接：https://pan.baidu.com/s/1-4vIztWULYucds9LmeICXQ
提取码：uf5e

---

1. 这是一个学习最小二乘法时候我参考的博客 ↩︎