| 查看: 541 | 回复: 1 | |||
[求助]
GM（1,1）的matlab程序
|
|
1. 这个程序只能运行一次，第二次就会出现维度不对的问题，改了很多遍也都这样，不知道怎么改才能让他正常运行 2. GM(1,1)的a，b这里没有输出，本人菜鸟一枚，不知道应该怎么编写这个程序，还望大神指点一二 3. 以上每一步计算出来的数据如果输出保存 很急，谢谢各位了~ h=xlsread('Microsoft.xlsx');%数据保存在excel里面调用 global k; for k=1:1:1254 %一行一行的计算 y=[h(k,: )]; n=length(y); yy=ones(n,1); yy(1)=y(1); clear for i=2:n yy(i)=yy(i-1)+y(i); end B=ones(n-1,2); for i=1: (n-1) B(i,1)=-(yy(i)+yy(i+1))/2; B(i,2)=1; end BT=B'; for j=1:n-1 YN(j)=y(j+1); end YN=YN'; A=inv(BT*B)*BT*YN; %出错的地方 a=A(1); u=A(2); t=u/a; t_test=1; i=1:t_test+n; yys(i+1)=(y(1)-t).*exp(-a.*i)+t; yys(1)=y(1); for j=n+t_test:-1:2 ys(j)=yys(j)-yys(j-1); end x=1:n; xs=2:n+t_test; yn=ys(2:n+t_test); det=0; for i=2:n det=det+abs(yn(i)-y(i)); end det=det/(n-1); disp(['百分绝对误差为：',num2str(det),'%']); disp(['预测值为： ',num2str(ys(n+1:n+t_test))]); end [ Last edited by OlivarD on 2013-5-12 at 01:41 ] |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何求这个matlab程序画出的legend曲线的法矢量？
已经有6人回复
- 急求，求帮忙写一个matlab程序！
已经有3人回复
- 最优化方法理论的MATLAB程序典型例子
已经有158人回复
- 想写matlab程序之前不妨看看matlab的fileexchange
已经有4人回复
- 算法程序用matlab，C语言，软件平台，硬件平台运行耗时之间的关系
已经有9人回复
- 求解相算法的相关matlab程序块
已经有3人回复
- 求窗口傅里叶变换的matlab程序
已经有6人回复
- 关于振动学方面的MATLAB程序怎么发核心？
已经有10人回复
- 中心差分求解二维热传导的matlab程序
已经有12人回复
- 如何知道一个MATLAB程序在计算中占的内存是多大？谢谢！
已经有3人回复
- matlab程序中调用m文件的问题
已经有13人回复
- 求助多输出支持向量机的matlab程序
已经有6人回复
- 模式搜索MATLAB程序
已经有9人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- William T. Freeman的基于例子超分辨率matlab程序跑不起来啊！
已经有6人回复
- 关于贝叶斯线性判别分析的matlab程序
已经有4人回复
- 同样的MATLAB程序，为何在别人的电脑上运行不了
已经有21人回复
- matlab 串口通信的问题
已经有3人回复
- 【求助】一组光谱数据的多元散射校正（MSC），求具体的MATLAB程序及用法
已经有6人回复
- 【求助】多目标pso算法的MATLAB程序
已经有3人回复
meatball1982
铜虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: -3.4
- 红花: 8
- 帖子: 297
- 在线: 108.3小时
- 虫号: 1530793
- 注册: 2011-12-09
- 性别: GG
- 专业: 理论和计算化学
2楼2013-05-15 14:59:44