北邮通原软件实验文档格式.docx

1、title（AM 信号频谱 ）xlabel（fylabel（S（f）axis（-25,25,0,max（abs（S1）;%xset（window,2）figure（2）plot（t,s1）AM 信号波形 ts（t）axis（-3,3,-3,3）;（ 4）实验结果AM 信 号 频 谱100908070605040302010-25 -20 -15 -10 -5 0 5 10 15 20 25f（kHz）2 DSB-SC信号（ 1 ）信号的产生和表达式c绘制时域波形及频谱3）源代码 %KHz %ms% kHz%DSB-SC 信号s= m.*cos（2*pi*fc*t）;plot（f,abs（S2）

2、DSB-SC 信号频谱 axis（-25,25,0,max（abs（S2）;figure（2）plot（t,s2）DSB-SC 信号波形 axis（-1,4,-3,3）;4）实验结果DSB-SC 信 号 波 形-1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4t（ms）DSB-SC 信 号 频 谱10 15 20 25-25 -20 -15 -10 -5DSC-SB 频谱3 SSB信号1 ）信号的产生和表达式s上 tsDSB SC AM t s下 tAcm t cos2 fct Ac m t sin 2 fctsSSB t Acm t cos2 fct mAc m t sin

3、2 fct2）流程图SSB信号s= m.*cos（2*pi*fc*t）-mh.*sin（2*pi*fc*t）;傅氏变换 S= t2f（s,fs）3）源代码：%SSB 信号的产生M= t2f（m,fs）;MH=-j*sign（f）.*M; %在频域进行希尔伯特变换mh= real（f2t（MH,fs）; %希尔伯特变换后的信号s= m.*cos（2*pi*fc*t）-mh.*sin（2*pi*fc*t）; %SSB signalplot（f,abs（S3）SSB 信号频谱 axis（-25,25,0,max（abs（S3）plot（t,s3）SSB 信号波形 axis（0,6,-3,3）321t

4、）（s 0-1-2-30123456SSB 信 号 频 谱 200 180 160 140 120 f）（S 100 80 60 40 20 0实验二假设基带信号为 m（t）=sin（2000 t）+2cos（1000 t）+4sin（500 t+ /3）, 载波频率为 40kHz，仿真产生 FM信号，观察波形与频谱，并与卡松公式作对照。 FM的频率偏移常数是 5kHz/V。（ 1 ）信号表达式（ 2）流程图 %kHzT= 16;f= -fs/2:Kf= 5; %kHz/Vfc= 40;m= sin（2*pi*fm*t）+2*cos（1*pi*fm*t）+4*sin（0.5*pi*fm*t+p

5、i/3）;phi= 2*pi*Kf*cumsum（m）*dt; %求相位s= cos（2*pi*fc*t+phi）; % s（t）plot（f,abs（S）.2）FM 信号功率谱 axis（-80,80,0,max（abs（S）.2）; % 功率谱密度为 |S|2plot（f,abs（S） 调制信号频谱 axis（-80,80,0,max（abs（S）;figure（3）plot（t,s）FM 信号波形 axis（0,3,-2,2）;FM 信 号 波 形调制信号频谱-80 -60 -40 -20 0 20 40 60 80fm 取 1kHz, 用卡松公式计算得到 FM 信号带宽 :Bfm=2*

6、（Kf*max（abs（m）+1） = 66.8325FM 频谱图比较，基本相等，说明实验 FM 信号带宽与理论值基本相符。FM 信号的特点。实验三25%、 50%、 75%以及 100%的单、双极性归零码波形及其clear all ;close all ;L=32; %每个码元间隔内的采样点数N=213; %总采样点数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（ 请输入占空比（

7、 01 ）： ）EP=zeros（1,N）;ch=input（ 请选择要观察的码型： 1- 单极性； 2- 双极性：for loop=1:1000 %1000 次样本函数取平均if ch=1a=（rand（1,M）0.5）; %生成单极性序列elsea=sign（rand（1,M）0.5）-0.5）; %生成双极性序列endtmp=zeros（L,M）; %一个码元的归零部分取零L1=L*L0; %占空比，求出一个码元不归零部分的采样点数tmp（1:L1,:）=ones（L1,1）*a; %将一个码元不归零部分的取样点置为s=tmp（:）;S=t2f（s,fs）;P=abs（S）.2/T; %

8、样本功率谱密度EP=EP*（1-1/loop）+P/loop; %随机部分的功率谱是各个样本功率谱的期望grid on title（ 时域图 ） xlabel（ ylabel（ S（t）axis（-3,3,-1.5,1.5）;figure（2） plot（f,abs（EP+eps） 功率谱图形 xlabel（ 功率 axis（-35,35,-5,max（EP+eps）;plot（f,10*log10（EP+eps）grid ontitle（ 功率谱图形 （dB）实验结果：（ 1 ） . 单极性）zHk/2V（ 率功功 率 谱 图 形 （dB）修改占空比可得到以下图形时域图v（ 率 功zkH率功

9、）zHk/Bd（）zHk/2v（ 率功-40 -30 -20 -10 0 1050%、 75%、 100%的波形图及功率谱密度图。从仿真单极性归零码的频谱主瓣宽度随占空比增加而减小， 且含有冲激。修改占空比后得到以下图形：率 功功率谱图形从上至下依次是占空比 50%、 75%、 100%。从仿真结果可以看出，随占空比增加，频谱主瓣宽度减小，且不含冲激。试验结论：单极性归零码和双极性归零码的图形由仿真得到， 其功率谱有一定特点， 单极性归零码的功率谱有支流分量，因为其均值不为零，双极性码均值为零，故没有直流分量。 占空比为100%时，相当于不归零码，功率谱符合部归零码的特点。实验四仿真测量滚降系

10、数为 =0.25 的根升余弦滚降系统的发送功率谱及眼图。（ 1 ）仿真模型：设定采样频率、采样点数、时间截短等设定滚降系数为 =0.25 的根升余弦滚降系统Hrcos产生 PAM信号 s2累计平均 EP 绘制发送功率谱信号通过理想高斯白噪基带信道得到r= s2+nw匹配滤波 绘制采样前信号的眼图（ 3）源代码clear allL=16;Rs=2;Ts=1/Rs;%生成升余弦alpha=0.25; %滚降系数hcos=zeros（1,N）;%升余弦表达式ii=find（abs（f）（1-alpha）/（2*Ts）&abs（f）=（1+alpha）/（2*Ts）;hcos（ii）=Ts/2*（1+

11、cos（pi*Ts/alpha*（abs（f（ii）-（1-alpha）/（2*Ts）;ii=find（abs（f）=（1-alpha）/（2*Ts）;hcos（ii）=Ts;%根升余弦hrcos=sqrt（hcos）;2000a=sign（randn（1,M）; %产生序列s1=zeros（1,N）;s1（1:L:N）=a*fs; %冲击序列S1=t2f（s1,fs）;S2=S1.*hrcos;s2=real（f2t（S2,fs）; %发送的 PAM信号P=abs（S2）.2/T; %累计平均%白高斯噪声end nw=sqrt（0.01*Bs）*randn（1,N）;r=s2+nw;R=t2

12、f（r,fs）;Y=R.*hrcos;y=real（f2t（Y,fs）; %采样前信号 figure（1）plot（f,EP）f（kHZ） 功率谱（ W/kHz ） axis（-1.5,1.5,0,max（EP）;eyediagram（y,3*L,3,9）edAm从发送功率谱的根升余弦功率谱可以看出，边缘比较陡峭，截止频率约为0 -1.51 .9 .8 .7 .6 .5 .4 .3 .2 00000000）ZHk/2V （谱 率功-0.50.51.5=0.25 的根升余弦发送功率谱Eye Diagram-1.5Time （ms）0 f（kHZ）接收眼图1.25 ，符合W=（ +1） *Rs/2 的公式，图形与理论基本相符。眼图噪声容限约为 1 ，张开较大，斜率较大，说明对定时误差的灵敏度较高，存在一定0，与理论相符。