改进SLM程序.doc

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); endenddd=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=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=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(indexx,tmp); fntmp1=f(index); fntmp2=conj(f_tmp(index); fntmp3=fntmp1.*fntmp2; sum1=sum(fntmp3); sum2=sum(abs(fntmp1).2); fenzi=real(sum1); fenmu=sum2; betaopt=fenzi./fenmu; % betaopt=abs(betaopt);%完成opt的求解 % T=abs(2-betaopt.*D).*X)-abs(betaopt.*D.*X); % T1,indextmp=sort(T,descend);%降序排列 T1,indextmp=sort(T);%升序排列 SS=sort(indextmp(1:Isize); alph=zeros(1,length(X); % for iii=1:length(alph) for jj=1:length(SS) if (iii=SS(jj) alph(iii)=2; end end end%产生SS元素 % C=alph.*X; C1=UtraSample(C,J); c=ifft(C1,J*N);%计算峰值抵消信号 x=x-c;%更新x X=X-C;%更新X Signal_Power = abs(x.2); Mean_Power = mean(Signal_Power,2); Peak_Power = max(Sig