%1-24622463-4443?時間序列
clear;?%清除變量
clc;

[Ydate]=xlsread（‘Y18j.xls‘）;

%p=sum（y）/length（y）;?%預處理
%td=std（y）;?%標準偏差
%for?i=1:length（y）
%???y（i）=（y（i）-p）/td;
%end

T=datenum（date）;

y=Y（1:2462）;
t=T（1:2462）;
y2=Y（2463:4443）;
t2=T（2463:4443）;

th=ivar（y8）;%最優輔助變量法建模
a=th2arx（th）;%轉換模型參數
data1=predict（yth）;%?對1-2462數據進行預測
data2=predict（Y（2463:4443）th）;%?對2463:4443數據進行預測
e?=pe（yth）;%計算預測誤差

%plot
scrsz?=?get（0‘ScreenSize‘）;
figure（‘Position‘[1?scrsz（4）/2?scrsz（3）?scrsz（4）/2.5]）??%[left?bottom?width?height]
plot（ty?‘k‘‘linewidth‘2）
hold?on
plot（?tdata1{11}‘:r‘‘linewidth‘1）%畫圖比較


legend（‘原始‘‘擬合‘）

datetick（‘x‘‘YYYY‘）;
set（gca‘YTick‘-10:1.5:5‘fontsize‘13）;

ylabel（‘位移（mm）‘）??????????%坐標標簽
title（‘18壩段正下徑向時間位移實測與擬合曲線‘）??????????????????%題目

%誤差

Q=0;
U=0;
Syy=0;
L=length（y）;
for?i=1:L
????Q=Q+（y（i）-data1{11}（i））^2;??%剩余（殘差）平方和
????U=U+（data1{11}（i）-mean（y））^2;
????Syy=Syy+（y（i）-mean（y））^2;

end
????


S=sqrt（Q/（L-10））;?%剩余標準差
R=sqrt（U/（U+Q））;?%復相關系數
RMSE=sqrt（Q/L）;
error=y-data1{11};?%殘差

scrsz?=?get（0‘ScreenSize‘）;
figure（‘Position‘[1?scrsz（4）/2?scrsz（3）?scrsz（4）/3]）??%[left?bottom?width?height]
plot（terror‘k‘‘linewidth‘2）;
legend（‘殘差‘）;?

datetick（‘x‘‘YYYY‘）;
set（gca‘YTick‘-2:0.5:2‘FontSize‘13）;



%legend（‘18壩段正下徑向‘）?????????????%曲線的‘外號’
%xlabel（‘年份‘）
ylabel（‘位移（mm）‘）??????????%坐標標簽
title（‘位移實測-擬合值殘差曲線‘‘FontSize‘16）??????????????????%題目

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%plot
scrsz?=?get（0‘ScreenSize‘）;
figure（‘Position‘[1?scrsz（4）/2?scrsz（3）?scrsz（4）/2.5]）??%[left?bottom?width?height]

plot（t2y2?‘k‘‘linewidth‘2）
hold?on
plot（?t2data2{11}‘:r‘‘linewidth‘1）%畫圖比較

legend（‘原始‘‘預測‘）

datetick（‘x‘‘YYYY‘）;
set（gca‘YTick‘-10:1.5:5‘fontsize‘13）;

ylabel（‘位移（mm）‘）??????????%坐標標簽
title（‘18壩段正下徑向時間位移實測與擬合曲線‘）??????????????????%題目

%誤差

Q2=0;
U2=0;
Syy2=0;
L2=length（y2）;
for?i=1:L2
????Q2=Q2+（y2（i）-data2{11}（i））^2;??%剩余（殘差）平方和
????U2=U2+（data2{11}（i）-mean（y2））^2;
????Syy2=Syy2+（y2（i）-mean（y2））^2;

end
????


S2=sqrt（Q2/（L2-10））;?%剩余標準差
R2=sqrt（U2/（U2+Q2））;?%復相關系數
RMSE2=sqrt（Q2/L2）;
error2=y2-data2{11};?%殘差

scrsz?=?get（0‘ScreenSize‘）;
figure（‘Position‘[1?scrsz（4）/2?scrsz（3）?scrsz（4）/3]）??%[left?bottom?width?height]
plot（t2error2‘k‘‘linewidth‘2）;
legend（‘殘差‘）;?

datetick（‘x‘‘YYYY‘）;
set（gca‘YTick‘-2:0.5:2‘FontSize‘16）;

%legend（‘18壩段正下徑向‘）?????????????%曲線的‘外號’
%xlabel（‘年份‘）
ylabel（‘位移（mm）‘）??????????%坐標標簽
title（‘位移實測-擬合值殘差曲線‘‘FontSize‘16）??????????????????%題目

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1008??2009-05-19?00:09??模型代碼+圖形+結果\時間序列模型\（matlab在時間序列分析中的應用）.asv

?????文件???????2866??2009-05-23?15:21??模型代碼+圖形+結果\時間序列模型\（matlab在時間序列分析中的應用）3.asv

?????文件???????2866??2009-05-23?15:24??模型代碼+圖形+結果\時間序列模型\（matlab在時間序列分析中的應用）3.m

?????文件?????107704??2009-05-23?14:58??模型代碼+圖形+結果\時間序列模型\1-2462.JPG

?????文件??????62842??2009-05-23?14:59??模型代碼+圖形+結果\時間序列模型\1-2462e.JPG

?????文件??????55631??2009-05-23?14:59??模型代碼+圖形+結果\時間序列模型\2463-4443.JPG

?????文件??????51282??2009-05-23?14:59??模型代碼+圖形+結果\時間序列模型\2463-e.JPG

?????文件??????51360??2009-05-23?14:00??模型代碼+圖形+結果\時間序列模型\acf.jpg

?????文件?????209328??2009-05-23?14:43??模型代碼+圖形+結果\時間序列模型\AIC+BIC.jpg

?????文件??????36657??2009-05-23?14:01??模型代碼+圖形+結果\時間序列模型\pacf.jpg

?????文件??????40448??2009-05-23?15:20??模型代碼+圖形+結果\時間序列模型\tsa.doc

?????文件?????210432??2009-05-17?19:43??模型代碼+圖形+結果\時間序列模型\Y18j.xls

?????文件????????956??2009-05-18?20:19??模型代碼+圖形+結果\時間序列模型\基于MATLAB的AR模型參數估計.asv

?????文件???????1037??2009-05-23?15:21??模型代碼+圖形+結果\時間序列模型\基于MATLAB的AR模型參數估計2.asv

?????文件???????1037??2009-05-23?15:24??模型代碼+圖形+結果\時間序列模型\基于MATLAB的AR模型參數估計2.m

?????文件????????347??2009-05-22?09:35??模型代碼+圖形+結果\時間序列模型\拖尾截尾時間序列.asv

?????文件????????416??2009-05-23?15:21??模型代碼+圖形+結果\時間序列模型\拖尾截尾時間序列1.asv

?????文件????????416??2009-05-23?15:24??模型代碼+圖形+結果\時間序列模型\拖尾截尾時間序列1.m

?????文件?????171394??2009-05-23?14:43??模型代碼+圖形+結果\時間序列模型\方差+FIE.jpg

?????文件?????160138??2009-05-23?21:21??模型代碼+圖形+結果\神經網絡模型\1.jpg

?????文件?????112750??2009-05-23?21:21??模型代碼+圖形+結果\神經網絡模型\1e.jpg

?????文件??????96580??2009-05-23?21:21??模型代碼+圖形+結果\神經網絡模型\2.jpg

?????文件??????73350??2009-05-23?21:21??模型代碼+圖形+結果\神經網絡模型\2e.jpg

?????文件?????172564??2009-05-23?21:22??模型代碼+圖形+結果\神經網絡模型\3.jpg

?????文件???????4004??2009-05-23?21:12??模型代碼+圖形+結果\神經網絡模型\3.m

?????文件?????107773??2009-05-23?21:22??模型代碼+圖形+結果\神經網絡模型\3e.jpg

?????文件?????125243??2009-05-06?16:56??模型代碼+圖形+結果\神經網絡模型\all.csv

?????文件??????46432??2009-05-23?19:03??模型代碼+圖形+結果\神經網絡模型\bp.jpg

?????文件???????4180??2009-05-23?20:57??模型代碼+圖形+結果\神經網絡模型\bp神經網絡.asv

?????文件???????4180??2009-05-23?20:57??模型代碼+圖形+結果\神經網絡模型\bp神經網絡.m

............此處省略45個文件信息