数字基带信号的功率谱计算

(+A)(-A) (+A)(-A) (+A*A)(+A*A)
0(0)
19
AMI码中出现an an+2组合的概率
|an| |an+1| |an+2| |an*an+2| p(anan+2)
00
0
0
0.125
00
A
0
0.125
0A
0
0
0.125
0A
A
0
0.125
A0
0
0
0.125
A0
A (-A*A) 0.125
• Na=4;
%示波器扫描宽度为4个码元
• SNR=30;
%信噪比为30dB
• ratio=0;
%误码率初始值设为0
• t=[-T/2+dt/2:dt:T/2];
%时域横坐标,－T/2＋dt的目的是对开分母为0的那个点
• f=[-Bs+df/2:df:Bs];
%频域横坐标
• EP=zeros(1,N);
S( k TS
)
2E2 (a) TS2| G(ຫໍສະໝຸດ k TS) |2(f
n TS
)
k是从负无穷到正无穷的整数。
当k=0，得到信号的直流成分。
10
例9-1 单极性二元码的功率谱计算。
假设单极性二元码中对应于输入信码0，1的 幅度取值为0，+A，输入信码为各态历经随机序列， 0，1的出现统计独立，则概率为1/2，即
C e j 2mfst m
m
其中Cm
1 TS
Ts
2 Ts
v(t )e
dt j 2mfst
2
27
其中Cm
1 TS
Ts
2 Ts
v(t )e
j
2mfst
dt
2
1
TS
Ts
2 Ts
2
e
j 2mfst
[
n
pg1 (t
nTS
)
(1
P)g2
(t
nTS
)]dt
fs
nTs Ts / 2
[
nTs Ts / 2
( f ) | G( f ) |2 A2[ 1 2( 1 ) cos(2f Ts)
TS
2
4
| G( f ) |2 A2 sin2 (f Ts) / Ts
22
*9.2.2 一般情况下的随机信号功率谱
一个二元序列，指数字“1”“0” 分别用两种不同的波形表示，如 2ASK，2FSK，2PSK，数字双相码 等。
(
f
)
1 TS 2
|
PG1(0)
(1
P)G2 (0)
|2
(
f
)
2 TS 2
m1
|
P
G1
(
m TS
)
(1
P)G2
(
m TS
)
|2
(f
m) TS
P(1 P) | G1() G2 () |2 (3)
TS
30
例9-3 双极性非归零二元码的功率谱计算。 假设双极性非归零码的幅度取值为+A，-A。 出现概率为1/2，即g1(t)=-g2(t)，P=1/2。
R(k )
E(anank )
A2 4
（查表）
13
*连续谱的功率谱密度函数为：（代入公式（1））
s (
f
)
1 TS
| G( f
) |2
{R(0)
E2[a] 2 (R(k)
k 1
E 2[a]cos(2fTS )}
1 TS
| G( f
) |2
{ A2 2
A2 4
2 ( A2 G1 4
A2 4
0 an
A
p1 2
Q1 2
11
解：先做出下表
二元码中an an+1组合的出现概率
an an+k an*an+k P(an an+k)
00
0
0.25
0A
0
0.25
A0
0
0.25
A A A*A
0.25
查表得：
11
A
E(a)
E(an )
*0 2
*A 2
2
12
R(0)
E(anan )
A2 2
，R（0）是交流功率。
2
TS 2
t TS
g2 (t) g1 (t)
32
代入公式(3)得
( f )
A2
(1
2P)2
n ,n0
2
n
2
(
f
n) TS
4
P(1
P)[
s
in 4 (f TS (f TS )2
)
]
A2TS
如果P=1/2，则离散线谱消失，
(
f
)
[
sin4 (fTS (fTS )2
)
]A2TS
33
• BaseBandSystem.m
%码速率是2Mb/s
• Ts=1/Rb;
%码元间隔
• TS=1/3;
%系统带宽为(1+alpha)/(2*TS)MHz=2.99MHz
• dt=Ts/L;
%时域采样间隔
• df=1/(N*dt);
%频域采样间隔
• T=N*dt;
%截短时间
• Bs=N*df/2;
%系统带宽
• alpha=0.33;
%升余弦滤波器滚降系数0.5
pg1
(t )
(1
P)g
2
(t ) ]e
j 2mfs
(t nTs ) dt
n
fs
[
pg1
(t )
(1
P)g2
(t ) ]e
j 2mfst
dt
f s [ pG1 (mfs ) (1 p)G2 (mfs )]
28
*傅氏级数的系数就是离散线谱的幅度,对于功 率谱有:
( f )
| Cm |2 ( f mfs )
m
| f s [ pG1 (mfs ) (1 p)G2 (mfs )] |2 ( f mfs )
m
1 TS 2
|
PG1 (0) (1 P)G2 (0) |2
(f
)
2 TS 2
m 1
|
P
G1
(
m TS
)
(1
P)G2
(
m TS
)
|2
(
f
m) TS
最后简化为：
（2）
29
纯随机二元序列的功率谱包括直流分量、 离散线谱和连续功率谱三项
k 1
E2[a]cos(2fTS )}
平均功率谱密度计算式（1）
*功率谱的连续部分与单个脉冲功率谱的平方成正比。 式中：G(f)是单个波形g(t)的频域特性。 E(a)是系数的均值。
E[a] E[an] an
R(k)是相关值。
R(k) E{an, ank} anank
9
*它的线谱部分计算式如：
•
%+-+-+-+-输入信号及其功率谱密度-+-+-+-%
• %global dt df t f T N L
%全局变量
• clear;close all ; clc;
%清理数据空间
• k=16;
%取样点幂次
• N=2^k;
%采样点数
• L=64;
%每码元的采样点数
• M=N/L;
%码元数
• Rb=2;
2、对于占空比50%的归零信号，脉冲时宽为 非归零信号的一半，带宽就为非归零信号的一倍。
16
归零信号的离散线谱中，除直流分量外，还有 奇次线谱，没有偶次线谱，由于有基频分量fs，可以 提取位定时信号。
17
例9-2 AMI码的功率谱计算。 假设AMI码的三种幅度取值为-A，0，+A。
输入信码为各态历经随机过程，0，1的出现概率 统计独立，概率各为1/2，由AMI编码规律可知：
现代通信原理
第九章 数字信号的基带传输
1
§9.2 数字基带信号的功率谱计算
数字基带信号是随机信号，只能计算功率 谱密度。
计算功率谱密度不是件容易的事，下面列 举几种方法。
2
首先补充单个脉冲波形的频谱。 1.矩形脉冲
3
2.半余弦形脉冲
4
3.升余弦脉冲
5
4. 三角形脉冲
6
9.2.1 相同波形随机序列的功率谱
• s(ii+[0:M-1]*L)=a;
• end
35
• %传输信道中加入高斯白噪声 • s_noise = awgn(s_in,SNR,'measured'); % Add white Gaussian noise. • Noise=T2F(s_noise,N,dt); • S_out=Noise.*S_RCos; • s_out=real(F2T(S_out,N,dt)); • temp=round(s_out); • Temp=T2F(temp,N,dt); • P=Temp.*conj(Temp)/T; • EP_out=(EP_out*(jj-1)+P)/jj; • [error_number,ratio]=biterr(s,round(s_out)); • ratio=(ratio*(jj-1)+ratio)/jj; • end • • %显示单极性不归零码序列 • disp('按任意键显示输入信号的相关信息');pause; • figure(1) • set(1,'position',[10,50,300,200]); %设定窗口位置及大小 • figure(2) • set(2,'position',[350,50,300,200]);%设定窗口位置及大小