%=============================================
%This?is?the?PAPR?simulator?for?SC-FDMA.
%=============================================
function?paprSCFDMA（）

dataType=‘16QAM‘;
totalSubcarriers=512;
numSymbols=16;
u=25;
A=0.03;
Q=totalSubcarriers/numSymbols;
Q_tilda=31;
subcarrierMapping=‘LFDMA‘;
pulseShaping=1;
filterType=‘rc‘;
rolloffFactor=0.0999999999;
%rolloffFactor=0.4;
Fs=5e6;
Ts=1/Fs;
Nos=4;

if?filterType==‘rc‘psFilter=rcPulse（TsNosrolloffFactor）;
elseif?filterType==‘rr‘psFilter=rrcPulse（TsNosrolloffFactor）;
end
numRuns=1e2;
papr=zeros（1numRuns）;
for?n=1:numRuns
????if?dataType==‘Q-PSK‘
????tmp=round（rand（numSymbols2））;
????????tmp=tmp*2-1;
????????data=（tmp（1:）+j*tmp（2:））/sqrt（2）;
????elseif?dataType==‘16QAM‘
????????dataSet=[-3+3i?-1+3i?1+3i?3+3i?...
????????????-3+i?-1+i?1+i?3+i?...
????????????-3-i?-1-i?1-i?3-i?...
????????????-3-3i?-1-3i?1-3i?3-3i];
????????dataSet=dataSet/sqrt（mean（abs（dataSet）.^2））;
????????tmp=ceil（rand（1numSymbols）*16）;
????????for?k=1:numSymbols
????????????if?tmp（k）==0
????????????????tmp（k）=1;
????????????end
????????????data（k）=dataSet（tmp（k））;
????????end
????????data=data.‘;
????end
????
????X=fft（data）;
????
????Y=zeros（1totalSubcarriers）;
????if?subcarrierMapping==‘IFDMA‘
????????Y（1:Q:totalSubcarriers）=X;
????elseif?subcarrierMapping==‘LFDMA‘
????????Y（1:numSymbols）=X;
????elseif?subcarrierMapping==‘DFDMA‘
????????Y（1:Q_tilda:Q_tilda*numSymbols）=X;
????end
????
????
????y=ifft（Y）;
????
????if?pulseShaping==1
?