EMPM.rar_DPD_EMPM_数字预失真_预失真_预失真 多项式

clear clc %预失真算法主函数 s=0; K=5; Q=2; fs=140*1.28e6; tin=load('10-12-00.mat'); %功放输入数据 tout=load('10-12-10.mat'); %功放输出数据 [x,y]=loop_delay(tin,tout); %环路延时算法--幅度差值最小法 x=x(500:2500); y=y(500:2500); %x1=500:2500; %figure; %stairs(x1,x,'r');hold on; %stairs(x1,y,'b');hold on; coef_pd=coef(x,y,s,K,Q); %计算预失真器系数 coef_pa=coef(y,x,s,K,Q); %计算功放模型的系数 z=out(x,coef_pd,s,K,Q); %输入信号经过预失真器的输出 pout=out(z,coef_pa,s,K,Q); %预失真器输出信号，再经过功放得到输出信号 %%%%%%%%%%%%%%%%%%%%%%功率谱密度输出%%%%%%%%%%%%%%%%%%%%%%%%%% [pxx1,ff1]=pwelch(x,hann(256),0,2048,fs,'twosided'); c1=max(pxx1); [pxx2,ff2]=pwelch(y,hann(256),0,2048,fs,'twosided');%未加预失真输出功率谱 c2=max(pxx2); [pxx3,ff3]=pwelch(pout,hann(256),0,2048,fs,'twosided');%加预失真输出功率谱 c3=max(pxx3); y1=10*log10(fftshift(pxx1/c1));%用fftshift将功率谱搬到中心 y2=10*log10(fftshift(pxx2/c2)); y3=10*log10(fftshift(pxx3/c3)); figure; plot(ff1-fs/2,y1,'r->','MarkerSize',2);hold on; plot(ff2-fs/2,y2,'g-<','MarkerSize',2);hold on; plot(ff3-fs/2,y3,'b-*','MarkerSize',2);hold on; axis([-fs/4,fs/4,-70,+5]); xlabel('频率(Hz)'); ylabel('功率谱密度(dB/Hz)'); %%%%—————————星座图输出———————————————%%%%%%% %scatterplot(x);hold on; %scatterplot(y);hold on; %scatterplot(pout);hold on; %%%%%%%————————AM-AM曲线输出—————————————————%%%%% figure; plot(abs(x),abs(y),'*r','MarkerSize',2);hold on; plot(abs(x),abs(pout),'+b','MarkerSize',2);hold on; title('预失真前后AM-AM曲线比较'); xlabel('|in|'); ylabel('|out|'); %%%%%%%————————AM-PM曲线输出—————————————————%%%%% figure; plot(abs(x),unwrap(unwrap(angle(y))-unwrap(angle(x))+2*pi),'*r','MarkerSize',2);hold on; plot(abs(x),unwrap(unwrap(angle(pout))-unwrap(angle(x))),'+b','MarkerSize',2);hold on; title('预失真前后AM-PM曲线比较'); xlabel('|in|'); ylabel('pahse of output');