MATLAB在数字调制解调中的应用.doc

MATLAB在数字调制解调中的应用系别：电信系 班级:XXXXXXXXX 姓名：XXX 学号：XXXXXXX XXXX 指导老师：XXX 一题目 ： 数字调制与解调二要求：用matlab模拟数字调制与解调三步骤：1基本知识：按数字调制的方法分类可以分为多进制幅度键控 (M-ASK)、正交幅度键控(QASK)、多进制频率键控(M-FSK)以及多进制相位键控(M-PSK)。数字调制包括数模转换和模拟调制两部分，如图6.25所示。 1M元幅度键控调制(M-ASK)M元幅度键控调制包含两部分，M-ASK映射和模拟幅度调制。M-ASK映射将输入的数字码符映射到区间-x, x，数字码符取值范围是 0, M-1区间内的整数。输出信号的幅值分别为-x, x的M-1等分点。 2M元正交键控调制(M-QASK) M-QASK是数字调制使用得最多的一种方法。它一般将输入的数字码符映射成为同相和正交的两个独立分量，然后用模拟QAM法对它们进行调制；在接收方，接收的信号被解调为同相和正交信号，从它们映射过程中恢复原始信号。有许多种方法可将输入信号映射成同相和正交分量。通信工具箱为M-QASK提供三种方案，即平面直角点集、圆点集和用户定义任意点集。3M 元频率键控调制(M-FSK) M 元频率键控调制是通过使用输入信号控制输出信号的频率来实现对数字信号的调制。M-FSK调制过程分成两部分，即映射和模拟调制。映射过程将输入信号反映成载波频率变化，模拟调制即FM。M-FSK的解调有两种方法，即相干法和非相干法。相干方法要求事先知道已调信号的相位，而非相干方法则不要求相位信息 , 它可以在解调过程中恢复已调信号的相位信息。4M-元相位键控调制(M-PSK)M-PSK通过改变已调信号的相位信息来实现对数字信号的调制。M-PSK设置不同的初 相移位以区别不同的数字码符。M-PSK调制器输入信号的取值区间0, M-1，数字 i对应的相位位移为 1 带通数字调制/解调1. 数字带通调制函数dmod( )格式：y=dmod(x, Fc, Fd, Fs, method )功能：用频率为Fc(Hz)的载波去调制数字信号，然后再输出采样频率为Fs(Hz)的仿真已调信号，Fs/Fd 必须是一个正整数。当输入x是一个矩阵时，x中每列做不相关处理，输出矩阵y的行数是Fs/Fd(x的行数)。若变量Fs是一个两维矢量。当它为矢量时，第一个元素是采样频率，第二个元素是调制载波信号的初相(弧度)，且默认值为0。要求FsFcFd。method 为调制方法，见表6.6。2数字带通解调函数ddemod( )格式：z=ddemod(y, Fc, Fd, Fs, method )功能：解调载波为Fc(Hz)的数字调制信号y，其抽样频率为Fd(Hz)，计算抽样频率为Fs(Hz)。该函数为dmod的逆过程，对频率与函数为dmod 中求要相同。当Fd为一个二维矢量时，Fd中的第二个值即为输出抽样点的时间偏移量(必须是整数，默认为0)。该函数测出接收码到编码中所有可能的码的距离，然后以离接收码最近的码作为输出数字码。除method =msk和method =fsk解调外，该函数可使用一个指定的低通滤波器，该滤波器的传输函数由分子num和分母den指定，滤波器的抽样时间为1/Fs，缺省滤波器时，在解调时无积分器。利用MATLAB对信号进行基带数字调制解调。 MATLAB程序如下：M=4; %设置M的数目Fd=1; %信号采样速率Fs=32; %采样频率SNRperBit=5; %信噪比adjSNR=SNRperBit-10*log10(Fs/Fd)+10*log10(log2(M);x=randint(5000,1,M); %源信号%正交FSK调制tone=0.5; %频隙间隙f=Fd/2randn(state,1945724);%设置RANDN产生器的状态w1=dmodce(x,Fd,Fs,fsk,M,tone); %调制y1=awgn(w1,adjSNR,measured,dB); % 对调制信号添加噪声z1=ddemodce(y1,Fd,Fs,fsk,M,tone);%调制ser1=symerr(x,z1) %输出符号误差率%非正交FSK调制tone=0.25;randn(state,1945724); %设置RANDN产生器的状态w2=dmodce(x,Fd,Fs,fsk,M,tone); %调制y2=awgn(w2,adjSNR,measured,dB); % 对调制信号添加噪声z2=ddemodce(y2,Fd,Fs,fsk,M,tone); %解调ser2=symerr(x,z2) %输出符号误差率2 .GUI界面的制作3.用到的程序M元正交调制解调M=16;Fc=10;Fd=1;Fs=50;x=randint(100,1,M);y=domd(x,Fc,Fd,Fs,qask,M);ynoisy=y+.01*randn(Fs/Fd*100,1);z=demod(ynoisy,Fc,Fd,Fs,qask,M);s=symerr(x,z);t=0.1:0.1:10;subplot(2,1,1);plot(t,x);title()subplot(2,1,2);plot(t,z);title()figure;modmap(qask,M)基带调制解调M=4;Fd=1;Fs=32;SURperBit=5;adjSNR=SNRperBit-10*log10(Fs/Fd)+10*log10(log2( M);x=randint(5000),1,M);tone=0.5;randn(state1945724);w1=dmodce(x,Fd,Fs,fskM,tone);y1=awgn(w1,adjSNR,measured,dB);z1=ddemodce(y1,Fd,FsfskM,tone);ser1=symeer(x,z1);tone=0.25;randn(state1945726);w2=dmodce(x,Fd,Fs,fskM,tone);y2=awgn(w2,adjSNR,measured,dB);z2=ddemodce(y2,Fd,FsfskM,tone);ser2=symeer(x,z2);模数转换a=4;f0=0.02;L=50;n=1:(L-1);x=a.*cos(2.*pi.*f0.*n);bit=5;range=25;plot(x);interval=(a*2)/(range-1);partition = -a:interval:a;interval2 = (a*2)/(range)codebook=-a:interval2:a;index quantized= quantiz(x,partition,codeb