?? ber_ofdm_64s_1e4_withclipping_snr20_differcr_1_2_3_nofilter.m
字號:
BitsTx = floor(rand(1,NumLoop*NumSubc)*2);%用rand產(chǎn)生0：1的隨機(jī)數(shù)，不會產(chǎn)生負(fù)值，如果用randn會產(chǎn)生負(fù)值；乘以2是為了使用floor函數(shù)
% --------------------------------------------- %
% QAM MODULATION %
% --------------------------------------------- %
% Modulate (Generates QAM symbols)
% input: BitsTx(1,NumLoop*NumSubc); output:SymQAM(NumLoop,NumSubc/2)
SymQAMtmp = reshape(BitsTx,2,NumLoop*NumSubc/2).';
SymQAMtmptmp = bi2de(SymQAMtmp,2,'left-msb');%/?????(bi2de)
%--------------------------------------------------------------------------
QAMTable = [-1-i -1+i 1-i 1+i];
SymQAM = QAMTable(SymQAMtmptmp+1);
% --------------------------------------------- %
% IFFT input: SymQAM(NumLoop,NumSubc/2); output: SymIFFT(NumSubc,NumLoop)
% --------------------------------------------- %
SymIFFT = zeros(NumSubc,NumLoop);
SymIFFTtmp = reshape(SymQAM,NumSubc/2,NumLoop);%（為什么載波數(shù)除以2）
SymIFFTtmptmp = zeros(NumSubc,NumLoop);
SymIFFTtmptmp(1,:) = real(SymIFFTtmp(1,:)); % 實(shí)數(shù)
SymIFFTtmptmp(NumSubc/2+1,:) = imag(SymIFFTtmp(1,:)); % 實(shí)數(shù)
SymIFFTtmptmp(2:NumSubc/2,:) = SymIFFTtmp(2:NumSubc/2,:);
SymIFFTtmptmp((NumSubc/2+2):NumSubc,:) = flipdim(conj(SymIFFTtmp(2:NumSubc/2,:)),1);
SymIFFT = ifft(SymIFFTtmp,NumSubc,1);%(后邊這個(gè)1是干什么的？)
% --------------------------------------------- %
% Add cyclic prefix %
% --------------------------------------------- %
% input: SymIFFT(NumSubc,NumLoop); output: SymCP(NumSubc + NumCP,NumLoop)
NumAddPrefix = NumSubc + NumCP;
SymCP = zeros(NumAddPrefix,NumLoop);
RowPrefix = (NumSubc - NumCP + 1):NumSubc;
SymCP = [SymIFFT(RowPrefix,:);SymIFFT];%(把保護(hù)間隔放在符號的前面)
% --------------------------------------------- %
% 限幅 %
% --------------------------------------------- %
CR = 0.8;
SymCP_Clip_last = zeros(NumAddPrefix,NumLoop);
SymCP_Clip_temp3 = zeros(NumAddPrefix,1);
%SymCP_Clip_temp2 = zeros()
for nSymbol = 1:NumLoop
Signal_Power = abs(SymCP(:,nSymbol).^2);
Peak_Power = max(Signal_Power); %Peak_Power = max(Signal_Power,[],2);
Mean_Power = mean(Signal_Power);%Mean_Power = mean(Signal_Power,2);
A = CR*sqrt(Mean_Power);
SymCP_Clip_temp1 = SymCP(:,nSymbol)';%把第nSymbol列的值轉(zhuǎn)換成行向量并付給一個(gè)向量
SymCP_Clip_temp2 = SymCP_Clip_temp1(Signal_Power>A);
SymCP_Clip = A*SymCP_Clip_temp2./abs(SymCP_Clip_temp2);
SymCP_Clip_temp1(Signal_Power>A) = SymCP_Clip;
SymCP_Clip_temp3 = SymCP_Clip_temp1'
SymCP_Clip_last(:,nSymbol) = SymCP_Clip_temp3;
end
% --------------------------------------------- %
% Go through the channel %
% --------------------------------------------- %
% input: SymCP(NumSubc + NumCP,NumLoop); output: SymCh(1,(NumSubc + NumCP)*NumLoop)
SymChtmp = zeros(1,(NumSubc + NumCP)*NumLoop);
SymChtmp = SymCP_Clip_last(:).'; % 進(jìn)行這個(gè)轉(zhuǎn)置操作之后就成了一個(gè)矢量; 并串變換
% 相當(dāng)于把矩陣的列向量依次排列 改變?yōu)橐粋€(gè)行向量
% % --------------------------------------------- %
% %濾波
% % --------------------------------------------- %
% Ch = [1 1/2 1/4];
% SymChtmptmp = filter(Ch,1,SymChtmp);%濾波
snrIndex3 = 0;
BerSnrTable_CR0P8 = zeros(7,2);%?????
for snr=0:2:13; % = SNR + 10*log10(log2(2));
snrIndex3=snrIndex3+1;
BerSnrTable_CR0P8(snrIndex3,1) = snr;
SymCh = awgn(SymChtmp,snr,'measured');
% --------------------------------------------- %
% Remove Guard Intervals %
% --------------------------------------------- %
% input: SymCh(1,(NumSubc + NumCP)*NumLoop); output: SymDeCP(NumSubc,NumLoop)
SymDeCP = zeros(NumSubc,NumLoop);
SymDeCPtmp = reshape(SymCh,NumSubc + NumCP,NumLoop);
SymDeCP = SymDeCPtmp((NumCP+1+SyncDelay):NumAddPrefix+SyncDelay,:);
% --------------------------------------------- %
% FFT %
% --------------------------------------------- %
% input: SymDeCP(NumSubc,NumLoop); output: SymFFT(NumSubc,NumLoop)
SymFFT = fft(SymDeCP,NumSubc,1);
% --------------------------------------------- %
% Make Decision(Include DeQAM) %
% --------------------------------------------- %
% SymFFT(NumSubc,NumLoop); output: SymDec(NumSubc,NumLoop)
SymDec = zeros(NumSubc,NumLoop);
SymEqtmp(1,:) = SymFFT(1,:)+i*SymFFT(NumSubc/2+1,:);
SymEqtmp(2:NumSubc/2,:) = SymFFT(2:NumSubc/2,:);%???????????(SymEqtmp是哪里來的？)
for m = 1:NumLoop
for n = 1:NumSubc/2
Real = real(SymEqtmp(n,m));
Imag = imag(SymEqtmp(n,m));
if( abs((Real -1)) < abs((Real +1)))%????????每個(gè)點(diǎn)到坐標(biāo)軸上1和-1的距離，哪個(gè)離1近就把1付給它，同理得-1
SymDec(2*n-1,m) = 1;
else
SymDec(2*n-1,m) = 0;%把實(shí)部放在奇數(shù)行
end
if( abs((Imag -1)) < abs((Imag +1 )) )
SymDec(2*n,m) = 1;
else
SymDec(2*n,m) = 0;%把虛部放在偶數(shù)行
end
end
end
%----------------------------------------------------------------------
% Bit Error
%----------------------------------------------------------------------
BitsRx = zeros(1,NumSubc*NumLoop);
BitsRx = SymDec(:).';
[Num,Ber] = symerr(BitsTx,BitsRx)%(BitsTx是原始數(shù)據(jù)符號)
%BerSnrTable(snr+1,2) = Num ;
BerSnrTable_CR0P8(snrIndex3,2) = Ber ;%BerSnrTable(snr+1,3) = Ber ;
end
BerSnrTable_CP0P8_plot = zeros(7,1);
BerSnrTable_CP0P8_plot = BerSnrTable_CR0P8(:,2);
save('CR_0P8','BerSnrTable_CP0P8_plot');
semilogy(BerSnrTable_CR0P8(:,1),BerSnrTable_CP0P8_plot,'*-g');
legend('CR=2');
% xlabel('SNR');
% ylabel('BER');
grid on;
hold on;
% ------------------------------------------------------------------------
time_of_sim = toc
%end;
%
%
%
%
%
%CR=0.6
% --------------------------------------------- %
% 參數(shù)定義 %
% --------------------------------------------- %
% Initialize the parameters
% NumLoop = 10000;
% NumSubc = 64;
% NumCP = 8;
% SyncDelay = 0;
%-----------------------------------------------
%-----------------------------------------------
% Generate the random binary stream for transmit test
%-----------------------------------------------
BitsTx = floor(rand(1,NumLoop*NumSubc)*2);%用rand產(chǎn)生0：1的隨機(jī)數(shù)，不會產(chǎn)生負(fù)值，如果用randn會產(chǎn)生負(fù)值；乘以2是為了使用floor函數(shù)
% --------------------------------------------- %
% QAM MODULATION %
% --------------------------------------------- %
% Modulate (Generates QAM symbols)
% input: BitsTx(1,NumLoop*NumSubc); output:SymQAM(NumLoop,NumSubc/2)
SymQAMtmp = reshape(BitsTx,2,NumLoop*NumSubc/2).';
SymQAMtmptmp = bi2de(SymQAMtmp,2,'left-msb');%/?????(bi2de)
%--------------------------------------------------------------------------
QAMTable = [-1-i -1+i 1-i 1+i];
SymQAM = QAMTable(SymQAMtmptmp+1);
% --------------------------------------------- %
% IFFT input: SymQAM(NumLoop,NumSubc/2); output: SymIFFT(NumSubc,NumLoop)
% --------------------------------------------- %
SymIFFT = zeros(NumSubc,NumLoop);
SymIFFTtmp = reshape(SymQAM,NumSubc/2,NumLoop);%（為什么載波數(shù)除以2）
SymIFFTtmptmp = zeros(NumSubc,NumLoop);
SymIFFTtmptmp(1,:) = real(SymIFFTtmp(1,:)); % 實(shí)數(shù)
SymIFFTtmptmp(NumSubc/2+1,:) = imag(SymIFFTtmp(1,:)); % 實(shí)數(shù)
SymIFFTtmptmp(2:NumSubc/2,:) = SymIFFTtmp(2:NumSubc/2,:);
SymIFFTtmptmp((NumSubc/2+2):NumSubc,:) = flipdim(conj(SymIFFTtmp(2:NumSubc/2,:)),1);
SymIFFT = ifft(SymIFFTtmp,NumSubc,1);%(后邊這個(gè)1是干什么的？)
% --------------------------------------------- %
% Add cyclic prefix %
% --------------------------------------------- %
% input: SymIFFT(NumSubc,NumLoop); output: SymCP(NumSubc + NumCP,NumLoop)
NumAddPrefix = NumSubc + NumCP;
SymCP = zeros(NumAddPrefix,NumLoop);
RowPrefix = (NumSubc - NumCP + 1):NumSubc;
SymCP = [SymIFFT(RowPrefix,:);SymIFFT];%(把保護(hù)間隔放在符號的前面)
% --------------------------------------------- %
% 限幅 %
% --------------------------------------------- %
CR = 0.6;
SymCP_Clip_last = zeros(NumAddPrefix,NumLoop);
SymCP_Clip_temp3 = zeros(NumAddPrefix,1);
%SymCP_Clip_temp2 = zeros()
for nSymbol = 1:NumLoop
Signal_Power = abs(SymCP(:,nSymbol).^2);
Peak_Power = max(Signal_Power); %Peak_Power = max(Signal_Power,[],2);
Mean_Power = mean(Signal_Power);%Mean_Power = mean(Signal_Power,2);
A = CR*sqrt(Mean_Power);
SymCP_Clip_temp1 = SymCP(:,nSymbol)';%把第nSymbol列的值轉(zhuǎn)換成行向量并付給一個(gè)向量
SymCP_Clip_temp2 = SymCP_Clip_temp1(Signal_Power>A);
SymCP_Clip = A*SymCP_Clip_temp2./abs(SymCP_Clip_temp2);
SymCP_Clip_temp1(Signal_Power>A) = SymCP_Clip;
SymCP_Clip_temp3 = SymCP_Clip_temp1'
SymCP_Clip_last(:,nSymbol) = SymCP_Clip_temp3;
end
% --------------------------------------------- %
% Go through the channel %
% --------------------------------------------- %
% input: SymCP(NumSubc + NumCP,NumLoop); output: SymCh(1,(NumSubc + NumCP)*NumLoop)
SymChtmp = zeros(1,(NumSubc + NumCP)*NumLoop);
SymChtmp = SymCP_Clip_last(:).'; % 進(jìn)行這個(gè)轉(zhuǎn)置操作之后就成了一個(gè)矢量; 并串變換
% 相當(dāng)于把矩陣的列向量依次排列 改變?yōu)橐粋€(gè)行向量
% % --------------------------------------------- %
% %濾波
% % --------------------------------------------- %
% Ch = [1 1/2 1/4];
% SymChtmptmp = filter(Ch,1,SymChtmp);%濾波
snrIndex2 = 0;
BerSnrTable_CR0P6 = zeros(7,2);%?????
for snr=0:2:13; % = SNR + 10*log10(log2(2));
snrIndex2=snrIndex2+1;
BerSnrTable_CR0P6(snrIndex2,1) = snr;
SymCh = awgn(SymChtmp,snr,'measured');
% --------------------------------------------- %
% Remove Guard Intervals %
% --------------------------------------------- %
% input: SymCh(1,(NumSubc + NumCP)*NumLoop); output: SymDeCP(NumSubc,NumLoop)
SymDeCP = zeros(NumSubc,NumLoop);
SymDeCPtmp = reshape(SymCh,NumSubc + NumCP,NumLoop);
SymDeCP = SymDeCPtmp((NumCP+1+SyncDelay):NumAddPrefix+SyncDelay,:);
% --------------------------------------------- %
% FFT %
% --------------------------------------------- %
% input: SymDeCP(NumSubc,NumLoop); output: SymFFT(NumSubc,NumLoop)
SymFFT = fft(SymDeCP,NumSubc,1);
% --------------------------------------------- %
% Make Decision(Include DeQAM) %
% --------------------------------------------- %
% SymFFT(NumSubc,NumLoop); output: SymDec(NumSubc,NumLoop)
SymDec = zeros(NumSubc,NumLoop);
SymEqtmp(1,:) = SymFFT(1,:)+i*SymFFT(NumSubc/2+1,:);
SymEqtmp(2:NumSubc/2,:) = SymFFT(2:NumSubc/2,:);%???????????(SymEqtmp是哪里來的？)
for m = 1:NumLoop
for n = 1:NumSubc/2
Real = real(SymEqtmp(n,m));
Imag = imag(SymEqtmp(n,m));
if( abs((Real -1)) < abs((Real +1)))%????????每個(gè)點(diǎn)到坐標(biāo)軸上1和-1的距離，哪個(gè)離1近就把1付給它，同理得-1
SymDec(2*n-1,m) = 1;
else
SymDec(2*n-1,m) = 0;%把實(shí)部放在奇數(shù)行
end
if( abs((Imag -1)) < abs((Imag +1 )) )
SymDec(2*n,m) = 1;
else
SymDec(2*n,m) = 0;%把虛部放在偶數(shù)行
end
end
end
%----------------------------------------------------------------------
% Bit Error
%----------------------------------------------------------------------
BitsRx = zeros(1,NumSubc*NumLoop);
BitsRx = SymDec(:).';
[Num,Ber] = symerr(BitsTx,BitsRx)%(BitsTx是原始數(shù)據(jù)符號)
%BerSnrTable(snr+1,2) = Num ;
BerSnrTable_CR0P6(snrIndex2,2) = Ber ;%BerSnrTable(snr+1,3) = Ber ;
end
BerSnrTable_CR0P6_plot = zeros(7,1);
BerSnrTable_CR0P6_plot = BerSnrTable_CR0P6(:,2);
save('CR_0P6','BerSnrTable_CR0P6_plot');
semilogy(BerSnrTable_CR0P6(:,1),BerSnrTable_CR0P6(:,2),'+:m');
legend('CR=1');
% xlabel('SNR');
% ylabel('BER');
grid on;
hold on;
% ------------------------------------------------------------------------
time_of_sim = toc
% end;
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -