>> clear;

%输入原始数据和预测长度

x0=[95.94,129.46,156.29,145.24,115.07,173.51,253.91];

x=[1,2,3,4,5,6,7,8,9,10];

n1=length(x0);

n2=length(x);

n1,n2,x0 %进行1-AGO生成

x1=zeros(1,7);

x1(1,1)=x0(1,1);

for k=2:7

x1(1,k)=x1(1,k-1)+x0(1,k);

end

%进行紧邻均值生产

z1=zeros(1,6);

for k=1:6

z1(1,k)=0.5.*[x1(1,k)+x1(1,k+1)];

end

%计算矩阵B和向量yn

B=zeros(6,2);

yn=zeros(6,1);

for k=1:6

B(k,1)=-z1(1,k);

B(k,2)=1;

yn(k,1)=x0(1,k+1);

end

%计算未知参数a和u

A=zeros(2,1);

A=inv(B.*B).* B.* yn;

a=A(1,1);

u=A(2,1);

%微分方程的离散化解

xx1=zeros(1,7);

xx1(1,1)=x0(1,1);

for k=1:6

xx1(1,k+1)=(x0(1,1)-u/a).*exp(-a.*k)+u/a;

end

%建立的GM(1,1)模型

xx0=zeros(1,7);

xx0(1,1)=x0(1,1);

for k=2:7

xx0(1,k)=xx1(1,k)-xx1(1,k-1);

end

xx0

%残差计算

e0=zeros(1,7);%绝对残差

for k=1:7

e0(k)=x0(1,k)-xx0(1,k);

end

e0

q=zeros(1,7);%相对残差

for k=1:7

q(k)=e0(k)/x0(1,k);

end

q

averq=zeros(1);%平均误差

averq=mean(abs(q));

%后验差计算

averx0=zeros(1);%x0平均值

averx0=mean(x0);

S1=zeros(1);%x0均方差

for k=1:7

S1=S1+(x0(k)-averx0)^2

end

S1=(S1/(6))^0.5;

avere0=zeros(1);%e0平均值

avere0=mean(abs(e0));

S2=zeros(1);%e0均方差

for k=1:7

S2=S2+(e0(k)-avere0)^2;

end

S2=(S2/(6))^0.5;

C=zeros(1);%后验差比值

C=S2/S1;

p=zeros(1);%小误差概率

count=zeros(1);

for k=1:7

if e0(k)-avere0<0.6745.*S1

count=count+1;

end

end

p=count/7;

%关联度计算

r=zeros(1);

for k=1:7

r=r+(min(abs(e0))+0.5.*max(abs(e0)))/(e0(k)+0.5.*max(abs(e0)))/7;

A=inv(B);

end

n1 =

7

n2 =

10

x0 =

95.9400  129.4600  156.2900  145.2400  115.0700  173.5100  253.9100

错误使用 inv

矩阵必须为方阵。

新手小白!求各位大神帮忙!!

矩阵方阵matlab,求助!!错误使用 inv 矩阵必须为方阵。相关推荐

  1. matlab矩阵 0,matlab zeros初始化为0矩阵

    zeros为创建一个值为零的数组: 如matrix1=zeros(4,5);%4*5的矩阵,矩阵中每个元素都为0 matrix2=zeros(4,5,3);%4*5*3的数组,数组中每个元素都为0 下 ...

  2. 矩阵正定 matlab,matlab 错误使用 chol 矩阵必须为正定矩阵

    警告: 矩阵为奇异值.接近奇异值或缩放错误.结果可能不准确.RCOND = NaN. 警告: 矩阵为奇异值.接近奇异值或缩放错误.结果可能不准确.RCOND = NaN. 警告: 矩阵为奇异值.接近奇 ...

  3. matlab 矩阵命令,matlab矩阵运算命令

    有哪位高人总结了matlab中的矩阵的基本运算命令?还有有关极限.积 1.1 矩阵的表示 1.2 矩阵运算 1.2.14 特殊运算 1.矩阵对角线元素的抽取 函数 diag 格式 X = diag(v ...

  4. 范德蒙德矩阵在MATLAB中怎么表示,Python 之 Python与MATLAB 矩阵操作总结

    Python 之 Python与MATLAB 矩阵操作小结 一.线形代数理论基础 线形代数(linear algebra)是数学的一个分支,研究矩阵理论.向量空间.线性变换和有限维线形方程组等内容. ...

  5. matlab 矩阵命令,matlab中的矩阵的基本运算命令

    matlab中的矩阵的基本运算命令 (2013-07-19 08:45:49) 1.1 矩阵的表示 1.2 矩阵运算 1.2.14 特殊运算 1.矩阵对角线元素的抽取 函数 diag 格式 X = d ...

  6. matlab怎么产生帕斯卡矩阵,MATLAB(一):矩阵基本操作

    MATLAB面向矩阵! MATLAB面向矩阵! MATLAB面向矩阵! 一些特殊矩阵 通用性特殊矩阵   如零矩阵,幺矩阵,单位矩阵等 用于专门学科的特殊矩阵   如魔方矩阵,范德蒙矩阵,希尔伯特矩阵 ...

  7. matlab 比较矩阵差异,Matlab矩阵操作

    1求特征值 [v,d]=eig(A); A为你的矩阵,V为特征向量矩阵,D为特征值矩阵 2除法运算 Matlab提供了两种除法运算:左除(\)和右除(/).一般情况下,x=a\b是方程a*x =b的解 ...

  8. matlab 调整矩阵形状,matlab对矩阵/向量的常用操作(拼接矩阵、向量逆序、改变矩阵形状、求行阶梯形矩阵、提取矩阵的一部分等)...

    几乎所有变量在matlab中都可以视为矩阵(1 x 1元素,1 x n向量,m x n矩阵等),matlab中对矩阵/向量的操作非常多,个人认为对矩阵的操作是体现matlab功底的地方:灵活搭配使用这 ...

  9. matlab 解方程组 矩阵,用MATLAB计算矩阵和解线性方程组.ppt

    用MATLAB计算矩阵和解线性方程组.ppt 第五章 矩阵与行列式 §5.6 用MATLAB计算 矩阵与行列式 用MATLAB计算矩阵与行列式 行列式的求值 矩阵的基本运算 矩阵的加.减 数与矩阵相乘 ...

最新文章

  1. Windows下通过Python 3.x的ctypes调用C接口
  2. 【Qt】通过QtCreator源码学习Qt(三):linux平台的信号、程序崩溃处理
  3. spring boot +spring data jpa +druid 多数据源配置
  4. 8.0 C++远征:对象数组
  5. 创建型、结构型、行为型模式(2)
  6. Python sorted功能
  7. UDS - 诊断协议和协议栈介绍
  8. HTML5新特性之标签实战应用:超星尔雅平台
  9. 小米盒子3增强版 android,小米盒子3增强版拆机评测 怎么样?好不好?
  10. C 关于unsigned int compzero = ~0;与unsigned int compzero = 0xFFFF; 的区别!
  11. 计算机组装diy,电脑diy,详细教您如何组装电脑
  12. 《老路用得上的商学课21—25》消费心理学(二)
  13. laravel maatwebsite/excel3.1 导入导出详解
  14. 微信小程序报错“对应的服务器证书无效”
  15. 计蒜客 - 猴子打字
  16. Java 添加Word项目符号、编号列表
  17. 电脑数据没有了怎么恢复?笔记本电脑数据丢失怎样恢复
  18. 关于IOS设备window onscroll滚动条滚动事件不触发的问题
  19. jdbc写入数据库乱码问题
  20. cdr轮廓字怎么把轮廓拆出来_用corelDraw怎么把文字做成轮廓?

热门文章

  1. 分布式事务中间件 Fescar - 全局写排它锁解读
  2. 使用Xshell连接Linux服务器
  3. springmvc DispatchServlet初始化九大加载策略(一)
  4. 「超全」工欲善其事必先利其器!
  5. 关于思科网际操作系统漏洞情况的通报
  6. MapReduce01
  7. Recurrent Neural Network系列2--利用Python,Theano实现RNN
  8. Java——容器(Set)
  9. ​Android中如何使用Intent在Activity之间传递对象[使用Serializable或者Parcelable]
  10. ADMM:交替方向乘子算法