ratio.m

OFDM信号的循环谱密度

%%%%% 循环谱检测OFDM信号 %%%%% 
clear;
clc;

%%% OFDM参数 %%%
%数据长度的1/8必须大于等于循环谱采样长度
%数据速率 6-BPSK, 12-QPSK
%%%%%%%%%%%%%%%%
TXVECTOR.LENGTH = 2000; % 数据长度
TXVECTOR.DATARATE = 6; % 数据速率
trst_rate = 20e6; % 信号发射速率,恒定
fc = 100e6;

%%% 循环谱检测参数 %%%
%检测带宽为 -fs/2 至 fs/2
%循环频率分辨率为 fs/N
%频率分辨率为 M*fs/N
%%%%%%%%%%%%%%%%%%%%%
fs = 300e6; % 采样频率
N = 2048; % 采样长度
M = 30; % 平滑点数

SNR = 7:0.5:15;
me = 100;
RT = zeros(length(SNR),me);
i = 1;
%%% 信道参数 %%%
for snr = SNR % 信噪比
    for k = 1:me
%%% 随机数据生成 %%%
PSDU = round(rand(1,8*TXVECTOR.LENGTH));

%%% OFDM信号生成 %%%
sig = transmitter(PSDU,TXVECTOR);

%%% 提升信号采样率 %%%
s_n = ceil(fs/trst_rate); % 检测采样率近似为OFDM信号速率的整数倍
sig = sig(ones(s_n,1),:); % 码元复制
sig = reshape(sig, 1, s_n*length(sig));

%%% 载波调制 %%%
sig_chnl = real(sig.*exp(j*2*pi*fc/fs*(0:length(sig)-1)));

%%% 高斯信道 %%%
sig_awgn = awgn(sig_chnl, snr);
%sig_awgn=sig_chnl;
%sig_awgn = awgn(zeros(1,N), snr);

%%% 循环谱检测 %%%
%cyclic_spectrum(sig_awgn, N, fs, M);

RT(i,k) = cs_threshold(sig_awgn, N, fs, M, fc, trst_rate);
    end   
i = i+1;    
end
rt = mean(RT,2)';
%normplot(RT(:));
save('204830', 'SNR', 'rt');
%plot(SNR, rt);