文章目录
- 通信系统一般模型
- 香农公式
- 内积、能量与许瓦兹不等式
- 周期信号的傅里叶级数:
- 傅里叶变换
- 能量信号的能量及其能量谱密度
- 功率信号的功率及其功率谱密度
- 确定信号的相关函数
- 确定信号通过线性系统
- 希尔伯特变换和解析信号
- 随机变量的统计特征
- 随机过程的统计特征
- 平稳随机过程
- 平稳随机过程的功率谱密度与自相关函数的关系
- 高斯过程的定义及性质
- 加性高斯白噪声(AWGN)
- 平稳随机过程通过线性系统
- 匹配滤波器
- 窄带高斯白噪声
- 窄带高斯白噪声 n c ( t ) n_c(t) nc(t)和 n s ( t ) n_s(t) ns(t)的统计特性
- 窄带高斯白噪声 a ( t ) a(t) a(t)和 ϕ ( t ) \phi(t) ϕ(t)的统计特性
- 循环平稳随机过程
- 幅度调制系统的原理
- Q函数和互补误差函数
Author🚹:CofCai
Email✉️:cai.dongjun@nexuslink.cn
QQ😙:1664866311
personal Page🏠
优快云 Page🏠
通信系统一般模型
模拟通信系统时:
- 发送机:调制
- 接收机:解调
数字通信系统时:
- 发送机:信源编码、信道编码、调制
- 接收机:解调、信道译码、信源译码
数字[模拟]通信系统的主要性能指标:
- 有效性[传输带宽]
- 可靠性[信噪比]
香农公式
C = B log 2 ( 1 + S N ) C=B\log_2^{(1+\frac{S}{N})} C=Blog2(1+NS)
内积、能量与许瓦兹不等式
令信号x(t),y(t)为任意信号,一般情况可以是时间的复函数,称
<
x
,
y
>
=
∫
−
∞
+
∞
x
∗
(
t
)
y
(
t
)
d
t
<x,y>=\int_{-\infin}^{+\infin}x^*(t)y(t)d_t
<x,y>=∫−∞+∞x∗(t)y(t)dt
为x(t),y(t)的内积、且
∫
−
∞
+
∞
x
∗
(
t
)
y
(
t
)
d
t
=
∫
−
∞
+
∞
X
∗
(
f
)
Y
(
f
)
d
f
\int_{-\infin}^{+\infin}x^*(t)y(t)d_t=\int_{-\infin}^{+\infin}X^*(f)Y(f)d_f
∫−∞+∞x∗(t)y(t)dt=∫−∞+∞X∗(f)Y(f)df
信号与其自身的内积就是能量:
E
x
=
∫
−
∞
+
∞
∣
x
(
t
)
∣
2
d
t
E_x=\int_{-\infin}^{+\infin}|x(t)|^2d_t
Ex=∫−∞+∞∣x(t)∣2dt
周期信号的傅里叶级数:
f T ( t ) = a 0 2 + ∑ n = 1 ∞ ( a n cos ( n ω 0 t ) + b n sin ( n ω 0 t ) ) n = 0 , 1 , 2 , . . . 其 中 a n = 2 T ∫ t 0 t 0 + T f T ( t ) cos ( n ω 0 t ) d t n = 0 , 1 , 2 , . . . b n = 2 T ∫ t 0 t 0 + T f T ( t ) sin ( n ω 0 t ) d t n = 0 , 1 , 2 , . . . f_T(t)=\frac{a_0}{2}+\sum_{n=1}^{\infin}(a_n{\cos}(n{\omega_0}t)+b_n{\sin}({n\omega_0}t)){\qquad}n=0,1,2,...\\ 其中\\ a_n=\frac{2}{T}\int_{t_0}^{t_0+T}f_T(t)\cos(n{\omega_0}t)d_t{\qquad}n=0,1,2,...\\ b_n=\frac{2}{T}\int_{t_0}^{t_0+T}f_T(t)\sin(n{\omega_0}t)d_t{\qquad}n=0,1,2,...\\ fT(t)=2a0+n=1∑∞(ancos(nω0t)+bnsin(nω0t))n=0,1,2,...其中an=T2∫t0t0+TfT(t)cos(nω0t)dtn=0,1,2,...bn=T2∫t0t0+TfT(t)sin(nω0t)dtn=0,1,2,...
或者写成指数形式:
f
(
t
)
=
1
2
π
∫
−
∞
+
∞
∫
−
∞
+
∞
f
(
t
)
e
−
j
ω
t
d
t
e
j
ω
t
d
ω
F
(
j
w
)
=
∫
−
∞
+
∞
f
(
t
)
e
−
j
ω
t
d
t
(
1
)
f
(
t
)
=
1
2
π
∫
−
∞
+
∞
F
(
j
w
)
e
j
ω
t
d
ω
(
2
)
f(t)=\frac{1}{2\pi}\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin}f(t)e^{-j{\omega}t}d_te^{j{\omega}t}d_{\omega}\\ F(jw)=\int_{-\infin}^{+\infin}f(t)e^{-j{\omega}t}d_t{\qquad}(1)\\ f(t)=\frac{1}{2\pi}\int_{-\infin}^{+\infin}F(jw)e^{j{\omega}t}d_{\omega}{\qquad}(2)
f(t)=2π1∫−∞+∞∫−∞+∞f(t)e−jωtdtejωtdωF(jw)=∫−∞+∞f(t)e−jωtdt(1)f(t)=2π1∫−∞+∞F(jw)ejωtdω(2)
傅里叶变换
f ( t ) = 1 2 π ∫ − ∞ + ∞ ∫ − ∞ + ∞ f ( t ) e − j ω t d t e j ω t d ω F ( j w ) = ∫ − ∞ + ∞ f ( t ) e − j ω t d t ( 1 ) f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( j w ) e j ω t d ω ( 2 ) f(t)=\frac{1}{2\pi}\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin}f(t)e^{-j{\omega}t}d_te^{j{\omega}t}d_{\omega}\\ F(jw)=\int_{-\infin}^{+\infin}f(t)e^{-j{\omega}t}d_t{\qquad}(1)\\ f(t)=\frac{1}{2\pi}\int_{-\infin}^{+\infin}F(jw)e^{j{\omega}t}d_{\omega}{\qquad}(2) f(t)=2π1∫−∞+∞∫−∞+∞f(t)e−jωtdtejωtdωF(jw)=∫−∞+∞f(t)e−jωtdt(1)f(t)=2π1∫−∞+∞F(jw)ejωtdω(2)
某个域中的面积等于另一个域,即:
x
(
0
)
=
∫
−
∞
+
∞
X
(
f
)
d
f
X
(
0
)
=
∫
−
∞
+
∞
x
(
t
)
d
t
x(0)=\int_{-\infin}^{+\infin}X(f)d_f\\ X(0)=\int_{-\infin}^{+\infin}x(t)d_t
x(0)=∫−∞+∞X(f)dfX(0)=∫−∞+∞x(t)dt
能量信号的能量及其能量谱密度
E = ∫ − ∞ + ∞ f 2 ( t ) d t = . . . = ∫ − ∞ + ∞ ∣ F ( f ) ∣ 2 d f E ( f ) = ∣ F ( f ) ∣ 2 E=\int_{-\infin}^{+\infin}f^2(t)d_t=...=\int_{-\infin}^{+\infin}|F(f)|^2d_f\\ E(f)=|F(f)|^2 E=∫−∞+∞f2(t)dt=...=∫−∞+∞∣F(f)∣2dfE(f)=∣F(f)∣2
功率信号的功率及其功率谱密度
P ( f ) = l i m T − > ∞ ∣ F ( f ) ∣ 2 T F ( f ) 是 f ( t ) 的 频 谱 函 数 P = ∫ − ∞ + ∞ P ( f ) d f P(f)=lim_{T->\infin}\frac{|F(f)|^2}{T}\\ F(f)是f(t)的频谱函数\\ P=\int_{-\infin}^{+\infin}P(f)d_f P(f)=limT−>∞T∣F(f)∣2F(f)是f(t)的频谱函数P=∫−∞+∞P(f)df
确定信号的相关函数
能量信号x(t)与其时延X(t+
τ
\tau
τ)的内积称为x(t)的自相关函数:
R
x
(
τ
)
=
∫
−
∞
+
∞
x
∗
(
t
)
x
(
t
+
τ
)
d
t
因
此
有
E
x
(
f
)
=
∣
X
(
f
)
∣
2
=
∫
−
∞
+
∞
R
x
(
τ
)
e
−
j
2
π
f
τ
d
τ
R_x(\tau)=\int_{-\infin}^{+\infin}x^*(t)x(t+\tau)d_t\\ 因此有{\qquad}E_x(f)=|X(f)|^2=\int_{-\infin}^{+\infin}R_x(\tau)e^{-j{2\pi}f\tau}d_{\tau}
Rx(τ)=∫−∞+∞x∗(t)x(t+τ)dt因此有Ex(f)=∣X(f)∣2=∫−∞+∞Rx(τ)e−j2πfτdτ
确定信号通过线性系统
y ( t ) = ∫ − ∞ + ∞ x ( u ) h ( t − u ) d u Y ( f ) = H ( f ) X ( f ) y(t)=\int_{-\infin}^{+\infin}x(u)h(t-u)d_u\\ Y(f)=H(f)X(f) y(t)=∫−∞+∞x(u)h(t−u)duY(f)=H(f)X(f)
希尔伯特变换和解析信号
希尔伯特变换就是经过一个冲激响应为
1
π
t
\frac{1}{{\pi}t}
πt1的线性系统,它本质是一个理想的宽带相移器。
解析信号就是
z
(
t
)
=
f
(
t
)
+
j
f
^
(
t
)
z(t)=f(t)+j\hat{f}(t)
z(t)=f(t)+jf^(t),其中
f
^
(
t
)
\hat{f}(t)
f^(t)是
f
(
t
)
f(t)
f(t)的希尔伯特变换,即虚部是实部的希尔伯特变换。
随机变量的统计特征
单个随机变量的统计特征:
- 分布函数、概率密度函数
- 数字特征
- 数学期望(统计平均值、均值)
- 方差(偏离程度)。且有 D [ X ] = E [ X − E [ X ] ] 2 = E ( X 2 ) − E 2 ( X ) = δ X 2 D[X]=E[X-E[X]]^2=E(X^2)-E^2(X)=\delta_X^2 D[X]=E[X−E[X]]2=E(X2)−E2(X)=δX2
多个随机变量的统计特征
- 联合分布函数: F ( x , y ) = P ( X ≤ x , Y ≤ y ) F(x,y)=P(X{\le}x,Y{\le}y) F(x,y)=P(X≤x,Y≤y),如果 α F ( x , y ) α x α y = p ( x , y ) \frac{\alpha{F(x,y)}}{\alpha{x}\alpha{y}}=p(x,y) αxαyαF(x,y)=p(x,y)存在,则称其为随机变量X和Y的联合概率密度函数。
- 数字特征
- 数学期望
- E ( X + Y ) = E ( X ) + E ( Y ) E(X+Y)=E(X)+E(Y) E(X+Y)=E(X)+E(Y)
- E ( X Y ) = E ( X ) E ( y ) E(XY)=E(X)E(y) E(XY)=E(X)E(y),前提是变量X和Y相互独立
- 方差: D ( X + Y ) = D ( X ) + D ( Y ) D(X+Y)=D(X)+D(Y) D(X+Y)=D(X)+D(Y),前提是变量X和Y相互独立
- 相关函数: R ( X , Y ) = E ( X Y ) = ∫ − ∞ + ∞ ∫ − ∞ + ∞ x y p ( x , y ) d x d y R(X,Y)=E(XY)=\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin}xyp(x,y)d_xd_y R(X,Y)=E(XY)=∫−∞+∞∫−∞+∞xyp(x,y)dxdy
- 协方差函数: C ( X , Y ) = E { [ X − E ( X ) ] [ Y − E ( Y ) ] } = R ( X , Y ) − m X m Y C(X,Y)=E\{[X-E(X)][Y-E(Y)]\}=R(X,Y)-m_Xm_Y C(X,Y)=E{[X−E(X)][Y−E(Y)]}=R(X,Y)−mXmY
- 归一化协方差函数(相关系数): ρ X Y = C ( X , Y ) δ X δ Y {\rho}_{XY}=\frac{C(X,Y)}{{\delta_X}{\delta_Y}} ρXY=δXδYC(X,Y)
- 数学期望
随机过程的统计特征
一个随机过程的统计特征
- 数学期望(统计平均值): E [ X ( t ) ] = ∫ − ∞ + ∞ x p ( x ; t ) d t = m X ( t ) E[X(t)]=\int_{-\infin}^{+\infin}xp(x;t)d_t=m_X(t) E[X(t)]=∫−∞+∞xp(x;t)dt=mX(t)
- 方差: D [ X ( t ) ] = E [ X 2 ( t ) ] − E 2 [ X ( t ) ] = ∫ − ∞ + ∞ x 2 p ( x ; t ) d x − m X 2 ( t ) D[X(t)]=E[X^2(t)]-E^2[X(t)]=\int_{-\infin}^{+\infin}x^2p(x;t)d_x-m_X^2(t) D[X(t)]=E[X2(t)]−E2[X(t)]=∫−∞+∞x2p(x;t)dx−mX2(t)
- 自相关函数: R X ( t 1 , t 2 ) = E [ X ( t 1 ) X ( t 2 ) ] = ∫ − ∞ + ∞ ∫ − ∞ + ∞ x 1 x 2 p ( x 1 , x 2 ; t 1 , t 2 ) d x 1 d x 2 R_X(t_1,t_2)=E[X(t_1)X(t_2)]=\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin}x_1x_2p(x_1,x_2;t_1,t_2)d_{x_1}d_{x_2} RX(t1,t2)=E[X(t1)X(t2)]=∫−∞+∞∫−∞+∞x1x2p(x1,x2;t1,t2)dx1dx2
- 自协方差: C X ( t 1 , t 2 ) = E [ X ( t 1 ) X ( t 2 ) ] − m X ( t 1 ) m X ( t 2 ) C_X(t_1,t_2)=E[X(t_1)X(t_2)]-m_X(t_1)m_X(t_2) CX(t1,t2)=E[X(t1)X(t2)]−mX(t1)mX(t2)
- 归一化协方差函数: ρ X ( t 1 , t 2 ) = C X ( t 1 , t 2 ) δ X ( t 1 ) δ X ( t 2 ) {\rho}_X(t_1,t_2)=\frac{C_X(t_1,t_2)}{{\delta}_X(t_1){\delta}_X(t_2)} ρX(t1,t2)=δX(t1)δX(t2)CX(t1,t2)
两个随机过程的统计特征
……
平稳随机过程
严平稳随机过程定义:随机过程X(t)的任意n维分布函数或概率密度函数与时间起点无关。
宽平稳随机过程的定义:
- 均值为常数
- 自相关函数仅与时间间隔相关
各态经历性:x(t)是平稳随机过程X(t)的任意一个样本函数,满足 E [ X ( t ) ] = x ( t ) ‾ , E [ X ( t ) X ( t + τ ) ] = x ( t ) x ( t + τ ) ‾ E[X(t)]=\overline{x(t)},E[X(t)X(t+\tau)]=\overline{x(t)x(t+\tau)} E[X(t)]=x(t),E[X(t)X(t+τ)]=x(t)x(t+τ)
平稳随机过程的功率谱密度与自相关函数的关系
D e f : P X ( f ) = l i m T − > ∞ E [ ∣ X T ( f ) ∣ 2 ] T R X ( τ ) = ∫ − ∞ + ∞ P X ( f ) e j 2 π f τ d f P X ( f ) = ∫ − ∞ + ∞ R X ( τ ) e − j 2 π f τ d τ Def:P_X(f)=lim_{T->\infin}\frac{E[|X_T(f)|^2]}{T}\\ R_X(\tau)=\int_{-\infin}^{+\infin}P_X(f)e^{j2{\pi}f\tau}d_f\\ P_X(f)=\int_{-\infin}^{+\infin}R_X(\tau)e^{-j2{\pi}f\tau}d_{\tau}\\ Def:PX(f)=limT−>∞TE[∣XT(f)∣2]RX(τ)=∫−∞+∞PX(f)ej2πfτdfPX(f)=∫−∞+∞RX(τ)e−j2πfτdτ
高斯过程的定义及性质
随机过程的任意n维分布都是正态分布,它完全取决于n个随机变量X(t1),X(t2),…,X(tn)的数学期望、方差和两两之间的归一化协方差函数。性质如下:
- 宽平稳和严平稳是等价的
- 任何两个时刻的随机变量统计独立和不相关是等价的
- 若干个高斯过程之和仍是高斯过程
- 经过线性系统后,输出仍是高斯过程
加性高斯白噪声(AWGN)
白噪声:功率谱密度均匀分布在整个频域范围内
- 其双边功率谱表示为: P n ( f ) = N 0 / 2 ( − ∞ < f < + ∞ ) P_n(f)=N_0/2{\qquad}(-\infin<f<+\infin) Pn(f)=N0/2(−∞<f<+∞)
- 单变功率谱表示为: P n ( f ) = N 0 ( 0 ≤ f < + ∞ ) P_n(f)=N_0{\qquad}(0{\le}f<+\infin) Pn(f)=N0(0≤f<+∞)
[加性]高斯白噪声:指即服从高斯分布、功率谱密度又是均匀分布的噪声。(一般都是指零均值加性高斯白噪声)
平稳随机过程通过线性系统
什么是输入?:平稳随机过程x(t)或X(t)
输出是什么?:
y
(
t
)
=
x
(
t
)
∗
h
(
t
)
=
∫
−
∞
+
∞
x
(
t
−
τ
)
h
(
τ
)
d
τ
o
r
Y
(
t
)
=
X
(
t
)
∗
h
(
t
)
=
∫
−
∞
+
∞
X
(
t
−
τ
)
h
(
τ
)
d
τ
y(t)=x(t)*h(t)=\int_{-\infin}^{+\infin}x(t-\tau)h(\tau)d_{\tau}\\ or\\ Y(t)=X(t)*h(t)=\int_{-\infin}^{+\infin}X(t-\tau)h(\tau)d_{\tau}
y(t)=x(t)∗h(t)=∫−∞+∞x(t−τ)h(τ)dτorY(t)=X(t)∗h(t)=∫−∞+∞X(t−τ)h(τ)dτ
随机过程**Y(t)的数学期望:
E
[
X
(
t
)
]
=
m
x
H
(
0
)
E[X(t)]=m_xH(0)
E[X(t)]=mxH(0)
随机过程Y(t)的自相关函数:
R
Y
(
t
1
,
t
1
+
τ
)
=
∫
−
∞
+
∞
∫
−
∞
+
∞
h
(
α
)
h
(
β
)
R
X
(
τ
+
α
−
β
)
d
α
d
β
R_Y(t_1,t_1+\tau)=\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin}h(\alpha)h(\beta)R_X(\tau+\alpha-\beta)d_{\alpha}d_{\beta}
RY(t1,t1+τ)=∫−∞+∞∫−∞+∞h(α)h(β)RX(τ+α−β)dαdβ
随机过程Y(t)**的功率谱密度:
P
Y
(
f
)
=
∣
H
(
f
)
∣
2
P
X
(
f
)
P_Y(f)=|H(f)|^2P_X(f)
PY(f)=∣H(f)∣2PX(f)
匹配滤波器
定义:采用一种线性滤波器,当信号和加性高斯白噪声通过它时,在某一特定抽样时刻,输出信噪比(信号瞬时功率与噪声平均功率之比)达到最大,此时这种线性滤波器就是匹配滤波器。
输入是什么?信号+噪声,即
s
(
t
)
+
n
W
(
t
)
s(t)+n_W(t)
s(t)+nW(t),且噪声是加性高斯白噪声。
目标是什么?找到输出信噪比(信号瞬时功率与噪声平均功率之比)最大的线性滤波器。
怎么找?如下:
s
o
(
t
)
=
∫
−
∞
+
∞
S
o
(
f
)
e
j
2
π
f
t
d
f
=
∫
−
∞
+
∞
S
(
f
)
H
(
F
)
e
j
2
π
f
t
d
f
P
n
o
(
f
)
=
P
i
o
(
f
)
∣
H
(
f
)
∣
2
=
N
0
2
∣
H
(
f
)
∣
2
m
i
n
S
N
R
o
=
∣
s
o
(
t
)
∣
2
P
n
o
(
f
)
g
o
t
H
(
f
)
=
K
S
∗
(
f
)
e
−
j
2
π
f
t
0
a
l
s
o
h
(
t
)
=
K
s
(
t
0
−
t
)
s_o(t)=\int_{-\infin}^{+\infin}S_o(f)e^{j2{\pi}ft}d_f=\int_{-\infin}^{+\infin}S(f)H(F)e^{j2{\pi}ft}d_f\\ P_{no}(f)=P_{io}(f)|H(f)|^2=\frac{N0}{2}|H(f)|^2\\ min{\quad}SNR_o=\frac{|s_o(t)|^2}{P_{no}(f)}\\ got{\quad}H(f)=KS^*(f)e^{-j2{\pi}ft0}\\ also{\quad}h(t)=Ks(t_0-t)
so(t)=∫−∞+∞So(f)ej2πftdf=∫−∞+∞S(f)H(F)ej2πftdfPno(f)=Pio(f)∣H(f)∣2=2N0∣H(f)∣2minSNRo=Pno(f)∣so(t)∣2gotH(f)=KS∗(f)e−j2πft0alsoh(t)=Ks(t0−t)
匹配滤波器(相关器)的输出:
s
o
(
t
)
=
K
R
(
t
−
t
0
)
s_o(t)=KR(t-t_0)
so(t)=KR(t−t0),有
s
o
(
t
0
)
=
K
R
(
0
)
=
K
E
s_o(t_0)=KR(0)=KE
so(t0)=KR(0)=KE
匹配滤波器输出在抽样时刻的一维概率密度函数:
输
出
信
号
瞬
时
值
:
s
o
(
t
0
)
=
K
R
(
0
)
=
K
E
输
出
噪
声
平
均
功
率
:
P
n
o
=
N
0
2
∫
−
∞
+
∞
∣
H
(
f
)
∣
2
d
f
=
K
2
N
0
E
2
输出信号瞬时值:s_o(t_0)=KR(0)=KE\\ 输出噪声平均功率:P_{no}=\frac{N_0}{2}\int_{-\infin}^{+\infin}|H(f)|^2d_f=\frac{K^2N_0E}{2}
输出信号瞬时值:so(t0)=KR(0)=KE输出噪声平均功率:Pno=2N0∫−∞+∞∣H(f)∣2df=2K2N0E
上面的我都明白,但是不明白为啥
K
E
KE
KE就是均值,
P
n
o
P_{no}
Pno就是方差。
窄带高斯白噪声
什么是窄带高斯白噪声:当高斯白噪声通过带通滤波器后,输出为带通高斯白噪声,又因为在通信系统中,带通噪声一般满足“窄带”的假设,所以带通噪声又称为窄带噪声。它的原理框图如下:
功率谱和时域波形如下:
此处请插入书上第58页的图!!!
可以看出,窄带高斯白噪声近似为正弦波,频率近似为
f
c
f_c
fc,包络缓慢变换,故我们可以设窄带高斯白噪声的时域表达式为:
n
(
t
)
=
a
(
t
)
cos
[
2
π
f
c
t
+
ϕ
(
t
)
]
a
(
t
)
≥
0
n
(
t
)
=
a
(
t
)
cos
ϕ
(
t
)
cos
2
π
f
c
t
−
a
(
t
)
sin
ϕ
(
t
)
sin
2
π
f
c
t
=
n
c
(
t
)
cos
2
π
f
c
t
−
n
s
(
t
)
sin
2
π
f
c
t
n(t)=a(t){\cos[2{\pi}f_ct+{\phi}(t)]}{\quad}a(t){\ge}0\\ n(t)=a(t){\cos{\phi(t)}{\cos{2{\pi}f_ct}}}-a(t){\sin{\phi(t)}{\sin{2{\pi}f_ct}}}\\ =n_c(t){\cos{2{\pi}f_ct}}-n_s(t){\sin{2{\pi}f_ct}}
n(t)=a(t)cos[2πfct+ϕ(t)]a(t)≥0n(t)=a(t)cosϕ(t)cos2πfct−a(t)sinϕ(t)sin2πfct=nc(t)cos2πfct−ns(t)sin2πfct
其中
n
c
(
t
)
=
a
(
t
)
cos
ϕ
(
t
)
n_c(t)=a(t)\cos{{\phi}(t)}
nc(t)=a(t)cosϕ(t)为同相分量,
n
s
(
t
)
=
a
(
t
)
sin
ϕ
(
t
)
n_s(t)=a(t)\sin{\phi(t)}
ns(t)=a(t)sinϕ(t)为正交分量,所以一个窄带高斯白噪声最重要的就是:包络、相位,从而可以得到同相分量、正交分量。
窄带高斯白噪声 n c ( t ) n_c(t) nc(t)和 n s ( t ) n_s(t) ns(t)的统计特性
通过证明,均值为0,方差为 δ n 2 \delta^2_n δn2的平稳窄带高斯白噪声 n ( t ) n(t) n(t)的同相分量和正交分量有如下性质:
-
同相分量和正交分量也是平稳高斯随机过程
-
同相分量和正交分量的均值也为0,方差(平均功率)均等于窄带高斯白噪声的平均功率(方差),即:
E [ n c ( t ) ] = E [ n s ( t ) ] = E [ n ( t ) ] = 0 D [ n c ( t ) ] = D [ n s ( t ) ] = D [ n ( t ) ] = δ n 2 E[n_c(t)]=E[n_s(t)]=E[n(t)]=0\\ D[n_c(t)]=D[n_s(t)]=D[n(t)]=\delta^2_n E[nc(t)]=E[ns(t)]=E[n(t)]=0D[nc(t)]=D[ns(t)]=D[n(t)]=δn2 -
同相分量和正交分量具有相同的自相关函数和功率谱函数
窄带高斯白噪声 a ( t ) a(t) a(t)和 ϕ ( t ) \phi(t) ϕ(t)的统计特性
可以证明,一个均值为0,方差为 δ n 2 \delta^2_n δn2的平稳窄带高斯白噪声 n ( t ) n(t) n(t),其包络 a ( t ) a(t) a(t)的一维概率密度服从瑞利分布;其相位 ϕ ( t ) \phi(t) ϕ(t)的一维概率密度函数服从均匀分布。
循环平稳随机过程
定义:随机过程s(t)的均值和自相关函数在时间上都是周期函数。
平稳随机过程X(t)和 cos ( 2 π f c t ) \cos(2{\pi}f_ct) cos(2πfct)同时通过乘法器,输出响应 Y ( t ) = X ( t ) cos ( 2 π f c t ) Y(t)=X(t)\cos(2{\pi}f_ct) Y(t)=X(t)cos(2πfct)是循环平稳随机过程,原因如下:
-
均值: E [ Y ( t ) ] = . . . = m X cos ( 2 π f c t ) E[Y(t)]=...=m_X\cos(2{\pi}f_ct) E[Y(t)]=...=mXcos(2πfct)
-
自相关函数:
R Y ( t , t + τ ) = R X ( τ ) 2 cos ( 2 π f c τ ) + R X ( τ ) 2 cos ( 4 π f c t + 2 π f c τ ) R_Y(t,t+\tau)=\frac{R_X(\tau)}{2}\cos(2{\pi}f_c\tau)+\frac{R_X(\tau)}{2}\cos(4{\pi}f_ct+2{\pi}f_c\tau) RY(t,t+τ)=2RX(τ)cos(2πfcτ)+2RX(τ)cos(4πfct+2πfcτ)
输出过程的功率谱密度:
P
Y
(
f
)
=
1
4
[
P
X
(
f
+
f
c
)
+
P
X
(
f
−
f
c
)
]
P_Y(f)=\frac{1}{4}[P_X(f+f_c)+P_X(f-f_c)]
PY(f)=41[PX(f+fc)+PX(f−fc)]
输出过程的平均功率为:
P
Y
=
.
.
.
=
1
2
R
X
(
0
)
P_Y=...=\frac{1}{2}R_X(0)
PY=...=21RX(0)
幅度调制系统的原理
调制谁?有用信号m(t),它是均值为0的基带信号。
用谁调制?载波信号
A
cos
(
2
π
f
c
t
)
A\cos(2{\pi}f_ct)
Acos(2πfct)
AM调制:
s
A
M
(
t
)
=
[
A
+
m
(
t
)
]
cos
(
2
π
f
c
t
)
s_{AM}(t)=[A+m(t)]\cos(2{\pi}f_ct)
sAM(t)=[A+m(t)]cos(2πfct)
解调:
-
包络检波器(利用电路检波)
-
相干解调
i n p u t : s A M ( t ) = [ A + m ( t ) ] cos ( 2 π f c t + ϕ 0 ) y ( t ) = [ A + m ( t ) ] cos ( 2 π f c t + ϕ 0 ) . 2 cos ( 2 π f c t + θ 0 ) = [ A + m ( t ) ] [ cos ( ϕ 0 − θ 0 ) + cos ( 4 π f c t + ϕ 0 + θ 0 ) ] s o ( t ) = [ A + m ( t ) ] [ cos ( ϕ 0 − θ 0 ) ] input:s_{AM}(t)=[A+m(t)]\cos(2{\pi}f_ct+\phi_0)\\ y(t)=[A+m(t)]\cos(2{\pi}f_ct+\phi_0).2\cos(2{\pi}f_ct+\theta_0)\\ =[A+m(t)][\cos({\phi_0}-{\theta_0})+\cos(4{\pi}f_ct+{\phi_0}+\theta_0)]\\ s_o(t)=[A+m(t)][\cos({\phi_0}-{\theta_0})] input:sAM(t)=[A+m(t)]cos(2πfct+ϕ0)y(t)=[A+m(t)]cos(2πfct+ϕ0).2cos(2πfct+θ0)=[A+m(t)][cos(ϕ0−θ0)+cos(4πfct+ϕ0+θ0)]so(t)=[A+m(t)][cos(ϕ0−θ0)]
频谱特性:
s
A
M
(
t
)
=
[
A
+
m
(
t
)
]
cos
(
2
π
f
c
t
)
=
A
cos
(
2
π
f
c
t
)
+
m
(
t
)
cos
(
2
π
f
c
t
)
S
A
M
(
f
)
=
A
2
[
δ
(
f
−
f
c
)
+
δ
(
f
+
f
c
)
]
+
1
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
P
A
M
=
A
2
2
+
m
2
(
t
)
ˉ
2
=
P
载
波
+
P
双
边
带
s_{AM}(t)=[A+m(t)]\cos(2{\pi}f_ct)=A\cos(2{\pi}f_ct)+m(t)\cos(2{\pi}f_ct)\\ S_{AM}(f)=\frac{A}{2}[\delta(f-f_c)+\delta(f+f_c)]+\frac{1}{2}[M(f-f_c)+M(f+f_c)]\\ P_{AM}=\frac{A^2}{2}+\frac{\bar{m^2(t)}}{2}=P_{载波}+P_{双边带}
sAM(t)=[A+m(t)]cos(2πfct)=Acos(2πfct)+m(t)cos(2πfct)SAM(f)=2A[δ(f−fc)+δ(f+fc)]+21[M(f−fc)+M(f+fc)]PAM=2A2+2m2(t)ˉ=P载波+P双边带
调制效率
η
A
M
=
P
双
边
带
/
P
A
M
\eta_{AM}=P_{双边带}/P_{AM}
ηAM=P双边带/PAM,效率并不高。主要原因是载波没有携带有用信息,但是它也会消耗功率。解决办法就是去掉载波的幅度A,也就是抑制载波的双边带调幅。
DSB调制: s D S B ( t ) = m ( t ) cos ( 2 π f c t ) s_{DSB}(t)=m(t)\cos(2{\pi}f_ct) sDSB(t)=m(t)cos(2πfct)
由于DSB的左右频谱是对称的,所以可以只传输一边的频带,即单边带调幅。
SSB调制: s D S B ( t ) = m ( t ) cos ( 2 π f c t ) s_{DSB}(t)=m(t)\cos(2{\pi}f_ct) sDSB(t)=m(t)cos(2πfct),只传输一半边的频谱。
除了AM调制,其它调制都需要用相干解调才能恢复原信号,因为对于其它调制方法来说,包络检波器只能得出是 ∣ m ( t ) ∣ |m(t)| ∣m(t)∣。
解调器输入噪声和输出噪声功率相等,即 P n o = P n i = N 0 B B P F P_{no}=P_{ni}=N_0B_{BPF} Pno=Pni=N0BBPF
调制方式 | 解调器输入信号 | 解调器输出信号 | 输入信号功率 | 输出信号功率 |
---|---|---|---|---|
AM | s A M ( t ) = [ A + m ( t ) ] cos ( 2 π f c t ) s_{AM}(t)=[A+m(t)]\cos(2{\pi}f_ct) sAM(t)=[A+m(t)]cos(2πfct) | m o ( t ) = m ( t ) m_o(t)=m(t) mo(t)=m(t) | S i = A 2 / 2 + m 2 ( t ) / 2 S_i=A^2/2+m^2(t)/2 Si=A2/2+m2(t)/2 | S o = m 2 ( t ) S_o=m^2(t) So=m2(t) |
DSB | s D S B ( t ) = m ( t ) cos ( 2 π f c t ) s_{DSB}(t)=m(t)\cos(2{\pi}f_ct) sDSB(t)=m(t)cos(2πfct) | m o ( t ) = m ( t ) m_o(t)=m(t) mo(t)=m(t) | S i = m 2 ( t ) / 2 S_i=m^2(t)/2 Si=m2(t)/2 | S o ( t ) = m 2 ( t ) S_o(t)=m^2(t) So(t)=m2(t) |
SSB | s S S B ( t ) = 1 2 [ m ( t ) cos ( 2 π f c t ) ± m ( t ) sin ( 2 π f c t ) ^ ] s_{SSB}(t)=\frac{1}{2}[m(t)\cos(2{\pi}f_ct){\pm}\hat{m(t)\sin(2{\pi}f_ct)}] sSSB(t)=21[m(t)cos(2πfct)±m(t)sin(2πfct)^] | m o ( t ) = m ( t ) / 2 m_o(t)=m(t)/2 mo(t)=m(t)/2 | S i = m 2 ( t ) / 4 S_i=m^2(t)/4 Si=m2(t)/4 | S o ( t ) = m 2 ( t ) / 4 S_o(t)=m^2(t)/4 So(t)=m2(t)/4 |
Q函数和互补误差函数
Q
(
x
)
=
P
(
Z
>
x
)
Q(x)=P(Z>x)
Q(x)=P(Z>x)表示
Z
>
x
Z>x
Z>x的概率**(前提是
Z
−
N
(
0
,
1
)
Z-N(0,1)
Z−N(0,1)分布)**,许多求解误码率的问题可归结为求概率
P
(
Z
>
a
)
P(Z>a)
P(Z>a)和
P
(
Z
<
−
a
)
P(Z<-a)
P(Z<−a)的问题,其中
a
>
0
a>0
a>0。对于均值为零的任意方差的高斯随机变量,令
Z
=
U
/
δ
Z=U/\delta
Z=U/δ,这时
Z
−
N
(
0
,
1
)
Z-N(0,1)
Z−N(0,1),有:
P
(
U
>
a
)
=
P
(
U
δ
>
a
δ
)
=
Q
(
a
δ
)
前
提
:
U
−
N
(
0
,
δ
2
)
P(U>a) = P(\frac{U}{\delta}>\frac{a}{\delta}) = Q(\frac{a}{\delta}) \\ 前提:\color{red}U-N(0,\delta^2)
P(U>a)=P(δU>δa)=Q(δa)前提:U−N(0,δ2)
对于任意方差的高斯随机变量
U
−
N
(
0
,
δ
2
)
U-N(0,\delta^2)
U−N(0,δ2),令
Z
=
U
/
(
2
δ
)
Z=U/(\sqrt{2}\delta)
Z=U/(2δ),这时
Z
−
N
(
0
,
1
/
2
)
Z-N(0,1/2)
Z−N(0,1/2),有:
P
(
U
>
a
)
=
P
(
U
2
δ
>
a
2
δ
)
=
1
2
e
r
f
c
(
a
2
δ
)
P(U>a) = P(\frac{U}{\sqrt{2}\delta}>\frac{a}{\sqrt{2}\delta}) = \frac{1}{2}erfc(\frac{a}{\sqrt{2}\delta})
P(U>a)=P(2δU>2δa)=21erfc(2δa)
两者的关系:
Q
(
x
)
=
1
2
e
r
f
c
(
x
2
)
=
1
−
Φ
(
x
)
Q(x) = \frac{1}{2}erfc(\frac{x}{\sqrt{2}}) = 1-\Phi(x)
Q(x)=21erfc(2x)=1−Φ(x)
Φ
(
x
)
\Phi(x)
Φ(x)是标准化正态分布函数,所以从上面看出,如果要使用Q函数或互补误差函数,则先要标准化自变量,即
U
−
N
(
0
,
1
)
U-N(0, 1)
U−N(0,1)。