2psk调制与解调的代码
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
max=10;
s=rand(1,max);%长度为max的随机二进制序列
s=round(s);
f=pi/4;
t = 0:2*pi/99:2*pi;
m1 = [];
c1 = [];
b1 = [];
for n = 1:length(s)
if s(n) == 0;
m = -ones(1,100);
b = zeros(1,100)
else s(n) == 1;
m = ones(1,100);
b = ones(1,100)
end
c = sin(f*t);
m1 = [m1 m];
c1 = [c1 c]
b1 = [b1 b];
end
psk = c1.*m1;
plot(b1);
title(['信号源:',num2str(s)])
axis([0 100*length(s) -0.2 1.1]);
% grid on;
jt=c1.*psk;
bn=fir1(101,[2/800 10/800]); %经过低通滤波器
jm=filter(bn,1,jt);
for m=1:100*length(s);
if jm(m)>-0.001;
jn(m)=1;
else jm(m)<-0.001;
jn(m)=0;
end
end
plot(psk)
title('2PSK信号');
plot(jt)
title ('乘本地载波后');
plot(jm)
title ('经低通滤波后');
plot(jn);
title ('PSK已调信号');
axis([0 100*10 -0.2 1.1]);