当前位置:首页 >> 其它课程 >> 数字信号处理 DSP 英文版课件5.1

数字信号处理 DSP 英文版课件5.1


In This Lecture:
5.4 Frequency Spectrum ·X ( z ) ~ X (? ) (z-domain view of frequency spectrum) 5.5 Inverse z-Transform · Partial fraction expansion method

5.4 Frequency Spec

trum
1. Relation of z-transform with spectrum

For a signal x(n) ,
X (? ) ?
n ? ??

x(n)e ? j?n ?
?

?

(DTFT)

(5.4.1)

X ( z) ?

n ? ??

x ( n) z ? n ?

(z-transform)

z ? re j?

?
?n ? j?n

X ( z) |z ?e j? ? ? x(n) z |z ?e j? ? ? x(n)e
n n

? X (?)

------- the evaluation of the z-transform on the unit circle

Illustration in Matlab:

z-plane

X ( z) ? ? x(n)r e
n

? n ? j?n

? DTFT[ x(n)r ]
z 变换扩大了能分析的 信号的范围

?n

X (?) ? X ( z) |z ?e j?
(Definition:

X (? ) ?

n ? ??

x(n)e ? j?n , ?

?

X ( z) ?

n ? ??

x ( n) z ? n ) ?

?

Frequency response of a linear system h(n) :

H (? ) ? DTFT[h(n)] ?

n???

h(n)e ? j?n ?

?

(5.4.3)

H (?) ? H ( z) | z?e j?

( H (z ) : Transfer function)

2. Meaning of the spectrum (a review)

2?f Digital frequency: ? ? fs
f : physical frequency (Hz) Nyquist interval:

(5.4.4)

?

fs

2

? f ?

fs

2

? ?? ? ? ? ? (5.4.5)

Inverse DTFT:

1 x ( n) ? 2?

??
?

?

X (? )e j?n d?

(5.4.6)

It expresses x(n) as a linear combination of discrete-time sinusoids e
j?n

of different frequencies.

The relative amplitudes and phases of these sinusoidal components are given by the DTFT X (? ) .
Sinusoidal frequency component

3. The spectrum of complex sinusoid

Double-sided complex sinusoid (双边复正弦):

x(n) ? e

j?0 n

, ??? n ? ?

?
?

Introducing ? (?) in frequency domain

X (?) ? 2?? (? ? ?0 ) +(Nyquist replica)
? 2?

m???

? ? (? ? ?

0

? 2?m)

Assuming

? ? ? ?0 ? ? ,

X (?) ? 2?? (? ? ?0 ) (?? ? ? ? ? )
Verification:

1 2?

2?? (? ? ?0 )e j?n d? ? e j?0n ? x(n) ???

?

or

X ( f ) ? ? ( f ? f0 )

( ? 0.5 ?

f ? 0.5 )

For a linear combination of two sinusoids:

x(n) ? A1e j?1n ? A2e j?2n
X (?) ? 2?A1? (? ? ?1 ) ? 2?A2? (? ? ?2 )

(?? ? ? ? ? )

(Assume both ?1 and interval)

?2 lie in the Nyquist

For real-valued cosine and sine signals:

cos(?0 n) ? ?? (? ? ?0 ) ? ?? (? ? ?0 ) sin(?0 n) ? ? j?? (? ? ?0 ) ? j?? (? ? ?0 )
j?0 n

( cos(?0 n) ?

e

?e 2

? j?0 n

, sin(?0 n) ?

e

j?0 n

?e 2j

? j?0 n



Causal, truncated complex sinusoid(单边复正弦):

x(n) ? e j?0n u (n)
X ( z) ? ? e
n ?0 ? j?0 n

(p.202, bottom)

z

?n

?

1 1? e
j?0

z

?1

(| z |? 1)

(The spectrum doesn’t exit)

The formal replacement of z by e

j?

yields:

X (? ) ?

1 1? e
j? 0

e

? j?

(

?0 : dominant frequency)

4. Geometric interpretation of frequency spectrum
频谱的几何确定法 ――通过X(z)零极点位置估计幅度、相位谱

A system with a zero z ? z1 and a pole

z ? p1 ,

1 ? z1z z ? z1 X ( z) ? ? ?1 z ? p1 1 ? p1z

?1

e j? ? z1 ? X (? ) ? j? e ? p1

?

? | e j? ? z1 | ?| X (? ) |? j? | e ? p1 | ? ?arg( X (? )) ? arg(e j? ? z ) ? arg(e j? ? p ) 1 1 ?

peak

dip

As

e

j?

moves around the unit circle, the denominator(分母)
j?

distance and numerator(分子) distance will vary.

e moves to the point with phase ? ? ?1 ? arg( p1 ) , it will cause a peak in | X (? ) | .
When When

angle

e

j?

moves

to

the

point

with .

phase

angle

? ? ?1 ? arg( z1 ) , it will cause a dip in | X (? ) |

Approximating the “spectrum” of

x ( n) ? e

j?0 n

u ( n) :

Suppose Pole:

x(n) ? r0 e
n

j?0n

u ( n)

( r0 ? 1 )

z= r0 e

j?0

Spectrum:

1 e X (? ) ? ? j? j?0 ? j? j?0 1 ? r0e e e ? r0e
| X (? ) |? | e j? 1 j?0 ? r0e |

j?

Magnitude spectrum:

When If

? ? ?0 ,

| e j? ? r0 e j?0 | decreases , | X (? ) | increases.

r0 ->1, then | X (?0 ) |? ?

提示: 几何确定法可用于定性估计已知系统的频响,也可根据频 响要求来设计系统。一般适用于零极点模式较简单的情况。

5. Properties of DTFT
◆ For real-valued signals x(n),

X (? ) ? X (?? )
*

(5.4.8)

?| X (? ) |?| X (??) | ?? ?arg X (?) ? ? arg X (??)

(5.4.9)

The magnitude response is even (偶) in response odd(奇)

? and the phase

For example, x =[1 1 1 1 1 1 1 1 1 1]

◆ y (n) ?

x ( n ) * h ( n ) ――> Y (? ) ? X (? ) H (? )
(5.4.10)

◆ Parseval’s equation

1 ? | x(n) | ? 2? n ? ??
2

?

??? | X (? ) | d?
2

?

(5.4.7)

5.5 Inverse z-Transform
提示:
求X(z)逆变换的数学依据是,函数X(z)在以它的极点所 在圆为边界的圆环域中有不同的唯一的罗朗级数展开式。 级数中z各幂次项的系数是反变换序列样值。

给定X(z)及某个ROC就有唯一的反变换x(n),如未指 定ROC,就可能有好几个x(n)。

求z反变换的公式是复平面上圆周积分,一般可借助于 留数方法求。但对较规则的X(z),可凑成级数展开式。

Partial fraction expansion method(部分分式展开法)

N ( z) N ( z) X ( z) ? ? D( z ) (1 ? p1 z ?1 )(1 ? p2 z ?1 )?(1 ? p M z ?1 )
A1 A2 AM ? ? ?? ?1 ?1 1 ? p1 z 1 ? p2 z 1 ? pM z ?1
(5.5.1)

pi (i ? 1 ~ M ) :zeros of the denominator polynomial D(z) ,
or the poles of X (z )

According to the degrees of N(z) and D(z), the representation differs slightly:
(1)The degree of N(z) is strictly less than M,

A1 A2 AM X ( z) ? ? ?? ?1 ?1 ?1 1 ? p1z 1 ? p2 z 1 ? pM z
Ai ? [(1 ? pi z ?1 ) X ( z )] | z ? pi (5.5.2) the coefficients

(2)The degree of N(z) is exactly equal to M,

A1 A2 AM X ( z ) ? A0 ? ? ?? ?1 ?1 1 ? p1 z 1 ? p2 z 1 ? pM z
(5.5.3)

A0 ? X ( z ) | z ?0
Ai (i ? 0)

(5.5.4)

are computed in the same way as (5.5.2)

(3)The degree of N(z) is strictly greater than M,

R( z ) X ( z ) ? Q( z ) ? D( z )
( , ) Q(z) :quotient(商式) R(z) :remainder(余式)

R( z ) Technique 1: ordinary partial fraction expansion for D( z)
Technique 2: Remove / Restore

1 X ( z) ? N ( z) ? ? N ( z )W ( z ) D( z )

Example 5.5.6: Determine the causal inverse z-transform

6? z X ( z) ? ?2 1 ? 0.25 z

?5

Solution: Technique 1: divide the denominator into the numerator
用长除法求出商和余(逐次消去幂次高的项,直至余式阶数小于 2) 。

6 ? 16 z ?1 X ( z ) ? ?16 z ?1 ? 4 z ?3 ? 1 ? 0.25 z ?2 19 13 ?1 ?3 ? ?16 z ? 4 z ? ? ?1 1 ? 0.5 z 1 ? 0.5 z ?1

? x(n) ? ?16? (n ?1) ? 4? (n ? 3)
? 19 ? 0.5 u(n) ? 13(?0.5) u(n)
n n

Technique 2:Remove / Restore

6 ? z ?5 X ( z) ? 1 ? 0.25 z ?2

1 0.5 0.5 W ( z) ? ? ? ?2 ?1 1 ? 0.25 z 1 ? 0.5 z 1 ? 0.5 z ?1

w(n) ? 0.5 ? 0.5 u (n) ? 0.5(?0.5) u (n)
n n

X ( z ) ? (6 ? z ?5 )W ( z ) ? 6W ( z ) ? z ?5W ( z )

? x(n) ? 6w(n) ? w(n ? 5)
? 3(0.5) n u(n) ? 3(?0.5) n u(n) ? 0.5(0.5) n?5 u(n ? 5) ? 0.5(?0.5) n?5 u(n ? 5)

For X(z) with real-valued coefficients, the poles are real or complex conjugate pairs(复共轭对).
* A1 A1 A2 X ( z) ? ? ? ?? ?1 * ?1 ?1 1 ? p1z 1 ? p1 z 1 ? p2 z

(p.211,bottom)

x(n) ? A p u (n) ? A p u (n) ? A2 p u (n) ? ...
n 1 1 * *n 1 1 n 2 n ? 2 Re[ A1 p1n ]u (n) ? A2 p2 u (n) ? ...

Problems: 5-8 (b) (c) (e)


更多相关文档:

科技英语 5数字信号处理器原文与翻译

搜 试试 7 帮助 全部 DOC PPT TXT PDF XLS ...科技英语 5数字信号处理器原文与翻译_英语学习_外语...In DSP, the term "core" refers to the section...

数字信号处理实验5

数字信号处理实验课件5 17页 免费 数字信号处理实验三(5) 1页 1财富值 数字信号处理 实验五 5页 2财富值 DSP实验5 4页 1财富值 数字信号处理(MTLAB版)实验...

数字信号处理英文翻译

搜 试试 7 帮助 全部 DOC PPT TXT PDF XLS ...5、 In troduction of DSP Today, DSP is w ...1数字滤波器介绍:模拟和数字滤波器 在信号处理、...

数字信号处理实验报告5

搜 试试 7 帮助 全部 DOC PPT TXT PDF XLS 百度文库 专业资料 自然科学 ...数字信号处理实验报告第... 4页 1下载券 DSP实验报告 (5) 3页 1下载券 哈...

数字信号处理 实验 (3)

关键词:数字信号处理dsp实验报告 同系列文档 数字信号处理 实验 (1) 数字信号处理...数字信号处理实验课件3 7页 免费 数字信号处理实验三 5页 1财富值喜欢...

数字信号处理 实验 (4)

数字信号处理实验课件4 20页 免费 dsp实验4 3页 1财富值喜欢此文档的还喜欢 ...3页 5财富值如要投诉违规内容,请到百度文库投诉中心;如要提出功能问题或意见建议...

数字信号处理(MTLAB版)实验5

数字信号处理(MTLAB版)实验5_工学_高等教育_教育专区。实验五今日...数字信号处理实验课件5 17页 免费 数字信号处理 实验五 5页 1下载券 DSP实验...

专业英语翻译之数字信号处理

搜 试试 7 帮助 全部 DOC PPT TXT PDF XLS ...5. Frequency domain Main article: Frequency domain...1.数字信号处理 数字信号处理(DSP)是关注的信号通过...

数字信号处理 复习 DSP_Review1_5

dsp数字信号处理复习课件... 153页 1下载券 数字信号处理器(DSP)复习... 9...(3.6.5) (3.6.4) and (3.6.5) is a called DTFT pair. Note that X ...

DSP 5

搜 试试 帮助 全部 DOC PPT TXT PDF XLS 百度文库 专业资料 IT/计算机 ...数字信号处理5 40页 1财富值 dsp 5-3 24页 2财富值喜欢此文档的还喜欢 DSP...
更多相关标签:
dsp数字信号处理器 | dsp数字信号处理 | 数字信号处理英文版 | 数字信号处理课件 | 数字信号处理课件下载 | 现代数字信号处理课件 | dsp课件 | dsp信号发生器 |
网站地图

文档资料共享网 nexoncn.com copyright ©right 2010-2020。
文档资料共享网内容来自网络,如有侵犯请联系客服。email:zhit325@126.com