matlab dpsk,Matlab之DPSK

% communication DPSK demodulation

% name:Chen Yu

% student ID:12353032

function Lab9

%

% 参数设置

%

Ts=1/100; %码元周期

Rs=1/Ts; %码元频率

n1=10; %载波频率系数

fc=Rs*n1; %载波频率

n2=64; %采样频率系数

fs=fc*n2; %采样频率，N_sample = 640

N_sample=fs*Ts; %每个码元的抽样点数=n1*n2,经测试小于128时，效果较差

dt=Ts/N_sample; %抽样时间间隔

N=10; %码元数

t=0:dt:((N+1)*N_sample-1)*dt;%传输序列持续时间,考虑了参考电平0

%

%基本码元g(t)=1，还有高斯白噪声

%

gt=cos(2*pi*fc*t(1:N_sample)+pi); %相位pi对应码元1

%噪声,均值=0,标准差=0.01

noise=normrnd(0,0.1,1,length(t));

s=zeros(1,(N+1)*N_sample);%采样点

absolute_code = randi([0,1],1,N); %生成绝对码

relative_code = zeros(1,N+1);

relative_code(1) = 0; %参考码元为0，为第一个码元

for k = 2:N+1

relative_code(k) = xor(absolute_code(k-1),relative_code(k-1));%生成相对码

end

%生成调制信号

s(1:N_sample) = -1*gt;

for i = 2:N+1

if relative_code(i) == 1 %得到波形图

s((i-1)*N_sample+1:i*N_sample) = gt; %码元为1

else

s((i-1)*N_sample+1:i*N_sample) = -1*gt; %码元为0

end

end

figure(1);

plot(t,s);

title('调制信号前两个码元');

axis([0,0.02,-1.5,1.5]);

grid on

%加入高斯白噪声

st = s + noise;

figure(2);

plot(t,st);

title('加入高斯白噪声后，调制信号前两个码元');

axis([0,0.02,-1.5,1.5]);

grid on

%带通滤波，去除AWGN

w1 = 2*(fc-Rs)/fs; %下截止频率

w2 = 2*(fc+Rs)/fs; %上截止频率

numerator= fir1(60,[w1,w2]); %FIR滤波器，numerator是FIR滤波器的系数

sfBandpass=filter(numerator,1,st); %st为时间信号

figure(3);

plot(t,sfBandpass);

title('通过带通滤波器后，调制信号前两个码元');

axis([0,0.02,-1.5,1.5]);

grid on

%进行相干解调

sfModulation = sfBandpass.*cos(2*pi*fc*t);

figure(4);

plot(t,sfModulation);

title('通过解调后，解调信号前两个码元');

axis([0,0.02,-1.5,1.5]);

grid on

%低通滤波

w1 = [];

w2 = 2*(Rs)/fs; %截止频率

numerator= fir1(60,[w1,w2]); %FIR滤波器，numerator是FIR滤波器的系数

sfLowpass=filter(numerator,1,sfModulation); %st为时间信号

figure(5);

plot(t,sfLowpass);

title('通过低通滤波器后的调制信号');

axis([0,0.1,-1.5,1.5]);

grid on

%抽样判决，我选取的是中间时刻

decide=zeros(1,N+1);

for k=1:N+1

if(sfLowpass((k-1)*N_sample+N_sample/2)>0)

decide(k)=1;

else

decide(k)=0;

end

end

%还原绝对码

get_absolute_code = zeros(1,N);

for j = 1:N

get_absolute_code(j) = xor(decide(j),decide(j+1));

end

%画出绝对码前后对比

n=1:N;

%输出原始码元

figure(6);

stem(n,absolute_code,'r','LineStyle','none');

%输出解调码元

hold on

stem(n,get_absolute_code,'*','b','LineStyle','none');

hleg=legend('原始码元','解调码元');

set(hleg,'Location','Northeast');