rcosflt_Jamming 1、根据干扰总功率要求

1、根据干扰总功率要求，在PSK调制信号带宽内产生单音干扰和多音干扰信号，并叠加到产生的信号源上。 2、将收到干扰的信号进行频率搬移到基带，然后通过匹配滤波器以及解调过程得到信息，统计误码率并与真实信息的解调误码率进行对比，评估音频干扰对PSK信号的干扰效果。

%%% 采用根升余弦滤波器

clear all;clc;close all;

r=4;%%% 过采样率
M=2; %% MPSK
N=100; %%% 符号个数
EbN0_v=[3:10]; 
a=0.35; %% 滚将滤波器系数
delay=3; %% 延迟
f=0.3; %% 载频
h=rcosfir（a delay r1‘sqrt‘）;%%% 采用根升余弦滤波器，使得匹配接收后为升余弦滤波器，从而避免ISI;
h=sqrt（r）*h/norm（h）; %%% 归一化滤波器系数
h_r=h（end:-1:1）; %%% 接收端匹配滤波器，由于线性相位对称，可不需要
K=1000; %% 仿真次数；
JSR=-10;%% 干信比
for m=1:length（EbN0_v）
    EbN0=EbN0_v（m）;
    EsN0=EbN0+10*log10（log2（M））;
    SNR=EsN0-10*log10（r）;
for k=1:K
    s=randsrc（N1[0:M-1]）;
    s_mod=pskmod（sM）;
    %scatterplot（s_mod）;
    s_base=rectpulse（s_modr）;
    s_base=zeros（rN）;
    s_base（1:）=s_mod.‘;
    s_base=s_base（:）;
    s_base=conv（h.‘s_base）;
    y_x=s_base.*exp（j*（2*pi*f*[0:length（s_base）-1].‘+rand*2*pi））;%Generate PSK modualted signal
    noise=sqrt（1/10^（SNR/10）/2）*（randn（size（s_base））+j*randn（size（s_base）））;%Generate noise signal
    %SNR10*log10（var（y_x）/var（noise））var（y_x）
    fj=f+（（1+a）*1/r）/8; %%% Frequency of jamming signal;
    Jam=sqrt（10^（JSR/10））*exp（j*（2*pi*fj*[0:length（s_base）-1].‘+rand*2*pi））;% Generate jamming signal
    %Jam=0;
    
    y=y_x+noise+Jam; %% 产生中频发射信号
   
    y_r=y.*exp（-j*（2*pi*f*[0:length（s_base）-1].‘））;  %%% 频谱搬移到基带
    s_r=conv（y_rh_r.‘）;  %% 匹配接收
    s_r=s_r（2*delay*r+1:r:2*delay*r+N*r）; %%% 抽取，无ISI
    
    s_est=pskdemod（s_r（2:end）./s_r（1:end-1）M）; %% 解调（点除消除相位模糊）
    s_true=pskdemod（s_mod（2:end）./s_mod（1:end-1）M）; 

    rate（mk）=length（find（s_est（1:length（s_true））~=s_true））/length（s_true）; %%% 误符号率
end
Rate（m）=mean（rate（m:））;

Rate_theo（m）=erfc（sqrt（2*（10^（EsN0/10）））*sin（pi/（2*M）））; %%% 理论误符号率
m
end

semilogy（EbN0_vRate）      %% 画图比较
hold onsemilogy（EbN0_vRate_theo‘r‘）
grid on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1856  2013-12-10 11:43  Communication_demo_zhg_rcosflt_Jamming.m