内容发布更新时间 : 2025/9/19 19:13:56星期一 下面是文章的全部内容请认真阅读。
axis([0 N -2 2]);xlabel('经低通滤波器后信号波形') %抽样判决
x=fun_panjue(x);%调用函数，进行抽样判决 figure(4);subplot(2,1,1);plot(x);grid on; xlabel('加噪后解调信号x(t)'); axis([0 N -2 2]); %消除延迟
x=fun_yanc(x); %调用函数，进行消除延迟 figure(4);subplot(2,1,2);plot(x);grid on; xlabel('加噪后去掉延迟的解调信号x(t)'); axis([0 N -2 2]); %误码率计算
Err1=length(find(x~=s)) %计算解调信号中错误码元个数 Pe_test1=Err1/N %计算实际误码率
Pe1=(1/2)*erfc(sqrt(snr)) %计算系统理论误码率 % 理论误码率曲线 Pe=[]; for SNR=1:10
am=0.7; %输入信号经信道后振幅由1衰减为0.7 E=am*am/2; snr=10^(SNR/10); N0=(am*am)/(2*snr);
no=N0/(2*200); %计算噪声功率
N0_db=10*log10(N0);%将噪声功率转换为dBW ni=wgn(1,N,N0_db);% 产生1行N列的高斯噪声
yi=e+ni; %BSK已调信号中加入白噪声，输入信噪比为SNR y=filter(b1,a1,yi);%对yi进行滤波(带通滤波器)，得到信号y x1=2*c.*y; %与恢复载波相乘 xx=filter(b2,a2,x1); %经低通滤波器滤波 xx=fun_panjue(xx);%抽样判决 xx=fun_yanc(xx); %消除延迟
20
snr=10^(SNR/10);
Pe=[Pe,(1/2)*erfc(sqrt(snr))]; %计算理论误码率 end Pe; figure; SNR=1:10;
semilogy(SNR,Pe,'b--');hold on %以log10(Pe)为纵坐标画图 grid on
子程序一：抽样判决
（由于存在“倒pi”现象，故以0为基准进行判决，大于0，判为0；小于0，判为1） function w=fun_panjue(w) N=length(w); if w(100)>0
w(1:100)=0; else
w(1:100)=1; end
for i=101:N if w(i)>0; w(i)=0; else
w(i)=1; end end
子函数二：消除延迟
（由于抽样判决后会产生延迟，故认为对解调信号进行时移，以减小甚至消除延迟，便于与基带信号进行比较，计算实际误码率） function m=fun_yanc(m) N=length(m); leng=0; if m(1)==0 for i=1:N
if m(i)==1 leng=i; break; end end else
for i=1:N
if m(i)==0 leng=i;
21
break; end end end
leng1=leng-(floor(leng/200))*200; for i=1:(N-leng1)
m(i)=m(i+leng1); end
for i=(N-leng1):N
m(i)=m(N-200+10); end
22