2psk调制通信系统

实用文案2psk调制通信统一，设计务与要求课程设计需要运用编程现调制解调过程输出其调制及解调过程中的波形，讨论其调制和解调效果。，实验本原理数字调制技术的两种方法：①利用模拟调制的方法去实现数式调制把数字调制看成是模拟调制的一个特例把数字基带信号当做模拟信号的特殊况处理。②利用数字信号的离散取值特点过开关键控载波而实现数字调制这种方法通常称为键控法，比如对载波的相位进行控，便可获得相移键控）基本的调制方式。相移键控是利用载波的相位变化来传递字信息，而振幅和频率保持不变。在2PSK中，通常用初始相位和π分别表示二进制1”和“2psk调制器可以采用相乘器，也可以采用相位选择器就模拟调制法而言，与产2ASK信号的方法比较是对s(t)要求不2PSK号可以看作是双极性基带信号作用下的DSB调幅信号。而就键控法来说，用字基带信号s(t)控制开关电路，选择不同相位的载波输出，这时s(t)为单极性NRZ或极性脉冲序列信号均可2PSK信号属于DSB信号。本次实验采用的的模拟相乘法即过载波和双极性不归零码的相乘得到2psk信号则2psk信号产生的调制原理框图和时域达式如下：_

,概,概率为-图1时域表达式图2调制原理框图2psk典型波形如下：标准文档

实用文案三，仿真方和参数设置参数设置如下所示：每码元采样点数Fn=500；码元数m=50；载波频率fc=2；码元速率Rm=1；加入的白噪声的信噪比snr分别10,30,50MATLAB产生2psk信号的程序框如下：标准文档

实用文案四，实验结与分析产生的双极性非归零码波形，2psk信号波形和频谱如;标准文档

////实用文案////

bds

分析如下：当二进制符号为“0”时，调制号相位差为0而当二进制符号为1”时，调制信号相位差为π。信号的频谱由连续和离散谱构成。当加入白噪声后，信号波和频谱如下：时调制信号

2psk调制信号频谱

10-1

543210123

456789

0-9-8-7-6-5-4-3

-2-1012345678910时调制信号

2psk调制信号频谱

10-1

543210123

456789

0-9-8-7-6-5-4-3

-2-1012345678910时调制信号

2psk调制信号频谱

10-1

54321012345

6789

0-9-8-7-6-5-4-3-2-101

2345678910由图可知加入白噪声后，2psk信发生了失真，随着信噪比的增加2psk失真越来越小。这是因为信噪比越大，2psk信的平均功率越大，信号越强，而噪声平均功率越小，噪声信号越弱因此失真越小。由于噪声没有频谱此无论信噪比如何变化，其最终的频谱都没发生任何变化，保持原样。五，源程序码标准文档

实用文案产生双极性非归零码和2psk信的子函数：function[dnrz,bpsk]=fbpsk(L,M,f)%产生双极型归零码%L为每个码元重复次数，则码间隔Ts%M为码元数目，f为载波频率t=0:1/(L-1):1;保载波与码元时域分辨率相同fc=[];n=1;赋值forii=1:M将进行重复拓展fc1=sin(2*pi*f*t);载fc=[fcfc1];按列方向合并为一个矩阵whilerem(n,2)~=0;保仿真时有偶数个1nrz=round(rand(1,M));随产生单极性非归零码n=length(find(nrz==1));endp=nrz==0;找中0的置nrz(p)=-1;将0为1dnrz=nrz(ones(1,L),:);将nrz进列拓展，并且每个元素重复L次dnrz=reshape(dnrz,1,L*M);重排成1*(L*M)矩阵endbpsk=dnrz.*fc;产信号End对信号进行FFT的子函数：function[X,f]=fftseq(x,fs)%本函数的功能为完成fft变换%x为输入信号%X为x的傅里叶变换L=length(x);%求出信x长度N=2^(nextpow2(L));%N为大等于n的的整数次幂的数中最小的一个按基二算法X1=fft(x,N);%计算x的N离散傅里叶变换X=fftshift(X1);X=X/fs;df=fs/N;%频域分率f=[0:df:df*(N-1)]-fs/2;主函数如下：clearall;closeall;标准文档

实用文案Fn=500;每码元采样点数为500m=50;%m码元数为50fc=2;%fc载波频率为Rm=1;%Rm码元速率为Tn=1/Rm;码元间隔tb=Tn/Fn;时域分辨率N=Fn*m;%N总样点数T=N*tb;计算时域上时间总长度t=linspace(-T/2,T/2,N);产生时间取样点矩阵，在长度内等分为N份[dnrz,bpsk]=fbpsk(Fn,m,fc);调用子函数，产生dnrz和dpsk信号[Sbpsk,f1]=fftseq(bpsk,Fn);%Fbpsk为信频谱Fbpsk=Sbpsk;figure(1);subplot(221);画双极性不归零码的波形plot(t,dnrz);title('双极性不归零码'xlabel('T/s');ylabel('幅度'axis([0,10,-1.5,1.5]);set(gca,'xtick',[0:0.5:10]);gridon;subplot(223);画2psk信的波形plot(t,bpsk);title('2psk调制信号');xlabel('T/s');ylabel('幅度'axis([0,10,-1.5,1.5])set(gca,'xtick',[0:0.5:10]);gridon;subplot(222);画2psk信的频谱plot(f1,Fbpsk);title('2psk调制信号频谱');xlabel('F/hz');ylabel('s/db');axis([-10,10,0,5.0]);gridon;snr1=10;增snr=10的噪声awgn1=awgn(bpsk,snr1);[Sawgn1,f2]=fftseq(awgn1,Fn);Fawgn1=Sawgn1;%Fawgn1为制信号的频谱标准文档

实用文案snr2=30;增snr=30的噪声awgn2=awgn(bpsk,snr2);[Sawgn2,f2]=fftseq(awgn2,Fn);Fawgn2=Sawgn2;%Fawgn2为制信号的频谱snr3=50;增snr=50的噪声awgn3=awgn(bpsk,snr3);[Sawgn3,f2]=fftseq(awgn3,Fn);Fawgn3=Sawgn3;%Fawgn3为制信号的频谱figure(2);subplot(321);画snr=10时2psk制信号的波形plot(t,awgn1);title('snr=10时2psk调制信');xlabel('T/s');ylabel('幅度'axis([0,10,-1.5,1.5]);gridon;subplot(322);画出snr=102psk调制信号的频谱plot(f2,Fawgn1);title('snr=10时2psk调制信频');xlabel('F/hz');ylabel('s/db');axis([-10,10,0,5.0]);set(gca,'xtick',[-10:1:10]);gridon;subplot(323);画出snr=302psk调制信号的波形plot(t,awgn2);title('snr=30时2psk调制信');xlabel('T/s');ylabel('幅度'axis([0,10,-1.5,1.5]);gridon;subplot(324);画出snr=302psk调制信号的频谱plot(f2,Fawgn2);title('snr=30时2psk调制信频');xlabel('F/hz');ylabel('s/db');axis([-10,10,0,5.0]);set(gca,'xtick',[-10:1:10]);gridon;subplot(325);画出snr=502psk调制信号的波形plot(t,awgn3);title('snr=50时2psk调制信');xlabel('T/s');ylabel('幅度'axis([0,10,-1.5,1.5]);标准文档

实用文案gridon;subplot(326);画出snr=50时2psk调信号的频谱plot(f2,Fawgn3);title('snr=50时2psk调制信频');xlabel('F/hz');ylabel