amam_ampm.rar

amam_ampm.rar

%%确定被模拟功放的模型参数，并画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