功放非线性特性MATLAB仿真

%%确定被模拟功放的模型参数，并画AM-AM,AM-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(1,SUM); 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(t,y) %plot(t,abs(s_bpsk)) %% 滤波器系数设置 N=50;%滤波器阶数 Wn1=[0.343255,0.443255];%1为fs的一半，10.23 滤波器的工作频率范围 注意滤波器的范围与采样率有关 Wn2=[0.05,0.15]; A=fir1(N,Wn1,'bandpass');%带通滤波器 B=fir1(N,Wn2,'bandpass'); s_bpsk=filter(B,1,s_bpsk);%带通滤波（A为分子，1为分母 % figure(2) % plot(t,abs(s_bpsk)) c1=[1.513+0.0904j,-0.007-0.001j,0.0289+0.0054j,-0.5-0.3j,0.2+0.2j,-0.0621-0.0932j,-0.9-0.6j,-0.2451-0.3735j,0.2+0.1j]; HPA_s1=volterra(c1,s_bpsk); h=spectrum.welch; hpsd=psd(h,HPA_s1,'fs',fs); %经标准功放模型后的信号功率谱密度 hold on; plot(hpsd) g=zeros(1,SUM); sita=zeros(1,SUM); sita1=zeros(1,SUM); sita2=zeros(1,SUM); 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(x,g); %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(x,y);