[交通运输]matlab仿真实验指导书-附程序.doc

matlab与通信仿真实验指导书(下)通信基础教研室上课时间： 学年第 学期 系 部： 班 级： 姓 名： 班内序号： 指导教师： 实验课程成绩： 目 录实验一 matlab基础实验1实验一成绩 实验二 绘图和确知信号分析实验8实验二成绩 实验三 随机信号与数字基带实验15实验三成绩 实验四 模拟调制实验24实验四成绩 实验五 模拟信号数字传输实验（一）32实验五成绩 实验六 模拟信号数字传输实验（二）41实验六成绩 实验七 数字频带传输系统实验47实验七成绩 实验八 通信系统仿真综合实验57实验八成绩 matlab与通信仿真实验一 matlab基础实验一、实验目的l 了解matlab 程序设计语言的基本特点，熟悉matlab软件运行环境l 掌握创建、保存、打开m文件及函数的方法l 掌握变量等有关概念，具备初步的将一般数学问题转化为对应的计算机模型并进行处理 的能力二、实验内容及步骤1.在command window里面计算；，计算：；，求；，输入复数矩阵；2.建立.m文件，用for循环语句生成1010的矩阵a：,将a矩阵进行水平和垂直翻转得到矩阵b和c。将a矩阵的前5行，5列变成0并赋值给d。3.建立.m文件，随机产生一个5050的矩阵，元素值为从0到255，要求用0和255对该矩阵进行标记，元素值大于128的标记为255，元素值小于128的标记为0。4.产生一个均值为2.4方差为0.2大小为34的随机矩阵。5.（选做）编写函数使用0.618搜索法（近似黄金分割法）求给定函数的极值：搜索法求解的基本过程：给出a,b，使得t*在a,b中。a,b称为搜索区间。迭代缩短a,b的长度。当a,b的长度小于某个预设的值，或者导数的绝对值小于某个预设的正数，则迭代终止。以函数，作为处理对象，其中搜索区间定为0,3,精度定为0.5。 四、实验报告：（调试好的程序，实验结果及分析）1.在command window里面计算；解：32；解：2.7384e-016，计算：；解：，求：解：，输入复数矩阵；2.建立.m文件，用for循环语句生成1010的矩阵a：,将a矩阵进行水平和垂直翻转得到矩阵b和c。将a矩阵的前5行，5列变成0并赋值给d。解：% 实验一的第二个实验 exp0102.mclear allx=1:10; % 生成1*50的矢量a=x; % 给a矩阵赋第一行元素for i=1:9 a=a;x+i; % 生成余下的49行a矩阵的元素，当前行元素是前一行的对应元素加1enda% 水平翻转ab=fliplr(a) % 本句和下一句都可实现水平翻转b=flipdim(a,2)% 垂直翻转ac=flipud(a) % 本句和下一句都可实现垂直翻转c=flipdim(a,1)% 将a矩阵的前10行，10列变成0并负值给da(1:5,1:5)=0d=a学生程序一：% 学生程序a=zeros(10);for i=1:10 for j=1:10 if i=j a(i,j)=2*i-1; else a(i,j)=i+j-1; end endenda学生程序二：% 学生程序a=zeros(10);for i=1:10 for j=1:10 a(i,j)= i+j-1; endenda学生程序三：% 学生程序b=rot90(a,3)%水平翻转b=rot90(a,1)%垂直翻转学生程序四：a=ones(10);n=linspace(1,10,10)；for m=1:10a(m,j=n+linspace(1,1,10);end% 学生程序b=rot90(a,3)%水平翻转b=rot90(a,1)%垂直翻转3.建立.m文件，随机产生一个1010的矩阵，元素值为从0到255，要求用0和255对该矩阵进行标记，元素值大于128的标记为255，元素值小于128的标记为0。解：% 实验一的第三个实验 exp0103.m% 建立.m文件，随机产生一个5050的矩阵，元素值为从0到255，要求用0和255对该矩阵进行标记，% 元素值小于128的标记为0。clear alln=5;%a=fix(90-10+1)*rand(n)+10)a=fix(255)*rand(5) % 产生产生一个5050的，元素值为从0到255随机矩阵a=128*ones(n,n) % 产生产生一个5050的，元素值全为128随机矩阵b=aa % 产生产生一个5050的矩阵，元素值元素值大于128的标记为1，元素值小于128的标记为0c=a128)% b=find(a128 a(i,j)=255; elseif a(i,j)128 a(i,j)=0; end endendaa = 51.7051 106.7556 128.2173 177.9641 168.3580 179.1987 249.8355 168.6679 57.6172 144.7963 50.6740 215.7865 180.9152 96.4851 87.2025 139.3756 69.2191 72.5242 147.8508 202.5237 153.9671 133.9139 109.3676 219.3030 73.8801 113.4445 64.3440 119.6522 193.8931 15.0916 69.4079 51.6751 77.6774 217.6821 87.0044 177.1146 223.3142 16.5192 135.1049 153.7316 50.6976 171.3951 48.3617 151.3585 136.1901 158.4341 188.0130 252.0254 163.3343 12.8185 3.8949 213.7202 49.3249 126.6209 185.4139 202.6794 34.8123 148.6119 53.3127 105.9206 190.4303 5.0081 173.9669 229.4411 78.8690 243.9951 2.9980 107.9915 96.8537 77.7747 113.4996 173.7257 77.2049 209.5154 213.8165 133.2605 227.9440 131.4555 199.7488 222.9636 237.6127 96.7677 138.1268 164.4521 144.8585 224.4363 50.7802 85.1576 173.6157 3.8274 118.8286 212.1080 38.4726 208.5835 94.4555 44.1038 76.1744 110.3912 117.5793 195.8273a = 0 0 255 255 255 255 255 255 0 255 0 255 255 0 0 255 0 0 255 255 255 255 0 255 0 0 0 0 255 0 0 0 0 255 0 255 255 0 255 255 0 255 0 255 255 255 255 255 255 0 0 255 0 0 255 255 0 255 0 0 255 0 255 255 0 255 0 0 0 0 0 255 0 255 255 255 255 255 255 255 255 0 255 255 255 255 0 0 255 0 0 255 0 255 0 0 0 0 0 2554.产生一个均值为2.4方差为0.2大小为34的随机矩阵。解：% 实验一的第四个实验 exp0104.m% 产生一个均值为2.4方差为0.2大小为34的随机矩阵。% clear allclose alla=2.4+sqrt(0.2)*randn(1000,4)mean(a)var(a)mean（a）= 2.3807 2.4204 2.4079 2.4118var(a)=0.1780 0.2127 0.2100 0.19655.（选做）编写函数使用0.618搜索法（近似黄金分割法）求给定函数的极值：搜索法求解的基本过程：给出a,b，使得t*在a,b中。a,b称为搜索区间。迭代缩短a,b的长度。当a,b的长度小于某个预设的值，或者导数的绝对值小于某个预设的正数，则迭代终止。% 学生程序function f=factor(a,b)t1=a+0.382*(b-a);t2=a+0.618*(b-a);if t13-2*t1+1=t23-2*t2+1 if t2-a=0.5 disp(t1); %break; else factor(a,t2); endelse if b-t1=0.5 disp(t2); %break; else factor(t1,b); endend三、实验问题：1.第一个实验的a矩阵是一个近似奇异的矩阵，因此导致其求逆不确定2.第二个实验，将50*50改为10*103.实验二 绘图和确知信号分析实验一、实验目的l 掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理l 理解周期信号的傅里叶级数展开的物理意义l 掌握信号的傅里叶变换及其反变换二、实验原理1.周期信号的傅里叶级数若一周期信号，其中为整数，成为信号的周期。若周期信号在一个周期内可积，则可通过傅里叶级数对该信号进行展开。其傅里叶展开式如下：，其中，为信号最小周期；为信号的基波；为傅里叶展开系数，其物理意义为频率分量的幅度和相位。2.信号的傅里叶变换及其反变换对于非周期信号，满足绝对可积的条件下，可利用傅里叶变换对其进行频域分析。，其中，称为信号傅里叶变换，表示了该信号的频谱特性。三、实验内容1.二维平面图形的绘制（任选3个） 假设n=12.对于m=4,5,7,10,在0n2n-1区间上画出,并添上适当标注。用plot和stem分别绘制该信号，并比较。 考虑信号,式中=2k/5.给出k=1,2,4,6,用stem画出每个信号在区间0n9内的图。利用subplot在同一幅图上用单独的坐标轴画出全部符号。 n=6，试画出;的图形。 在0n31内画出下面每一个信号： 。 用stem画出信号定义：y1n=xn-2, y2n=xn+1,y3n=x-n, y4n=x-n+1, 用stem分别画出y1y4，并用legend 命令给出图例。2.设周期信号一个周期的波形为，求该信号傅里叶级数展开式，并用matlab画出傅里叶技术展开后的波形，并通过展开式项数的变化考察其对的逼近程度，考察其物理意义。3.设非周期信号，求该信号的傅里叶变换，matlab画出傅里叶变换后的频谱，并对频谱进行反变换，画出的波形。四、实验程序，实验结果及分析1.二维平面图形的绘制（任选3个） 假设n=12.对于m=4,5,7,10,在0n2n-1区间上画出,并添上适当标注。用plot和stem分别绘制该信号，并比较。% 程序m=4,5,7,10;n=12;n=0:2*n-1;for i=1:length(m) xn=sin(2*pi*m(i)*n/n); subplot(length(m),1,i) plot(n,xn);hold on s=int2str(m(i) title(m= s); stem(n,xn,r);end结果： 考虑信号,式中=2k/5.给出k=1,2,4,6,用stem画出每个信号在区间0n9内的图。利用subplot在同一幅图上用单独的坐标轴画出全部符号。程序：k=1,2,3,6;n=0:9;figurefor i=1:length(k) xk=sin(2*pi*k(i)*n/5); subplot(length(k),1,i) % plot(n,xk);hold on s=int2str(k(i) title(k= s); stem(n,xk,r);end结果： n=6，试画出;的图形。程序：n=6;n=0:100;xn1=cos(2*pi*n/n)+2*cos(3*pi*n/n);xn2=cos(2*n/n)+2*cos(3*n/n);xn3=cos(2*pi*n/n)+3*sin(5*pi*n/n);figuresubplot(3,1,1)stem(n,xn1);holdplot(n,xn1,r)subplot(3,1,2)stem(n,xn2);holdplot(n,xn2,r)subplot(3,1,3)stem(n,xn3);holdplot(n,xn3,r)结果： 在0n31内画出下面每一个信号： 。程序：n=0:31;xn1=cos(pi*n/4).*sin(pi*n/4);xn2=cos(pi*n/4).2;xn3=cos(pi*n/8).*sin(pi*n/4);figuresubplot(3,1,1)stem(n,xn1);holdplot(n,xn1,r)subplot(3,1,2)stem(n,xn2);holdplot(n,xn2,r)subplot(3,1,3)stem(n,xn3);holdplot(n,xn3,r)结果： 用stem画出信号定义：y1n=xn-2, y2n=xn+1,y3n=x-n, y4n=x-n+1, 用stem分别画出y1y4，并用legend 命令给出图例。程序：function x=xn(t)%方法一%x1=2*ones(1,length(t);%x2=ones(1,length(t);%x3=-1*ones(1,length(t);%x4=3*ones(1,length(t);%x=x1.*(t=0)+x2.*(t=2)+x3.*(t=3)+x4.*(t=4);% 方法二x=zeros(1,length(t);x=2*(t=0)+1*(t=2)+(-1)*(t=3)+3*(t=4);%主程序t= -10 : 10 ;y0=xn(t) ;y1=xn(t-2) ;y2=xn(t+1) ;y3=xn(-t) ;y4=xn(-t+1) ;subplot(511);p=plot(t,y0,k);set(p,linewidth,2);hold on;stem(t,y0,r)xlabel(t);ylabel(x(t);title(原始信号 x(t);grid;subplot(512);p=plot(t,y1,k);set(p,linewidth,2);hold on;stem(t,y1,r)xlabel(t);ylabel(x(t-2);title(第一个变换);grid;subplot(513);p=plot(t,y2,k);set(p,linewidth,2);hold on;stem(t,y2,r)xlabel(t);ylabel(x(t+1);title(第二个变换);grid;subplot(514);p=plot(t,y3,k);set(p,linewidth,2);hold on;stem(t,y3,r)xlabel(t);ylabel(x(-t);title(第三个变换);grid;subplot(515);p=plot(t,y4,k);set(p,linewidth,2);hold on;stem(t,y4,r)xlabel(t);ylabel(x(-t+1);title(第四个变换);grid;2.设周期信号一个周期的波形为，求该信号傅里叶级数展开式，并用matlab画出傅里叶技术展开后的波形，并通过展开式项数的变化考察其对的逼近程度，考察其物理意义。程序：% 实验二 exp0302.m% 设周期信号一个周期的波形为f(t)=1,|t|tao/2,f(t)=0,其他，求该信号傅里叶级数展开式，% 并用matlab画出傅里叶技术展开后的波形，并通过展开式项数的变化考察其对% 的逼近程度，考察其物理意义。clear all;n=100;% 取展开式的项数为2n+1项t=1;% 周期信号的周期tao=0.5; % 脉宽tao为0.5fs=1/t;% 基频n_sample=128;% 为了画波形，设置每个周期的采样点数dt=t/ n_sample;% 时间分辨率t=0:dt:10*t-dt;% 为10个周期的波形，准备时间轴n=-n:n;% 展开项的自变量n，从-n到n，共2n+1项fn=sinc(tao*n/t)/t; % 求傅立叶系数fn(n+1)=0;% 大于n的项f(n+1)为0ft=zeros(1,length(t);% 取一个长度与t相同的全0矢量for m=-n:n% 求(2-1)式共2n+1项的和，来近似f(t)ft=ft+fn(m+n+1)*exp(j*2*pi*m*fs*t) ;% fn是一个数组，其序号是从1开始的，到% 2n1结束， % 故该语句中为fn(m+n+1) % 而当n=0时，fn=0，在数组中的位置% 为第n+1个元素，故令fn(n+1)=0endplot(t,ft)3.设非周期信号，求该信号的傅里叶变换，matlab画出傅里叶变换后的频谱，并对频谱进行反变换，画出的波形。程序：close allclear all % 删除所有工作空间的变量tao=1; % 脉宽tao_sample=512; % 为了画波形，在脉宽内的采样点数%n_sample=30; % 为了画波形，在时间窗内的采样点数dt=tao/ tao_sample; % 时间分辨率add_zeros=1000;t=(-tao/2-add_zeros*dt):dt:tao/2+(add_zeros-1)*dt; % 所有时间采样值n_sample=length(t); % 为了画波形，在时间窗内的采样点数t=t(end)-t(1); % 时间窗的大小st=zeros(1,add_zeros) ones(1,(n_sample-2*add_zeros) zeros(1,add_zeros);% 依据t将信号离散化subplot(411); % 设置3*1的1号窗plot(t,st); % 在1号窗中画时间波形axis(-1 1 0 2); % 设置1号窗坐标轴的范围xlabel(t);ylabel(s(t); % 标出横轴、纵轴变量subplot(412) ; % 设置3*1的2号窗f,sf=t2f(t,st) ; % 求st的fftplot(f,abs(sf) ; % 在2窗中画st频谱sf波形的绝对值axis(-10 10 0 2); % 设置2号窗坐标轴的范围xlabel(f);ylabel(|s(f)|); % 标出横轴、纵轴变量%sff=tao*sinc(f*tao) ; %例6信号谱sff=tao*(sin(tao*pi*f)./(tao*pi*f) ; %例6信号谱subplot(413) ; % 设置3*1的2号窗plot(f,abs(sff),r) % 在2窗中画st频谱sf波形的绝对值，用红色axis(-10 10 0 2); % 设置2号窗坐标轴的范围tt,st= f2t_3(f,sf);% 进行离散傅立叶反变换，求原始信号subplot(414) ; % 设置3*1的3号窗xlabel(t);ylabel(恢复的s(t); % 标出横轴、纵轴变量plot(tt,st) ;hold off ;% 在3窗中画sf回复的st波形，关闭共画开关axis(-1 1 0 2);% 标出横轴、纵轴变量结果：实验三 随机信号与数字基带实验一、实验目的l 掌握库函数产生随机数方法l 基带信号波形生成和其功率谱密度l 理解采用蒙特卡罗算法仿真的思想二、实验原理（一）.库函数产生随机数 均匀分布的随机数利用matlab库函数rand产生。rand函数产生（0，1）内均匀分布的随机数，使用方法如下：1）x=rand(m)；产生一个mm的矩阵，所含元素取值均为在（0，1）内均匀分布的随机数。2）x=rand(m，n)；产生一个mn的矩阵，所含元素取值均为在（0，1）内均匀分布的随机数。3）x=rand；产生一个随机数。4）调用rand(1,n)给出均值为0.5，功率为1/12=0.083的白噪声。5）调用rand(1,n)给出均值为0.5，功率为的白噪声，现在要均值为0，功率为0.01。 调整均值：用u(n)减上均值即可；调整功率：令希望的功率为p，则需要求出常数a，用a乘u(n)。 高斯分布的随机数randn函数产生均值为0，方差为1的高斯分布的随机数，使用方法如下：1）x=randn(m)；产生一个mm的矩阵，所含元素都是均值为0，方差为1的高斯分布的随机数。2）x=randn(m，n)；产生一个mn的矩阵，所含元素都是均值为0，方差为1的高斯分布的随机数。3）x=randn；产生一个均值为0，方差为1的高斯分布的随机数。4）调用randn，生成均值为0，方差为1，服从高斯分布的白噪声信号u(n)。调整均值：用u(n)加上均值即可；调整功率：令希望的功率为p，则需要求出常数a，用a乘u(n)。（二）.基带信号波形生成和信号的功率谱密度1).要画出完整的基带信号波形，每一个码元要采 n个样。2）信号的功率谱密度为：。（三）.蒙特卡罗算法蒙特卡罗估计是指通过随机实验估计系统参数值的过程。蒙特卡罗算法的基本思想：由概率论可知，随机实验中实验的结果是无法预测的，只能用统计的方法来描述。故需进行大量的随机实验，如果实验次数为，以表示事件发生的次数。若将发生的概率近似为相对频率，定义为。这样，在相对频率的意义下，事件发生的概率可以通过重复无限多次随机实验来求得，即：在二进制数字通信系统中，若是发送端发送的总码元数，是差错发生的次数，则总误码率可通过蒙特卡罗算法计算。三、实验内容1.产生随机数利用rand函数产生54的（0，2）内均匀分布的随机数；利用randn函数产生54的均值为0，方差为2的高斯分布的随机数。2. 利用随机数产生单极性基带信号，每一个码元要采 8个样值，并画出其波形和功率谱密。3.利用蒙特卡罗算法仿真二进制基带通信系统的误码率假定通信系统满足以下条件： 信源输出的数据符号是相互独立和等概的双极性基带信号 发送端没有发送滤波器，接收端没有接收滤波器 信道是加性高斯白噪声信道数字基带信号传输系统模型如图1所示：图1 数字基带信号传输系统模型当时：最佳判决门限：，误码率：，利用，可得，故可用q函数表示误码率。抽样判决器输入信噪比： 抽样判决器输入信号为：为判决器输入有用信号电压，为信道输入的均值为0，方差为高斯噪声。通信系统的蒙特卡罗仿真模型如图2所示。编程实现二进制基带通信系统的误码率的蒙特卡罗仿真，并和理论误码率比较。图2 通信系统的蒙特卡罗仿真模型四、预习报告：（原理，流程图，编程思想，基本程序）五、实验报告：（调试好的程序，实验结果及分析）实验四 模拟调制实验一、实验目的l 掌握线性模拟调制信号的波形及产生方法；l 掌握线性模拟调制信号的频谱特点；l 掌握线性模拟调制信号的解调方法；l 掌握线性模拟调制系统的matlab仿真实现。二、实验原理1. am调制am信号的时域表示式：频谱：调制器模型如图3所示：图3 am调制器模型am的时域波形和频谱图如图4所示：时域 频域图4 am调制时、频域波形am的调制效率最高为1/3；其优点：可以采用包络检波进行解调，实现起来容易；缺点：调制效率不高。2. dsb-sc调制dsb信号的时域表示式频谱：dsb的时域波形和频谱图如下：时域 频域图5 dsb调制时、频域波形dsb的相干解调模型图：图6 dsb调制器模型dsb的调制效率为100；其优点：节省了载波功率；缺点：不能用包络检波，需用相干检波，较复杂3. ssb调制ssb信号的时域表示式频谱：分为上边带和下边带，均为双边带的一半。ssb的调制效率为100；其优点：节省发射功率，而且它所占用的频带宽度比am、dsb减少了一半；缺点：调制实现起来困难，解调需用相干检波，较复杂。三、实验内容（任选一）1.用matlab产生一个频率为1hz，振幅为1的余弦信源，设载波频率为10hz，a=2。am系统中，当相干解调器输入信噪比为15db时，系统各点的波形及其频谱，并对调制信号和解调信号进行比较分析。分析比较不同信噪比下的解调信号。2.用matlab产生一个频率为1hz，振幅为1的余弦信源，设载波频率为10hz。dsb-sc系统中，当相干解调器输入信噪比为15db时，系统各点的波形及其频谱，并对调制信号和解调信号进行比较分析。分析比较不同信噪比下的解调信号。3.用matlab产生一个频率为1hz，振幅为1的余弦信源，设载波频率为10hz。ssb系统中，当相干解调器输入信噪比为15db时，系统各点的波形及其频谱，并对调制信号和解调信号进行比较分析。分析比较不同信噪比下的解调信号。四、预习报告：（原理，流程图，编程思想，基本程序）五、实验报告：（调试好的程序，实验结果及分析）实验五 模拟信号数字传输实验（一）一、实验目的l 掌握低通信号抽样定理l 理解13折线a率逐次比较型pcm编码仿真的思想二、实验原理1.低通信号的抽样定理一个频带为0, fh的低通信号，可以无失真地被抽样速率的抽样序列所恢复。2.验证低通信号的抽样定理抽样的过程是将输入的模拟信号与抽样信号相乘而得，通常抽样信号是一个周期为ts的周期脉冲信号，抽样后得到的信号称为抽样序列。理想抽样信号定义如下：其中，称为抽样速率。因此抽样后的信号为经带宽为低通滤波器后可恢复抽样信号。3.pcm的基本原理框图如下图所示：抽样是对模拟信号进行周期性的扫描， 把时间上连续的信号变成时间上离散的信号。我们要求经过抽样的信号应包含原信号的所有信息， 即能无失真地恢复出原模拟信号， 抽样速率的下限由抽样定理确定。 量化是把经抽样得到的瞬时值进行幅度离散，即指定q规定的电平，把抽样值用最接近的电平表示。 编码是用二进制码组表示有固定电平的量化值。实际上量化是在编码过程中同时完成的。4.逐次比较型pcm编码的产生原因：均匀量化时其量化信噪比随信号电平的减小而下降。产生这一现象的原因就是均匀量化时的量化级间隔为固定值，而量化误差不管输入信号的大小均在(-/2, /2)内变化。故大信号时量化信噪比大，小信号时量化信噪比小。对于语音信号来说，小信号出现的概率要大于大信号出现的概率，这就使平均信噪比下降。同时，为了满足一定的信噪比输出要求，输入信号应有一定范围(即动态范围)， 由于小信号信噪比明显下降，也使输入信号范围减小。要改善小信号量化信噪比，可以采用量化间隔非均匀的方法，即非均匀量化。压扩特性曲线有a律压扩特性曲线，其表达式为工程上使用13折线来近似a律压扩特性曲线，目的是为了使用逐次比较算法来提高运算速度。5.逐次比较型编码的编码方法码位安排：极性码 段落码 段内码m1 m2m3m4 m5m6m7m8其中： 第一位m1 表示量化值的极性正负。m1=1 代表信号极性为正 m1=0 代表信号极性为负；后面7位分为段落码和段内码两部分,用于表示量化值的绝对值。 m2m3m4为段落码，分别对应表1段落编码m5m6m7m8为段内码，分别对应表1pcm编码的后四位.下表为a律13折线的编码规则，输入的模拟信号经过抽样、量化后，每个抽样值编码成8个比特的二进制码组，如表1所示。表1 a律pcm编码三、实验内容1.利用matlab软件验证低通抽样定理若低通信号为，则完成画出该低通信号的波形；画出抽样速率为的抽样序列，画出抽样后的信号波形；画出经低通滤波器恢复的波形。2.设输入一个样值x-2048，+2048，对x进行a律pcm编码。要求编写成函数，该函数输入变量为样值，输出变量为a律13折线逐次比较的八位pcm编码，设码元宽度为1，画出其波形。四、预习报告：（原理，流程图，编程思想，基本程序）五、实验报告：（调试好的程序，实验结果及分析）实验六 模拟信号数字传输实验（二）一、实验目的l 掌握13折线a律逐次比较型pcm编，译码原理二、实验原理参见实验五模拟信号的数字传输仿真（一）实验原理部分三、实验内容输入信号，对该信号进行采样，要求符合低通采样定理，然后将采样信号进行非均匀量化，对每个样值的量化结果利用模拟信号数字传输的仿真（一）所编写的编码函数进行pcma律13折线编，译码。（1） 画出原始模拟信号的波形（2） 画出编码结果的波形（3） 试编写程序将编码结果进行译码并通过低通滤波器来恢复模拟信号，对恢复的信号与原始模拟信号进行比较分析。四、预习报告：（原理，流程图，编程思想，基本程序）五、实验报告：（调试好的程序，实验结果及分析）实验七 数字频带传输系统实验一、实验目的l 掌握数字频带传输系统调制解调的仿真过程l 掌握数字频带传输系统误码率仿真分析方法二、实验原理数字频带信号通常也称为数字调制信号，其信号频谱通常是带通型的，适合于在带通型信道中传输。数字调制是将基带数字信号变换成适合带通型信道传输的一种信号处理方式，正如模拟通信一样，可以通过对基带信号的频谱搬移来适应信道特性，也可以采用频率调制、相位调制的方式来达到同样的目的。1.调制过程1）2ask如果将二进制码元“0”对应信号0，“1”对应信号，则2ask信号可以写成如下表达式：，。可