GM(1,N)动态模型的建立,用灰色系统理论的matlab怎么样编写程序 clc,clearx0=[692.5 649.6 669.8 596.0 607.9 617.9635.2 652.2 667.8 682.8 675.1 695.7];n=length(x0);lamda=x0(1:n-1)./x0(2:n)range=minmax(lamda)x1=cumsum(x0)for i=2:nz(i)=0.5*(x1(i)+x1(i-1));endB=[-z(2:n)',ones(n-1,1)];Y=x0(2:n)';u=B\\Yx=dsolve('Dx+a*x=b','x(0)=x0');x=subs(x,{'a','b','x0'},{u(1),u(2),x1(1)});yuce1=subs(x,'t',[0:n-1]);digits(6),y=vpa(x)%为提高预测精度,先计算预测值,再显示微分方程的解yuce=[x0(1),diff(yuce1)]epsilon=x0-yuce%计算残差delta=abs(epsilon./x0)%计算相对误差rho=1-(1-0.5*u(1))/(1+0.5*u(1))*lamda%计算级比偏差值改变x0的输入值就行了灰色预测matlab代码怎么写 这是我曾经写过的一个灰色预测的程序:第一个文件为函数,需要在调用时输入原始数据x0和预测周期T,第二个文件用于计算灰色关联度,使用时直接修改相关参数和原始数据。第一个文件(用于灰色建模):grymdl.mfunction GM=grymdl(x0,T)输入原始数据x0T为从最后一个历史数据算起的第T时点x1=zeros(1,length(x0));B=zeros(length(x0)-1,2);yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T);Hatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T);epsilon=zeros(length(x0),1);omega=zeros(length(x0),1);for i=1:length(x0)for j=1:ix1(i)=x1(i)+x0(j);endendfor i=1:length(x0)-1B(i,1)=(-1/2)*(x1(i)+x1(i+1));B(i,2)=1;yn(i)=x0(i+1);endHatA=(inv(B'*B))*B'*yn;GM(1,1)模型参数估计for k=1:length(x0)+THatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1);endHatx0(1)=Hatx1(1);for k=2:length(x0)+THatx0(k)=Hatx1(k)-Hatx1(k-1);累计还原得到历史数据的模拟值endfor i=1:length(x0)%开始模型检验epsilon(i)=x0(i)-Hatx0(i);omega(i)=(epsilon(i)/x0(i))*100;endx0;HatA;Hatx0;epsilon;omega;c=std(epsilon)/std(x0。SPSS可以做灰色关联度的分析么 可以的,但是最好用matlabmatlab7.0工具箱中有没有灰色关联度的函数 没有谁能教我matlab 求灰色关联度,或者帮我分析下我的数据,QQ聊下 我想起来当时有上过灰色系统这门课,那本书上面有算法的,你吧他整到matlab 里面就是的,不过当时是上的选修课,也就考试的时候拿计算器算过一次,现在都不记得了,你找找灰色系统那本书谁提供一个可以用的灰色关联度分析软件?? EXCEL就行 用灰色系统建模软件,括号内给你地址了(http://igss.nuaa.edu.cn/institute/),一键式的,可以计算灰色关联度,而且我觉得你应该用matlab来写,而不是sas如何用matlab实现灰色关联度的检验 这要有数据的,根据不同要求,求出相关系数。可把数据发给我QQ2674716548,我写个程序,求出灰色关联度。怎样用matlab做灰色关联度分析方法? 1、function f=grayrelated(X,Y)这里X是标准化后的参考序列,Y是评价矩阵Y=71.8 90.1 0.57 0.45 051 40.2 0.38 0.55 10.552 25 0.22 0.52 1268 90 0.38 0.38 2128 40 0.32 0.3 18.551 45 0.15 0.3 576 95 0.7 0.55 1287 95 0.7 0.5 9.876 90 0.57 0.5 1150 35 0.32 0.35 2068 90 0.57 0.35 18.582 95 0.7 0.35 0100 200 1 1 097.5 180 0.94 0.95 1.395 160 0.88 0.9 2.586.3 105 0.68 0.75 6.382.5 90 0.6 0.7 7.578.8 75 0.53 0.65 8.875 60 0.45 0.7 7.568.8 52.5 0.41 0.55 13.862.5 45 0.38 0.5 17.556.3 37.5 0.34 0.45 21.343.8 26.3 0.28 0.35 50.650 30 0.3 0.4 2537.5 22.5 0.25 0.3 7531.3 18.8 0.23 0.25 10018.8 11.3 0.15 0.15 168.825 15 0.2 0.2 12512.5 7.5 0.1 0.1 212.56.3 0.8 0.05 0.05 256.32、];输入评价矩阵YX=[1 1 1 1 1];X为参考序列,均为1,个数就是指标个数,情形不同要修改个数Len=size(Y,2);取Y矩阵的列数,也就是指标的个数Wen=size(Y,1);取行数,就是目标个数for i=1:LenY(:,i)=(Y(:,i)-mean(Y(:,i)))/sqrt(var(Y(:,i)));将Y矩阵用统计方法标准化标准化,endfor i=1:Len-1S(:,i)=(Y(:,i)-。
随机阅读
- 水稻土游离氧化铁 水稻土的划分
- 中山横栏去中山影视城远吗 请问,丛中山影视城站到横栏镇六沙村怎么去?
- 镇巴永乐镇领导班 石泉自驾游攻略?
- 柳岩越整越丑 同样是整容,为什么主动承认自己整容的明星反而会获得观众好感?
- 瓦房店西杨乡渤海村住蓿 大连海洋大学的宿舍条件如何?校区内有哪些生活设施?
- 余杭农村商业银行本本贷 余杭农村商业银行市民卡丰收借记卡是什么卡?
- 广州东站怎么去三水? 广州东站到三水客运站
- 主要成果,贡献怎么填(教师) 教师贡献值分
- 求几本好看的书? 最美购物袋
- 请问各个城市2017年度中国质量协会质量专业人员6西格玛绿带考试考点分别在哪? 天行健大厦南门
- 激光电子琴通过按键 键盘是什么
- 小学校长简历 翠微小学校长张彦祥的履历
- 广州白云区四海一家自助餐 广州有哪些可以吃自助餐的地方?
- 德拉诺救出纳兹格雷尔 问一下奥妮克希亚有人型吗。
- 五峰山国际青少年营地 大家觉得最意难平的CP是哪一对呢?
- 乡村爱情谢永强出轨是第几部第几集 乡村爱情9谢永强跳河第几集
- 谁的扣篮霸气 NBA里,谁的扣篮最霸气?
- 大学之道在明明德。“明明德”怎么理解? 明德一路德邦快递地址
- 有哪些充满正能量的电影? 女生节奥莉维亚
- 公共卫生服务责任区域要求 个人开办社区卫生服务站都需要什么?