利用Matlab simulink实现DPSK通信
通信原理DPSKFSKASKPSK检波及波形实现的MATLAB程序
通信原理DPSKFSKASKPSK检波及波形实现的MATLAB程序十、附录程序程序一2ASK调制解调程序及注释clear all close all i=10;%10个码元j=5000; t=linspace(0,5,j);%0-5之间产生5000个点行矢量,即分成5000份fc=10;%载波频率fm=i/5;%码元速率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号x=(rand(1,i))%rand函数产生在0-1之间随机数,共1-10个figure(2) plot(x) a=round(x);%随机序列,round取最接近小数的整数figure(3) stem(a)%火柴梗状图st=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st(m)=0; end else for m=j/i*(n-1)+1:j/i*n st(m)=1; end end end figure(1); subplot(421); plot(t,st); axis([0,5,-1,2]); title('基带信号st'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波s1=cos(2*pi*fc*t); subplot(422); plot(s1); title('载波信号s1'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制e_2ask=st.*s1; subplot(423); plot(t,e_2ask); title('已调信号'); noise =rand(1,j); e_2ask=e_2ask+noise;%加入噪声subplot(424);plot(t,e_2ask); title('加入噪声的信号'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调at=e_2ask.*cos(2*pi*fc*t); at=at-mean(at);%因为是单极性波形,还有直流分量,应去掉subplot(425); plot(t,at); title('与载波相乘后信号'); [f,af] = T2F(t,at);%通过低通滤波器[t,at] = lpf(f,af,2*fm); subplot(426); plot(t,at);title('相干解调后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1; if at(1,m*500+250)+0.5 for j=m*500+1:(m+1)*500; at(1,j)=0; end else for j=m*500+1:(m+1)*500; at(1,j)=1; end end end subplot(427); plot(t,at); axis([0,5,-1,2]); title('抽样判决后波形')程序二2FSK调制解调程序及注释clear all close all i=10;%基带信号码元数j=5000; a=round(rand(1,i));%产生随机序列t=linspace(0,5,j); f1=10;%载波1频率f2=5;%载波2频率fm=i/5;%基带信号频率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号st1=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end st2=t; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%基带信号求反for n=1:j; if st1(n) st2(n)=0; else st2(n)=1; end end; figure(1);subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]); subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=cos(2*pi*f1*t) s2=cos(2*pi*f2*t) subplot(413),plot(s1); title('载波信号s1'); subplot(414),plot(s2); title('载波信号s2'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制F1=st1.*s1;%加入载波 1 F2=st2.*s2;%加入载波 2 figure(2); subplot(411); plot(t,F1); title('F1=s1*st1'); subplot(412); plot(t,F2); title('F2=s2*st2'); e_fsk=F1+F2; subplot(413); plot(t,e_fsk); title('2FSK信号')%键控法产生的信号在相邻码元之间相位不一定连续nosie=rand(1,j); fsk=e_fsk+nosie; subplot(414); plot(t,fsk); title('加噪声后信号') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调st1=fsk.*s1;%与载波1相乘[f,sf1] = T2F(t,st1);%通过低通滤波器[t,st1] = lpf(f,sf1,2*fm); figure(3); subplot(311); plot(t,st1); title('与s1相乘后波形'); st2=fsk.*s2;%与载波2相乘[f,sf2] = T2F(t,st2);%通过低通滤波器[t,st2] = lpf(f,sf2,2*fm); subplot(312); plot(t,st2); title('与s2相乘后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1; if st1(1,m*500+250)st2(1,m*500+250); for j=m*500+1:(m+1)*500; at(1,j)=0; end else forj=m*500+1:(m+1)*500; at(1,j)=1; end end end; subplot(313); plot(t,at); axis([0,5,-1,2]); title('抽样判决后波形')程序三2PSK调制解调程序及注释clear all close all i=10; j=5000; fc=4;%载波频率fm=i/5;%码元速率B=2*fm; t=linspace(0,5,j); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号a=round(rand(1,i));%随机序列,基带信号figure(3); stem(a); st1=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end figure(1); subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%基带信号求反%由于PSK中的是双极性信号,因此对上面所求单极性信号取反来与之一起构成双极性码st2=t; for k=1:j; if st1(k) st2(k)=0; else st2(k)=1; end end; subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]); st3=st1-st2; subplot(413); plot(t,st3); title('双极性基带信号st3'); axis([0,5,-2,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=sin(2*pi*fc*t); subplot(414); plot(s1); title('载波信号s1'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制e_psk=st3.*s1; figure(2); subplot(511); plot(t,e_psk);title('e_2psk'); noise=rand(1,j); psk=e_psk+noise;%加入噪声subplot(512); plot(t,psk); title('加噪后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调psk=psk.*s1;%与载波相乘subplot(513);plot(t,psk); title('与载波s1相乘后波形'); [f,af] = T2F(t,psk);%%%%%%%%%%%通过低通滤波器[t,psk] = lpf(f,af,B); subplot(514); plot(t,psk); title('低通滤波后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1; if psk(1,m*500+250) for j=m*500+1:(m+1)*500; psk(1,j)=0; end else for j=m*500+1:(m+1)*500; psk(1,j)=1; end end end subplot(515); plot(t,psk); axis([0,5,-1,2]); title('抽样判决后波形')程序四2DPSK调制解调程序及注释clear all close all i=10; j=5000; fc=4;%载波频率fm=i/5;%码元速率B=2*fm; t=linspace(0,5,j); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号a=round(rand(1,i)); figure(4); stem(a); st1=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end figure(1); subplot(321); plot(t,st1); title('绝对码'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%差分变换b=zeros(1,i);%%%%%%%%全零矩阵b(1)=a(1); for n=2:10 if a(n)if b(n-1)=1 b(n)=0; elseb(n)=1; end else b(n)=b(n-1); end end st1=t; for n=1:10 if b(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end subplot(323); plot(t,st1); title('相对码st1'); axis([0,5,-1,2]); st2=t; for k=1:j; if st1(k) st2(k)=0; else st2(k)=1; end end; subplot(324); plot(t,st2); title('相对码反码st2'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=sin(2*pi*fc*t); subplot(325); plot(s1); title('载波信号s1'); s2=sin(2*pi*fc*t+pi); subplot(326); plot(s2); title('载波信号s2'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%信号调制d1=st1.*s1; d2=st2.*s2; figure(2); subplot(411); plot(t,d1); title('st1*s1'); subplot(412); plot(t,d2); title('st2*s2'); e_dpsk=d1+d2; subplot(413); plot(t,e_dpsk); title('调制后波形'); noise=rand(1,j); dpsk=e_dpsk+noise;%%%%%%%%%%加入噪声subplot(414); plot(t,dpsk); title('加噪声后信号');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调dpsk=dpsk.*s1;%%%%与载波s1相乘figure(3); subplot(411); plot(t,dpsk); title('与载波相乘后波形'); [f,af]=T2F(t,dpsk);%%%%通过低通滤波器[t,dpsk]=lpf(f,af,B); subplot(412); plot(t,dpsk); title('低通滤波后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决st=zeros(1,i);%%全零矩阵for m=0:i-1; if dpsk(1,m*500+250)st(m+1)=0; for j=m*500+1:(m+1)*500; dpsk(1,j)=0; end else for j=m*500+1:(m+1)*500; st(m+1)=1; dpsk(1,j)=1; end end end subplot(413); plot(t,dpsk); axis([0,5,-1,2]); title('抽样判决后波形') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%码反变换dt=zeros(1,i);%%全零矩阵dt(1)=st(1); for n=2:10; if (st(n)-st(n-1))=0(st(n)-st(n-1)) dt(n)=0; else dt(n)=1; end end st=t; for n=1:10 if dt(n) for m=j/i*(n-1)+1:j/i*n st(m)=0; end else for m=j/i*(n-1)+1:j/i*n st(m)=1; end end end subplot(414); plot(t,st); axis([0,5,-1,2]); title('码反变换后波形');。
BPSK(DPSK)调制解调实验指导书
电子科技大学通信学院《二相BPSK(DPSK)调制解调实验指导书》二相BPSK(DPSK)调制解调实验班级学生学号教师二相BPSK(DPSK)调制解调实验指导书二相BPSK(DPSK)调制解调实验一、实验目的1、掌握二相BPSK(DPSK)调制解调的工作原理。
2、掌握二相绝对码与相对码的变换方法。
3、熟悉BPSK(DPSK)调制解调过程中各个环节的输入与输出波形。
4、了解载波同步锁相环的原理与构成,观察锁相环各部分工作波形。
5、了解码间串扰现象产生的原因与解决方法,能够从时域和频域上分析经过升余弦滚降滤波器前后的信号。
6、掌握Matlab软件的基本使用方法,学会Simulink环境的基本操作与应用。
二、实验原理数字信号载波调制有三种基本的调制方式:幅移键控(ASK),频移键控(FSK)和相移键控(PSK)。
它们分别是用数字基带信号控制高频载波的参数如振幅、频率和相位,得到数字带通信号。
PSK调制在数字通信系统中是一种极重要的调制方式,它的抗干扰噪声性能及通频带的利用率均优于ASK幅移键控和FSK频移键控。
由于PSK调制具有恒包络特性,频带利用率比FSK高,并在相同的信噪比条件下误码率比FSK低。
同时PSK调制的实现也比较简单。
因此,PSK技术在中、高数据传输中得到了十分广泛的应用。
BPSK是利用载波相位的变化来传递数字信息,而振幅和频率保持不变。
在BPSK中,通常用初始相位0和π分别表示二进制“1”和“0”。
其调制原理框图如图1所示,解调原理框图如图2所示。
图1 BPSK的模拟调制方式由于在BPSK 信号的载波恢复过程中存在着载波相位0 和180 的不确定性反向,所以在实际的BPSK 通信系统设计中,往往采用差分编解码的方法克服这个问题。
差分编解码是利用前后信号相位的跳变来承载信息码元,不再是以载波的绝对相位传输码元信息。
差分编解码的原理可用下式描述。
1n n n d b d -=⊕ 1ˆˆˆn n n b d d -=⊕ 其中第一个公式为差分编码原理,第二个公式为差分解码原理。
实验四 2DPSK系统的仿真实验
图1
Communications Blockset / Source Coding 库下的 Differential Encoder 模块
Communications Blockset / Utility Functions 库下的 Unipolar to Bipolar Converter 模块
Communications Bll Passband Modulation/PM 库下的
码元间隔与数字信号 采样周期要一致
图4
进制数 要一致
图5 仿真相对相移键控方法产生 2DPSK 时,可以用 M-DPSK Modulator Passband 模块完成 2DPSK 的调制功能,该模块的参数设置如图 7 所示,注意各参数的设置方法;M-DPSK Demodulator Passband 模块完成 2DPSK 的解调功能,该模块的参数设置如图 8 所示,注意解 调模块的参数设置必须与调制模块的相关参数设置一致。
图6 46
2DPSK 信号的频谱如图 9 所示。2DPSK 信号时域波形如图 10 所示。
图7
四、实验内容
1. 进一步熟悉并掌握 Matlab/Simulink 基本库、通信库和 DSP 库中较为重要的一些功能 模块的作用以及相应功能参数的物理意义与设置方法。
2. 搭建 2DPSK 模拟法仿真模型如图 1 所示。设置系统参数并调试,同时观测并记录 A~D 各点的时域波形以及 D 点的频谱。
数字信号 的进制数
数字信号 采样周期
图3 43
数设置如图 4 所示,注意该模块差分的性质;单双极性变换 Unipolar to Bipolar Converter 模 块的参数设置如图 5 所示,注意进制数的设置值必需与 Random Integer Generator 模块一致; Random Number 模块产生一个高斯型分布的随机噪声,该模块的参数设置如图 6 所示。
通信网络基础实验报告 基于MATLAB-SIMULINK设计ASK、PSK、FSK通信仿真系统以及Simulink编程的优点和不足
通信网络基础实验实验报告姓名:学号:班级:通信网络系统仿真设计与实现实验目的:1、学习MATLAB软件,掌握MATLAB-SIMULINK模块化编程。
2、理解并掌握通信网络与通信系统的基本组成及其工作方式。
实验环境:1、软件环境:Windows2000/XP2、硬件环境:IBM-PC或兼容机实验学时:4学时、必做、综合实验实验内容:1、ASK调制解调的通信仿真系统;2、PSK调制解调的通信仿真系统:3、FSK调制解调的通信仿真系统。
实验要求:1、基于MATLAB-SIMULINK分别设计一套ASK、PSK、FSK通信系统。
2、比较各种调制的误码率情况,讨论其调制效果。
实验步骤:独立自主完成分析思考:通信系统中滤波器的参数你是如何设计选择的,为什么?Simulink编程的优点和不足是什么?1、对于ASK调制解调的通信系统调制:仿真结果显示如下:上图中CH1表示的是调制前的信号频谱,CH2表示的是ASK调制后的信号波形。
上图中第一张图是幅度调制前原始基带信号的波形,第二张图是幅度调制后通带信号的波形,第三张是解调信号的波形图。
有图可看出信道有一定的延迟。
由于在解调过程中没有信道和噪声,所以误码率相对较小,一般是由于码间串扰或是参数设置的问题,此系统的误码率为0.3636。
2、对于FSK调制解调的通信系统调制:仿真结果如下:2FSK基带调制信号频谱图CH1表示的是基带信号的功率谱,而CH2表示2FSK调制后通带信号的功率谱。
2FSK信号解调各点时间波形经过系统的仿真可以观察出系统的误码率为0.7273,如下图所示:3、对于PSK的调制解调通信系统调制:在二进制数字调制中,当正弦载波的相位随二进制数字基带信号离散变化时,则产生二进制移相键控(2PSK)信号. 在此用已调信号载波的0°和180°分别表示二进制数字基带信号的1 和0.用两个反相的载波信号进行调制,其方框图如下:2PSK信号调制的simulink的模型图其中Sin wave和Sin wave1是反相的载波,正弦脉冲作为信号源。
通信原理matlab课程设计2ASK、2FSK、2PSK、2DPSK调制解调matlab仿真
通信原理matlab课程设计--2ASK、2FSK、2PSK、2DPSK调制解调matlab仿真南昌大学通信原理课程设计报告题目: 2ASK、2FSK、2PSK、2DPSK调制解调matlab仿真姓名:学院:信工学院专业:指导教师:完成日期:2013 年5 月5日一、设计要求课程设计需要运用MATLAB 编程实现2ASK,2FSK,2PSK ,2DPSK 调制解调过程,并且输出其源码,调制后码元以及解调后码元的波形。
二、基本原理二进制数字调制技术原理数字信号的传输方式分为基带传输和带通传输,在实际应用中,大多数信道具有带通特性而不能直接传输基带信号。
为了使数字信号在带通信道中传输,必须使用数字基带信号对载波进行调制,以使信号与信道的特性相匹配。
这种用数字基带信号控制载波,把数字基带信号变换为数字带通信号的过程称为数字调制。
通常使用键控法来实现数字调制,比如对载波的振幅、频率和相位进行键控。
(1)振幅键控是利用载波的幅度变化来传递数字信息,而其频率和相位保持不变,在2ASK 中,载波的幅度只有两种变化状态,分别对应二进制信息‘0’和‘1’。
OOK (通-断键控)是一种常用的二进制振幅键控式模拟调制器法 键控法包络检波法)开关电路2e2e同步检测法(2) 一个2FSK 信号可以看成是两个不同载波的2ASK 信号的叠加。
其解调和解调方法和ASK 差不多。
2FSK 信号的频谱可以看成是f1和f2的两个2ASK 频谱的组合。
2FSK 信号的产生方法采用模拟调频电路来实现:信号在相邻码元之间的相位是连续变化的。
采用键控法来实现:相邻码元之间的相位不一定连续。
2FSK 信号的解调方法相干解调2e FSK2e FSK非相干解调(3) 2PSK 以载波的相位变化作为参考基准的,当基带信号为0时相位相对于初始相位为0, 当基带信号为1时相对于初始相位为180°。
调制器原理方框图如下:检控法2PSK 信号的解调器原理方框图(4) 2DPSK 是利用前后相邻码元的载波相对相位变化传递数字信息,所以又称相对相移键控。
利用MATLAB实现DPSK调制及解调
利用MATLAB实现DPSK调制及解调
DPSK(Differential Phase Shift Keying)调制是一种数字信息传输调制方式。
它采用相位差的改变来表示数字信息,具有抗噪声和波动的能力,因此在数字通信领域得到了广泛的应用。
MATLAB是一种适合数字信号处理的工具,可以有效地实现DPSK调制及解调。
以下是具体的实现步骤:
DPSK调制
1. 生成数字信息比特流,转换为1和-1形式。
2. 将比特流进行差分编码得到差分比特流。
3. 将差分比特流分组,每组2个比特。
4. 根据相邻两个比特的差异,确定相位差。
差分比特流为00或11时,相位差为0;差分比特流为01或10时,相位差为π。
5. 根据相位差,生成相位进行调制得到调制信号。
可以使用sinc函数或高斯函数对信号进行脉冲整形。
DPSK解调
1. 将DPSK调制后的信号送入相干解调器。
2. 使用带通滤波器去除高频噪声。
3. 再次进行相干解调,得到调制信号。
4. 对调制信号进行差分解码还原差分比特流。
5. 对差分比特流进行译码得到数字信息比特流。
利用MATLAB实现DPSK调制及解调的代码可在Matlab官网上找到并学习使用。
通信原理课程设计(基于MATLAB的PSK,DPSK仿真)
通讯原理课程设计报告题目鉴于Matlab的2PSK,2DPSK仿真学院电子信息工程学院专业学生姓名学号年级指导教师职称讲师2013年12月20日设计报成功绩(依据优、良、中、及格、不及格评定)指导教师考语:指导教师(署名)年月日说明:指导教师评分后,设计报告交院实验室保留。
鉴于 Matlab 的 2PSK,2DPSK仿真专业:学号:学生:指导老师:纲要:现代通讯系统要求通讯距离远、通讯容量大、传输质量好,作为其重点技术之一的调制技术向来是研究的一个重要方向。
本设计主要表达了数字信号的调制方式,介绍了2PSK数字调制方式的基来源理,功率谱密度,并运用MATLAB软件对数字调制方式2PSK进行了编程仿真切现,在MATLAB 平台上成立2PSK和 2DPSK调制技术的仿真模型。
进一步学习了MATLAB编程软件,将 MATLAB与通讯系统中数字调制知识联系起来,为此后在通讯领域学习和研究打下了基础在计算机上,运用MATLAB 软件来实现对数字信号调制技术的仿真。
重点词:数字调制与解调;MA TLAB ; 2PSK; 2DPSK ;目录第 1 章绪论 (1)1.1 调制方式 (1)1.2 设计要求 (1)设计内容 (1)设计仪器 (1)第 2 章 2PSK,2DPSK原理 (2)2.1 2PSK 原理 (2)2PSK 基来源理 . (2)2PSK 调制原理 . (2)2PSK 解调原理 . (3)2.2 2DPSK 原理 (4)2DPSK 基来源理 . (4)2DPSK 调制原理 . (5)2DPSK 解调原理 . (6)第 3 章实验过程 (8)3.1 2PSK 仿真部分 (8)2PSK 仿真图 . (8)2PSK 模块的参数设置: . (8)3.2 2DPSK 仿真部分 (9)2DPSK 仿真图 . (9)2DPSK 模块的参数设置: . (10)第 4 章仿真结果 (15)4.1 2PSK 仿真结果 (15)4.2 2DPSK 仿真结果 (15)总结 . (16)参照文件 . (17)道谢 . (18)第1章绪论1.1调制方式数字通讯系统 ,按调制方式能够分为基带传输和带通传输。
实验四--基于simulink的2PSK、2DPSK数字调制与解调的仿真
3、锻炼学生分析问题和解决问题的能力
2、实验原理
1.12PSK调制原理
数字调相:如果两个频率相同的载波同时开始振荡,这两个频率同时达到正最大值,同时达到零值,同时达到负最大值,它们应处于"同相"状态;如果其中一个开始得迟了一点,就可能不相同了。如果一个达到正最大值时,另一个达到负最大值,则称为"反相"。一般把信号振荡一次(一周)作为360度。如果一个波比另一个波相差半个周期,我们说两个波的相位差180度,也就是反相。当传输数字信号时,"1"码控制发0度相位,"0"码控制发180度相位。载波的初始相位就有了移动,也就带上了信息。
模拟调相法原理方框图如下图所示,极性变器将输入的二进制单极性码转换成双极性不归零码,然后与载波直接相乘,以实现2PSK
模拟调相法
键控法原理方框图如下图所示,用数字基带信号s(t)控制开关电路,以选择不同相位的载波输出。此时s(t)通常是单极性的,当s(t)=0时,输出e2PSK(t)=cosωct;当s(t)=1时,输出e2PSK(t)=-cosωct。
码型变化器参数设置:采用1变0不变调制,故极性设为“Negative”.
多路选择器参数设置:当二进制序列大于0时,输出第一路信号;当二进制序列小于0时,输出第二路信号。
高斯噪声的设置:信噪比90
带通滤波器参数:载波频率为4HZ,而基带号带宽为1HZ,考滤到滤波器的边沿缓降,故设置为2~7HZ。
低通滤波器参数设置:截止频率为1HZ,二进制序列的带宽为1HZ,故取1HZ。
差分编码为 (2)
差分译码为 (3)
再将相对码 进行2PSK调制,则所得到的即是2DPSK已调信号,其在任一码元时间 内的表达式为
毕业设计(最终版)-2dpsk调制解调仿真系统设计[管理资料]
论文题目:2DPSK调制解调仿真系统设计专业:通信技术学号: 0807040211学生:苏瑜签名:_____指导教师:张鸣签名:_____摘要二进制差分相移键控简称二相相对调相,记作2DPSK。
它是数据通信中最常用的一种调制方式,这种方式的优点是简单,易于实现。
与2PSK的波形不同,2DPSK 波形的同一相位并不对应相同的数字信息符号,而前后码元的相对相位才唯一确定信息符号。
为了使数字信号在带通信道中传输,必须用数字基带信号对载波进行调制,以使信号与信道特性相匹配。
调制解调技术是实现现代通信的重要手段,研究数字通信调制解调理论,提供有效的调制方式,有着重要意义。
Matlab中的通信工具箱可以用来进行通信领域的研究、开发、系统设计和仿真。
本文主要研究了利用Matlab Simulink软件对2DPSK调制解调系统的设计。
首先介绍了2DPSK调制解调的基本原理和Matlab Simulink仿真软件。
然后根据2DPSK的调制原理及解调原理进行了仿真模型搭建,并分析了仿真结果。
最后利用Matlab Simulink 中的DBPSK Modulator Baseband和DBPSK Demodulator Baseband模块实现了2DPSK调制解调系统设计。
【关键词】Matlab Simulink 2DPSK 调制解调【论文类型】设计型Title: 2DPSK of Modulation and Demodulation Simulation System Design Major: Communication Technology Number: 0807040211 Name:Su Yu Signature:Supervisor:Zhang Ming Signature:ABSTRACTBinary differential phase shift keying modulation is relatively short two-phase, denoted by 2DPSK. It is one of the most commonly used modulation in data communications, this approach has the advantage of simple, easy to implement. 2DPSK waveform is different from 2PSK waveform. 2DPSK waveform does not correspond to the same phase of the same number of information symbols, and the relative phase of the before and after symbol determine the only information symbols. In order to make the digital signal transmitted in a communication channel,we must use digital baseband signal to modulate the carriert, which can make the signal match the channel characteristics. Modulation and demodulation technology is an important means to achieve modern communications, modulation and demodulation of digital communications theory,and it can provide an effective way of great significance. Matlab toolbox of communication in the field of communications can be used for research, development, system design and simulationThis paper studies the use of Matlab Simulink software to design 2DPSK modem system. First it introduces the basic principles of 2DPSK modulation and demodulation and Matlab Simulink simulation software. Then according to the principle of 2DPSK modulation and demodulation, it builds the simulation model and analysis of simulation results. Finally, in the 2DPSK modulator and 2DPSK demodulator module of Matlab Simulink,it implements the design of 2DPSK modem system.【Key words】Matlab Simulink 2DPSK Modulate Demodulate 【Type of Thesis】Design前言本文主要是熟悉Matlab的运用,了解2DPSK的原理及其原理图,利用Matlab/Simulink仿真软件,设计一个2DPSK调制与解调系统,并得出调制后的信号,用示波器对调制前与解调后的波形进行比较,最后进行性能分析。
2DPSK差分相干解调器Simulink仿真
课程设计班级:通信09—2班*名:**学号:**********指导教师:***成绩:电子与信息工程学院通信工程系目录(contents)摘要-----------------------------------------------------------02 第一章:绪论---------------------------------------------------02 第二章:2DPSK信号调制原理--------------------- ----------------02 第三章:2DPSK调制和解调原理------------------------------------03 3.1 2DPSK信号调制原理------------------------------------------04 3.2 2DPSK信号解调原理------------------------------------------04 第四章:Simulink模型建立---------------------------------------05 4.1 2DPSK调制模块-----------------------------------------------05 4.2 相干解调----------------------------------------------------07 4.3差分相干解调-------------------------------------------------09 第五章:心得体会-------------------------------------------------11 参考文献--------------------------------------------------------11基于simulink的2DPSK差分相干解调器设计Abstract:using communication principle in the knowledge and understanding of matlab simulink, use simulin function design 2 DPSK difference coherent modem and simulation.摘要:利用在通信原理中所学知识和对matlab中simulink的了解,利用simulin功能设计2DPSK差分相干解调器和进行仿真。
基于Simulink的FSK、DPSK体系系统仿真
基于Simulink的FSK、DPSK系统仿真
一、实验目的
理解数字调制解调系统基本原理,并熟悉Matlab:Simulink通信软件包的使用。
二、实验内容
Simulink通信软件包模拟FSK、DPSK系统的调制和解调。
三、实验波形显示
FSK:
系统Sample Time:0.01s
输入信号
上为转化为二进制的信号,下为量化的信号
上为数字信号,两种频率不同的调制信号,w1= 1000rad/s w2=10000 rad/s
信号。
判决前信号及输入数字信号。
DPSK:系统Sample Time 0.5s
上为输入信号,下为0差分编码的结果
调制信号,ω=200*pi rad/s 两调制信号相位差为pi
调制信号和延时后信号
上为解调信号,下为调制信号
上为输入数字信号,下为解调后判决获得的信号。
四、实验结果分析
实验结果与书上原理相符,获得了预期的结果。
可以比较发现,DPSK系统更为稳定,FSK
要求调制频率较高,这样减少误码率。
DPSK相位差的方法,ω=200*pi的效果已经很好了五、实验小结
此实验加强了对通信系统调制解调的理解,当然还有许多没有搞清的问题,例如采样量化的方面,有待以后解决。
六、引用资料
1、张志涌、徐彦琴.《Matlab 教程—基于6.X版本》北京:北京航空航天出版社
2、孙屹、吴磊.《Simulink通信仿真开发手册》湖南:国防工业出版社
3、孙学军、王秉钧.《通信原理》北京:电子工业出版社。
基于MATLAB的DPSK通信系统仿真设计
〈通信原理》课程实验报告基于MATLAB的DPSK通信系统仿真院系: 班级: 学号: 姓名: 教师:2012年12月7日目录1设计题目....................... 3 .................................................. 2设计要求....................... 3 .................................................. 3设计原理....................... 3 ..................................................3.1 DPSK 调希9 ....................................... 3 ....................................................3.2 DPSK 解调......................... 4 .................................................... 4设计流程....................... 5 .................................................. 5代码实现....................... 5 ...........................................................5.1基带信号的获取......................... 5 ...................................................5.2 差分编码.......................... 6 ......................................................5.3 DPSK 调希9 ..................................... 7 ......................................................5.4 AWGN 信道.......................... 8 ...................................................5.5 限噪处理.......................... 8 ......................................................5.6差分相干解调........................... 9 ...................................................5.4 AWGN 信道......................... 11 ...................................................5.7功率谱密度图的绘制6测试结果...................... •........6.1调制解调正确性的测试 (12)6.2系统抗噪性能的测试........................ 14 -7总结........................ 17 ................... 8心得体会...................... 17 ............... 附:完整代码................... 18 ••-1设计题目用MATLAB仿真DPSK通信系统2设计要求基本参数如下:1) Fc=1800Hz2) Rb=1200bps3) 考虑加入AWGN信道,Eb/N0=20dB4)自定义一小段信息序列结果要求:1)绘制各阶段信号波形2)绘制信号的PSD3设计原理二进制差分相移键控常简称为二相相对调相,记为2DPSK。
关于BPSK、MPSK、MDPSK及MATLAB的操作说明
关于BPSK、MPSK、MDPSK及MATLAB的操作说明一、示范程序的运行方法本课程提供的示范程序包括三个方面,分别是BPSK、M(M=4、8)进制PSK和DPSK 通信系统的仿真。
前者由c++实现,后两者用Matlab语言编写。
示范程序已经上载于我院的ftp服务器上,同学们可以从“ftp://202.112.108.122/”的“通信原理与计算机仿真”文件夹中下载到本地机。
1、BPSK示范程序的运行直接运行文件夹下的可执行文件“A wgnbpsk.exe”,观察高斯白噪声的BPSK仿真;直接运行文件夹下的可执行文件“Rayleigh.exe”,观察噪声为瑞利分布的BPSK仿真。
2、M(M=4、8)进制PSK和DPSK示范程序的运行一个完整的matlab仿真程序一般由多个文件组成(即主函数和子函数)。
另外,Matlab 强调了“当前工作目录”,执行matlab仿真程序的过程,简单的说,就是在“Matlab的当前工作目录为该程序的主函数所在文件夹”的前提下,键入主函数名,然后回车运行。
一系列的操作类似于DOS命令。
举例来说,我们要运行四进制PSK仿真程序,首先把示范程序中的“M-psk”文件夹拷至Matlab的根目录下(假设其路径为“c:/matlab5.3/M-psk”),观察四进制PSK仿真程序,它的主函数为“psk.m”,路径是“c:/matlab5.3/M-psk/4psk/psk.m”。
运行此仿真程序的过程为:a.用“pwd”命令显示当前目录,其默认值一般为“c:/matlab5.3/work”;b.用“cd”命令改变当前工作目录,使其指向“c:/matlab5.3/M-psk/4psk”,即主函数“psk.m”所在文件夹,命令为“cd c:/matlab5.3/M-psk/4psk”;c.键入“psk”运行仿真程序。
以此类推,如果我们要运行八进制DPSK仿真程序(假设事先已将文件夹“M-dpsk”复制到路径“c:/matlab5.3/work/M-dpsk”下,则此时主函数“dpsk8.m”所在路径是“c:/matlab5.3/work/M-dpsk/8dpsk”),故只需键入命令“cd c:/matlab5.3/work/M-dpsk /8dpsk”,回车后运行“dpsk8”即可。
基于MATLAB的2DPSK调制和解调仿真
基于MATLAB的2DPSK调制和解调仿真1. MATLAB/SIMULINK简介美国Mathworks公司于1967年推出了矩阵实验室“Matrix Laboratory”(缩写为Matlab)这就是Matlab最早的雏形。
开发的最早的目的是帮助学校的老师和学生更好的授课和学习。
从Matlab诞生开始,由于其高度的集成性及应用的方便性,在高校中受到了极大的欢迎。
由于它使用方便,能非常快的实现科研人员的设想,极大的节约了科研人员的时间,受到了大多数科研人员的支持,经过一代代人的努力,目前已发展到了7.X版本。
Simulink是MATLAB中的一种可视化仿真工具,也是目前在动态系统的建模和仿真等方面应用最广泛的工具之一。
确切的说,Simulink是一个用来对动态系统进行建模、仿真和分析的软件包,它支持线性和非线性系统,连续、离散时间模型,或者是两者的混合。
系统还可以使多种采样频率的系统,而且系统可以是多进程的。
在Simulink环境中,它为用户提供了方框图进行建模的图形接口,采用这种结构画模型图就如同用手在纸上画模型一样自如、方便,故用户只需进行简单的点击和拖动就能完成建模,并可直接进行系统的仿真,快速的得到仿真结果。
2.二进制差分相位键控(2DPSK)在2PSK信号中,信号相位的变化是以未调正弦载波的相位作为参考,用载波相位的绝对数值表示数字信息的,所以称为绝对移相.由图 2 - 14 所示2PSK 信号的解调波形可以看出, 由于相干载波恢复中载波相位的180°相位模糊,导致解调出的二进制基带信号出现反向现象,从而难以实际应用. 为了解决2PSK信号解调过程的反向工作问题, 提出了二进制差分相位键控(2DPSK). 2DPSK方式是用前后相邻码元的载波相对相位变化来表示数字信息.假设前后相邻码元的载波相位差为Δφ,可定义一种数字信息与Δφ之间的关系为则一组二进制数字信息与其对应的 2DPSK 信号的载波相位关系如下所示:二进制数字信息: 1 1 0 1 0 0 1 1 1 02DPSK信号相位: 0 π 0 0 πππ 0 π 0 0或π 0 ππ 0 0 0 π 0 ππ数字信息与Δφ之间的关系也可以定义为可以看出,2DPSK信号的实现方法可以采用:首先对二进制数字基带信号进行差分编码,将绝对码表示二进制信息变换为用相对码表示二进制信息,然后再进行绝对调相,从而产生二进制差分相位键控信号。
基于MATLAB实现的PSK和DPSK信号系统仿真设计
文档从互联网中收集,已重新修正排版,word格式支持编辑,如有帮助欢迎下载支持。
目录摘要 ................................................................................................................... 错误!未定义书签。
Abstract............................................................................................................. 错误!未定义书签。
1.设计目的与要求 ............................................................................................ 错误!未定义书签。
2.方案的选择 .................................................................................................... 错误!未定义书签。
2.1调制部分............................................................................................. 错误!未定义书签。
2.2解调部分............................................................................................. 错误!未定义书签。
3.单元电路原理和设计 .................................................................................... 错误!未定义书签。
基于MATLAB的DPSK系统的仿真
呼伦贝尔学院论文开题报告专业:电子信息工程论文题目:基于MATLAB/Simulink的DPSK系统的仿真学生姓名:郭东芹学号:201208102089指导教师:刘丽开题报告日期:2015年11月13日一、选题意义、国内外发展现状与发展趋势1.课题的意义为了使数字信号在带通信道中传输,必须用数字基带信号对载波进行调制,以使信号与信道特性相匹配。
调制解调技术是实现现代通信的重要手段,研究数字通信调制解调理论,提供有效的调制方式,有着重要意义。
二进制差分相移键控又称相对相移键控,记作2DPSK。
它不是利用载波相位的绝对数值传送数字信息,而是用前后码元的相对载波相位值传送数字信息。
所谓相对载波相位是指本码元初相与前一码元初相之差。
2DPSK信号的解调有两种方式,一种是差分相干解调,另一种是相干解调-码变换法,后者又称为极性比较-码变换法。
与2PSK的波形不同,2DPSK波形的同一相位并不对应相同的数字信息符号,而前后码元的相对相位才唯一确定信息符号。
这说明解调2DPSK信号时,并不依赖于某一固定的载波相位参考值,只要前后码元的相对相位关系不破坏,则鉴别这个相位关系就可正确恢复数字信息。
这就避免了 2PSK方式中的“倒π”现象发生。
2.国内外现状20世纪70年代,美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担,用FORTRAN编写了最早的MATLAB。
1984年由Little、Moler、Steve Bangert合作成立了的MathWorks公司正式把MATLAB推向市场。
到20世纪90年代,MATLAB已成为国际控制界的标准计算软件。
随着计算机技术的发展MATLAB已经被应用到各个领域。
在许多诸如控制论、时间序列分析、系统仿真、图像信号处理等方面产生了大量的矩阵及其相应的计算问题.自己去编写大量的繁复的计算程序,不仅会消耗大量的时间和精力,减缓工作进程,而且往往质量不高。
matlab设计2DPSK信号调制 (1)
2DPSK调制与解调系统的仿真1、 2DPSK基本原理1.1 2DPSK信号原理2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。
现假设用Φ表示本码元初相与前一码元初相之差,并规定:Φ=0表示0码,Φ=π表示1码。
则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的,在接收端只能采用相干解调,它的时域波形图如图2.1所示。
图1.1 2DPSK信号在这种绝对移相方式中,发送端是采用某一个相位作为基准,所以在系统接收端也必须采用相同的基准相位。
如果基准相位发生变化,则在接收端回复的信号将与发送的数字信息完全相反。
所以在实际过程中一般不采用绝对移相方式,而采用相对移相方式。
定义∆Φ为本码元初相与前一码元初相之差,假设:∆Φ=0→数字信息“0”;∆Φ=π→数字信息“1”。
则数字信息序列与2DPSK信号的码元相位关系可举例表示如下:数字信息: 1 0 1 1 0 1 1 1 0 1DPSK信号相位:π 0 0 π 0 0 π 0 π π 0或:0 π π 0 π π 0 π 0 0 π1.2 2DPSK信号的调制原理一般来说,2DPSK信号有两种调试方法,即模拟调制法和键控法。
2DPSK信号的的模拟调制法框图如图1.2.1所示,其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。
图1.2.1 模拟调制法2DPSK 信号的的键控调制法框图如图1.2.2所示,其中码变换的过程为将输入的基带信号差分,即变为它的相对码。
选相开关作用为当输入为数字信息“0” 时接相位0,当输入数字信息为“1”时接pi 。
图1.2.2 键控法调制原理图1.3 2DPSK 信号的解调原理2DPSK 信号最常用的解调方法有两种,一种是极性比较和码变换法,另一种是差分相干解调法。
1.3.1 2DPSK 信号解调的极性比较法它的原理是2DPSK 信号先经过带通滤波器,去除调制信号频带以外的在信道中混入的噪声,再与本地载波相乘,去掉调制信号中的载波成分,再经过低通滤波器去除高频成分,得到包含基带信号的低频信号,将其送入抽样判决器中进行抽样判决的到基带信号的差分码,再经过逆差分器,就得到了基带信号。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
利用Matlab simulink实现DPSK通信刘易 2011001010010工作简介此次工作完成了以下几个任务:1.Simulink的DPSK实现:利用MATLAB simulink模块实现了自定义的DPSK信号的调制和解调,其中通信是利用simulink的UDP send和UDP receive来实现的,将两台电脑直接用网线连接起来,便可在两边分别运行UDP send和UDP receive实现通信,其中目的IP用“255.255.255.255”以进行广播。
2.考虑噪声的通信:在考虑噪声的情况下传输自定义信号,并成功复原原信号。
3.语音信号的传输:在实现任意信号的通信后,用一段语音信号进行DPSK通信。
先对语音信号进行量化,然后通过DPSK调制,再用相干解调得到原传输信号,最后复原语音信号并进行播放。
信号的传输是在有噪声的情况下进行的。
下面分别介绍上述工作:Simulink的DPSK实现1.1 首先展示simulink的发送方的模块:1.2 定义需要传输的信号:{1 0 1 0 1 0 1 0 1 1}1.3 生成需要传输的信号:题目中要求:Tc=1800HZ,Rb=1200bps,但是为了模拟显示得更方便直观,我仅保持了这两个量的比例,将两个量改为了:Tc=3HZ,Rb=2bps。
(消息信号生成的代码参考附录1:产生消息信号)1.4 将消息信号与载波信号相乘得到需要传输的信号:需要指出,我只定义了10个需要传输的信号,即需要传输5秒,但我都是以6秒来发送后的,后面的5-6秒的两个信号用0填充了,其原因是后面接收方的低通滤波器有延迟,需要留足够的空间来显示完整信号。
下面给出传输信号的功率谱密度,这里的计算是利用离散信号做FFT再除以信号长度得到的(由于主要集中在某一范围因此仅给出部分范围的频谱):1.5然后UDP send模块需要将小心信号发送到接收端。
2.1 接收方的Simulink模块2.2接收方的接收信号(由scope6显示):从接收模块可以看到,在直接接收的信号到scope6之间有一个我自己定义的模块。
其原因是在直接接收到的信号中,有许多不相关的信号,且一般都比需要的信号要小。
如果将直接接收的信号存下来观察,可以得到的是255×1×6001的信号,但我们所需要的只是1×6001的信号,因此我添加了一个模块来滤去其他无关信号。
(该模块代码参考附录2:接收信号的提取)下面是处理过后的接收信号:2.3 将接收信号再乘一个与载波信号同频率的正弦信号(由scope8显示):2.4 将信号通过低通滤波器,得到如下信号(由scope5显示):可以看到这个低通滤波器在时间上延迟了大概0.75秒,这就是为什么一开始要以6秒的长度来传输数据。
2.5 接着通过一个自定义的模块,进行抽样判决,大于0记为1,小于0为0,这样便复原了原信号:(抽样判决模块代码参考附录3:抽样判决)考虑时延后,抽样判决输出的序列为{1 0 1 0 1 0 1 1},和原序列相同。
以上便完成了DPSK调制通信的全部过程,成功地复原了原信号。
其实这个系统可以跑在两台电脑上,通过网线连接,用UDP通信,但在实行的时候出现了些问题,我在报告最后再来说说两台电脑的通信问题。
以上虽然成功地完成了DPSK的通信,但是很重要的是没有考虑噪声,在实际生活中是不可能不考虑噪声的,因此下面在考虑噪声的情况下实现DPSK通信。
考虑噪声的通信只需在发送方的模块中添加一个产生噪声的模块即可模拟有噪声的通信:下面依次给出接收方收到的信号以及处理之后的信号波形:1.接收的信号:2.乘cos(2πfct)后的信号:3.经过低通滤波器后的信号:4.抽样判决输出:看到这个信号,我的确有点吃惊,滤波出来的波形居然会这么好,但仔细想想,我发送的信号的波特率很低,也就很靠近0频,而带宽也很小。
而且我把滤波器的通带调得很小,因此噪声的大部分能量都被滤去了。
因此结果能够达到这么好。
语音信号的传输语音信号这个部分我就没用使用simulink来模拟了,主要原因是语音信号较长,而且一面simulink模块中,运行低通滤波器需要很大消耗,因此我考虑编程来模拟语音信号的量化,调制,以及解调。
其中的滤波器我需要重新设计,以减少程序运行时间。
下面是一段自己录制的语音信号(大概1秒左右,内容是“你好”):当前语音信号的最大值为0.1290,最小值为-0.1333。
因此我将其进行从-0.15到0.15的均匀量化,并采用格雷码进行编码(量化代码参考附录4:语音信号的量化,以及格雷码的生成代码,参考附录5:格雷码的产生)。
之所以采用格雷码进行编码,是因为如果接收产生误码,可以不至于声音差得太多。
下面是量化后的语音信号,是采用8位二进制数进行量化,即将-0.15到0.15均匀分成了256份。
试听以后效果不错,声音基本没有失真。
通过以上处理可以得到一串语音的二进制码,采用与上文中相同的DPSK传输方法,将语音信号调制并传送出去(代码参考附录6:语音信号的调制)。
由于信号量化后序列较长,因此仅给出一小部分波形:给信号添加噪声,以模拟真实的噪声信道,添加噪声后信号的部分波形如下:给信号乘上与载波同频的正弦信号,得到如下波形:这里自行设计了滤波器,采用切比雪夫2型滤波器,阶数为9,阻带为80dB,阻带频为0.07(代码参考附录7:切比雪夫低通滤波器),通过低通滤波器后得到如下波形(对应序列为000000100):由于有噪声影响,滤除的效果并不算特别好,但仍然能进行较准确的抽样判决了。
接下来就通过抽样判决的函数输出原量化后的语音信号的序列(代码参考附录8:抽样判决),接着就根据格雷码对信号进行复原了(代码参考附录9:语音信号的复原)。
下面给出抽样判决后的部分序列,和格雷码复原后的信号波形:通过上面的操作,便完成了语音信号的量化、调制、加噪、解调,成功的,模拟了语音信号利用DPSK的通信过程,最后播放得到的语音信号,声音辨识度很高,基本没有失真。
最后的一点话:其实本次工作一开始我是希望能够在两台电脑上实现通信的,由于没有串口线,我想我们都有网线,就可以直接把两台电脑用网线相连,利用运输层的协议就可以实现通信了,如下图一样:之后我就分别在两台电脑上运行下面两个Simulink模块:可以说,和我预想的一样,可以接收到发送的信号,但是接收到的信号却让人琢磨不透,如下图(是多维的信号,第一个模块接收到的信号,其余三个是部分信号的部分波形):和在本机模拟的一样,接收到的都是255×1×6001的信号,但是不同的是本机接收到的信号中,只需通过能量最大就可以找出需要的信号,而在两台电脑上通信时,接收方接收到的就是一些有规律但完全看不出什么的信号,我在这个问题上困扰了很久,最终也没有解决。
后来又想通过互联网来传送,但是在simulink的UDP模块中设置好目的IP,和目的端口号后,接收方却没有收到任何信号,也没有解决这个问题。
于是最后只好折中,采用本机模拟两电脑基于UDP的通信。
但好在的是,除了上面的问题,其他方面虽然遇到的问题,都被一一解决了,比如滤波器始终不理想的问题。
总的来说我觉得这次工作对我帮助还是挺大的,至少我在simulink仿真方面有了很大进步。
程序附录:附录1:产生消息信号(simulink模块)function y = fcn2(t)a=[1 -1 1 -1 1 -1 1 -1 1 1 0 0];if t==6y=a(10);return;endy=a(floor(t/0.5)+1);附录2:接收信号的提取(simulink模块)function y=fun2(u)[n,m,c]=size(u);sum1=0;index=1;for i=1:nsum2=sum(abs(u(i,:,:)));if sum2>sum1sum1=sum2;index=i;endendy=zeros(1,c);y(1:c)=u(index,1,1:c);附录3:抽样判决(simulink模块)function y=funpanjue(u)if u>=0y=1;elsey=0;end附录4:语音信号的量化function [out_voice_seq,out_voice_signal]=voice_signal_quantization(voice_seq) n=length(voice_seq);m=8;out_voice_seq=zeros(n*m,1);out_voice_signal=zeros(n,1);interval=0.3/(2^m);code1=grey_code(m);for i=1:nfor j=1:2^mif voice_seq(i)<=(-0.15+j*interval)out_voice_seq((i-1)*m+1:i*m)=code1(:,j);out_voice_signal(i)=-0.15+j*interval-interval/2;break;endendend附录5:格雷码的产生function code=grey_code(m)code=zeros(m,2^m);code(1,2)=1;for i=2:mcode(i,2^(i-1)+1:2^i)=1;for j=1:2^(i-1)code(1:i-1,2^(i-1)+j)=code(1:i-1,2^(i-1)-j+1);endendend附录6:语音信号的调制function [modulated_signal,cos_sample]=voice_signal_modulation(voice_signal) n=length(voice_signal);modulated_signal=zeros(n*100,1);cos_sample=zeros(n*100,1);for i=1:nfor j=1:100t=(i-1)+(j-1)*1/100;cos_sample((i-1)*100+j)=cos(3*pi*t);modulated_signal((i-1)*100+j)=cos(3*pi*t)*voice_signal(i);endendend附录7:切比雪夫低通滤波器function output_signal=lowpassdesign(input_signal,p1,p2,p3) [A,B]=cheby2(p1,p2,p3);output_signal=filter(A,B,input_signal);end附录8:抽样判决function output_seq=sample_judge(input)n=length(input);m=n/100;output_seq=zeros(m,1);for i=1:mif input(i*100)>=0output_seq(i)=1;elseoutput_seq(i)=0;endend附录9:语音信号的复原function output=signal_restoration(input)n=length(input);m=n/8;greycode=grey_code(8);output=zeros(m,1);for i=1:mfor j=1:2^8if sum(input((i-1)*8+1:i*8)==greycode(:,j))==8output(i)=-0.15+j*0.3/(2^8)-0.3/(2^9);break;endendendend。