add_pbj_jamming.m

好东西

function [out,jam]=add_pbj_jamming(input_ofdm_symbol,jpow,partial_fraction,jfre_begin)
%在输入信号上加部分频带干扰，input_ofdm_symbol为待加扰信号，jpow为干扰总功率
%delta_f为子载波间隔，t_interval为采样时间间隔，partial_fraction为部分频带干扰占通信带宽的比例
%jfre_begin部分频带干扰起始频率,
 [N_carrier,N_ofdm]=size(input_ofdm_symbol);
 q=round(partial_fraction*N_carrier);%部分频带干扰所占子载波个数,N_carrier-N_cp子波数
jam=zeros(q,N_ofdm);
sgma=sqrt(jpow/2);
for i=1:q
    for j=1:N_ofdm
        %jam(i,j)=normrnd(0,sgma);
        jam(i,j)=normrnd(0,sgma)+normrnd(0,sgma)*sqrt(-1);
    end
end
jamming=zeros(N_carrier,N_ofdm);
jamming(jfre_begin:jfre_begin+q-1,:)=jam;
jam=jamming;
out=input_ofdm_symbol+jamming;