%*******************modulation****************
D=[1 0 1 1 0 0 1 0]; %产生随机数
N=length(D); %基带信号符号数目
t=0:0.01:7.99;
y=cos(2*pi*t); %载波
d0=D;
d1=2*D-1; %单极性信号变双极性 1 -1 1 1 -1 -1 1 -1
x=ones(1,100); %ones(1,100)返回一个值全为1的1*100维矩阵
D0=[]; %空矩阵
for k=1:N
D0=[D0,d1(k)*x];%在矩阵D0中增加
end %开关信号,为便于观察
z=y.*D0; %bpsk的表达式
%*****************demodulation***************
SNR=30; % 信噪比,以dB为单位,一般不应该低于70dB,高保真音箱的信噪比应达到110dB以上
z=awgn(z,SNR,'measured'); %在信号z中加入高斯白噪声,信噪比为SNR,‘measured’表示函数将在加入噪声之前测定信号强度
recos=z.*y;
meg=find(recos>0); % find(recos>0):寻找recos中大于0的指标
ch(meg)=1; % >0的置1
meg=find(recos<0); % find(recos<0):寻找recos中小于0的指标
ch(meg)=-1; % <0的置0
demo=(ch+1)/2;
%***************drawing**********************
x=ones(1,100); %ones(1,100)返回一个值全为1的1*100维矩阵
D0=[]; %空矩阵
for k=1:N
D0=[D0,d0(k)*x];%在矩阵D0中增加
end %开关信号,为便于观察
subplot(3,1,1);plot(D0);grid ;title('original signal');axis([0,800,-2,2]);
subplot(3,1,2);plot(z);grid ;title('BPSK modulation');axis([0,800,-2,2]);
subplot(3,1,3);plot(demo);grid ;title('BPSK demodulation');axis([0,800,-2,2]);
% xlabel('1 0 1 1 0 0 1 0');
BPSK.rar_OFDMQPSK
版权申诉
102 浏览量
2022-09-23
10:25:04
上传
评论
收藏 815B RAR 举报
JaniceLu
- 粉丝: 83
- 资源: 1万+
评论0