【科研基础】matlab绘制PARR

作者在使用matlab编程时，常忘记一些函数的使用规则，因此在这里记录，先起二级标题，如果能够记录达到一定规模，就归纳一级标题。

matlab官方文档

配色

fillcolor1=[0.85, 0.33, 0.10]、fillcolor2=[0.93, 0.69, 0.13]、fillcolor3=[0.00, 0.45, 0.74]

绘制PARR

关于如何求PAPR的CCDF以及ecdf函数的使用

% Code to simuulate PAPR for OFDMA
clc
clear

M = 16; % M-ary QAM
totalSubcarriers = 512; % Number of total subcarriers.
numSymbols = 32; % Input Data block size.
Q = totalSubcarriers/numSymbols; % Bandwidth spreading factor of IFDMA.
Q_tilda = 15; % Bandwidth spreading factor of DFDMA. Q_tilda < Q.
numRuns = 1e5; % Number of runs.
papr = zeros(1,numRuns); % Initialize the PAPR results. (LFDMA)
papr1 = zeros(1,numRuns); % Initialize the PAPR results. (IFDMA)
papr2 = zeros(1,numRuns); % Initialize the PAPR results. (DFDMA)
for n = 1:numRuns
    % Generate random data.
    inputData = randi([0 M-1], 1, numSymbols);
    inputSymbols = qammod(inputData, M, 'UnitAveragePower', false);
    % Initialize the subcarriers.
    Y = zeros(totalSubcarriers,1); %For LFDMA
    Y1 = zeros(totalSubcarriers,1); %For IFDMA
    Y2 = zeros(totalSubcarriers,1); %For DFDMA
    % Subcarrier mapping. LFDMA
    Y(1:numSymbols) = inputSymbols;
    % Subcarrier mapping. IFDMA
    Y1(1:Q:totalSubcarriers) = inputSymbols;
    % Subcarrier mapping. DFDMA
    Y2(1:Q_tilda:Q_tilda*numSymbols) = inputSymbols;
    % IDFT
    y = ifft(Y);
    y1 = ifft(Y1);
    y2 = ifft(Y2);
    % Calculate PAPR.
    papr(n) = 10*log10(max(abs(y).^2)/mean(abs(y).^2));
    papr1(n) = 10*log10(max(abs(y1).^2)/mean(abs(y1).^2));
    papr2(n) = 10*log10(max(abs(y2).^2)/mean(abs(y2).^2));
end

% Plot CCDF.
%X是横轴
%N是每个bin的值
[N,X] = histcounts(papr, 200,'Normalization','count');
[N1,X1] = histcounts(papr1, 200,'Normalization','count');
[N2,X2] = histcounts(papr2, 200,'Normalization','count');
figure
 semilogy(X(1:end-1),1-cumsum(N)/max(cumsum(N)),'-r') % For LFDMA
 hold on
 semilogy(X1(1:end-1),1-cumsum(N1)/max(cumsum(N1)),'-g') % For IFDMA
 hold on
 semilogy(X2(1:end-1),1-cumsum(N2)/max(cumsum(N2)),'-b') % For DFDMA
 hold on
 legend('LFDMA','IFDMA','DFDMA')
 grid on
 xlabel('PAPR')
 ylabel('CCDF')
 title('OFDMA PAPR')

纵坐标表示的是，$\rho$>PAPR（PAPR0是横坐标）的概率,
$$\rho =\frac{\max {\left|x_{\mathrm{RF}}(t)\right|}^2}{\mathbb{E}\left[{\left|x_{\mathrm{RF}}(t)\right|}^2\right]}$$

histcounts

matlab计算概率密度函数pdf，使用ksdensity histcounts区别
用于计算概率和概率密度的相关函数

clear;
clc;
close all
%X = randn(n) 返回由正态分布的随机数组成的 n×n 矩阵
x = randn(100000,1);%生成服从正态分布的数
width = 0.2;
xi = -4:width:4;
edge = [xi-0.5*width,xi(end)+0.5*width];%   -4-0.5*width:width:4+0.5*width;
%直方图的根数（nbin）
%直方图的宽度设置（edge，即每根的左右边界）），得到的值与nbin一致，比edge少一个
[dens,edge,bin] = histcounts(x,edge,'normalization','pdf');%可以计数，计算概率密度，概率密度函数，累积概率密度函数等，通过‘normalization'来选择
[yi,xi]=ksdensity(x,xi);
figure
hold on 
%把x的数值用edge个直方图表示
histogram(x,edge,'normalization','pdf')
plot(xi,dens,'o-')
plot(xi,yi,'*-')
legend('histogram','histcounts','ksdensity')
 
id = ismember(bin,1); %bin矩阵中是否包含1，返回0或1
x(id);
id = ismember(bin,20);
a = categorical(bin);
summary(a)

semilogy

clc
clear
x = 0:0.1:10;
y = exp(x);
figure
semilogy(x,10.^x) %如何理解对数坐标：我们可以保持对数值等距离摆放，这就是对数坐标系。