通原软件实验报告

目录

实验一调幅信号波形频谱仿真 (1)

一、实验要求 (1)

二、基本原理 (1)

1、AM调制原理 (1)

2、DSB-SC调制原理 (1)

3、SSB调制原理 (1)

三、程序框图 (2)

四、源代码 (2)

五、实验结果及分析 (5)

实验二调频信号波形频谱仿真 (9)

一、实验题目 (9)

二、基本原理 (10)

三、程序框图 (10)

四、仿真源代码 (10)

六、实验结果及分析 (12)

实验三单双极性归零码波形及功率谱仿真 (14)

一、实验题目 (14)

二、基本原理 (14)

1、单极性归零码 (14)

2、双极性归零码 (14)

3、各种码的比较 (14)

三、程序框图 (14)

四、仿真源代码 (15)

五、实验结果及分析 (22)

实验四根升余弦滚降功率谱密度及眼图仿真 (31)

一、实验题目 (31)

二、基本原理 (31)

1、升余弦滚降 (31)

2、眼图 (31)

三、程序框图 (32)

四、仿真源代码 (32)

五、实验结果及分析 (33)

选做题数字基带系统仿真 (37)

一、实验题目 (37)

二、基本原理 (37)

1、PAM信号的产生 (37)

2、常用码型 (38)

3、信道建模 (38)

4、眼图特性 (39)

5、取样判决 (39)

三、仿真思路 (40)

四、程序框图 (41)

五、仿真源码 (42)

六、实验结果及分析 (44)

附：代码中用到的变换函数代码 (48)

实验心得体会 (49)

实验一 调幅信号波形频谱仿真

一、实验要求

假设基带信号为()sin(2000)2cos(1000)m t t t ππ=+，载波频率为20kHz ，请仿真出AM 、DSB-SC 、SSB 信号，观察已调信号的波形及频谱。

二、基本原理

1、AM 调制原理 对于单音频信号

()sin(2)m m m t A f t π=

进行AM 调制的结果为

()(sin(2))sin 2(1sin(2))sin 2AM c m m c c m c s t A A A f t f t A A a f t f t ππππ=+=+

其中调幅系数m

A a A

=

，要求1a ≤以免过调引起包络失真。 由max A 和min A 分别表示AM 信号波形包络最大值和最小值，则AM 信号的调幅系数为

max min

max min

A A a A A -=

+

2、DSB-SC 调制原理

DSB 信号的时域表达式为

()()cos DSB c s t m t t ω=

频域表达式为

1

()[()()]2

DSB c c S M M ωωωωω=-++

3、SSB 调制原理

SSB 信号只发送单边带，比DSB 节省一半带宽，其表达式为：

1

1

()cos cos sin 2

2

ssB m m m m c s t A t w t A t w t =

三、程序框图

四、源代码

AM 信号：

global dt df t f T N fs =800;%khz T =200;%ms N =T *fs ; dt =1/fs ;

t =[-T /2:dt :T /2]; df =1/T ;

f =[-fs /2:df :fs /2];

fm1=1;%khz fm2=0.5;%khz fc =20;%khz

m =sin (2*pi *fm1*t )+2*cos (2*pi *fm2*t );%基带信号 a =cos (2*pi *fc *t );%载波

s =(1+m ).*a ;%已调信号 S =t2f (s ,fs );%已调信号频谱

figure (1) plot (t ,s )

axis ([0,3.5,-5,5])

xlabel (('t(ms)'));%横坐标标签 ylabel ('s(t)(V)');%纵坐标标签 figure (2) plot (f ,abs (S ))

axis ([-25,25,min (abs (S )),max (abs (S ))]) xlabel ('f(kHZ)'); %横坐标标签

产生载波和调制信号m(t)

SSB 、DSB 、AM 信号表达式

FFT 变换得各调制信号频谱

作图

ylabel('S(f)(V/kHZ)');%纵坐标标签

grid on;

DSB-SC

global dt df t f T N

fs=800;%khz

T=200;%ms

N=T*fs;

dt=1/fs;

t=[-T/2:dt:T/2];

df=1/T;

f=[-fs/2:df:fs/2];

fm1=1;%khz

fm2=0.5;%khz

fc=20;%khz

m=sin(2*pi*fm1*t)+2*cos(2*pi*fm2*t);%基带信号a=cos(2*pi*fc*t);

s=m.*a;%已调信号

S=t2f(s,fs);%已调信号频谱

figure(1)

plot(t,s)

axis([0,3.5,-5,5])

xlabel(('t(ms)'));%横坐标标签

ylabel('s(t)(V)');%纵坐标标签

figure(2)

plot(f,abs(S))

axis([-25,25,0,max(abs(S))])

xlabel('f(kHZ)');%横坐标标签

ylabel('S(f)(V/kHZ)');%纵坐标标签

SSB信号

global dt df t f T N

fs=800;%khz

T=200;%ms

N=T*fs;