%%確定被模擬功放的模型參數，并畫AM-AMAM-PM特性曲線
clc;
close?all;??%關閉畫圖窗口
format?long;
%BPSK信號
Gold_Wave=load（‘E:\paper\graduation\Program\sl_simulation\gold_wave.mat‘）;?
Gold_Wave=Gold_Wave.Gold_Wave;?
SUM=20000;
fs=100*10^6;%采樣速率
fc=20*10^6;%載波
fcode=10*10^6;?%碼速率
fc1=60*10^6;???%??功放范圍
ts=1/fs;
tcode=1/fcode;
count=0;
n1=0;
f_in=zeros（1SUM）;
code=Gold_Wave;
for??n=0:SUM-1
?????n1=floor（n*ts/tcode）;
?????f_in（n+1）=code（n1+1）;??%BPSK信號信號的幅度不能超過1（多項式模型決定的，否則輸出隨信號幅度增大迅速增大）
end???
s_bpsk_cos=0.7*f_in.*cos（2*pi*fc*（0:SUM-1）/fs）;
s_bpsk_sin=0.7*f_in.*sin（2*pi*fc*（0:SUM-1）/fs）;
s_bpsk=s_bpsk_cos+j*s_bpsk_sin;
t=1:length（s_bpsk）;
%?y=fft（s_bpsk）;
%?figure（1）;
%?plot（ty）
%plot（tabs（s_bpsk））

%%?濾波器系數設置
N=50;%濾波器階數
Wn1=[0.3432550.443255];%1為fs的一半，10.23??濾波器的工作頻率范圍???注意濾波器的范圍與采樣率有關
Wn2=[0.050.15];

A=fir1（NWn1‘bandpass‘）;%帶通濾波器
B=fir1（NWn2‘bandpass‘）;


s_bpsk=filter（B1s_bpsk）;%帶通濾波（A為分子，1為分母
%?figure（2）
%?plot（tabs（s_bpsk））
c1=[1.513+0.0904j-0.007-0.001j0.0289+0.0054j-0.5-0.3j0.2+0.2j-0.0621-0.0932j-0.9-0.6j-0.2451-0.3735j0.2+0.1j];
HPA_s1=volterra（c1s_bpsk）;
h=spectrum.welch;
hpsd=psd（hHPA_s1‘fs‘fs）;??%經標準功放模型后的信號功率譜密度
hold?on;
plot（hpsd）
g=zeros（1SUM）;
sita=zeros（1SUM）;
sita1=zeros（1SUM）;
sita2=zeros（1SUM）;

for?n=1:SUM
%????g（n）=10*log10（y（n）/x（n））;
???sita1（n）=atan（imag（HPA_s1（n））/real（HPA_s1（n）））;
???sita2（n）=atan（imag（s_bpsk（n））/real（s_bpsk（n）））;

end
figure（4）
plot（abs（s_bpsk）abs（（HPA_s1）））;

?%?plot（10*log10（x）g）;
??%plot（xg）;
???%plot（10*log10（x）10*log10（y））;

??figure（5）
%??plot（abs（s_carri_b）sita‘r‘）
%??hold?on
??plot（abs（s_bpsk）sita1-sita2）

%?plot（xy）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1875??2012-04-23?11:15??amam_ampm.m

?????文件????????413??2012-04-16?16:03??volterra.m

-----------?---------??----------?-----??----

?????????????????2288????????????????????2