【图像传输】QPSK调制解调数字带通图像传输（含误码率）【含Matlab源码 3602期】

最新推荐文章于 2024-06-24 13:15:19 发布

Matlab领域

最新推荐文章于 2024-06-24 13:15:19 发布

阅读量1.7k

点赞数 22

分类专栏： Matlab图像处理（高阶版）文章标签： matlab

本文链接：https://blog.csdn.net/weixin_63266434/article/details/134912987

版权

Matlab图像处理（高阶版）专栏收录该内容

297 篇文章 18 订阅

订阅专栏

在这里插入图片描述

⛄一、QPSK调制解调数字带通图像传输

软件无线电在无线通信领域被称为是自模拟通信过渡到数字通信之后的又一次革命，在军用和民用方面都有着广阔的应用。它是一种新的无线通信技术，基于通用的可编程的硬件平台，把尽可能多的通信功能用软件实现，从而使系统的改进和升级都非常方便，容易实现不同通信系统之间的兼容。调制与解调是整个数字通信系统基带处理的基本组成部分，也基本上构建了数字通信系统的总体性能框架。其中，QPSK调制与解调是一种性能较好，易于实现且已广泛应用于实践的调制解调方式。

正交相移键控(Quadrature Phase Shift Keying，QPSK)是一种数字调制方式。它分为绝对相移和相对相移两种。由于绝对相移方式存在相位模糊问题，所以在实际中主要采用相对移相方式DQPSK。目前已经广泛应用于无线通信中，成为现代通信中一种十分重要的调制解调方式。

正交相移键控（Quadrature Phase Shift Keying，QPSK）是一种数字调制方式。它分为绝对相移和相对相移两种。由于绝对相移方式存在相位模糊问题，所以在实际中主要采用相对移相方式DQPSK。QPSK是一种四进制相位调制，具有良好的抗噪特性和频带利用率，广泛应用于卫星链路、数字集群等通信业务。本课题通过使用QPSK建立如下的图像传输流程：
在这里插入图片描述
采用QPSK 调制方式，同时保证了信号传输的效率和误码性能。

一般的QPSK的实施，也表明高阶PSK的实施。在星座图中的正弦和余弦波用来传输方面的书面符号：

这就产生了四个阶段π/ 4，3π/ 4，5π/ 4和7π/ 4需要。

这个结果与单位的基础上功能在一个两维的信号空间被用作信号的同相分量和正交分量信号的第二首的基础功能。

而扩频通信，其特点是传输信息所用的带宽远大于信息本身带宽。在这种信息传输方式中，发送端以扩频编码（一般是伪随机码）进行扩频调制。在接收端则用同样的码进行相关同步接收、解扩及恢复所传信息数据。

⛄二、部分源代码

clear
N1=2400;
[BIT,data,data2,a1,a2,a3] =PICTURE_EXCHANGE(‘my_picture3.jpg’,‘603_1063_3.txt’);%读取高清图片
figure(1);
imshow(uint8(data2));
title(‘发送的图片’);
L1=length(BIT);
fs=8000;
alpha=0.2;
RB=2400;
Tb=1/RB;
fc=1800;
Rb=2RB;
M=4;
SNR=10;
ERROR_COUNT=0;
Redundancy=8;
PE_COUNT=0;
SNR_OLD=0;
SNR_NEW=0;
SNR_NEW_OUT=0;
My_Redundancy=[0 0 0 0];
[BIT,N] = ADD_ZERO(BIT,Rb-Redundancy);
RECEIVER_1= zeros(1,length(BIT));
h2=my_fitter(RB,fs,fc,Tb,alpha,‘RECEIVE_COS’);
draw_filter(h2,fs,3,2,8,‘COS的带通滤波器’);
h3=my_fitter(RB,fs,fc,Tb,alpha,‘RECEIVE_SIN’);
draw_filter(h3,fs,3,3,8,‘SIN的带通滤波器’);
h1=my_bpcosfir(fs,0,Rb2,N1,alpha,‘COS’);
h4=match_filter(2401,Tb,0.2);
draw_filter(h4,fs,3,1,8,‘发送滤波器’);

LPF_FLAG=‘Enable_FIR’;
BPF_FLAG=‘Enable_FIR’;

% 分帧处理，一次发送2400个码元，即处理4800个bit
for i=1:1:length(BIT)/(Rb-Redundancy)
i
% 接收缓冲区
Temp_1=BIT((Rb-Redundancy)*(i-1)+1:(Rb-Redundancy)*i);
%加冗余
Temp=[My_Redundancy’; Temp_1 ;My_Redundancy’];
% GENERATED_QPSK为发射机发送的QPSK信号，h1为发送滤波器响应,LPF_FLAG为是否使用发送滤波器
[GENERATED_QPSK,u4_send]=GENERATE_QPSK(Temp,fc,fs,h1,LPF_FLAG);
% 画频谱图
draw_fft(9,‘加噪前’,3,1,GENERATED_QPSK,fs);
% 加噪声，SNR为信噪比，GENERATED_QPSK为加噪后的信号
% GENERATED_QPSK_NOISE=awgn(GENERATED_QPSK,SNR,‘measured’);
GENERATED_QPSK_NOISE=GENERATED_QPSK;
% [GENERATED_QPSK,NOISE] = add_noise(GENERATED_QPSK,SNR);
SNR_OLD=SNR_singlech(GENERATED_QPSK,GENERATED_QPSK_NOISE);
% 画频谱图
draw_fft(9,‘加噪后’,3,2,GENERATED_QPSK_NOISE,fs);
% 载波同步
% x_sync = sync_two_signals( x_resampled,x_training_wave,idx )
% 相干解调，h2与h3为两个带通滤波器，BPF_FLAG为是否使用接收带通滤波器
[Demodulated_QPSK,COS_SIN] = Demodulation_QPSK(GENERATED_QPSK_NOISE,fc,fs,h2,h3,BPF_FLAG);
% draw_fft(9,‘带通后’,3,3,COS_SIN,fs);
SNR_NEW=SNR_singlech(GENERATED_QPSK(12000:36000),real(COS_SIN(12000:36000)));

SNR_NEW_OUT=(SNR_NEW+SNR_NEW_OUT);
% 计算误码率并赋值
for k=1:1:Rb-Redundancy
if( Demodulated_QPSK(k)~= Temp_1(k))
ERROR_COUNT=ERROR_COUNT+1;
end
if(mod(k,2)==1)
if(Demodulated_QPSK(k)~= Temp_1(k)||Demodulated_QPSK(k+1)~= Temp_1(k+1))
PE_COUNT=PE_COUNT+1;
end

end
    
    RECEIVER_1((Rb-Redundancy)*(i-1)+k)=Demodulated_QPSK(k);

end
end
SNR_NEW_OUT=SNR_NEW_OUT/(length(BIT)/(Rb-Redundancy));
%将二进制数据变为图片
Picture_Show(RECEIVER_1,L1,a1,a2,a3);
FIR_FLAG(LPF_FLAG,BPF_FLAG);
% fprintf(1,‘当前带通滤波前信噪比为：%dDB\n’,SNR);
% fprintf(1,‘滤波后信噪比为：%fDB\n’,SNR_NEW_OUT);
fprintf(1,‘输出信号的误码率为：%d\n’,PE_COUNT/L1)
fprintf(1,‘理论误码率为：%d\n’,1-(1-0.5*erfc(sqrt((10^{(SNR_NEW_OUT/10))/2)))}2)