技术标签： matlab 弹出提示
%灰色预测
function gm(x0,m)
%定义函数 gm(x0,m),x0为原始数据，m为预测个数
n=length(x0);
x1=zeros(1,n);
x1(1)=x0(1);
for i=2:n
x1(i)=x1(i-1)+x0(i); %计算累加序列x1
end
%对原始数列平行移位并赋给 y
y(1:n-1)=x0(2:n);
y=y'; %将 y变成列向量
B=ones(n-1,2);
for i=1:n-1;
%计算数据矩阵B的第一列数据
B(i,1)=-0.5*(x1(i)+x1(i+1));
end
au=inv(B'*B)*B'*y;%计算参数 a,u 矩阵
a=au(1);
u=au(2);
t=u/a;
i=1:m+n;
ago(i+1)=(x0(1)-t).*exp(-a.*i)+t;
ago(1)=x0(1);
for j=n+m:-1:2
yc(j)=ago(j)-ago(j-1);
end
for i=2:n;
error(i)=yc(i)-x0(i); %计算残差值
end
c=std(error)/std(x0);
%计算后验差比
p=0;
for i=2:n
if(abs(error(i)-mean(error))<0.6745*std(x0))
p=p+1;
end
end
p=p/(n-1);
w1=min(abs(error));
w2=max(abs(error));
for i=1:n;
%计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
end
w=sum(w)/(n-1); %
%% 绘图%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x=1:n;
xs=2:n+m;
yn=yc(2:n+m);
plot(x,x0,'^r',xs,yn,'*-b');
% axis([2009+m,2014+m,1,10]); %x，y,(z)坐标范围
% plot([2009:n+2008],x0,'+',[2009:n+2008+m],yc,'*'); %
% grid on;
% xlabel('年份');
% ylabel('销量(mm)');
% title('销售量灰色模型预测拟和曲线');
% legend('实测值','预测值',4);
% fprintf('a,u值:')
% fprintf('%g ',au) %输出参数 a,u 的值
% fprintf('\n累加数列:\n')
% fprintf('%g ',ago) %输出累加数列 ago 的值
% fprintf('\n原始序列:\n')
% fprintf('%g ',x0) %输出原始序列值
% fprintf('\n预测:\n')
% fprintf('%g ',yc) %输出预测的值
% fprintf('\n残差:\n')
% fprintf('%g ',error) %输出残差的值
% fprintf('\n后验差比:\n')
% fprintf('%g',c) %输出后验差比的值
% fprintf('\n小误差概率:\n')
% fprintf('%g',p) %输出小误差概率的值
% fprintf('\n关联度:\n')
% fprintf('%g\n',w) %输出关联度 w
取消
评论
文章浏览阅读188次。文章目录Chapter 4 查找和引用数据的高手LOOKUP函数，查询中的王者INDIRECT函数，将文本转为引用OFFSET函数，求符合条件的区域Chapter 4 查找和引用数据的高手LOOKUP函数，查询中的王者LOOKUP函数，与MATCH函数，第3参数为“1”类似= MATCH(A2,C2:C9,1) 参数1代表函数将查找小于等于查找值的最大值，并返回该数值在查找区域的..._function中单元格引用
文章浏览阅读192次。1081 子段求和题目链接：https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1081题目给出一个长度为N的数组，进行Q次查询，查询从第i个元素开始长度为l的子段所有元素之和。例如，1 3 7 9 -1，查询第2个元素开始长度为3的子段和，1 {3 7 9} -1。3 + 7 + 9 = 19，输出19_51nod 1081
文章浏览阅读319次。https://ftp.gnu.org/gnu/gdb/ 下载gdb-11.2.tar.gz，版本最好同local交叉编译器的对齐。2.拷贝install目录下的bin/gdb和test以及main.c拷贝到u盘。1.板子上运行gdbserver指定端口12345。3.gdb server会显示已经连接。4.server gdb开始debug。2.服务器上运行gdb。_交叉编译gdb 和 gdb server 教程
文章浏览阅读1.3k次。介绍线程的状态与状态间的转换_java thread.sleep占不占用cpu资源
文章浏览阅读657次。一、解读CDH1、CM服务（1）状态（2）Cloudera Management Service实例: 5个进程（3）配置（4）单个服务的配置服务: HDFS Kafka角色/实例: 所属的服务的进程以HDFS为例：2、自定义监控图表（1）查询图标库选定一个服务，选择图表库（2）打开图表（3）复制查询语句SELE..._cdh 添加内存监控
文章浏览阅读6.6k次，点赞4次，收藏25次。一、控制外设时钟的寄存器有以下几个：PCLKCR0（外设时钟寄存器0）:作用于上图标黄外设，所有时钟默认不开启注：这里需要注意TBCLKSYNC这个寄存器，他是让ePWM模块时钟同步，且当检测到TBCLK的上升沿时，同步开启。另一个需要注意的是，及时这个寄存器没有使能ePWM,如果PCLKCR1使能了一样有效。对于使能ePWM推荐做法：• Enable ePWM module clocks in the PCLKCR1 register.• Set TBCLKSYNC to 0.• Confi._tbclksync
文章浏览阅读8.5w次，点赞80次，收藏480次。为什么需要数据仓库？ 传统的数据库中，存放的数据都是一些定制性数据较多，表是二维的，一张表可以有很多字段，字段一字排开，对应的数据就一行一行写入表中，特点就是利用二维表表现多维关系。 但这种表现关系的上限和下限就定死了，比如QQ的用户信息，直接通过查询info表，对应的username、introduce等信息即可，而此时我想知道这个用户在哪个时间段购买了什么？修改..._数据仓库
文章浏览阅读61次。转：http://www.php7.site/tw/book/php7/23.html#4_php7 关于对象的新特性
文章浏览阅读2.5k次，点赞5次，收藏20次。NXP的官方开发板很便宜，比如下图这个MPC5744p的开发板，才300块RMB，而且还提供免费的Simulink自动代码生成工具NXP MBDT（Model-based design toolbox），这点非常友好。但是也存在着如下的问题：1、板载了OpenSDA调试器，在Simulink中编译后可直接下载进开发板中，但是要想把这个OpenSDA调试器移植到自己设计的控制板中，那是相当难的——据我所知，面向PowerPC的OpenSDA固件并不免费提供。2、那我们自己设计控制板，又想要使用NXP _model-based design toolbox
文章浏览阅读771次。效果图：实现代码：html： var list = ['red', 'orange', 'yellow', 'green', 'blue', 'purple','black']; var html = ''; for(var i = 0; i < list.length; i++){ var ul = document.createElement('ul'_sass生成色板
文章浏览阅读114次。引言嵌入式Linux是指对Linux进行剪裁后，将其固化在单片机或者存储器中，应用于特定场合的专用Linux系统。嵌入式系统要求实时性能高，但Linux为分时系统设计的操作系统，尽管最新的内核在实时性能方面有所提高，但它仍然不是一个实时系统，在很多场合不能满足实时性要求。一般地，通过改造Linux的内核以提高其实时性能有2种策略：一种是采用底层编程的方法对Linux内核进行修改(如调度算法、时钟修..._apic 框图
文章浏览阅读5.1k次。Windows server 2012R2 设置文件共享目录报错：无法连接到C$管理共享已验证文件夹xxx在计算机xxx上是否存在问题描述：Windows server 2012R2 机器，在安装完文件共享服务后（安装步骤参考了https://jingyan.baidu.com/article/5553fa821b3a1d25a3393401.html），设置文件共享目录时报错：无法连接到C$管理共享已验证文件夹xxx在计算机xxx上是否存在。解决方案windows server 2012r2 在_无法连接到管理共享以验证文件夹是否存在