clc
clear
close?all



%?產生訓練樣本與測試樣本
%?訓練樣本
P1=load（‘P11.txt‘）;
T1=load（‘T11.txt‘）;

???M=200;??
???
??A2=zeros（M1）;?
??A=zeros（M2）;
?for?i=1:M
?????????A2（i）=-T1（i）.*（P1（i））.^2;
?end
?A（:1）=P1;
?A（:2）=A2;
????????
?a=A\T1?;???%?a為AM-AM失真的兩個參數
?????
?????
?%?測試樣本
P2=load（‘P21.txt‘）;
T2=load（‘T21.txt‘）;????

M=100;
Y2=zeros（M1）;
for?i=1:M
???Y2（i）=a（1）*P2（i1）./（1+a（2）*P2（i1）.^2）;??
end?


E2=Y2-T2;
SUMAY2=0;
SUME2=0;
M=100;
for?i=1:M
????SUMAY2=SUMAY2+（Y2（i））^2;
????SUME2=SUME2+（E2（i））^2;
end
%?RMSE_E2=sqrt（SUME2/M）
%?RMSE_AY2=sqrt（SUMAY2/M）;
%?NRMS2=10*log10（RMSE_E2/RMSE_AY2）
RMSE_E2=