Matlab已调带通功率谱密度,MATLAB通信原理仿真五章作业.doc

通信原理仿真作业

第五章 模拟调制

1. AM、DSB调制及解调

用matlab产生一个频率为1Hz，功率为1的余弦信源，设载波频率，,试画出：

l AM及DSB调制信号的时域波形；

l 采用相干解调后的AM及DSB信号波形；

l AM及DSB已调信号的功率谱；

l 调整载波频率及m0，观察分的AM的过调与DSB反相点现象。

l 在接收端带通后加上窄带高斯噪声，单边功率谱密度，重新解调。

2. SSB调制及解调

用matlab产生一个频率为1Hz，功率为1的余弦信源，设载波频率，,试画出：

l SSB调制信号的时域波形；

l 采用相干解调后的SSB信号波形；

l SSB已调信号的功率谱；

l 在接收端带通后加上窄带高斯噪声，单边功率谱密度，重新解调。

3. FM调制及解调

设输入信号为，载波中心频率为，VCO的压控振荡系数为，载波平均功率为1W。试画出：

l 已调信号的时域波形；

l 已调信号的振幅谱；

l 用鉴频器解调该信号，并与输入信号比较。

加入相同功率的高斯白噪声信号，

DSB

信号解调较

AM

信号解调误差较大。

二代码

clc;

closeall;

clearall;

fm=1;

fc=10;

T=5;

t=linspace(-5,5,1024);%

采样频率

a=sqrt(2);

mt=a*cos(2*pi*fm*t);%

带通滤波器

A=2;

s_db=mt.*cos(2*pi*fc*t);

s_am=(A+mt).*cos(2*pi*fc*t);

figure(1);

subplot(2,1,1);

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

-1

0

1

t/s

m

(

t

)

AM

相

干

解

调

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

-1

0

1

t/s

m

(

t

)

DSB

相

干

解

调

plot(t,s_am);

holdon;

plot(t,A+mt,r--);

title(AM

调制信号

);

xlabel(t/s);

ylabel(

幅度

);

subplot(2,1,2);

plot(t,s_db);

title(DSB

调制信号

);

xlabel(t/s);

ylabel(

幅度

);

holdon;

plot(t,mt,r--);

sp=s_am.*cos(2*pi*fc*t);

sp2=s_db.*cos(2*pi*fc*t);

fs1=0.3;fp1=0.5;fpu=1.8;fsu=2.0;Fs=50;Rp=1;Rs=42;

ws1=2*pi*fs1/Fs;wp1=2*pi*fp1/Fs;wpu=2*pi*fpu/Fs;wsu=2

*pi*fsu/Fs;

wd1=0.5*(wp1+ws1);wd2=0.5*(wpu+wsu);

bt=min(abs(ws1-wp1),abs(wsu-wpu));

n=ceil(6.6*pi/bt);

wn=[wd1/pi,wd2/pi];

hn=fir1(n-1,wn,bandpass,hamming(n));

y=conv(sp,hn);

y2=conv(sp2,hn);

y=2*y;

y2=2*y2;

figure(2);

subplot(2,1,1);

plot(t,y(1:1024));

holdon;

plot(t,mt,r--);

axis([05-1.51.5]);

xlabel(t/s),ylabel(m(t)),title(AM

相干解调

);

subplot(2,1,2);

plot(t,y2(1:1024));

holdon;

plot(t,mt,r--);

axis([05-1.51.5]);

xlabel(t/s),ylabel(m(t)),title(DSB

相干解调

AM

、

DSB

调制及解调

用

matlab

产生一个频率为

1Hz

，

功率为

1

的余弦信源

(

)

m

t

，

设载波频率

10

c

Hz

w

=

，

0

2

m

=

,

试画出：

l

AM

及

DSB

调制信号的时域波形；

l

采用相干解调后的

AM

及

DSB

信号波形；

l

AM

及

DSB

已调信号的功率谱；

l

调整载波频率及

m0

，观察分的

AM

的过调与

DSB

反相点现象。

l

在接收端带通后加上窄带高斯噪声，

单边功率谱密度

0

0.1

n

=

，

重新解调。

一图表输出及结果分析

1

AM

及

DSB

调制信号的时域波形

AM

信号的包络波形与其调制信号的波形完全一样。

DSB

调制信号在载波处反向。

2

采用相干解调后的

AM

及

DSB

信号波形

AM

信号可用包络检波器检波，滤除直流后可恢复出原信号。与

AM

相比，

DSB

信

号由于不存在载波分量，全部功率都用来传输信息。

-5

-4

-3

-2

-1

0

1

2

3

4

5

-5

0

5

AM

调

制

信

号

t/s

幅

度

-5

-4

-3

-2

-1

0

1

2

3

4

5

-2

0

2

DSB

调

制

信

号

t/s

幅

度

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

-1

0

1

t/s

m

(

t

)

AM

相

干

解

调

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

-1

0

1

t/s

m

(

t

)

DSB

相

干

解

调

3

AM

及

DSB

已调信号的功率谱

AM

信号所需的传输带宽为调制信号带宽的两倍。

DSB

节省了载波功率，但所需

的传输带宽仍为调制信号带宽的两倍。

4

调整载波频率及

m0

，观察分的

AM

的过调与

DSB

反相点现象。

AM

信号用于包络检波器检波时，当调制信号的最大值大于直流信号时，会出现

“过调制”现象，此时会发生失真。

DSB

相位在

0.2

π

角频率处频率处明显存在反相点现象。

5

在接收端带通后加上窄带高斯噪声，单边功率谱密度

0

0.1

n

=

，重新解调。

0

20

40

60

80

100

120

0

500

1000

AM

已

调

信

号

功

率

谱

f/hz

w

0

20

40

60

80

100

120

0

100

200

DSB

已

调

信

号

功

率

谱

f/hz

w

-5

-4

-3

-2

-1

0

1

2

3

4

5

-5

0

5

AM

不

过

调

t/s

幅

度

-5

-4

-3

-2

-1

0

1

2

3

4

5

-5

0

5

AM

过

调

t/s

幅

度

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0

200

400

DSB

调

幅

波

w/(2*pi)

幅

度

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

-5

0

5

w/(2*pi)

相

位

加入相同功率的高斯白噪声信号，

DSB

信号解调较

AM

信号解调误差较大。

二代码

clc;

closeall;

clearall;

fm=1;

fc=10;

T=5;

t=linspace(-5,5,1024);%

采样频率

a=sqrt(2);

mt=a*cos(2*pi*fm*t);%

带通滤波器

A=2;

s_db=mt.*cos(2*pi*fc*t);

s_am=(A+mt).*cos(2*pi*fc*t);

figure(1);

subplot(2,1,1);

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

-1

0

1

t/s

m

(

t

)

AM

相

干

解

调

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

-1

0

1

t/s

m

(

t

)

DSB

相

干

解

调

plot(t,s_am);

holdon;

plot(t,A+mt,r--);

title(AM

调制信号

);

xlabel(t/s);

ylabel(

幅度

);

subplot(2,1,2);

plot(t,s_db);

title(DSB

调制信号

);

xlabel(t/s);

ylabel(

幅度

);

holdon;

plot(t,mt,r--);

sp=s_am.*cos(2*pi*fc*t);

sp2=s_db.*cos(2*pi*fc*t);

fs1=0.3;fp1=0.5;fpu=1.8;fsu=2.0;Fs=50;Rp=1;Rs=42;

ws1=2*pi*fs1/Fs;wp1=2*pi*fp1/Fs;wpu=2*pi*fpu/Fs;wsu=2

*pi*fsu/Fs;

wd1=0.5*(wp1+ws1);wd2=0.5*(wpu+wsu);

bt=min(abs(ws1-wp1),abs(wsu-wpu));

n=ceil(6.6*pi/bt);

wn=[wd1/pi,wd2/pi];

hn=fir1(n-1,wn,bandpass,hamming(n));

y=conv(sp,hn);

y2=conv(sp2,hn);

y=2*y;

y2=2*y2;

figure(2);

subplot(2,1,1);

plot(t,y(1:1024));

holdon;

plot(t,mt,r--);

axis([05-1.51.5]);

xlabel(t/s),ylabel(m(t)),title(AM

相干解调

);

subplot(2,1,2);

plot(t,y2(1:1024));

holdon;

plot(t,mt,r--);

axis([05-1.51.5]);

xlabel(t/s),ylabel(m(t)),title(DSB

相干解调

);

S=fft(s_am,1024);

p=S.*conj(S)/1024;

S2=fft(s_db,1024);

16

展开阅读全文