北邮通信原理软件实验报告

1、通信原理软件实验报告学院：信息与通信工程学院班级：12 / 31通信原理Matlab仿真实验实验八一、实验内容假设基带信号为m(t)=sin(2000*pi*t)+2cos(1000*pi*t),载波频率为20kHz,请仿真出AM、DSB-SC、SSB信号，观察已调信号的波形和频谱。二、实验原理1、具有离散大载波的双边带幅度调制信号AM该幅度调制是由DSB-SCAM信号加上离散的大载波分量得到，其表达式及时间波形图为：s=.cos2/7+cos2/77c二,41+股口包络应当注意的是，m的绝对值必须小于等于1,否则会出现下图的过调制:I 5.1JW I，）AM信号的频谱特性如下图所示:$(/)

2、=五可/-勿+可"勿+叫(/-勿+因江+m£由图可以发现，AM信号的频谱是双边带抑制载波调幅信号的频谱加上离散的大载波分量。2、双边带抑制载波调幅(DSB-SCAM)信号的产生双边带抑制载波调幅信号s(t)是利用均值为0的模拟基带信号m和正弦载波c相乘得到，如图所示：Sn筋。)=神(')/8其”，十%)m和正弦载波s的信号波形如图所示:若调制信号m是确定的，其相应的傅立叶频谱为M，载波信号c的傅立叶频谱是C，调制信号s(t)的傅立叶频谱S由M和C相卷积得到，因此经过调制之后，基带信号的频谱被搬移到了载频fc处，若模拟基带信号带宽为W,则调制信号带宽为2W,并且频谱中

3、不含有离散的载频分量，只是由于模拟基带信号的频谱成分中不含离散的直流分量。3、单边带条幅SSB信号双边带抑制载波调幅彳S号要求信道带宽B=2W,其中W是模拟基带信号带宽。从信息论关点开看，此双边带是有剩余度的，因而只要利用双边带中的任一边带来传输，仍能在接收机解调出原基带信号，这样可减少传送已调信号的信道带宽。单边带条幅SSBAM信号的其表达式：与和(')=4期Wco$2笈+4碗("sin2/ifct或$5(，)=4用COS-4而(f)sin17rfct其频谱图为：上边惜下边橹上也带颠谪/V叫LSBLSB三、仿真设计1、流程图:2、实验结果盼析讨论实验仿真结果FigureNo

4、.1nxFigureNo.2Filt：Edi+ViewInser+Toolsyind口wHelpLDsclLL.100500-2bp-H-1-H-H-1-ih-5-20-15-10-50510152025nx从上至下依次是AM信号、DSB信号、SSB信号。从仿真结果看，AM调制信号包络清晰，可利用包络检波恢复原信号，接收设备较为简单。其频谱含有离散大载波，从理论分析可知，此载波占用了较多发送功率，使得发送设备功耗较大。DSB-SC信号波形和频谱，其时域波形有相位翻转，频谱不含离散大载波。SSB信号比DSB信号节省一半带宽，适合于语声信号的调制，因为其没有直流分量，也没有很低频的成分。3、结果分

5、析：根据通原理论课的知识可知，信号的AM调制比较容易实现，但其功率谱中有相当大一部分是载频信号，效率非常低；DSB-SC调制解决了AM信号效率低下的问题，但仍然存在的问题是调制信号的带宽为基带信号的两倍，频谱利用率较低；SSB调制方式在频谱利用上又做出了改进，为原先的一半，但其可靠性降低了，总之，可靠性与有效性是难以两全其美的，为一对矛盾体。四、程序代码fs=800;%kHzT=200;N=T*fs;t=linspace(-T/2,T/2,N);f=linspace(-fs/2,fs/2,N);fm1=1;%kHzfm2=0.5;%kHzfc=20;%kHzmt=sin(2*pi*fm1*t)

6、+2*cos(2*pi*fm2*t);%am!M制信号波形和频谱A=2;a=0.3;AM=(1+a*mt).*cos(2*pi*fc*t);figure(1);subplot(2,1,1);plot(t,AM);xlabel('时间t(ms)');ylabel('AM');axis(-3,3,-3,3);title('AM波形');gridon;FAM=12f(AM,fs);subplot(2,1,2);plot(f,abs(FAM);xlabel('频率f(kHz)');ylabel('FAM');axis(-2

7、5,25,0,150);title('AM幅频特性');gridon;%DSB-SC信号波形和频谱DSB=mt.*cos(2*pi*fc*t);figure(2);subplot(2,1,1);plot(t,DSB);xlabel('时间t(ms)');ylabel('DSB');axis(-3,3,-3,3);title('DSB波形')；gridon;FDSB=12f(DSB,fs);subplot(2,1,2);plot(f,abs(FDSB);xlabel('频率f(kHz)');ylabel('F

8、DSB');axis(-25,25,0,150);title('DSB幅频特性');gridon;%SSB言号波形和频谱M=12f(mt,fs);MH=-j*sign(f).*M;mh=real(f2t(MH,fs);SSB=mt.*cos(2*pi*fc*t)-mh.*sin(2*pi*fc*t);FSSB=12f(SSB,fs);figure(3);subplot(2,1,1);plot(t,SSB);xlabel('时间t(ms)');ylabel('SSB');axis(-3,3,-5,5);title('SSB波形

9、9;)；gridon;subplot(2,1,2);plot(f,abs(FSSB);xlabel('频率f(kHz)');ylabel('FSSB');axis(-25,25,0,150);title('SSB幅频特性);gridon;实验九一、实验内容假设基带信号为m(t)=sin(2000兀t)+2cos(1000兀t)+4sin(500兀t+兀/3),载波频率为40kHz,仿真产生FM信号，观察波形与频谱。FM的频偏常数为5kHz/V。二、实验原理1、调频信号表达式为：中卜区侬"-/93=2符/用()号一频率偏移常数2、调频信号的频谱特

10、征：单频调制信号包括无穷多频率分量，在实际应用中，我们运用等效带宽这个概念。等效带宽的定义为：包含98%或者99%的已调信号总功率的带宽，根据卡松公式有：Bfm定2(Kfm(t)max+fm)三、仿真设计1、流程图基带信号m(t)J利用scilab中的、cumsum函数对基带信号m(t)进行求累积和，相当于对L其进行积分将累计求和的结果作为瞬时相位偏移，得到已调信号s(t)的表送式，即FM信号傅立叶变换广作图(时域图&频域图)<2、实验结果盼析讨论3、结果分析由于max(m(t)=6.4833,信号带宽fm=1KHz,Kf=5KHz/V,由卡松公式Bfm 制 2(Kf|m(t)|

11、max +fm)，计算得带宽为66.8320KHZ,与频域仿真结果基本相符。（事实上，由于调制信号不是单频信号，卡松公式不满足）四、程序代码fs=800;%kHzT=200;N=T*fs;dt=1/fs;t=linspace(-T/2,T/2,N);f=linspace(-fs/2,fs/2,N);fm1=1;%kHzfm2=0.5;%kHzfm3=0.25;%kHzfc=40;%kHzk=5;%kHzmt=sin(2*pi*fm1*t)+2*cos(2*pi*fm2*t)+4*sin(2*pi*fm3*t+pi/3);phi=2*pi*k*cumsum(mt)*dt;FM=cos(2*pi*

12、fc*t+phi);subplot(2,1,1);plot(t,FM);xlabel('时间t(ms)');ylabel('FM');axis(-2,2,-3,3);title('FM波形');gridon;FFM=12f(FM,fs);subplot(2,1,2);plot(f,abs(FFM);xlabel('频率f(kHz)');ylabel('FFM');axis(-80,80,0,20);title('FM幅频特性');gridon;实验十一实验内容通过仿真测量占空比为50%、75%以及1

13、00%的单、双极性归零码波形以及其功率谱，分析不同占空比对仿真结果的影响。二、实验原理1、时域特性1.1 二进制单极性归零码：17 / 311.2 二进制双极性归零码：b=0,1oa=-L+lfrfr(t1gT(t)=Arect-2、频域特性（功率谱）：数字基带信号s的功率谱密度为:=J：K("*=%=j-pa(/)iGr(/)|：数字基带信号s的功率谱密度与随机序列“”的功率谱特性，（/）以及发送滤波器的频率特性Gr（/）有关。在实随机序列卜”的各符号互不相关时，s的功率谱为上""会认加+分停）H?一三）2.1单极性归零码:三、仿真设计1、流程图：利用设置占空比

14、利用公式作图randn(1,（分别为P=abs(S1).A2/T（波形M）以及25%、50%、,求样本信号的&功率逻辑判75%、100%)功率谱密度，然k谱）断，产生后对各个样本单（双）的功率谱密度极性数求数学期望，得据到功率谱密度Jk,JLJk2、实验结果&分析讨论2.1 单极性从上至下依次是占空比为50%、75%、100%。从仿真结果可以看出，单极性归零码的频谱主瓣宽度随占空比增加而减小，且含有冲激2.2 双极性从上至下依次是占空比50%、75%、100%。从仿真结果可以看出，随占空比增加，频谱主瓣宽度减小，且不含冲激。3、结果分析：3.1 单极性由于单极性码含有直流分量，

15、所以表现在频域内为在直流处奇级次谐波处有一个冲激，而其功率谱主瓣宽度随着占空比的不同而不同，对于50%、75%、100%的占空比的带宽分别对应为码元速率B的4倍、2倍、1.33倍和1倍。3.2 双极性由于双极性码不含直流分量，所以没有单极性码所具有的直流处和奇次谐波处的冲激，但在带宽上与单极性码表现一致。3.3 对比因为很多传输线路不能传送直流分量，而且在误码特性上双极性马要比单极性码好很多，故此，一般都优先选用双极性码。四、程序代码%通过仿真测量占空比为50%、75%以及100%的单双极性归零码波形及其功率谱，分析不同占空比对仿真结果的影响clearall;xdel(winsid();/关闭

16、所有图形窗口exec12f.sci;execf2t.sci;L=32;/每个码元间隔内的采样点数N=2A13;/总采样点数M=N/L;/总码元数Rb=2;/码元速率Ts=1/Rb;/比特间隔fs=L/Ts;/采样速率T=N/fs;/截短时间Bs=fs/2;/系统带宽t=-T/2+0:N-1/fs;/时域采样点f=-Bs+0:N-1/T;/频域采样点L0=input('请输入占空比(01):')EP=zeros(1,N);ch=input('请选择要观察的码型：1-单极性；其他-双极性')forloop=1:1000/1000次样本函数取平均ifch=1a=(ra

17、nd(1,M)>0.5)+0;/生成单极性序列elsea=sign(rand(1,M)>0.5)-0.5);/生成双极性序列endtmp=zeros(L,M);/一个码元的归零部分取零L1=L*L0;/占空比，求出一个码元不归零部分的取样点数tmp(1:L1,:)=ones(L1,1)*a;/将一个码元不归零部分的取样点值置为1s=tmp(:)'S=t2f(s,fs);/傅里叶变化P=abs(S)A2/T;/样本的功率谱密度EP=EP*(1-1/loop)+P/loop;/随机过程的功率谱是各个样本的功率谱的数学期望endxset("window",1)

18、plot(t,s)set(gca(),"grid",1,1)title('时域图')xlabel('t')ylabel('S(t)')mtlb_axis(-3,3,-1.5,1.5);xset("window",2)plot(f,abs(EP+%eps)set(gca(),"grid",1,1)title('功率谱图形')xlabel('f')ylabel('功率')mtlb_axis(-35,35,-5,max(EP+%eps);xset

19、("window",3)plot(f,10*log10(EP+%eps)16/31set(gca(),"grid",1,1)title('功率谱图形(dB)')xlabel('f')ylabel('功率')实验十二一、实验内容仿真测量滚降系数为a=0.25的根升余弦滚降系统的发送功率谱密度及眼图二、实验原理1.理想限带情况下的最佳基带传输:发送滤波信道接收滤波白味声取样判决nw(t)a的存在，使得传输频带增加，降低了频谱利用率。A愈小，波形振荡起伏愈大，传输频带扩展小；反之，a愈大，波形振荡起伏愈小，频带扩

20、展增大。实际中一般选取a不小于0.2。2、眼图实际通信系统中，数字信号经过非理想的传输系统产生畸变，总是在不同程度上存在码间干扰的，系统性能很难进行定量的分析，常常甚至得不到近似结果。而眼图可以直观地估价系统码问干扰和噪声的影响，是常用的测试手段。眼图分析中常用结论：(1)最佳取样时刻应选择在眼睛张开最大的时刻；(2)眼睛闭合的速率，即眼图斜边的斜率，表示系统对定时误差灵敏的程度，斜边愈陡，对定位误差愈敏感；(3)在取样时刻上，阴影区的垂直宽度表示最大信号失真量；(4)在取样时刻上，上下两阴影区的间隔垂直距离之半是最小噪声容限，噪声瞬时值超过它就有可能发生错误判决；(5)阴影区与横轴相交的区间

21、表示零点位置变动范围，它对于从信号平均零点位置提取定时信息的解调器有重要影响。三、仿真设计1、流程图产生M个0、2等概随机码4定义升余弦滚降系统系统函数序列经过升余弦滚降滤波器傅里叶反变换得时域波形作出功率谱和眼图一2、实验结果&分析讨论发送功率谱密度：实验结果：仿真结果与计算公式WuRsO+a"2=2(1+0.25)/2=1.25基本相等。眼图:Figure：No.2五、程序代码clearallglobaldttdfNcloseallN=2A14;L=32;M=N/L%采样点数%每码元的采样点数%码元数Rb=2;Ts=1/Rb;dt=Ts/L;df=1/(N*dt)T=N*

22、dt%码速率是2Mb/s%码元间隔%时域采样间隔%频域采样间隔%时域截短时间Bs=N*df/2Na=4;alpha=0.25;t=-T/2+dt/2:dt:T/2;f=-Bs+df/2:df:Bs;%系统带宽%设示波器扫描宽度为4个码元%时域横坐标%频域横坐标g1=sin(pi*t/Ts)./(pi*t/Ts);g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).A2);g=g1.*g2；G=t2f(g,df);%升余弦脉冲波形%升余弦脉冲波形的傅氏正变换figure(1)%set(1,'Position',100,50,300,200)设定窗口位

23、置及大小%set(2,'Position',350,50,300,200); % figure(2) hold on%grid xlabel('t in us') ylabel('s(t) in V')EP=zeros(size(f)+eps;%for ii=1:2%a=sign(randn(1,M);imp=zeros(1,N);%imp(L/2:L:N)=a/dt; S=t2f(imp,df).*G;%s=f2t(t2f(imp,df).*G,df);%s=real(s); P=S.*conj(S)/T;%EP=(EP*(ii-1)+P+ep

24、s)/ii;figure(1) plot(f,30+10*log10(EP),'g'); grid axis(-3,+3,-200,200) xlabel('f (MHz)') ylabel('Ps(dBm/MHz)') figure(2)tt=0:dt:Na*L*dt; % for jj=1:Na*L:N-Na*L %Na*L %最后一点的起始点。plot(tt,s(jj:jj+Na*L); % endend%设定窗口位置及大小图形叠加，即重复扫描力Deps是为了躲零设累计平均的循环产生冲激序列升余弦信号的傅氏变换升余弦信号的时域波形升余弦信号

25、的功率谱%累计平均设定画眼图的坐标的位置是每次画眼图的间隔步长，N-Na*L是画眼图的画升余弦信号的眼图对应累计平均的循环实现数字基带传输系统实验目的：利用编程实现数字基带传输系统1）画出发送端输入码序列波形和功率谱、发送滤波器输出波形和功率谱2）画出接收端采样判决后码序列波形和功率谱、接收滤波器输出波形和功率谱3）画出接收滤波器输出信号眼图（在升余弦滚将系数分别为0、0.5、1的情况下）4）分别画出升余弦滚将系数为0、0.5、1,采样判决点在眼图最大处的系统的实际误码曲线（Pes/n0曲线），并在同坐标系中画出理论误码曲线5）改变采样点重复1） 4）。作出发端功率谱图和眼图作出收端功率谱、眼

26、图、判决输出海形程序代码:主程序：clearallexec12f.sci;execf2t.sci;execeyes.sci;k=13;N=2”;采样点数L=16;/每码元的采样点数M=N/L;/码元数Rb=2;/码速率是2Mb/sTs=1/Rb;码元间隔dt=Ts/L;时域采样间隔fs=1/dt;采样速率df=1/(N*dt);/频域采样间隔T=N*dt;截短时间Bs=N*df/2;系统带宽f=-Bs+df/2:df:Bs;/频域横坐标t=-T/2+dt/2:dt:T/2;/时域横坐标alpha=1;/升余弦滚降系数Hcos=zeros(1,N);ii=find(abs>(1-alpha

27、)/(2*Ts)&abs(f)<=(1+alpha)/(2*Ts);Hcos(ii尸Ts/2*(1+cos(%pi*Ts/alpha*(abs(f(ii)-(1-alpha)/(2*Ts);ii=find(abs<=(1-alpha)/(2*Ts);Hcos(ii)=Ts;GT=sqrt(Hcos);GR=GT;最佳系统的发送接收滤波器的傅式变换EPC=zeros(1,N);EPB=zeros(1,N);EP1=zeros(1,N);EP2=zeros(1,N);RECT=ones(L,1);for11=1:20;Eb_N0(11)=(11-1);信噪比eb_n0(11)=

28、10A(Eb_N0(11)/10);Eb=1;n0=Eb/eb_n0(11);/信道的噪声谱密度sita=n0*Bs;信道中噪声功率n_err=0;误码计数for1oop=1:100a=sign(rand(1,M,"norma1");tmp1=RECT*a;seq_send=tmp1(:)'发送端产生序列！S1=t2f(seq_send,fs);P1=abs(S1)A2/T;样本信号的功率谱密度EP1=EP1*(1-1/1oop)+P1/1oop;/随机过程的功率谱是各个样本的功率谱的数学期望EP11=20+10*1og10(EP1+%eps);sp=zeros(1

29、,N);产生冲激序列sp(1:L:N尸a/dt;SP=t2f(sp,fs);SH=SP.*GT;sh=real(f2t(SH,fs);通过发送滤波器后的信号！PB=abs(SH).A2/T;/过升余弦后能量谱密度EPB=EPB*(1-1/loop)+PB/loop;/求平均EPBB=20+10*log10(EPB+%eps);nw=sqrt(sita)*rand(1,N,"normal");信道噪声SY=SH.*GR;sy=real(f2t(SY,fs);r=sy+nw;接收信号R=t2f(r,fs);PC=abs(R).A2/T;EPC=EPC*(1-1/loop)+PC

30、/loop;收端EPCC=20+10*log10(EPC+%eps);sd=r(L/4:L:N);/*取样！sdp=sign(sd);判决tmp2=RECT*sdp;seq_reci=tmp2(:)'接收端信号！S2=t2f(seq_reci,fs);P2=abs(S2)A2/T;样本信号的功率谱密度EP2=EP2*(1-1/loop)+P2/loop;EP22=20+10*log10(EP2+%eps);n_err=n_err+length(find(sdp=a);/错误累计endPe(l1)=n_err/(M*loop);xset("window",5)plot

31、(Eb_N0,log10(Pe+%eps),'g');/PeEb/N0曲线画图xlabel('Eb/N0');ylabel('Pe');title("PeEb/N0曲线)eb_n0=10.A(Eb_N0/10);set(gca(),"auto_clear","off")k=log10(0.5*erfc(sqrt(eb_n0);plot(Eb_N0,k)mtlb_axis(0,15,-3.5,0);xlabel('Eb/N0')ylabel('Pe')legend(&

32、#39;实际的,'理论的');set(gca(),"auto_clear","on")endxset("window",1)subplot(2,2,1)plot(t,seq_send)title("最初产生序列输出波形")mtlb_axis(0,L,-1.5,1.5)xgridsubplot(2,2,2)plot(t,sh)title("通过发送滤波器后的信号波形”)xlabel("t(ms)")ylabel("s2(t)(V)")mtlb_axis

33、(0,L,-2.5,2.5)xgridsubplot(2,2,3)plot(t,sy)title("采样前的彳t号波形")xlabel("t(ms)")ylabel("y(t)(V)")mtlb_axis(0,L,-2,2)xgridsubplot(2,2,4)plot(t,seq_reci)title("判决后输出序列”)mtlb_axis(0,L,-1.5,1.5)xgridxset("window",2)subplot(2,2,1)plot(f,EP11)xlabel("f/kHz")ylabel("|FMf|/(V/Hz)”)mtlb_axis(-10,10,-20,max(EP11)title("发送序列频率谱”)xgridsubplot(2,2,2)plot(f,EPBB)title("发送滤波器输出功率谱”)xlabel("f(kHz)”)ylabel("功率谱(W/kHz)")mtlb_axis(-10,10,0,max(EPBB)xgridsubplot(2,2,3)plot(f,EPCC)title("判决后输出功率谱&quo