%%%使用BPTT算法完成Elman網(wǎng)絡(luò)的訓(xùn)練
%%?程序說(shuō)明
%??1、數(shù)據(jù)為7天，四個(gè)時(shí)間點(diǎn)的空調(diào)功耗，用前三個(gè)推測(cè)第四個(gè)訓(xùn)練，依次類(lèi)推。第七天作為檢驗(yàn)
%??2、Elman網(wǎng)絡(luò)輸入結(jié)點(diǎn)為12，輸出結(jié)點(diǎn)為4個(gè)，隱藏結(jié)點(diǎn)18個(gè)
clear?all;
clc;
%%?加載數(shù)據(jù)并進(jìn)行歸一化處理
[train_datatest_data]=Elman_data_process（）;
data_length=size（train_data1）;????????????%每個(gè)樣本的長(zhǎng)度
data_num=size（train_data2）;???????????????%樣本數(shù)目?
%%?網(wǎng)絡(luò)參數(shù)初始化
%結(jié)點(diǎn)數(shù)目設(shè)置
input_layer_num=data_length;
hidden_layer_num=18;
output_layer_num=size（test_data1）;
%權(quán)重初始化
weight_input_hidden=rand（hidden_layer_numinput_layer_num）/1000;
weight_prehidden_hidden=rand（hidden_layer_numhidden_layer_num）/1000;
weight_hidden_output=rand（output_layer_numhidden_layer_num）/1000;
output_state=zeros（output_layer_num3）;
%%?網(wǎng)絡(luò)訓(xùn)練學(xué)習(xí)部分??采用BPTT算法
yita=0.001;??????????????%每次迭代的步進(jìn)
for?num=1:4000
for?t=1:3???????%%時(shí)間步，使用三推一
????%從輸入到隱層
????if?（t==1）
????????pre_hidden_state=weight_input_hidden*train_data（:t）;
????else
???????pre_hidden_state=weight_input_hidden*train_data（:t）+weight_prehidden_hidden*hidden_state（:t-1）;
????end

????%從隱層到輸出
????for?n=1:hidden_layer_num
????????hidden_state（nt）=1/（1+exp（-pre_hidden_state（n:）））;??????????????%%通過(guò)sigmoid函數(shù)
????end
????output_state（:t）=weight_hidden_output*hidden_state（:t）;

????%誤差計(jì)算
????Error=output_state（:t）-test_data（:t）;
????Error_cost（1num）=sum（（output_state（:t）-test_data（:t））.^2）;
????if（Error_cost（1num）<1e-4）
????????break;
????end
????%權(quán)值更新
????[weight_input_hiddenweight_prehidden_hiddenweight_hidden_output]=updata_weight（tyitaErrortrain_datahidden_stateweight_input_hiddenweight_prehidden_hiddenweight_hidden_output）;

end
????if（Error_cost（1num）<1e-4）
????????break;
????end
end
%%?繪制COST曲線
for?n=1:1:num
????text（nError_cost（1n）‘*‘）;
????axis（[0num01]）;
????xlabel（‘迭代次數(shù)‘）;
????ylabel（‘誤差‘）;
end
%%?結(jié)果檢驗(yàn)
%%由前1-3天=>第4天
pre_hidden_state=weight_input_hidden*train_data（:1）;
for?n=1:hidden_layer_num
????hidden_state（n1）=1/（1+exp（-pre_hidden_state（n:）））;??????????????%%通過(guò)sigmoid函數(shù)
end
output_state（:1）=weight_hidden_output*hidden_state（:1）;
for?i=4:7
fprintf（‘第%d天前面3天推測(cè)第%d天的結(jié)果%d?\n‘iioutput_state（i-33））;
fprintf（‘真實(shí)第%d天結(jié)果%d?\n‘itest_data（i-34））;
end
%使用第七組檢驗(yàn)
test_final=[0.4557?0.4790?0.7019?0.8211?0.4601?0.4811?0.7101?0.8298?0.4612?0.4845?0.7188?0.8312]‘;
test_final=test_final/sqrt（sum（test_final.^2））;
test_output=test_data（:4）;
pre_hidden_state=weight_input_hidden*test_final+weight_prehidden_hidden*hidden_state（:3）;
for?n=1:hidden_layer_num
????hidden_state1（n1）=1/（1+exp（-pre_hidden_state（n:）））;??????????????%%通過(guò)sigmoid函數(shù)
end
t_output=weight_hidden_output*hidden_state1（:1）;
test_output;
%計(jì)算誤差度
wucha=sum（abs（test_output-t_output）./t_output）/4

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件?????1084464??2019-08-16?00:32??賈偉偉_RNN算法推導(dǎo)過(guò)程及代碼.pdf
?????文件????????2899??2019-08-15?17:21??運(yùn)行代碼\Elman_BPTT.m
?????文件????????1304??2019-08-15?16:42??運(yùn)行代碼\Elman_data_process.m
?????文件????????1324??2019-08-15?16:42??運(yùn)行代碼\updata_weight.m
?????文件??????594299??2019-08-16?00:30??RNN算法推導(dǎo)過(guò)程及代碼.docx