function gm(x) %定义函数gm(x) %清屏,以使计算结果独立显示
format long; %设置计算精度
if length(x(:,1))==1 %对输入矩阵进行判断,如不是一维列矩阵,进行转置变换
x=x';
end;
n=length(x); %取输入数据的样本量
z=0;
for i=1:n %计算累加值,并将值赋与矩阵be
z=z+x(i,:);
be(i,:)=z;
end
for i=2:n %对原始数列平行移位
y(i-1,:)=x(i,:);
end
for i=1:n-1 %计算数据矩阵B的第一列数据
c(i,:)=-0.5*(be(i,:)+be(i+1,:));
end
for j=1:n-1 %计算数据矩阵B的第二列数据
e(j,:)=1;
end
for i=1:n-1 %构造数据矩阵B
B(i,1)=c(i,:);
B(i,2)=e(i,:);
end
alpha=inv(B.'*B)*B.'*y; %计算参数α、μ矩阵
for i=1:n+3 %计算数据估计值的累加数列,如改n+1为n+m可预测后m-1个值
ago(i,:)=(x(1,:)-alpha(2,:)/alpha(1,:))*exp(-alpha(1,:)*(i-1))+alpha(2,:)/alpha(1,:);
end
var(1,:)=ago(1,:);
for i=1:n+2 %如改n为n+m-1,可预测后m-1个值
var(i+1,:)=ago(i+1,:)-ago(i,:); %估计值的累加数列的还原,并计算出下一预测值
end
for i=1:n
error(i,:)=var(i,:)-x(i,:); %计算残差
end
c=std(error)/std(x); %调用统计工具箱的标准差函数计算后验差的比值c
ago %显示输出预测值的累加数列
alpha %显示输出参数α、μ数列
var %显示输出预测值
error %显示输出误差
c %显示后验差的比值c
GM.rar_GM11_matlab代入_灰色
版权申诉
14 浏览量
2022-09-22
22:32:22
上传
评论
收藏 1KB RAR 举报
![avatar](https://profile-avatar.csdnimg.cn/76d4ee81f5f6465b8156d97495319ccf_weixin_42659252.jpg!1)
JaniceLu
- 粉丝: 85
- 资源: 1万+
最新资源
- 省 市 县 三级联动 数据字典
- 1D-moc-to-CKJ.py
- VMware vCenter Server 7.0U3r安装ISO的度盘地址,由Patch界面提供,产品界面仍然为7.0U3p
- python脚本-批量修改文件后缀
- VMware vCenter Server 8.0U2d安装ISO的度盘地址,由Patch界面提供,下载界面仍然为8.0U2c
- bge-reranker-base本地模型
- 前端学习笔记,做一个简单的网站-学生成绩.HTML
- go语言实现的简单搜索引擎demo,使用了redis,colly,gin-gonic等技术
- 基于时空图卷积(ST-GCN)的骨骼动作识别(python源码+毕业设计).zip
- JAVAEE开发+MyBatis+Spring+SpringMVC+个人博客管理系统+用于发布日志和博客
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback-tip](https://img-home.csdnimg.cn/images/20220527035111.png)