矩阵方阵matlab,求助!!错误使用 inv 矩阵必须为方阵。
>> 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 矩阵必须为方阵。相关推荐
- matlab矩阵 0,matlab zeros初始化为0矩阵
zeros为创建一个值为零的数组: 如matrix1=zeros(4,5);%4*5的矩阵,矩阵中每个元素都为0 matrix2=zeros(4,5,3);%4*5*3的数组,数组中每个元素都为0 下 ...
- 矩阵正定 matlab,matlab 错误使用 chol 矩阵必须为正定矩阵
警告: 矩阵为奇异值.接近奇异值或缩放错误.结果可能不准确.RCOND = NaN. 警告: 矩阵为奇异值.接近奇异值或缩放错误.结果可能不准确.RCOND = NaN. 警告: 矩阵为奇异值.接近奇 ...
- matlab 矩阵命令,matlab矩阵运算命令
有哪位高人总结了matlab中的矩阵的基本运算命令?还有有关极限.积 1.1 矩阵的表示 1.2 矩阵运算 1.2.14 特殊运算 1.矩阵对角线元素的抽取 函数 diag 格式 X = diag(v ...
- 范德蒙德矩阵在MATLAB中怎么表示,Python 之 Python与MATLAB 矩阵操作总结
Python 之 Python与MATLAB 矩阵操作小结 一.线形代数理论基础 线形代数(linear algebra)是数学的一个分支,研究矩阵理论.向量空间.线性变换和有限维线形方程组等内容. ...
- matlab 矩阵命令,matlab中的矩阵的基本运算命令
matlab中的矩阵的基本运算命令 (2013-07-19 08:45:49) 1.1 矩阵的表示 1.2 矩阵运算 1.2.14 特殊运算 1.矩阵对角线元素的抽取 函数 diag 格式 X = d ...
- matlab怎么产生帕斯卡矩阵,MATLAB(一):矩阵基本操作
MATLAB面向矩阵! MATLAB面向矩阵! MATLAB面向矩阵! 一些特殊矩阵 通用性特殊矩阵 如零矩阵,幺矩阵,单位矩阵等 用于专门学科的特殊矩阵 如魔方矩阵,范德蒙矩阵,希尔伯特矩阵 ...
- matlab 比较矩阵差异,Matlab矩阵操作
1求特征值 [v,d]=eig(A); A为你的矩阵,V为特征向量矩阵,D为特征值矩阵 2除法运算 Matlab提供了两种除法运算:左除(\)和右除(/).一般情况下,x=a\b是方程a*x =b的解 ...
- matlab 调整矩阵形状,matlab对矩阵/向量的常用操作(拼接矩阵、向量逆序、改变矩阵形状、求行阶梯形矩阵、提取矩阵的一部分等)...
几乎所有变量在matlab中都可以视为矩阵(1 x 1元素,1 x n向量,m x n矩阵等),matlab中对矩阵/向量的操作非常多,个人认为对矩阵的操作是体现matlab功底的地方:灵活搭配使用这 ...
- matlab 解方程组 矩阵,用MATLAB计算矩阵和解线性方程组.ppt
用MATLAB计算矩阵和解线性方程组.ppt 第五章 矩阵与行列式 §5.6 用MATLAB计算 矩阵与行列式 用MATLAB计算矩阵与行列式 行列式的求值 矩阵的基本运算 矩阵的加.减 数与矩阵相乘 ...
最新文章
- Windows下通过Python 3.x的ctypes调用C接口
- 【Qt】通过QtCreator源码学习Qt(三):linux平台的信号、程序崩溃处理
- spring boot +spring data jpa +druid 多数据源配置
- 8.0 C++远征:对象数组
- 创建型、结构型、行为型模式(2)
- Python sorted功能
- UDS - 诊断协议和协议栈介绍
- HTML5新特性之标签实战应用:超星尔雅平台
- 小米盒子3增强版 android,小米盒子3增强版拆机评测 怎么样?好不好?
- C 关于unsigned int compzero = ~0;与unsigned int compzero = 0xFFFF; 的区别!
- 计算机组装diy,电脑diy,详细教您如何组装电脑
- 《老路用得上的商学课21—25》消费心理学(二)
- laravel maatwebsite/excel3.1 导入导出详解
- 微信小程序报错“对应的服务器证书无效”
- 计蒜客 - 猴子打字
- Java 添加Word项目符号、编号列表
- 电脑数据没有了怎么恢复?笔记本电脑数据丢失怎样恢复
- 关于IOS设备window onscroll滚动条滚动事件不触发的问题
- jdbc写入数据库乱码问题
- cdr轮廓字怎么把轮廓拆出来_用corelDraw怎么把文字做成轮廓?
热门文章
- 分布式事务中间件 Fescar - 全局写排它锁解读
- 使用Xshell连接Linux服务器
- springmvc DispatchServlet初始化九大加载策略(一)
- 「超全」工欲善其事必先利其器!
- 关于思科网际操作系统漏洞情况的通报
- MapReduce01
- Recurrent Neural Network系列2--利用Python,Theano实现RNN
- Java——容器(Set)
- ​Android中如何使用Intent在Activity之间传递对象[使用Serializable或者Parcelable]
- ADMM:交替方向乘子算法