%This?program?is?useful?to?forecaste?the?ridership?of?Beijing?
%public?transportation


%?Get?the?input?sampledata?from?the?.txt?file?“sample.txt“
fid=fopen（‘sample.txt‘‘rt‘）;
%fopend?打開文件函數(shù)
originalData=fscanf（fid‘%f‘?[2014]）;
%fscanf?將文件轉(zhuǎn)成矩陣
status=fclose（fid）;
%fclose?關(guān)閉文件函數(shù)
%?transport?sampledata?into?input?training?sampledata
inputSampledata=originalData‘;
%‘矩陣轉(zhuǎn)置

%?Get?the?output?sampledata?from?the?.txt?file?“goal.txt“
fid=fopen（‘goal.txt‘‘rt‘）;
outputData=fscanf（fid‘%f‘）;
status=fclose（fid）;
%?transport?sampledata?into?input?training?sampledata
outputSampledata=outputData‘;

%?creating?neural?network?and?setting?trainging?parameters
gwwnet=newff（minmax（inputSampledata）[41]{‘tansig‘‘purelin‘}‘traingdm‘）;
%newff?建立神經(jīng)網(wǎng)絡(luò)函數(shù)
%minmax?得到最小值和最大值
%traingdm?下降法

gwwnet.trainParam.show?=?50;
%每50循環(huán)顯示一次
gwwnet.trainParam.lr?=?0.05;
%步長
gwwnet.trainParam.epochs?=?50000;
%50000次停止學(xué)習(xí)
gwwnet.trainParam.goal?=?1e-3;
%le-3停止學(xué)習(xí)

%data?scaling?（converting?the?network?input?and?output?data?to?the?intervel?[-11]）
[inputmininputmaxinputoutputminoutputmaxoutput]?=?premnmx（inputSampledataoutputSampledata）;
%premnmx?-1到1之間的數(shù)

%training
[gwwnettr]=train（gwwnetinputoutput）;
%訓(xùn)練
y=sim（gwwnetinput）;
%仿真

%data?offset?（converting?the?network?output?data?to?it?original?unit）
nnoutput?=?postmnmx（yminoutputmaxoutput）;
%postmnmx?從-1到1返回實(shí)際值

%plot?畫圖
time=1978:1:1997;
plot（timeoutputSampledata‘-‘timennoutput‘o‘）;
%legend（‘a(chǎn)ctual?output‘‘NN?output‘）;
xlabel（‘time‘）;ylabel（‘Learning?fitting?curve‘）;


%scenario1?forecasting?process?預(yù)測輸入
column=10;
for?i=1:column;????
??????SceInput（1i）=inputSampledata（120）*（1.0464^i）;
??????SceInput（2i）=inputSampledata（220）*（1.0631^i）;
??????SceInput（3i）=inputSampledata（320）*（1.0872^i）;
??????SceInput（4i）=inputSampledata（420）*（1.2044^i）;
??????SceInput（5i）=inputSampledata（520）*（1.2326^i）;
??????SceInput（6i）=inputSampledata（620）*（1.0605^i）;
??????SceInput（7i）=2*（1.01^i）;
??????SceInput（8i）=42*（1.02^i）;
??????SceInput（9i）=inputSampledata（920）*（1.1426^i）;
??????SceInput（10i）=inputSampledata（1020）*（1.017^i）;
??????SceInput（11i）=inputSampledata（1120）*（1.0205^i）;
??????SceInput（12i）=inputSampledata（1220）*（1.1336^i）;
??????SceInput（13i）=inputSampledata（1320）*（1.1599^i）;
??????SceInput（14i）=inputSampledata（1420）*（1.1783^i）;
?end

for?j=1:20;
???for?i=1:14;
??????recalldata（ij）=inputSampledata（ij）;
???end
end

for?j=21:30;
???for?i=1:14;
??????recalldata（ij）=SceInput（ij-20）??
??end
end

[alterinputmininputmaxinput]?=?premnmx（recalldata）;

%alterinput和input發(fā)生變化
%可以制定一個最大值和最小值


%training
fvalue=sim（gwwnetalterinput）;

%data?offset?（converting?the?network?output?data?to?it?original?unit）
forecastvalue?=?postmnmx（fvalueminoutputmaxoutput）;

%plot

waitforbuttonpress;
clf;
time=1978:1:2007;
time1=1978:1:1997;
plot（timeforecastvalue‘o‘time1outputSampledata‘-‘）;
%legend（‘預(yù)測曲線‘‘實(shí)際曲線‘

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

?????文件???????3183??2007-03-21?17:43??matlab\a-publictran1.m

?????文件????????161??2003-03-29?14:56??matlab\ahp.m

?????文件????????161??2003-04-01?14:59??matlab\ahpAB.m

?????文件????????274??2003-04-01?15:05??matlab\ahpB1C.m

?????文件????????218??2003-04-01?15:15??matlab\ahpB2C.m

?????文件????????229??2003-04-01?15:25??matlab\ahpB3C.m

?????文件????????318??2003-03-01?18:56??matlab\forecast.m

?????文件???????2832??2003-03-10?09:29??matlab\gln.m

?????文件????????178??1997-04-27?11:20??matlab\goal.txt

?????文件?????????34??2008-03-11?21:02??matlab\gwfile.txt

?????文件???????2831??2009-04-07?11:32??matlab\publictran.m

?????文件???????3137??2007-03-21?17:39??matlab\publictran1.asv

?????文件????????725??2003-03-25?15:05??matlab\ridership.m

?????文件???????1788??1997-04-27?11:19??matlab\sample.txt

?????文件????????455??2003-03-09?21:23??matlab\sample1.m

?????文件????????327??2003-03-09?00:28??matlab\sample2.m

?????文件????????970??2007-03-21?17:48??matlab\smooth.m

?????文件???????2838??2006-01-10?19:56??matlab\student2.m

?????文件????????427??2006-01-10?19:19??matlab\students.m

?????文件????????320??2003-02-14?12:12??matlab\XORNET.M

?????文件????????401??2009-05-04?22:56??matlab\yiyuanxianxinghuigui.m

?????文件??????73216??2003-03-09?11:32??matlab\歷史統(tǒng)計數(shù)據(jù)表.doc

?????文件????????229??2003-04-01?15:25??matlab\復(fù)件?ahpB3C.m

?????文件????????821??2011-06-09?10:04??matlab\徐榮.txt

?????文件??????34816??2009-04-07?11:44??matlab\程序說明.doc

?????文件?????189952??2014-03-14?20:45??matlab\系統(tǒng)預(yù)測上機(jī)實(shí)踐2010.doc

?????目錄??????????0??2014-03-14?21:05??matlab

-----------?---------??----------?-----??----

???????????????320841????????????????????27



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