| 查看: 390 | 回复: 0 | ||
wodeai87金虫 (正式写手)
|
[求助]
大家帮我看看这个代码错在哪里？（无偏灰色预测模型代码）万谢
|
|
function gm=GM(x0) %灰色系统GM（1，1）预测 x0=[1.01 1.6387 2.7283 4.4717 7.3991 12.1725]' T=input('T='); x1=zeros(1,length(x0)); B=zeros(length(x0)-1,2); yn=zeros(length(x0)-1,1); hatx0=zeros(1,length(x0)+T); hatx00=zeros(1,length(x0)); hatx1=zeros(1,length(x0)+T); epsilon=zeros(length(x0),1); omega=zeros(length(x0),1); for i=1:length(x0) for j=1:i x1(i)=x1(i)+x0(j);%累加生成 end end x1 for i=1:length(x0)-1 B(i,1)=(-1/2)*(x1(i)+x1(i+1)); B(i,2)=1; yn(i)=x0(i+1); end hatA=(inv(B'*B))*B'*yn %求a，u for k=1:length(x0)+T hatx1(k)=(x0(1)-hatA(2)/hatA(1))*exp(-hatA(1)*(k-1))+hatA(2)/hatA(1); end hatA(1)=log((2-hatA(1))/(2+hatA(1))); hatA(2)=2*hatA(2)/(2+hatA(1)); hatx0(1)=hatx1(1); for k=2:length(x0)+T hatx0(k)=hatA(2)*exp(hatA(1)*(k-1));%累减还原 end hatx0 %历史数据的模拟值 for i=1:length(x0)%开始检验 epsilon(i)=x0(i)-hatx0(i); omega(i)=(epsilon(i)/x0(i))*100; end c=std(epsilon)/std(x0); p=0; for i=1:length(x0) if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) p=p+1; end end p=p/length(x0) if p>0.95&c<0.35 disp('预测精度好，预测值为:') disp(hatx0(length(x0)+T)) elseif p>0.85&c<0.5 disp('预测精度合格，预测值为:') disp(hatx0(length(x0)+T)) elseif p>0.70&c<0.65 disp('预测精度勉强合格，预测值为:') disp(hatx0(length(x0)+T)) elseif p<=0.7&c>0.65 disp('预测精度不合格') end for i=1:length(x0) hatx00(i)=hatx0(i); end z=1:length(x0); plot(z,x0,'-',z,hatx00,':') text(2,x0(2),'历史数据：实线') text(length(x0)/2,hatx00(length(x0))/2,'模拟数据：虚线') end 本人菜鸟，实在找不到错误，还希望了解的人帮忙！！！ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
求助求定位，大家看看我的这个情况，该怎么办呢？万分感谢了
已经有14人回复
- 请大家帮我看看编辑的信，最好是懂植物学的朋友们，谢谢。
已经有11人回复
- 请大家帮我看看这个计算结果，多谢！
已经有7人回复
- 大家帮我看看图片中的黄色的支撑哪有卖的，万谢
已经有5人回复
- 大家帮我看看图片中的黄色的支撑哪有卖的，万谢
已经有4人回复
- 大家帮我看看这个图是啥意思谢谢啊
已经有4人回复
- 请各位帮忙看看，我这种情况该如何选择公派，如何操作？万分感谢
已经有8人回复
- 【求助】帮忙再看看这个化合物，是白芷香豆素类的化合物吗？万分感谢 ，
已经有6人回复
- 大家帮我看看这句话是什么意思，谢谢了
已经有9人回复
- 【请教】大家帮我看看这种薄膜样品的截面SEM照片是怎么获得的，谢谢了
已经有18人回复
- 【求助】大家帮我看看反应的第一步第二步是什么机理，谢谢了。
已经有17人回复
- 【求助/交流】帮我看看这两种植物学名叫什么，谢谢。
已经有12人回复
- 几乎要崩溃了！大家帮我看看是否有挽回的余地？谢谢了
已经有112人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
