function?main（）
InDim?=?4;??????????????????????????????%輸入維數
OutDim?=?1;?????????????????????????????%輸出維數???????????????????
SamNum?=?3000;??????????????????????????%訓練樣本數
TestSamNum?=?500;???????????????????????%測試樣本數
SP?=?0.6;???????????????????????????????%隱節點擴展常數
ErrorLimit?=?0.9;???????????????????????%目標誤差
x（1）?=?1.2;?????????????????????????????%初始值從1開始
%訓練起始點
StartSamNum?=?202;
%?訓練終止點
EndSamNum?=?StartSamNum?+?SamNum?-1;
%測試樣本起始點
TestStartNum?=?5002;
%測試樣本終止點
EndTestNum?=?TestStartNum?+?TestSamNum?-?1;
%Mackey-Glass?序列終止點
EndMG?=?5600;
StartSamOutNum?=?StartSamNum?+?85;?????????????????????????%期望輸出起始點
StartTestSamOutNum?=?TestStartNum?+85;?????????????????????%測試期望輸出起始點
EndSamOutNum?=?EndSamNum?+?85;?????????????????????????????%期望輸出終止點
EndTestSamOutNum?=?EndTestNum+85;??????????????????????????%測試期望輸出終止點
for?t?=?1:EndMG
if?t?<=?17
????x（t+1）?=?x（t）?-0.1*x（t）;
else
???x（t+1）?=?x（t）?+?（0.2?*?x（t-17）/（1+x（t-17）^10）-0.1*x（t））;
end
end
%figure
%plot（1:t+1x）
%SumIn?=[x（t-18）x（t-12）x（t-6）x（t）]?;?
%TestSumIn?=?x（:5002:5501）;
SamIn?=?[];
for?t?=?StartSamNum:EndSamNum
?????SamInOne?=?[x（t-18）x（t-12）x（t-6）x（t）]‘;
?????SamIn?=?[SamIn?SamInOne];?????????????????????????????????????????%訓練樣本輸入4*3000；
end
TestSamIn=[];
for?t?=?TestStartNum:EndTestNum
????TestSamInOne?=??[x（t-18）x（t-12）x（t-6）x（t）]‘;
????TestSamIn?=?[TestSamIn?TestSamInOne]?;?????????????????????????????%測試樣本輸入4*500
end
rand（‘state‘sum（100*clock））
NoiseVar?=?0.1;
Noise?=?NoiseVar?*?randn（1SamNum）;
SamOutNNoise?=?[];
for?i?=?StartSamOutNum:EndSamOutNum??????????????????????????????????????%無噪聲期望輸出
????SamOutNNoise?=?[SamOutNNoise?x（i）];
end
SamOut?=??SamOutNNoise?+?Noise;????????????????????????????????????%添加噪聲的期望輸出
TestSamOut?=?[];
for?i