查看: 420|回復: 0|關注: 0
打印 上一主題 下一主題

[未答復] 求助,第六章根據灰色預測GM(1,1)的代碼怎么改成非均勻步長的灰色預測

[復制鏈接]

新手

5 麥片

財富積分


050


1

主題

2

帖子

0

最佳答案
跳轉到指定樓層
1#
發表于 2020-1-30 13:42:19 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
求助,第六章根據灰色預測GM(1,1)的代碼怎么改成非均勻步長的灰色預測?比如說用第六章的灰色預測解決第四章的一元回歸問題?
  1. % 灰色預測
  2. % 《MATLAB數學建模方法與實踐》(《MATLAB在數學建模中的應用》升級版),北航出版社,卓金武、王鴻鈞編著.
  3. clear
  4. syms a b;
  5. c=[a b]';
  6. A=[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670];
  7. B=cumsum(A);  % 原始數據累加
  8. n=length(A);
  9. for i=1:(n-1)
  10.     C(i)=(B(i)+B(i+1))/2;  % 生成累加矩陣
  11. end
  12. % 計算待定參數的值
  13. D=A;D(1)=[];
  14. D=D';
  15. E=[-C;ones(1,n-1)];
  16. c=inv(E*E')*E*D;
  17. c=c';
  18. a=c(1);b=c(2);
  19. % 預測后續數據
  20. F=[];F(1)=A(1);
  21. for i=2:(n+10)
  22.     F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a ;
  23. end
  24. G=[];G(1)=A(1);
  25. for i=2:(n+10)
  26.     G(i)=F(i)-F(i-1); %得到預測出來的數據
  27. end
  28. t1=1999:2008;
  29. t2=1999:2018;
  30. G
  31. plot(t1,A,'ko', 'LineWidth',2)
  32. hold on
  33. plot(t2,G,'k', 'LineWidth',2)
  34. xlabel('年份', 'fontsize',12)
  35. ylabel('利潤/(元/年)','fontsize',12)
  36. set(gca,  'LineWidth',2);
復制代碼
這是第六章列舉的灰色預測代碼。



x=[23.80,27.60,31.60,32.40,33.70,34.90,43.20,52.80,63.80,73.40];
y=[41.4,51.8,61.70,67.90,68.70,77.50,95.90,137.40,155.0,175.0];       

這是第四章的數據。

非常感謝!
回復主題 已獲打賞: 0 積分

舉報

您需要登錄后才可以回帖 登錄 | 注冊

本版積分規則

關閉

站長推薦上一條 /4 下一條

快速回復 返回頂部 返回列表
哪一款德州扑克还能玩 南粤36选7 12107 上市公司发行股票 至尊棋牌app官方 捕鱼王者现金版 2017 星星武汉麻将胡牌技巧 双色球弧形码 山西麻将扣点点代理 互联网营销 大众麻将是什么麻将 15选5专家预测号 福州全民麻将下载安装 在家赚钱的方法 今日大盘股票 网上赚钱团队真实吗 股票投资的收益等于 pk10计划软件冠军五码