原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
灰色关联度模型建立和求解
灰色关联度模型的建立与求解游客人数第一步：级比检验首先为了保证建模方法的可行性，需要对已知数据列做必要的检验处理。建立游客人数平均数据时间序列如下：x(0)=(x(0)(1),x(0)(2),…,x(0)(7))=（8438.73,9398.53,9959.17,10949.99,11145.92,11800,12700)求级比 λ(k)λ=(λ(2),λ(3),…λ(7))=（0.8979,0.9437,0.9095,0.9824,0.9446,0.9291）级比判断由于所有的 λ(k)∈[0.8979,0.9824], k=2,3,…,7,故可以用x(0)作满意的GM（1,1）建模。第二步：GM（1,1）建模(1) 对原始数据x(0)做一次累加，即x(1)=（8438.73,17800.15,28088.80,39310.32,51470.10,64573.59,78626.25,93633.56,109601.1). 构造数据矩阵及数据向量，计算 由最小二乘法，求得使达到最小值的 = 于是a=-0.0001，b=8.7076(4)建立模型,相应的白化微分方程为 求解得：求生成数列值及模型还原值.令k=1,2,3,4,5,6,由上面的时间响应函数可算得，其中取由得（8439,9466,10028,.10624,11256,11925,12634)第三步：模型检验模型的各种检验指标值的计算结果见表表2.2 游客人数GM（1,1）模型预测及分析经验证，该模型的精度较高，可进行预测和预报。程序如下：clc,clear x0=[8438.73 9398.53 9959.17 10949.99 11145.92 11800 12700]; n=length(x0); lamda=x0(1:n-1)./x0(2:n) range=minmax(lamda) x1=cumsum(x0) for i=2:n z(i)=0.5*(x1(i)+x1(i-1)); end B=[-z(2:n),ones(n-1,1)]; Y=x0(2:n); u=B\Y x=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 %计算级比偏差值
文档评论（0）