i=10;%基带信号码元数
j=5000;
t=linspace(0,10,j);%0-10之间产生5000个点行矢量，即将[0,5]分成5000份
f1=10;%载波1频率
f2=5;%载波2频率
fm=i/5;%基带信号频率 码元数是10，而时域长度是5，也就是一个单位2个码元
a=round(rand(1,i));%产生随机序列
%产生基带信号
st1=t;
for n=1:10
if a(n)<1
for m=j/i*(n-1)+1:j/i*n
st1(m)=0;
end
else
for m=j/i*(n-1)+1:j/i*n
st1(m)=1;
end
end
end
figure(1);
subplot(511);
plot(t,st1);
title('基带信号st1');
axis([0,10,-0.2,1.2]);
%基带信号求反
st2=t;
for n=1:j
if st1(n)==1
st2(n)=0;
else
st2(n)=1;
end
end
%载波信号
s1=cos(2*pi*f1*t);
s2=cos(2*pi*f2*t);
subplot(512),plot(s1);
title('载波信号s1');
subplot(513),plot(s2);
title('载波信号s2');
%调制
F1=st1.*s1;%加入载波1
F2=st2.*s2;%加入载波2
e_fsk=F1+F2;
subplot(514);
plot(t,e_fsk);
title('2FSK信号');%键控法产生的信号在相邻码元之间相位不一定连续
%加噪
nosie=rand(1,j);
fsk=e_fsk+nosie;
%相干解调
st1=fsk.*s1; %与载波1相乘
[f,sf1] = T2F(t,st1);%傅里叶变换
[t,st1] = lpf(f,sf1,2*fm);%通过低通滤波器
st2=fsk.*s2;%与载波2相乘
[f,sf2] = T2F(t,st2);%通过低通滤波器
[t,st2] = lpf(f,sf2,2*fm);
%抽样判决
for m=0:i-1
if st1(1,m*500+250)>st2(1,m*500+250)
for j=m*500+1:(m+1)*500
at(1,j)=1;
end
else
for j=m*500+1:(m+1)*500
at(1,j)=0;
end
end
end
subplot(515);
plot(t,at);
axis([0,10,-0.2,1.2]);
title('抽样判决后波形')
评论0