BPSK_bpsk调制解调MATLAB仿真_

clear; clc; NSym=10; x=randint(1,NSym); xDualPole=2*x-1; fb=1000; fc=4000; fs=16000; OverSamp=fs/fb; Delay=5; alpha=0.25; h_sqrt=rcosine(1,OverSamp,'fir/sqrt',alpha,Delay); SendSignal_OverSample=kron(xDualPole,[1 zeros(1,OverSamp-1)]); SendShaped=conv(SendSignal_OverSample,h_sqrt); figure; subplot(2,1,1);plot(SendShaped);title('脉冲成型后的时域波形') subplot(2,1,2);plot(abs(fft(SendShaped)));title('脉冲成型后的频域波形') N=0:length(SendShaped)-1; CarrierWave=sin(2*pi*fc*N/fs); ModemWave=SendShaped.*CarrierWave; figure; subplot(2,1,1);plot(ModemWave);title('调制后的时域波形') subplot(2,1,2);plot(abs(fft(ModemWave)));title('调制后的频域波形') DemodWave=ModemWave.*CarrierWave; figure; subplot(2,1,1);plot(DemodWave);title('解调后的时域波形') subplot(2,1,2);plot(abs(fft(DemodWave)));title('解调后的频域波形') RcvMatched=conv(DemodWave,conj(h_sqrt)); figure; subplot(2,1,1);plot(RcvMatched);title('匹配接收后的时域波形') subplot(2,1,2);plot(abs(fft(RcvMatched)));title('匹配接收后的频域波形') SynPosi=Delay*OverSamp*2; SymPosi=SynPosi+(0:OverSamp:(NSym-1)*OverSamp); RcvSignal=RcvMatched(SymPosi); for i=1:NSym if(RcvSignal(i)>0) RcvBit(i)=1; else RcvBit(i)=-1; end end figure; subplot(2,1,1);stem(xDualPole);title('发送的信号波形') subplot(2,1,2);stem(RcvBit);title('接收的信号波形')