Clc;clear
%产生比特信号
t=0:0.01:9.99;
a=randint(1,10);
m=a(ceil(t+0.01));
figure(1)
plot(t,m)
title('产生随机十比特二进制比特序列');
axis([0,10,-2,2]);
%载波信号%
fc=4800;
fs=100000;
ts=0:1/fs:(1000-1)/fs;
carry=cos(2*pi*fc*ts/2.5);
figure(2)
plot(t,carry)
title('载波信号')
axis([0,10,-2,2]);
%2psk 调制%
e_2psk=cos(2*pi*fc*ts/2.5+pi*m);
figure(3)
plot(t,e_2psk)
title('2psk 调制信号');
axis([0,10,-2,2]);
%产生高斯白噪声
a=0.1;
noise=a*wgn(1,1000,10);
%产生 1 行 1000 列的高斯白噪声为 10dBW
figure(4)
plot(t,noise)
title('高斯白噪声');
axis([0,10,-2,2]);
%调制信号叠加高斯白噪声
e2psk=e_2psk+noise;
figure(5)
plot(t,e2psk)
title('2psk 叠加白噪声波形')
axis([0,10,-2,2]);