2DPSK_原创精品文档.pdfVIP

clear;clc;

closeall;

%系统仿真参数

A=1;%载波振幅

fc=2;%载波频率（Hz）

snr=15;%信噪比dB

N_Sample=8;%基带信号中每个码元的的采样点数

N=10000;%码元数

Ts=1;%码元宽度

df=0.01;%频率分辨率

B=1/Ts;%传码率

f_start=fc-B;%带通信号的下截止频率

f_cutoff=fc+B;%带通信号的上截止频率

fs=fc*N_Sample;%系统采样频率，即考虑载波后一个码元内的采样点数

ts=Ts/fs;%系统采样间隔

t=0:ts:N*Ts-ts;%时间序列

Lt=length(t);%采样点数

%画出调制信号波形及功率谱

%产生二进制信源

d=sign(randn(1,N));%产生双极性二进制序列

s=(d+1)/2;%转化为单极性

%差分编码

b(1)=abs(1-s(1));%1为参考

fori=2:N%差分编码

b(i)=abs(s(i)-b(i-1));

end

dd=sigexpand(s,fc*N_Sample);

gt=ones(1,fc*N_Sample);

d_NRZ=conv(dd,gt);

d_NRZ1=d_NRZ(1:Lt);

d_sjx=2*d_NRZ-1;

d_sjx1=d_sjx(1:Lt);

figure(1);

subplot(321);

plot(t,d_sjx1);

axis([min(t),10,min(d_sjx1)-0.2,max(d_sjx)+0.2]);

xlabel(t);title(原始双极性信号);

subplot(322);

[d_sjx1f,d_sjx1,df1,f]=T2F_new(d_sjx1,ts,df,fs);

plot(f,10*log10(abs(fftshift(d_sjx1f).^2/length(f))));

axis([-5*B,5*B,-50,0]);

xlabel(f);title(原始双极性信号PDF);

dd=sigexpand(b,fc*N_Sample);

gt=ones(1,fc*N_Sample);

d_NRZ=conv(dd,gt);

d_NRZ1=d_NRZ(1:Lt);

d_sjx=2*d_NRZ-1;

d_sjx1=d_sjx(1:Lt);

subplot(323);

plot(t,d_sjx1);

axis([min(t),10,min(d_sjx1)-0.2,max(d_sjx1)+0.2]);

xlabel(t);title(差分编码后的双极性信号);

subplot(324);

[d_sjx1f,d_sjx1,df1,f]=T2F_new(d_sjx1,ts,df,fs);

plot(f,10*log10(abs(fftshift(d_sjx1f).^2/length(f))));

axis([-5*B,5*B,-50,0]);

xlabel(f);title(双极性信号PDF);

ht=A*sin(2*pi*fc*t);%载波

%对数字基带信号进行2DPSK调制

s_2dpsk=d_sjx(1:Lt).*ht;%生成2DPSK信号

subplot(325);

plot(t,s_2dpsk);

axis([min(t),10,min(s_2dpsk)-0.2,max(s_2dpsk)+0.2]);

xlabel(t);title(2DPSK信号);

subplot(326);

[s_2dpskf,s_2dpsk,df1,f]=T2F_new(s_2dpsk,ts,df,fs);

plot(f,10*log10(abs(fftshift(s_2dpskf).^2/length(f))));

axis([-fc-3*B