北邮通原软件实验

实验一

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

1．AM信号：

（1）信号的表达式

（3）流程图

AM信号

s= (1+0.3*m).*cos(2*pi*fc*t);

傅氏变换S= t2f(s,fs)

绘制时域波形及频谱

（2）源代码

%AM信号的产生

fs= 800; %采样频率KHz

T= 200; %截短时间ms

N= T*fs; %采样点数

dt= 1/fs;

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

df= 1/T;

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

fm= 1; % kHz

fc= 20; % kHz

m= sin(2*pi*fm*t)+2*cos(1*fm*pi*t);

s= (1+0.3*m).*cos(2*pi*fc*t); %AM 信号

S= t2f(s,fs);

figure(1)

plot(f,abs(S1))

title('AM信号频谱')

xlabel('f')

ylabel('S(f)')

axis([-25,25,0,max(abs(S1))]);

%xset('window',2)figure(2)

plot(t,s1)

title('AM信号波形')

xlabel('t')

ylabel('s(t)')

axis([-3,3,-3,3]);

（4）实验结果

-3

-2-1

0123

-3-2

-1

1

2

3

AM 信号波形

t(ms)

s (t )

-25

-20

-15

-10

-5

05

10

15

20

25

0102030405060708090

100AM 信号频谱

f(kHz)

S (f )

2．DSB-SC信号

（1）信号的产生和表达式

（2）流程图

DSC-SB信号

s= m.*cos(2*pi*fc*t);

傅氏变换S= t2f(s,fs)

绘制时域波形及频谱

（3）源代码

fs= 800; %KHz

T= 200; %ms

N= T*fs;

dt= 1/fs;

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

df= 1/T;

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

fm= 1; % kHz

fc= 20; % kHz

m= sin(2*pi*fm*t)+2*cos(1*fm*pi*t);

s= m.*cos(2*pi*fc*t); %DSB-SC 信号S= t2f(s,fs);

figure(1)

plot(f,abs(S2))

title('DSB-SC信号频谱')

xlabel('f')

ylabel('S(f)')

axis([-25,25,0,max(abs(S2))]);

figure(2)

plot(t,s2)

title('DSB-SC信号波形')

xlabel('t')

ylabel('s(t)')

axis([-1,4,-3,3]);

（4）实验结果

()()()()()

cos

c c c

s t m t c t m t A t

ωϕ

==+

-1

-0.500.51

1.52

2.53

3.54

-3-2

-1

1

2

3

DSB-SC 信号波形

t(ms)

s (t )

-25

-20

-15

-10

-5

05

10

15

20

25

DSB-SC 信号频谱

f(kHz)

S (f )

2 DSC-SB 频谱

3．SSB 信号

（1）信号的产生和表达式

（2）流程图

（3）源代码：

%SSB 信号的产生 fs= 800; %KHz T= 200; %ms N= T*fs; dt= 1/fs;

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

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

fm= 1; % kHz fc= 20; % kHz m= sin(2*pi*fm*t)+2*cos(1*fm*pi*t); 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 signal S= t2f(s,fs); figure(1)

plot(f,abs(S3))

title('SSB 信号频谱') xlabel('f') ylabel('S(f)')

axis([-25,25,0,max(abs(S3))]) figure(2) plot(t,s3)

title('SSB 信号波形') xlabel('t') ylabel('s(t)') axis([0,6,-3,3]) ()()

()

()()()

()

()cos 2sin 2cos 2sin 2DSB SC AM c c c c SSB c c c c s t s t s t A m t f t A m t f t

s t A m t f t A m t f t

ππππ--=-=-=下上