大数实验2与实验3.pptVIP

综合实验2、3讲解 以模型（2）为例： 步骤一：方程求解 输入命令：dsolve(DW=(u*W)*exp(-r*t),t) 运行结果为：ans =C1*exp(-u/r*exp(-r*t)) 相关程序：（一）首先定义两个函数 （二）编写程序，作相关计算 程序编写： * （拟选用的模型） Step 1、计算如下问题： 作出散点图； 应用matlab求模型的解（解析解或数值解）； 进行非线性模型拟合，并根据表中的数据求出模型的参数及相关系数； 作出模型相应的实测点和预测点的比较图； 计算模型的相关特征点：即当t何时使得dW/dt=0（称拐点周龄），并求出在拐点周龄处的体重（称拐点体重）。 Step 2、在相关指标计算及图形比较的基础上，分析四个模型在应用上的优劣。 解 function z=ff(k,t) %定义函数 %z=k(1)./(1-k(2)*exp(-k(3)*t)); % Logistic模型 %z=k(1)*exp(-k(2)./(k(3)*exp(-k(3)*t))); % 模型2 z=k(1)*(1-k(2)*exp(-k(3)*t)).^3; % Bertlanffy模型 %z=k(1)*exp(-k(2)*exp(-k(3)*t)); % Gompertz模型 function y=df2(k,x) %