变异系数:v=std(A)./abs(mean(A)) k阶原点矩、中心距: ak=mean(A.^k) bk=mean((A-mean(A)).^k) %中心距系统命令bk=moment(A,k) 3.极差、四分位极差(上、下分位数之差) R=rangr(A) R1=iqr(A) 4.异常点判别(截断点) XJ=parctile(A,[25])-1.5*R1 SJ=parctile(A,[75])+1.5*R1 ...
functionT=application_QAR1(yy,ww,valfa,vb,nombrearchivo,s1,s2,s3,s4)%% 指定一个QAR(1)模型%yy=时间序列 nyx1%ww=独立时间序列 nyx1%valfa=阿尔法向量(分位数水平)%vb=自举子样本大小的向量%rechcvmsub= nbx3 拒绝矩阵ny=length(yy);nw=length(ww);na=length(valfa);nb=length(vb);%%rechcvmsub...
生成的对角矩阵的Matlab调用格式为:A=diag(v,k) 生成第 k 个对角线由向量 v 组成的对角阵, k 可以是正数,零或负数。当 k 是零时指主对角线(可简记为A= diag(v) ), k 是负数时 v 在对角线下, k 是正数时 v 在对角线上。 (6)幻方矩阵(魔阵) 即每行每列和对角线上元素之和相等 命令: magic...
引入高斯核函数K(x, y) = exp(-‖x-y‖^2/2σ^2),其中x和y表示原始数据点,σ是高斯核函数的宽度参数。 通过高斯核函数将原始数据映射到高维特征空间,得到新的数据表示。 在特征空间中使用L2约束的最小二乘法进行拟合,即求解如下优化问题:min‖w‖^2 + C∑(yi - w^Tφ(xi))^2,其中w是拟合参数,...
令P(k)=(p1(k),p2(k),pn(k),表示第k层元素对第k-1层个元素的排序,则第k层元素对于总目标的总排序为:w(k)=(w1(k),w2(k),wn(k)T= p(k) w( 18、k-1)或 I=1,2,n同样,也需要对总排序结果进行一致性检验。假定已经算出针对第k-1层第j个元素为准则的C.I.j(k)、R.I.j(k)和C.R...
由于判断矩阵构造过程只是在做两两比较,看不到这些传递过来的关系,故不可能做到完全的一致性,所以需要对判断矩阵进行一致性检验:即保证一致性的偏差不能太大。 定义2.1若正互反矩阵A满足: a_{i j} a_{j k}=a_{i k}, \quad i, j, k=1, \cdots, n \\ ...
D-W= ∑(Et-Et-1)^2/∑Et^2,Et是第t期的残差,Et-1是第t-1期的残差,∑是对t从第2期到第t期求和,^2表示平方.在D-W小于等于2时,D-W检验法则规定：如D-W＞d U,认为ei无自相关；有自相关.如4-D-W＜dL,认为ei存在负自相关；如dL＜4-D-W＜dU,不能确定是否 有自相关.例如：...
w 白色. 实心小点 hold on 添加图形 r 红色 x 叉号形状 grid on 添加网格 k 黑色 s 方形 - 实线 d 菱形 -- 虚线 ^ 向上箭头 3、可以用subplot(3,3,1)表示将绘图区域分为三行三列,目前使用第一区域。此时如要画不同的图形在一个窗口里,需要hold on。 附录1 函数分类查询 ...
p=0.0001;i=2;k=abs(m(2)-m(1));while k>pi=i+1;x(:,i)=A*y(:,i-1);m(i)=max(x(:,i));y(:,i)=x(:,i)/m(i);k=abs(m(i)-m(i-1));enda=sum(y(:,i));w=y(:,i)/a;t=m(i);disp('权向量');disp(w);disp('最大特征值');disp(t);%以下是一致性检验CI=(...
没有问题的，带回去A的行列式是0的，不过是由于是数值计算，结果和0是差一点点，但基本上都是0，这是数值计算带来的误差。