时空分布预测模程序.docVIP

下载本文档

1
0
约3.32万字
约 33页
2019-05-29 发布于江苏
举报
版权申诉

时空分布预测模程序.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

一、做散点图： clear all; x1=xlsread(18-24号时空分布.xls,1); x2=xlsread(18-24号时空分布.xls,2); x3=xlsread(18-24号时空分布.xls,3); x4=xlsread(18-24号时空分布.xls,4); x5=xlsread(18-24号时空分布.xls,5); x6=xlsread(18-24号时空分布.xls,6); x7=xlsread(18-24号时空分布.xls,7); for m=1:64 X=[x1(m),x2(m),x3(m),x4(m),x5(m),x6(m),x7(m)];%a-a交通量一周数据序列 m plot(X,k-) %hold on figure end xlabel(时间 /天) ylabel(各小区之间交通量 /辆) 二：用灰色数列预测：G（1,1）模型预测：由于时空分布数据不呈现出一定的规律，比较适合使用灰色预测模型进行预测。 clear all; x1=xlsread(18-24号时空分布.xls,1); x2=xlsread(18-24号时空分布.xls,2); x3=xlsread(18-24号时空分布.xls,3); x4=xlsread(18-24号时空分布.xls,4); x5=xlsread(18-24号时空分布.xls,5); x6=xlsread(18-24号时空分布.xls,6); x7=xlsread(18-24号时空分布.xls,7); for m=1:64 m X=[x1(m),x2(m),x3(m),x4(m),x5(m),x6(m),x7(m)];%a-a交通量一周数据序列 %对X数据序列作一次累加生成X1 X1=zeros(1,7); X1(1)=X(1); for i=2:7 X1(i)=X1(i-1)+X(i); end B=zeros(6,2); for j=1:6 B(j,1)=-0.5*(X1(j+1)+X1(j)); end B=[B(:,1),ones(6,1)]; Yn=[X(2),X(3),X(4),X(5),X(6),X(7)]; a=inv((B*B))*B*Yn ;%最小二乘法求解系数a %X_k+1=(X(0)-a(2)/a(1))*e^(-a(1)*k+a(2)/a(1); X_y=zeros(1,9); k=1:8; X_1=(X(1)-a(2)/a(1))*exp(-a(1)*k)+a(2)/a(1) %建立GM(1,1)模型 X_1=[X1(1),X_1]; %预测 for k=1:8 X_y(k+1)=X_1(k+1)-X_1(k); end X_y(1)=X_1(1) plot(X,k-) hold on plot(X_y, r-) figure e=(X-X_y(1:7))./X %真值与预测值得绝对误差 delt_e=sum(e)/7 %平均相对误差 delt_e0.01则精度为一级 ss=0; ss1=0; for k=2:6 ss=ss+(X(k)-X(1)); ss1=ss1+(X_y(k)-X_y(1)); end s=abs(ss+0.5*(X(7)-X(1))); s1= abs(ss1+0.5*(X_y(7)-X_y(1))); delt_s=abs(ss+0.5*(X(7)-X(1))- ss1-0.5*(X_y(7)-X_y(1))) w=(1+s+s1)/(1+s+s1+delt_s) %w为关联度,w0.90关联度为一级。 C=(mean((e-delt_e).^2) /mean((X-mean(X)).^2))^0.5 %均方差比 C0.35为一级 end m = 1 X_1 = 1.0e+004 * Columns 1 through 6 1.3931 2.0409 2.7020 3.3767 4.0653 4.7681 Columns 7 through 8 5.4854 6.2175 X_y = 1.0e+003 * Columns 1 through 6 7.5850 6.3465 6.4773 6.6108 6.7471 6.8861 Columns 7 through 9 7.0281 7.1729 7.3208 e = Columns 1 through 6 0 0.1085 -0.1189 -0.0985 0.0500 0.0066 Column 7 0.0155 delt_e