BPSK的matlab调制解调程序

clear %产生比特信号 
t=0:0.01:9.99; 
a=randint（110）;
snr=10;                %信噪比
m=a（ceil（t+0.01））;     %产生方波
figure（1）
plot（tm）  
title（‘产生随机十比特二进制比特序列‘）; 
axis（[010-22]）; %载波信号% 
fc=4800; 
fs=100000; 
ts=0:1/fs:（1000-1）/fs; 
carry=cos（2*pi*fc*ts/2.5）;  %因为信息速率2400b/s，载频 4800Hz，一个Ts两个载波信号，不除2.5是五个载波  
figure（2） 
plot（tcarry） 
title（‘载波信号‘）
axis（[010-22]）; %2psk调制% 
e_2psk=cos（2*pi*fc*ts/2.5+pi*m）;
figure（3）
plot（te_2psk） 
title（‘2psk调制信号‘）; 


ff= linspace（-fs/2 fs/2 length（e_2psk））;%频率坐标，单位Hz
figure（4）
plot（ff/1000fftshift（abs（fft（e_2psk））））;
title（‘未加噪声前的2psk信号的频谱‘）;
xlabel（‘频率（KHz）‘）;

%%
% axis（[010-22]）; %产生高斯白噪声 
% a=0.1;
% noise=a*wgn（1100010）;  %产生1行1000列的高斯白噪声为10dBW
%噪声
P = sum（abs（e_2psk）.^2）/length（e_2psk）;
% 计算噪声功率
Var_noise = P / 10^（snr/10）;
% 生成噪声并归一化
Noise = randn（1length（e_2psk））; 
Noise = Noise - mean（Noise）;
Noise = Noise / std（Noise）;%std函数是用来计算标准偏差的一个函数
% 混合
noise = Noise*sqrt（Var_noise）;%//噪声

figure（5）
plot（tnoise）
title（‘高斯白噪声‘）; 
%axis（[010-22]）;

%调制信号叠加高斯白噪声 
e2psk=e_2psk+noise;
figure（6）
plot（te2psk） 
title（‘2psk叠加白噪声波形‘）
figure（7）
plot（ff/1000fftshift（abs（fft（e2psk））））;
title（‘加噪声后的2psk信号的频谱‘）;
xlabel（‘频率（KHz）‘）;
%%

% %axis（[010-22]）; %带通滤波器的设计%
% [BB1AB1]=ellip（50.560[24006400]*2/100000）;      %带通滤波器通带上、下界频率分别为2400Hz6400Hz通带最大衰减为0.5dB阶数为5，计算带通滤波器系统函数分子分母多项式系数向量BB1和AB1  
% e_2psk1=filter（BB1AB1e2psk）;                       %带通滤波器滤除 带外噪声  
fp=5000;%通带截止
Fs= 1750;%阻带截止
%Fs= 1000;%采样频率
%阻带attn= 40dB
ws=Fs*2/fs
wp=fp*2/fs
[N Wp] = ellipord（wpws140）;
[ba]=ellip（N140Wp）;
[Hw] = freqz（ba）;
plot（w/2/pi*fs20*log10（H））; grid;

e_2psk1=filter（bae2psk）;%用设计的滤波器进行滤波
figure（8） 
plot（te_2psk1） 
title（‘调制信号通过带通滤波器‘） %相干解调与载波相乘 
e_2psk2=e_2psk1.*carry*2; 
figure（9） 
plot（te_2psk2）  
title（‘调制信号与载波相乘进行相干解调‘） 
%axis（[010-22]）;
%%

%设计低通滤波器  
[BB2AB2]=ellip（50.560700*2/100000）;  %通滤波器通带截止频率分别为700Hz通带最大衰减为0.5dB阶数为5，计算低通滤波器系统函数分子分母多项式系数向量BB2和AB2 
e_2psk3=filter（BB2AB2e_2psk2）; 
figure（10）
plot（te_2psk3） 
title（‘解调信号通过低通滤波器‘） %抽样判决 
for i=0:9    
    if（e_2psk3（（i+1）*100）>0）      
        e_2psk4（i*100+1:（i+1）*100）=zeros（1100）;   
    else         e_2psk4（i*100+1:（i+1）*100）=ones（1100）; 
    end
end



 figure（11） 
 plot（te_2ps