BPSK.zip_BPSK调制解调_bpsk调制

%BPSK %公共参数设置 clc; %清屏幕 close all; clear all; %清内存 code_len=50;%码元数 fre=1; %码元速率 T=1/fre; %码元周期 fs=100*fre;%抽样频率 Ts=1/fs; %抽样周期 SNR=0 for i=1:code_len data(i)= 2*round(rand)-1;%产生0、1等概的随机数据并进行双极性化处理 end exdata=[];%把数据扩展成方波 for i=1:length(data) for rep=1:100 exdata= [exdata data(i)];%离散化处理,因为matlab不能作连续的波形，一个周期采50个点，共采样2500个点 end end %调制 t=Ts:Ts:(code_len*T); carrier=cos(4*pi*fre*t);%离散化处理,与随机数对应，抽样2500个点； bpsk=carrier.*exdata; bpsk1=bpsk; %加噪声 Eb=sum((abs(bpsk)).^2)/(code_len);%计算信号功率 No=Eb/(10^(SNR/10)); %计算噪声功率，SNR=Eb/No noise=randn(size(bpsk)).*sqrt(No/2); bpsk=bpsk+noise; %解调 uncarry=bpsk.*carrier;%信号乘以本地载波，便是主要的解调步骤 % 低通滤波器 f=1/100; % 截止频率 qifs=1; %采样率 [a,b]=butter(4,f/(qifs/2)); bpsk2=filter(a,b,uncarry); dec1=[]; for inc=1:100:length(uncarry) dec=trapz(inc:inc+99,uncarry(inc:inc+99));%计算uncarry关于inc的积分 dec1=[dec1 dec]; end %门限判别 demod=[]; for i=1:length(dec1) if dec1(i)>0 demod=[demod 1];%大于0就取为1 else demod=[demod -1];%小于0就取为-1 end end %计算误码率 error=0; for i=1:length(demod) if data(i)~=demod(i) error=error+1; end end error ber=error/code_len %画图 figure(1);%比较发送码与接收码 subplot(2,1,1); stem(data,'r'); title('(1,-1)code'); subplot(2,1,2); stem(demod); title('(1,-1)demod code'); hold off; figure(2);%比较调制波形、加噪波形、解调波形、经过低通后的波形 subplot(4,1,1); plot(bpsk1,'r'); title('bpsk'); subplot(4,1,2); plot(bpsk); title('bpsk+noise'); subplot(4,1,3); plot(uncarry); title('demod bpsk'); subplot(4,1,4); stairs(bpsk2); title('low pass filter'); hold off;