clear?%產(chǎn)生比特信號?
t=0:0.01:9.99;?
a=randint（110）;
snr=10;????????????????%信噪比
m=a（ceil（t+0.01））;?????%產(chǎn)生方波
figure（1）
plot（tm）??
title（‘產(chǎn)生隨機(jī)十比特二進(jìn)制比特序列‘）;?
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調(diào)制%?
e_2psk=cos（2*pi*fc*ts/2.5+pi*m）;
figure（3）
plot（te_2psk）?
title（‘2psk調(diào)制信號‘）;?


ff=?linspace（-fs/2?fs/2?length（e_2psk））;%頻率坐標(biāo)，單位Hz
figure（4）
plot（ff/1000fftshift（abs（fft（e_2psk））））;
title（‘未加噪聲前的2psk信號的頻譜‘）;
xlabel（‘頻率（KHz）‘）;

%%
%?axis（[010-22]）;?%產(chǎn)生高斯白噪聲?
%?a=0.1;
%?noise=a*wgn（1100010）;??%產(chǎn)生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函數(shù)是用來計算標(biāo)準(zhǔn)偏差的一個函數(shù)
%?混合
noise?=?Noise*sqrt（Var_noise）;%//噪聲

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

%調(diào)制信號疊加高斯白噪聲?
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]）;?%帶通濾波器的設(shè)計%
%?[BB1AB1]=ellip（50.560[24006400]*2/100000）;??????%帶通濾波器通帶上、下界頻率分別為2400Hz6400Hz通帶最大衰減為0.5dB階數(shù)為5，計算帶通濾波器系統(tǒng)函數(shù)分子分母多項式系數(shù)向量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）;%用設(shè)計的濾波器進(jìn)行濾波
figure（8）?
plot（te_2psk1）?
title（‘調(diào)制信號通過帶通濾波器‘）?%相干解調(diào)與載波相乘?
e_2psk2=e_2psk1.*carry*2;?
figure（9）?
plot（te_2psk2）??
title（‘調(diào)制信號與載波相乘進(jìn)行相干解調(diào)‘）?
%axis（[010-22]）;
%%

%設(shè)計低通濾波器??
[BB2AB2]=ellip（50.560700*2/100000）;??%通濾波器通帶截止頻率分別為700Hz通帶最大衰減為0.5dB階數(shù)為5，計算低通濾波器系統(tǒng)函數(shù)分子分母多項式系數(shù)向量BB2和AB2?
e_2psk3=filter（BB2AB2e_2psk2）;?
figure（10）
plot（te_2psk3）?
title（‘解調(diào)信號通過低通濾波器‘）?%抽樣判決?
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