%%%?采用根升余弦濾波器

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