clear?all;?clc;
%?%?%?SamNum?=?100;
%?%?%?rand（‘state‘sum（100*clock））
%?%?%?Noise=（0.01+0.01*1j）*rand（1SamNum）;
%?%?%?SamIn?=0.01*pi:0.02*pi:2*pi;
%?%?%?SamInputData=sin（SamIn）.‘+j*sin（SamIn）.‘+Noise.‘;
%?%?%?SamOut=SamInputData;
%?%?%?
%?%?%?TestIn?=0.01*pi:0.02*pi:2*pi;
%?%?%?TestInputData=sin（TestIn）.‘+j*sin（TestIn）.‘;
%?%?%?TestOut=TestInputData;
SamNum?=?315;?%?訓練樣本數
%?根據目標函數獲得訓練樣本輸入輸出
rand（‘state‘sum（100*clock））
Noise=（0.01+0.01*1j）*rand（1SamNum）;
SamIn?=0:0.02:2*pi;
SamInputData=exp（1j*SamIn）.‘;
SamOut=1+0.02*exp（1+2*SamInputData）.*（1+0.2*SamInputData+SamInputData.^2）+Noise‘;%SamOut形式必須與TestOut相同

%?根據目標函數獲得測試樣本輸入輸出
TestSamIn?=?0:0.02:2*pi;
TestInputData=exp（1j*TestSamIn）.‘;
TestOut?=1+0.02*exp（1+2*TestInputData）.*（1+0.2*TestInputData+TestInputData.^2）;%TestOut形式必須與SamOut相同

[InDimMaxUnitNum]=size（SamIn）;??%樣本輸入維數和最大允許隱節點數

d=SamOut‘;%期望輸出賦給d
NumSam?=?size（SamInputData1）;?%訓練樣本數
center=SamInputData;
P=[];
??%--------求回歸矩陣P-----------%
??for?icenter=1:MaxUnitNum?%MaxUnitNum為隱層神經元數
??