%%?第一種神經(jīng)網(wǎng)絡(luò)：BP神經(jīng)網(wǎng)絡(luò)--幅度/相位模型
%%?清空環(huán)境變量
clc
clear
set（0‘defaultfigurecolor‘‘w‘）;
%%?訓(xùn)練數(shù)據(jù)預(yù)測(cè)數(shù)據(jù)提取
%下載輸入輸出數(shù)據(jù)
load（‘AMAM9040_1.txt‘）
load（‘PMPM9040_1.txt‘）
%將數(shù)據(jù)分類(lèi)存儲(chǔ)
A?=?AMAM9040_1（:2）;
P?=?PMPM9040_1（:2）;
T?=?AMAM9040_1（:1）;
T0?=?T（1:300）‘;
Ain?=?A（1:300）;
Aout?=?A（301:600）;
Pin?=?P（1:300）;
Pout?=?P（301:600）;

%顯示原始數(shù)據(jù)圖形
plot（T0Ain）;
hold?on
plot（T0Aout）;
title（‘9040\_原始時(shí)域輸入輸出幅度‘）
xlabel（‘時(shí)間/S‘）;
ylabel（‘幅度/V‘）;

grid?on?
figure
plot（T0Pin）;
hold?on
plot（T0Pout）;
title（‘9040\_原始時(shí)域輸入輸出相位‘）
xlabel（‘時(shí)間/S‘）;
ylabel（‘相位/角度‘）;
grid?on?

Pin?=?Pin*pi/180;
Pout?=?Pout*pi/180;
for?i?=?1:300
????if?Pout（i）<-0.81
????????Pout（i）=Pout（i）+2*pi;
????end
end

figure
plot（AinAout‘r.‘）;
title（‘9040\_原始輸入輸出幅度關(guān)系‘）
xlabel（‘輸入幅度/V‘）;
ylabel（‘輸出幅度/V‘）;
legend（‘原數(shù)據(jù)‘‘預(yù)測(cè)數(shù)據(jù)‘）;
grid?on?
figure
plot（PinPout‘r.‘）;
title（‘9040\_原始輸入輸出相位關(guān)系‘）
xlabel（‘輸入相位/弧度‘）;
ylabel（‘輸入相位/弧度‘）;
legend（‘輸入相位‘‘輸出相位‘）;
grid?on?

Ain_train?=?Ain（1:200）‘;
Pin_train?=?Pin（1:200）‘;
Aout_train?=?Aout（1:200）‘;
Pout_train?=?Pout（1:200）‘;
Vinput_train?=?[Ain_train;Pin_train];
Voutput_train?=?[Aout_train;Pout_train];

Ain_test?=?Ain（201:300）‘;
Pin_test?=?Pin（201:300）‘;
Aout_test?=?Aout（201:300）‘;
Pout_test?=?Pout（201:300）‘;
Vinput_test?=?[Ain_test;Pin_test];
Voutput_test?=?[Aout_test;Pout_test];


%%?BP網(wǎng)絡(luò)訓(xùn)練
%?%初始化網(wǎng)絡(luò)結(jié)構(gòu)
net=newff（Vinput_trainVoutput_train5）;

net.trainParam.epochs=100;
net.trainParam.lr=0.1;
net.trainParam.goal=0.00004;

%網(wǎng)絡(luò)訓(xùn)練
net=train（netVinput_trainVoutput_train）;
%netBPAP780_5挺好的
save（‘D:\Study\Graduation_project\Code\Mycode\netBPAP9040_1.mat‘‘net‘）;
%load???netBPAP9040???%net為已保存的網(wǎng)絡(luò)
%%?BP網(wǎng)絡(luò)預(yù)測(cè)

%網(wǎng)絡(luò)驗(yàn)證輸出
BPoutput_train?=?sim（netVinput_train）;

%網(wǎng)絡(luò)預(yù)測(cè)輸出
BPoutput_test?=?sim（netVinput_test）;
%?預(yù)測(cè)輸出整理
A_BPoutput_train?=?BPoutput_train（1:）;
P_BPoutput_train?=?BPoutput_train（2:）;
A_BPoutput_test?=?BPoutput_test（1:）;
P_BPoutput_test?=?BPoutput_test（2:）;
%%?結(jié)果分析
%?我們使用圖像來(lái)看網(wǎng)絡(luò)對(duì)非線性函數(shù)的擬合效果

figure
plot（Ain_trainAout_train‘r.‘）;
hold?on
plot（Ain_trainA_BPoutput_train‘bo‘）;
title（‘9040\_BP網(wǎng)絡(luò)驗(yàn)證輸出幅度‘）
xlabel（‘Ainput‘）
ylabel（‘Aoutput‘）
grid?on?

figure
plot（Pin_trainPout_train‘r.‘）;
hold?on
plot（Pin_trainP_BPoutput_train‘bo‘）;
title（‘9040\_BP網(wǎng)絡(luò)驗(yàn)證輸出相位‘）
xlabel（‘Pinput‘）
ylabel（‘Poutput‘）
grid?on?

figure
plot（Ain_testAout_test‘r.‘）;
hold?on
plot（Ain_testA_BPoutput_test‘bo‘）;
title（‘9040\_BP網(wǎng)絡(luò)預(yù)測(cè)輸出幅度‘）
xlabel（‘Ainput/V‘）
ylabel（‘Aoutput/V‘）
legend（‘原數(shù)據(jù)‘‘預(yù)測(cè)數(shù)據(jù)‘）;
grid?on?

figure
plot（Pin_testPout_test‘r.‘）;
hold?on
plot（Pin_testP_BPoutput_test‘bo‘）;
title（‘9040\_BP網(wǎng)絡(luò)預(yù)測(cè)輸出相位‘）
xlabel（‘Pinput/rad‘）
ylabel（‘Poutput/rad‘）
legend（‘原數(shù)據(jù)‘‘預(yù)測(cè)數(shù)據(jù)‘）;
grid?on?

%%?觀察誤差
A_BPerroproportion?=?abs（（A_BPoutput_test-Aout_test）./Aout_test）;
P_BPerroproportion?=?abs（（P_BPoutput_test-Pout_test）./Pout_test）;
A_BPnum1_10?=0;
P_BPnum1_10?=0;
for?i?=?1:100
????if?A_BPerroproportion（i）?<=?0.1
????????A_BPnum1_10=A_BPnum1_10+1;
????end?
????if?P_BPerropro

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????27660??2018-05-30?21:20??AMAM9040_1.txt
?????文件????????4061??2018-06-10?15:47??BP_AP_9040_1.m
?????文件????????3624??2018-06-10?17:09??Elman_AP_9040_1.m
?????文件????????9911??2018-06-10?16:01??netBPAP9040_1.mat
?????文件?????????362??2018-05-29?20:44??NRMS.m
?????文件???????27917??2018-05-30?21:23??PMPM9040_1.txt
?????文件????????3308??2018-06-03?17:38??RBF_AP_9040_1.m
?????文件??????196329??2018-05-16?15:26??Rd_Nu.mat
?????文件??????????51??2018-05-16?15:26??Rd_Num.m
?????文件?????????750??2019-02-28?18:18??readme.txt
?????文件?????????267??2018-05-29?20:43??RMSE.m