改进SLM程序.docVIP

clc; clf; clear all; close all; A = 4; mod_mode = 4; N = 128; M = sqrt(mod_mode)/2; J = 4; L = 8; MAX_SYMBOLS = 10e4; PAPR_Orignal = zeros(1,MAX_SYMBOLS); PAPR_SLM = zeros(1,MAX_SYMBOLS); PAPR_tmp = zeros(L,MAX_SYMBOLS); for m=1:M for n=1:M s1(m,n)=sqrt((m-1).^+(n-1).^2); end end dd=4/mod_mode*sqrt(6/(mod_mode-1))*sum(sum(s1)); dd=sqrt(N)/(2*dd); %%% tic; for nSymbol=1:MAX_SYMBOLS bitsource=randint(1,N*sqrt(mod_mode)); X=modulation(bitsource,mod_mode); X1=UtraSample(X,J); x=ifft(X1,J*N); Signal_Power = abs(x.^2); Mean_Power = mean(Signal_Power,2); Peak_Power = max(Signal_Power,[],2); PAPR_temp = 10*log10(Peak_Power./Mean_Power); PAPR_Orignal(nSymbol)=PAPR_temp; eita=10*log10(PAPR_temp); threshold=10^(A/10)*Mean_Power;%%%关键 if (eita=A) PAPR_SLM(nSymbol)=PAPR_temp; else x_iter = zeros(J*N,L); for iter=1:L SS=[]; f=zeros(size(x)); for ii=1:length(x); if (Signal_Power(ii)=threshold); f(ii)=x(ii); else f(ii)=0; end end F=fft(f,J*N); F=[F(1:N/2),F((J*N-N/2+1):(J*N))]; D=zeros(1,N); D=real(F.*conj(X)./((abs(X)).^2)); if (all(D=0)) flagg=0; break;%%%不对 else flagg=1; end IdB=eita-A; Isize=ceil(dd*IdB./2); FFF=D.*X; FFF=UtraSample(FFF,J); f_tmp=ifft(FFF,J*N); indexx=find(Signal_Power=threshold); tmp=[]; iii=1; ftmp=[0 f 0]; for jj=1:length(f) if ((abs(ftmp(jj+1))abs(ftmp(jj)))(abs(ftmp(jj+1))=abs(ftmp(jj+2)))) tmp(iii)=jj; iii=iii+1; end end %Sp %%% index=intersect(