要更改调制参数，请打开 .m 文件。

在瑞利衰落信道中使用 PSK 调制的 LDPC 编码大多数函数都内置在 comm 中的 matlab 函数中。 工具箱。 这些也可以在代码中轻松配置以满足您的要求。 在给定的代码中，设置 SNR、codeRate、调制阶数 (psk) 和要模拟的帧数。 请随时与我联系，让我知道您的意见。

在此代码的输出中可以看到星座图和消息信号，发射信号和接收信号的波形。

包含带通滤波器、低通滤波器、基带信号生成、波形表示等子函数。使用相干解调实现2psk的功能。里面的子函数可以单独调用。

摘要 Psk调制是通信系统中最为重要的环节之一，Psk调制技术的改进也是通信系统性能提高的重要途径。本文首先分析了数字调制系统的基本调制解调方法，然后，运用Matlab及附带的图形仿真工具 ——Simulink设计了这几种数字调制方法的仿真模型。通过仿真，观察了调制解调过程中各环节时域和频域的波形，并结合这几种调制方法的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。最后，在仿真的基础上分析比较了各种调制方法的性能，并通过比较仿真模型与理论计算的性能，证明了仿真模型的可行性。另外，本文还利用Matlab的图形用户界面（GUI）功能为仿真系统设计了一个便于操作的人机交互界面，使仿真系统更加完整，操作更加方便。 关键词：数字调制；分析与仿真；Matlab；Simulink；GUI图形界面

相移键控是一种数字调制方案，它通过改变或调制参考信号（载波）的相位来传送数据。 在代码中，载波是正弦波，按照发射器和接收器组件的框图，获得输出。

数字信号对载波相位调制称为相移键控（即相位键控）PSK（ Phase－Shift Keying ）。 数字相位调制（相位键控）是用数字基带信号控制载波的相位，使载波的相位发生跳变的一种调制方式。 二进制相位键控用同一个载波的两种相位来代表数字信号。 由于PSK系统抗噪声性能优于ASK和FSK，而且频带利用率较高，所以，在中、高速数字通信中被广泛采用。

matlab开发-ldpc在一个GNChannel中使用psk调制进行编码。在awgn通道中使用psk执行lpdc

用于生成4psk调制的小程序，很好用，试试看

clear all; close all; fs=8e5;%抽样频率 fm=20e3;%基带频率 n=2*(6*fs/fm); final=(1/fs)*(n-1); fc=2e5; % 载波频率 t=0:1/fs:(final); Fn=fs/2;%耐奎斯特频率 %用正弦波产生方波 %========================================== twopi_fc_t=2*pi*fm*t; A=1; phi=0; x = A * cos(twopi_fc_t + phi); % 方波 am=1; x(x>0)=am; x(x<0)=-1; figure(1) subplot(321); plot(t,x); axis([0 2e-4 -2 2]); title('基带信号'); grid on car=sin(2*pi*fc*t);%载波 ask=x.*car;%载波调制 subplot(322); plot(t,ask); axis([0 200e-6 -2 2]); title('PSK信号'); grid on; %===================================================== vn=0.1; noise=vn*(randn(size(t)));%产生噪音 subplot(323); plot(t,noise); grid on; title('噪音信号'); axis([0 .2e-3 -1 1]); askn=(ask+noise);%调制后加噪 subplot(324); plot(t,askn); axis([0 200e-6 -2 2]); title('加噪后信号'); grid on; %带通滤波 %====================================================================== fBW=40e3; f=[0:3e3:4e5]; w=2*pi*f/fs; z=exp(w*j); BW=2*pi*fBW/fs; a=.8547;%BW=2(1-a)/sqrt(a) p=(j^2*a^2); gain=.135; Hz=gain*(z+1).*(z-1)./(z.^2-(p)); subplot(325); plot(f,abs(Hz)); title('带通滤波器'); grid on; Hz(Hz==0)=10^(8);%avoid log(0) subplot(326); plot(f,20*log10(abs(Hz))); grid on; title('Receiver -3dB Filter Response'); axis([1e5 3e5 -3 1]); %滤波器系数 a=[1 0 0.7305];%[1 0 p] b=[0.135 0 -0.135];%gain*[1 0 -1] faskn=filter(b,a,askn); figure(2) subplot(321); plot(t,faskn); axis([0 100e-6 -2 2]); title('通过带通滤波后输出'); grid on; cm=faskn.*car;%解调 subplot(322); plot(t,cm); axis([0 100e-6 -2 2]); grid on; title('通过相乘器后输出'); %低通滤波器 %================================================================== p=0.72; gain1=0.14;%gain=(1-p)/2 Hz1=gain1*(z+1)./(z-(p)); subplot(323); Hz1(Hz1==0)=10^(-8);%avoid log(0) plot(f,20*log10(abs(Hz1))); grid on; title('LPF -3dB response'); axis([0 5e4 -3 1]); %滤波器系数 a1=[1 -0.72];%(z-(p)) b1=[0.14 0.14];%gain*[1 1] so=filter(b1,a1,cm); so=so*10;%add gain so=so-mean(so);%removes DC component subplot(324);