bpsk_ber.m

来自「Digital Bandpass Transmission」· M 代码 · 共 19 行

%N=no of bits detected in errors
%M=total number of bits transmitted

function [N,M]=bpsk_ber(Rxout,Txdata)
Rxdata=Rxout(1:16:length(Rxout));
Rxdata=Rxdata(11:length(Rxdata));
figure
stem(Rxdata)
hold on
plot(Txdata(1:length(Rxdata)),'ro')
mean(Rxdata.^2)
Rxdata=sign(Rxdata);
figure
errors=Rxdata-Txdata(1:length(Rxdata));
stem(Rxdata-Txdata(1:length(Rxdata)))
N=sum(abs(errors));
M=length(errors);

PER = N/M