协作图clc
clear all
num=10;%码元个数
ek=round(rand(1,num,1)); %产生载波
t=0:2*pi/999:2*pi; %周期
num1=1000; %个数
len=length(t);
c=-sin(t); %调制
ektemp=[];
ctemp=[];
for n=1:num
if ek(n)==0
ctemp=[ctemp,-c];
a=zeros(1,len);
ektemp=[ektemp,a];
elseif ek(n)==1
ctemp=[ctemp,c];
a=ones(1,len);
ektemp=[ektemp,a];
end
end
subplot(2,1,1); %分别画出原信号、已调信号示意 plot(ektemp,'LineWidth',1.5);
title('原码元');
grid on;
axis([0 num1*num -2.5 2.5]);
subplot(2,1,2);
plot(ctemp,'LineWidth',1.5);
大渡河造林局
title('已调信号');
grid on;
axis([0, num1*num, -2.5, 2.5]);%相干解调 %本地载波
loca=[];
for n=1:num
loca=[loca,c];%此处c的正负决定解调后是否倒pi
end
locb=loca;
%相乘
locc=2*ctemp.*locb;
%低通滤波器
[b,a] = butter(3,2*pi*0.0003,'low'); %信号频率为0.001
locd=filter(b,a,locc);%低通滤波后的信号
figure;
subplot(2,1,1); %相乘后的波形
plot(locc,'LineWidth',1.5);
title('相乘后的信号');
grid on;
axis([0 num1*num -2.5 2.5 ]);
dd=locd+randn(1,10000);
%低通滤波后的信号波形图
subplot(2,1,2);
plot(dd,'LineWidth',1.5);
title('低通滤波后的信号');
grid on;
axis([0 num1*num -2.5 2.5 ]);
%抽样判决
for i=1:num1*num
if(locd(i)>0) %判决,得到解调结果
locd(i)=1;
else
locd(i)=-1;
end
end
figure;
subplot(2,1,1); %分别画出原信号、解调后信号示意
plot(ektemp,'LineWidth',1.5);
title('原码元')
grid on;
axis([0 num1*num -2.5 2.5]);
聚氧化乙烯subplot(2,1,2);
plot(locd,'LineWidth',1.5);
title('解调后的信号');
grid on;
axis([0 num1*num -2.5 2.5 ]);
基于单片机的信号发生器
p=symerr(ektemp,locd)/10 %误码率 %误码率曲线
figure(4)
嘌呤代谢紊乱
r=0:2:10;
湖南农业大学综合信息服务系统
rr=10.^(r/10);
pe1=1/2*exp(-rr);%相干解调的误码率曲线
hold on
plot(r,pe1,'r');grid on;
pe2=(1-1/2*erfc(sqrt(rr))).*erfc(sqrt(rr));%差分相干解调的误码率曲线
plot(r,pe2,'b');xlabel('bpsk,dpsk误码率曲线');
set(gca,'XTick',-6:3:18);