# MATLAB代做-基于matlab的PTS算法的仿真

### 时间：2018-12-15 0:29:56 点击：

核心提示：MATLAB代做-基于matlab的PTS算法的仿真...

% Original PTS--PAPR Simulation

%
% Reference: S. H. Müller and J. B. Huber, “OFDM with Reduced
% Peak–to–Average Power Ratio by Optimum Combination
% of Partial Transmit Sequences,” Elect. Lett., vol. 33,
% no. 5, Feb. 1997, pp. 368–69.%
% Author: Bin Jiang
% National Mobile Communication Research Laboratory
% Southeast University, 210096 Nanjing, China
% Email: bjiang@seu.edu.cn

clear all; clc; close all;

K = 128;                                                                    % SIZE OF FFT
V = 4;                                                                      % NUMBER OF SELECTIONS
QPSK_Set  = [1 -1 j -j];
Phase_Set = [1 -1];                                                         % W = 2
% 产生的序列个数 = W ^(V-1) = 8
% 针对QPSK共有16种相位变化组合
Choose = [1 1 1 1; 1 1 1 2; 1 1 2 1; 1 2 1 1; 2 1 1 1;...
1 1 2 2; 1 2 1 2; 1 2 2 1; 2 2 1 1; 2 1 2 1; 2 1 1 2;...
2 2 2 1; 2 2 1 2; 2 1 2 2; 1 2 2 2; 2 2 2 2];
Choose_Len = 16;

MAX_SYMBOLS  = 1e5;
PAPR_Orignal = zeros(1,MAX_SYMBOLS);
PAPR_PTS     = zeros(1,MAX_SYMBOLS);

for nSymbol=1:MAX_SYMBOLS
nSymbol
Index = randint(1,K,length(QPSK_Set))+1;
X = QPSK_Set(Index(1,:));                                               % Orignal Frequency domain signal
x = ifft(X,[],2);                                                       % Time domain signal
Signal_Power = abs(x.^2);
Peak_Power   = max(Signal_Power,[],2);
Mean_Power   = mean(Signal_Power,2);
PAPR_Orignal(nSymbol) = 10*log10(Peak_Power./Mean_Power);

% PTS (交织分割/相邻分割)
A = zeros(V,K);
for v=1:V
%        A(v,v:V:K) = X(v:V:K);
A(v,(1+(v-1)*K/4):1:v*K/4) = X((1+(v-1)*K/4):1:v*K/4);
end
a = ifft(A,[],2);
% 门限()
min_value = 10;
% 此处为遍历寻找最优辅助信息，可以用迭代法寻求次优辅助信息，以减小计算量，以V=4为例，可以由16次减少到V次迭代；参见Paper
for n=1:Choose_Len
temp_phase = Phase_Set(Choose(n,:)).';
%        temp_max = max(abs(sum(a.*repmat(temp_phase,1,K))));
b = repmat(temp_phase,1,K);
c = a.*b;
d = sum(c);
e = abs(d);
temp_max = max(e);
if temp_max<min_value
min_value = temp_max;
Best_n = n;
end
end
aa = sum(a.*repmat(Phase_Set(Choose(Best_n,:)).',1,K));

Signal_Power = abs(aa.^2);
Peak_Power   = max(Signal_Power,[],2);
Mean_Power   = mean(Signal_Power,2);
PAPR_PTS(nSymbol) = 10*log10(Peak_Power./Mean_Power);

end

[cdf1, PAPR1] = ecdf(PAPR_Orignal);
[cdf2, PAPR2] = ecdf(PAPR_PTS);

%--------------------------------------------------------------------------
semilogy(PAPR1(1:40:end),1-cdf1(1:40:end),'-bo',PAPR2(1:40:end),1-cdf2(1:40:end),'-rx')
legend('Orignal','O-PTS')
title('V=4')
xlabel('PAPR0 [dB]');
ylabel('CCDF (Pr[PAPR>PAPR0])');
grid on

QQ ：1224848052

Tags:PTS_PAPR

• 百度搜索
• 查阅资料过程中
• 论坛发现
• 百度贴吧发现
• 朋友介绍
• 没有相关文章
• 大名：
• 内容：