%該函數實現功放的輸入輸出功能
%其中，A為功放的門限;x是一個列向量
%此函數實現硬clipping操作
function y=pa_nonlinear(x,A)
len = length(x);

for it = 1 : len
  rad = abs( x(it) );
  ang = angle( x(it) );
  
  if rad >= A
    x(it) = A * exp( j * ang );
  end
end

y = x;
end

%此函數的功放模型是固態功放模型
% function y = pa_nonlinear( x, A )
% p = 3;
% rad = abs(x);
% ang = angle(x);
% rad = rad ./ ( 1 + ( rad ./ A ) .^ (2 * p) ) .^ ( 1 / ( 2 * p ) );
% ang = ang;
% y = rad .* exp( j * ang );
% end