数学建模①

解： (1) 设商品重量 ( 质量 ) 为 w, 表面积为 S, 那么商品价格 C=aw+bS+f，其中 a, b, f 为常数

 作图

a=1;k=1;f=1;
w=1:10;
c=a+k*w.^(-1/3)+f./w;
plot(w,c)

注意点：①必须要点乘   ②格式是w^-1，实际写/w

求一阶导数：

diff(c)

注意这里一直报错但是不是常见原因

这里和烂大街的例子不同的是它是点乘

答案是直接手动求，同时k,f都进行了保留，不赋值

 

 解：设鱼的形状相似，那么

 这里分析身长L和胸围C的比例，即L/C=常数，那么直接简化

L=[36.8, 31.8 43.8 36.8 32.1 45.1 35.9 32.1];
W=[765 482 1162 737 482 1389 652 454];
C=[24.8 21.3 27.9 24.8 21.6 31.8 22.9 21.6];
%模型 1（公式(3)）
X1=L.^3;
fun1=@(m,x)m*x;
m=lsqcurvefit(fun1, 0.1, X1, W) %计算得到 m=0.0146
W1=fun1(m,X1);%拟合值
norm(W-W1)%拟合误差 98.48

3、考虑新冠疫情传播预测 Logistic 模型，讨论下列问题，并就上海参数作定量分析。

（1）若防护措施使得日接触数λ降低一半，其他参数不变，新增病人的高峰时刻和相应人

数会有怎样的变化？

（2）若隔离、疫苗等措施使得易感人群下降到总人数的 70%，其他参数不变，新增病人的

高峰时刻和相应人数会有怎样的变化？

（3）考虑感染者大约 7 天后不再具有传染性，模型应作怎样修正？若其他参数不变，新增

病人的高峰时刻和相应人数会有怎样的变化？

解：（1） 新增病人的高峰时刻推迟为原来的两倍时间；

                    新增病人的新增人数速度变为原来的一半。

（2） 新增病人的新增人数速度降低；

 高峰时刻往后推迟。

clear;close;
x=[0.023 0.031 0.022 0.047 0.064 0.083 0.1 0.131 0.164];
for t=9:40
 x(t+1)=x(t)+0.28*(x(t)-x(t-7))*(1-x(t));
end
dx=diff(x);
subplot(1,2,1);plot(1:41,x);title('累计感染')
subplot(1,2,2);plot(2:41,dx);title('新增感染')

结果：

4、根据中国疾病预防控制中心提供的全国报告人群新型冠状病毒抗原检测阳性率数据， 对本轮疫情全国感染数作推算，并作出合理的解释。