实验四调制解调实验报告.docx

文档介绍：实验四调制解调(BPSK,QPSK,信噪比)
本章目标
掌握数字频带传输系统调制解调的仿真过程
掌握数字频带传输系统误码率仿真分析方法
调制解调原理
本章以BPSK为例,仿真说明数字频带传输的整个过程
假定:信道为加性高斯白噪声信道,其均值为0、方差为,采用矩形成形;发射端BPSK调制信号为:


经信道传输,接收端输入信号为:
经相干解调,匹配滤波,定时恢复后输出:
当1,0独立等概出现时,BPSK系统的最佳判决门限电平。故判决规则为在取
样时刻的判决值大于0,判1,小于0,判0。
实验要求
1) 分别编写BPSK与QPSK调制解调系统的Matlab仿真程序,要求:
发送滤波器与接收滤波器均为根升余弦滚降滤波器;
信道噪声为加性高斯白噪声

2) 绘制BPSK与QPSK调制下的误码率与信噪比曲线图,并与理论曲线进行对比
:
调制解调:
编程思想:
1. 先用radint函数产生二进制信源。
2. 用sym[-1;1], A=sym(b+1)进行符号映射。
3. 进行采样,用rcosine函数产生升余弦滤波器,信号通过滤波器
4. 调制,X=S. cos(2*pi*fc*t);
5. 用awgn函数加高斯白噪声
,r=x. cos(2*pi*fc*t)
7. 采样,补领
8. 求误码率
主程序:
%BPSK snr=5
clear all;
close all;
%------------------
Snr=5;
fs=16;
N=1000;
Ts=1;
fd=1;
ts=1/fs;
t=0:ts:N*Ts-ts;
t=t.';
N_sample=Ts*fs ;
%-----------»
%
B=randint(1,N);
sym=[-1;1];
A=sym(B+1);

H=rcosine(fd, fs, 'sqrt');
S1=upsample([A], fs/fd);
S=filter(H, 1, S1);

X1=S.*cos(2*pi*fc*t);
figure(1)

plot(t(1 :200),X1(1 :200));

X=awgn(X1,Snr-10*log10(fs/(2*fd)),'measured');

r=X.*cos(2*pi*fc*t);


S2=filter(H, 1, r);
S3=S2(5/fd*fs+1:end);
Sr=zeros(N-5,1);
for k=1:(N-5)
M=S3(k*Ts*fs);
if(M>0)
Sr(k)=1;
else
Sr(k)=-1;
end
end
A2=A(1:end-5);
biterror=sum(A2~=Sr)/(N-5);
Ber=biterror


运行程序,得到调制信号图和误比特率
Ber =
调制解调误码率曲线。
编程思想:
讲上述的bpsk调制系统的snr=5改为 snr=-10::12 用for 循环得到每个snr下的误码率,再用semilogy函数画出误码率曲线,同时画出理论误码率曲线。
主程序:
%BPSK
clear all;
close all;
%------------------
fc=2;
fs=16;
N=100000;
Ts=1;
fd=1;
ts=1/fs;
t=0:ts:N*Ts-ts;
t=t.';

N_sample=Ts*fs ;


snr=-10::12;


Ber=size(snr);
ber0=size(snr);
for j=1:length(snr)

B=randint(1,N);
sym=[-1;1];
A=sym(B+1);

Snr=snr(j);

H=rcosine(fd, fs, 'sqrt');
S1=upsample([A], fs/fd);
S=filter(H, 1, S1);

X1=S.*cos(2*pi*fc*t);

X=awgn(X1,Snr-10*log10(fs/(2*fd)),'measured');
r=X.*cos(2*pi*fc*t);


S2=filter(H, 1, r);
S3=S2(5/fd*fs+1:end);
Sr=zeros(N-5,1);
for k=1:(N-5)