saleh的仿真实现

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=