灰色verhulst模型
- 格式:docx
- 大小:17.29 KB
- 文档页数:3
森林可燃物灰色Verhulst模型动态预测
居恩德;何忠秋
【期刊名称】《森林防火》
【年(卷),期】1994(000)002
【摘要】森林可燃物灰色verhulst模型动态预测居恩德,何忠秋,刘艳红,张宗林(东北林业大学)(黑龙江省森林经营研究所)(东北林业大学)(讷河县林业局)1数据来源与调查方法数据来自1988年周再知对塔河林业局三种可燃物类型的调查资料。
可燃物类型包括杜鹃一落...
【总页数】4页(P44-46,52)
【作者】居恩德;何忠秋
【作者单位】不详;不详
【正文语种】中文
【中图分类】S718.55
【相关文献】
1.基于挣得进度的灰色Verhulst模型对工程工期的动态预测 [J], 贾铭钰;杨秀芸;朱美洁;李招弟
2.基于灰色Verhulst动态模型的人口预测——以天津为例 [J], 陈立新;王琳
3.基于灰色Verhulst的互联网上网人数动态预测模型 [J], 朱苗苗;牛国锋;乐德广
4.灰色马尔科夫Verhulst动态模型在滑坡形变预测中的应用 [J], 邓洪高;姚鹏远;孙希延;纪元法;严素清
5.灰色Verhulst动态新陈代谢模型在产品价格预测与需求预测中的应用 [J], 曹萃文;顾幸生
因版权原因,仅展示原文概要,查看原文内容请购买。
灰色verhulst模型
close;
clear;
x1=[0.025 0.023 0.029 0.044 0.084 0.164 0.332 0.521 0.97 1.6 2.45 3.11 3.57 3.76 3.96 4 4.46 4.4 4.49 4.76 5.01];
n=length(x1);
year=0:n-1;
figure(1);
plot(year,x1,'b+');
x0=diff(x1);
x0=[x1(1),x0];
fori=2:n
z1(i)=0.5*(x1(i)+x1(i-1));
end
z1;
B=[-z1(2:end)',z1(2:end)'.^2];
Y=x0(2:end)';
abvalue=B\Y;
x=dsolve('Dx+a*x=b*x^2','x(0)=x0');
x=subs(x,{'a','b','x0'},{abvalue(1),abvalue(2),x(1)}); forecast=subs(x,'t',0:n-1);
digits(6);x=vpa(x);
forecast
hold on;
plot(year,forecast,'g-.','linewidth',4);
xlabel('时间均匀采样/5h');
ylabel('细菌培养液吸光度/OD600');
legend('实际数量','预测数量');
title('大肠杆菌s型整张曲线');
axis tight;
灰色verhulst模型
nian=0:20;%原始数据
x1=[0.025 0.023 0.029 0.044 0.084 0.164 0.332 0.521 0.97 1.6 2.45 3.11 3.57 3.76 3.96 4 4.46 4.4 4.49 4.76 5.01];%原始数据
n=length(x1);
plot(nian,x1,'o-');
x0=diff(x1);%作累减生成
x0=[x1(1),x0];
z1=0.5*(x1(2:n)+x1(1:n-1)); %求紧邻均值生成序列
B=[-z1',z1'.^2];
Y=x0(2:end)';
ab_hat=B\Y; %估计参数a,b的值
x=dsolve('Dx+a*x=b*x^2','x(0)=x0'); %求解常微分方程
x=subs(x,{'a','b','x0'},{ab_hat(1),ab_hat(2),x1(1)});%代入
参数值
yuce=subs(x,'t',0:20); %计算预测值
%下面显示微分方程的解,为了提高计算精度,把该语句放在计算预测值之后
x=vpa(x,6);
epsilon=x1_all-yuce; %计算残差
delta=abs(epsilon./x1_all); %计算相对误差
delta_mean=mean(delta); %计算平均相对误差
x1_all_0=x1_all-x1_all(1); %观测值数据列的始点零化像yuce_0=yuce-yuce(1); %预测值数据列的始点零化像
s0=abs(sum(x1_all_0(2:end-1))+0.5*x1_all_0(end));
s1=abs(sum(yuce_0(2:end-1))+0.5*yuce_0(end));
tt=yuce_0-x1_all_0;
s1_s0=abs(sum(tt(2:end-1))+0.5*tt(end));
absdegree=(1+s0+s1)/(1+s0+s1+s1_s0);%计算灰色绝对关联度
c=std(epsilon,1)/std(x1_all,1); %计算标准差比值
s1_s0=abs(sum(tt(2:end-1))+0.5*tt(end));
absdegree=(1+s0+s1)/(1+s0+s1+s1_s0)%计算灰色绝对关联度c=std(epsilon,1)/std(x1_all,1) %计算标准差比值。