查看: 8722|回复: 8

[建模教程] Verhulst 预测模型

字体大小: 正常放大

20 主题	2 听众	72 积分

升级 70.53%

该用户从未签到

电梯直达

1^#

发表于 2009-9-8 22:04 |只看该作者 |正序浏览

|招呼Ta 关注Ta

clc,clear
x1=[4.93 5.33 5.87 6.35 6.63 7.15 7.37 7.39 7.81 8.35
9.39 10.59 10.94 10.44];
n=length(x1);
nian=1990:2003;
plot(nian,x1,'o-');
x0=diff(x1);
x0=[x1(1),x0]
for i=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)'
abhat=B\Y %估计参数a,b 的值
x=dsolve('Dx+a*x=b*x^2','x(0)=x0'); %求解常微分方程
x=subs(x,{'a','b','x0'},{abhat(1),abhat(2),x1(1)}); %代入参数值
yuce=subs(x,'t',0:14) %计算预测值
digits(6); x=vpa(x) %显示微分方程的解，为了提高计算精度，把该语句放在计算预测值
之后，或者不使用该语句
yuce(16)=yuce(15);
x1_all=[x1,9.92,10.71];
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(1:end-1))+0.5*x1_all_0(end));
s1=abs(sum(yuce_0(1:end-1))+0.5*yuce_0(end));
tt=yuce_0-x1_all_0;
s1_s0=abs(sum(tt(1:end-1))+0.5*tt(end));
absdegree=(1+s0+s1)/(1+s0+s1+s1_s0) %计算灰色绝对关联度
c=std(epsilon,1)/std(x1_all,1) %计算标准差比值
请高手帮忙一下，里面哪里错误了，看了一个晚上，老是不对！！！！
这是Verhulst 预测模型 GM（1，1）改进过来的

zan

转播0 淘帖0 分享0 收藏2 支持0 反对0 微信

使用道具举报

liu123456cn

6 主题	4 听众	492 积分

升级 64%

TA的每日心情

	奋斗 2015-2-6 13:42

签到天数: 112 天

[LV.6]常住居民II

群组: 中国矿业大学数学建模协会

群组: 全国大学生数学建模竞

群组: 数学建摸协会

9^#

发表于 2010-8-1 10:42 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

高手帮忙啊！~~~~~~~~~求助

我也说一句

发表

使用道具举报

gssdzc

0 主题	2 听众	941 积分

升级 85.25%

该用户从未签到

群组: 兰州大学数学建模协会

8^#

发表于 2010-6-17 22:04 |只看该作者 |招呼Ta 关注Ta

是灰色理论的。谢谢

我也说一句

发表

使用道具举报

liuxn_2005

1 主题	6 听众	693 积分

升级 23.25%

TA的每日心情

	开心 2021-2-3 08:59

签到天数: 9 天

[LV.3]偶尔看看II

群组: 2017himcm交流群组

7^#

发表于 2010-4-13 10:01 |只看该作者 |招呼Ta 关注Ta

matlab 程序。。。。。。。。。。。。。。。。。。。。。。

我也说一句

发表

使用道具举报

d815856116

0 主题	3 听众	106 积分

升级 3%

该用户从未签到

6^#

发表于 2009-9-10 10:10 |只看该作者 |招呼Ta 关注Ta

我也不会看不懂、、、、、、、

我也说一句

发表

使用道具举报

xinzhiyong

20 主题	2 听众	72 积分

升级 70.53%

该用户从未签到

5^#

发表于 2009-9-9 22:37 |只看该作者 |招呼Ta 关注Ta

xiexie,juejie

我也说一句

发表

使用道具举报

杨帆

4 主题	4 听众	626 积分

升级 6.5%

TA的每日心情

	奋斗 2012-12-4 13:38

签到天数: 314 天

[LV.8]以坛为家I

群组: 数学趣味、游戏、IQ等

群组: 数学建模

群组: 我们一定会赢

4^#

发表于 2009-9-8 23:09 |只看该作者 |招呼Ta 关注Ta

clc,clear
x1=[4.93 5.33 5.87 6.35 6.63 7.15 7.37 7.39 7.81 8.35 9.39 10.59 10.94 10.44];
n=length(x1);
nian=1990:2003;
plot(nian,x1,'o-');
x0=diff(x1);
x0=[x1(1),x0]
for i=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)'
abhat=B\Y %估计参数a,b 的值
x=dsolve('Dx+a*x=b*x^2','x(0)=x0'); %求解常微分方程
x=subs(x,{'a','b','x0'},{abhat(1),abhat(2),x1(1)}); %代入参数值
yuce=subs(x,'t',0:14) %计算预测值
digits(6); x=vpa(x) %显示微分方程的解，为了提高计算精度，把该语句放在计算预测值之后，或者不使用该语句
yuce(16)=yuce(15);
x1_all=[x1,9.92,10.71];
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(1:end-1))+0.5*x1_all_0(end));
s1=abs(sum(yuce_0(1:end-1))+0.5*yuce_0(end));
tt=yuce_0-x1_all_0;
s1_s0=abs(sum(tt(1:end-1))+0.5*tt(end));
absdegree=(1+s0+s1)/(1+s0+s1+s1_s0) %计算灰色绝对关联度
c=std(epsilon,1)/std(x1_all,1) %计算标准差比值

你数据输入有问题