Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Vychislitelny_praktikum

.pdf
Скачиваний:
26
Добавлен:
16.03.2015
Размер:
7 Мб
Скачать
☆

 

 

 

 

 

 

Продолжение приложения 8

Частотные характеристики ортогональных фильтров Якоби (α = 0, β = 1)

 

k = 0, 2, 4

 

 

 

k = 1, 3, 5

 

 

 

0

 

 

 

 

 

2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

0.8

 

 

 

 

 

1.2

 

 

 

 

 

Im(W11(k,ω,γ))

 

 

 

 

 

Im(W11(k,ω,γ))

 

 

 

 

 

1.2

 

 

 

 

 

2.8

 

 

 

 

 

1.6

 

 

 

 

 

4.4

 

 

 

 

 

2 1

0.2

0.6

1.4

2.2

3

6 4

2

0

2

4

6

 

 

Re(W11(k,ω,γ))

 

 

 

 

Re(W11(k,ω,γ))

 

 

10

 

 

 

 

 

10

 

 

 

 

 

6

 

 

 

 

 

6

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

Im(W11(k,ω,γ))

 

 

 

 

 

Im(W11(k,ω,γ))

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

6

 

 

 

 

 

6

 

 

 

 

 

10 10

6

2

2

6

10

10 10

6

2

2

6

10

 

 

Re(W11(k,ω,γ))

 

 

 

 

Re(W11(k,ω,γ))

 

 

10

 

 

 

 

 

20

 

 

 

 

 

4

 

 

 

 

 

12

 

 

 

 

 

2

 

 

 

 

 

4

 

 

 

 

 

Im(W11(k,ω,γ))

 

 

 

 

 

Im(W11(k,ω,γ))

 

 

 

 

 

8

 

 

 

 

 

4

 

 

 

 

 

14

 

 

 

 

 

12

 

 

 

 

 

20 10

5

0

5

10

15

20 20

12

4

4

12

20

 

 

Re(W11(k,ω,γ))

 

 

 

 

Re(W11(k,ω,γ))

 

 

203

Продолжение приложения 8

Частотные характеристики ортогональных фильтров Якоби (α =0, β = 2)

k = 0, 2, 4

 

 

 

 

k = 1, 3, 5

 

 

 

1

 

 

 

 

 

5

 

 

 

 

 

0.2

 

 

 

 

 

2

 

 

 

 

 

0.6

 

 

 

 

 

1

 

 

 

 

 

 

Im(W12(k,ω,γ))

 

 

 

 

 

Im(W12(k,ω,γ))

 

 

 

 

 

 

1.4

 

 

 

 

 

4

 

 

 

 

 

 

2.2

 

 

 

 

 

7

 

 

 

 

 

 

3 1

0

1

2

3

4

10 10

6

2

2

6

10

 

 

Re(W12(k,ω,γ))

 

 

 

 

 

Re(W12(k,ω,γ))

 

 

20

 

 

 

 

 

20

 

 

 

 

 

12

 

 

 

 

 

8

 

 

 

 

 

4

 

 

 

 

 

4

 

 

 

 

 

 

Im(W12(k,ω,γ))

 

 

 

 

 

Im(W12(k,ω,γ))

 

 

 

 

 

 

4

 

 

 

 

 

16

 

 

 

 

 

 

12

 

 

 

 

 

28

 

 

 

 

 

 

20 20

12

4

4

12

20

40

40

24

8

8

24

40

 

 

Re(W12(k,ω,γ))

 

 

 

 

 

Re(W12(k,ω,γ))

 

 

40

 

 

 

 

 

50

 

 

 

 

 

24

 

 

 

 

 

30

 

 

 

 

 

8

 

 

 

 

 

10

 

 

 

 

 

Im(W12(k,ω,γ))

 

 

 

 

 

Im(W12(k,ω,γ))

 

 

 

 

 

 

8

 

 

 

 

 

10

 

 

 

 

 

 

24

 

 

 

 

 

30

 

 

 

 

 

 

40 40

24

8

8

24

40

50

50

30

10

10

30

50

 

 

Re(W12(k,ω,γ))

 

 

 

 

 

Re(W12(k,ω,γ))

 

 

204

Продолжение приложения 8

Частотные характеристики ортогональных фильтров Якоби (α = −0,5, β = 0)

 

k = 0, 2, 4

 

 

 

 

 

k = 1, 3, 5

 

 

 

0

 

 

 

 

 

1

 

 

 

 

 

0.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

0.2

 

 

 

 

 

 

Im(W6(k,ω,γ))

 

 

 

 

 

Im(W6(k,ω,γ))

 

 

 

 

 

 

0.6

 

 

 

 

 

0.8

 

 

 

 

 

 

0.8

 

 

 

 

 

1.4

 

 

 

 

 

 

1 0

0.4

0.8

1.2

1.6

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W6(k,ω,γ))

 

 

 

 

 

Re(W6(k,ω,γ))

 

 

2

 

 

 

 

 

2

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

Im(W6(k,ω,γ))

 

 

 

 

 

Im(W6(k,ω,γ))

 

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

 

2 2

1.2

0.4

0.4

1.2

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W6(k,ω,γ))

 

 

 

 

 

Re(W6(k,ω,γ))

 

 

2

 

 

 

 

 

2

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

Im(W6(k,ω,γ))

 

 

 

 

 

Im(W6(k,ω,γ))

 

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

 

2 2

1.2

0.4

0.4

1.2

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W6(k,ω,γ))

 

 

 

 

 

Re(W6(k,ω,γ))

 

 

205

Продолжение приложения 8

Частотные характеристики ортогональных фильтров Якоби (α = 0,5, β = 0)

 

k = 0, 2, 4

 

 

 

 

 

k = 1, 3, 5

 

 

 

0

 

 

 

 

 

1

 

 

 

 

 

0.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

0.2

 

 

 

 

 

 

Im(W7(k,ω,γ))

 

 

 

 

 

Im(W7(k,ω,γ))

 

 

 

 

 

 

0.6

 

 

 

 

 

0.8

 

 

 

 

 

 

0.8

 

 

 

 

 

1.4

 

 

 

 

 

 

1 0

0.4

0.8

1.2

1.6

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W7(k,ω,γ))

 

 

 

 

 

Re(W7(k,ω,γ))

 

 

2

 

 

 

 

 

2

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

Im(W7(k,ω,γ))

 

 

 

 

 

Im(W7(k,ω,γ))

 

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

 

2 2

1.2

0.4

0.4

1.2

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W7(k,ω,γ))

 

 

 

 

 

Re(W7(k,ω,γ))

 

 

2

 

 

 

 

 

2

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

Im(W7(k,ω,γ))

 

 

 

 

 

Im(W7(k,ω,γ))

 

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

 

2 2

1.2

0.4

0.4

1.2

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W7(k,ω,γ))

 

 

 

 

 

Re(W7(k,ω,γ))

 

 

206

Продолжение приложения 8

Частотные характеристики ортогональных фильтров Дирихле

 

k = 0, 2, 4

 

 

 

 

 

k = 1, 3, 5

 

 

 

0

 

 

 

 

 

2

 

 

 

 

 

0.2

 

 

 

 

 

1.4

 

 

 

 

 

0.4

 

 

 

 

 

0.8

 

 

 

 

 

Im(W5(k,ω,γ))

 

 

 

 

 

Im(W5(k,ω,γ))

 

 

 

 

 

 

0.6

 

 

 

 

 

0.2

 

 

 

 

 

0.8

 

 

 

 

 

0.4

 

 

 

 

 

 

1 0

0.4

0.8

1.2

1.6

2

1

2

1.2

0.4

0.4

1.2

2

 

 

Re(W5(k,ω,γ))

 

 

 

 

 

Re(W5(k,ω,γ))

 

 

2

 

 

 

 

 

2

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

Im(W5(k,ω,γ))

 

 

 

 

 

Im(W5(k,ω,γ))

 

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

 

2 2

1.2

0.4

0.4

1.2

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W5(k,ω,γ))

 

 

 

 

 

Re(W5(k,ω,γ))

 

 

2

 

 

 

 

 

2

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

Im(W5(k,ω,γ))

 

 

 

 

 

Im(W5(k,ω,γ))

 

 

 

 

 

 

0.4

 

 

 

 

 

0.4

 

 

 

 

 

 

1.2

 

 

 

 

 

1.2

 

 

 

 

 

 

2 2

1.2

0.4

0.4

1.2

2

2

2

1.2

0.4

0.4

1.2

2

 

 

Re(W5(k,ω,γ))

 

 

 

 

 

Re(W5(k,ω,γ))

 

 

207

Приложение 9

ВЫЧИСЛИТЕЛЬНЫЙ ПРАКТИКУМ 4. «ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ ОРТОГОНАЛЬНЫХ ФУНКЦИЙ»

Цель работы: исследование частотных характеристик ортогональных функций.

1. Построить частотные характеристики, соответствующие выражениям (4.6) - (4.9); m = 0 – 4,

γ = const.

m := 0.. 4

γ := 0.87

ω := 0,0 + 0.01.. 50

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

φ(k,ω) := atan

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ (k + 1)

 

 

 

 

 

 

 

 

 

 

ReW(m,ω) :=

 

 

1

 

 

cos (φ(m,ω)) cos φ(m,ω) +

 

0 if m

 

 

 

0

 

 

 

 

 

 

 

 

 

 

γ (m +

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

m−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ∑

φ(s ,ω)

otherwise

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s = 0

 

 

 

 

 

ImW(m,ω) := −

1

 

 

 

cos (φ(m,ω)) sin φ(m,ω) +

 

0 if m

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ (m + 1)

 

 

 

m−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ∑

 

φ(s ,ω)

otherwise

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s = 0

 

 

 

 

 

MoW (m,ω) :=

cos (φ(m,ω))

 

 

 

 

 

 

 

 

 

 

 

 

 

γ (m + 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m−1

 

 

 

 

 

 

 

 

 

φ(m,ω) + 2 ∑

φ(s ,ω)

 

 

m−1

 

 

 

 

 

 

 

 

φ (m,ω) + 2 ∑

 

+ trunc

 

 

s = 0

 

 

 

ΦW(m,ω) := −

φ (s ,ω)

 

 

π

 

+ 0.5

π

 

s = 0

 

 

 

 

 

 

 

 

208

 

 

1.5

 

 

 

 

 

ReW(0,ω) 1.1

 

 

 

 

 

ReW(1,ω)

 

 

 

 

 

 

ReW(2,ω) 0.7

 

 

 

 

 

ReW(3,ω) 0.3

 

 

 

 

 

ReW(4,ω)

 

 

 

 

 

 

 

− 0.1

 

 

 

 

 

 

− 0.50

4

8

12

16

20

 

 

 

 

 

ω

 

 

 

 

0.4

 

 

 

 

 

ImW(0

,ω)

0.2

 

 

 

 

 

ImW(1

,ω)

 

 

 

 

 

 

ImW(2

,ω)

0

 

 

 

 

 

 

 

 

 

 

 

ImW(3

,ω)− 0.2

 

 

 

 

 

ImW(4,ω)

− 0.4

− 0.6

 

 

 

 

 

 

 

0

4

8

12

16

20

 

 

 

 

ω

 

 

 

1.5

MoW(0,ω)1.2

MoW(1,ω)

MoW(2,ω)0.9

MoW(3,ω)0.6

MoW(4,ω)

0.3

0

 

 

 

 

 

 

 

0

3

6

9

12

15

 

 

 

 

ω

 

 

 

209

 

2

π

 

 

ΦW(0,ω)

1.2

2

 

ΦW(1,ω)

ΦW(2,ω) 0.4

ΦW(3,ω)− 0.4

ΦW(4,ω)

− 1.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−

 

 

 

− 2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

10

20

30

40

 

 

 

 

50

ω

2. Построить зависимость квадрата модуля частотной характеристики ортогонального фильтра от порядка m; m = 0 - 4, γ = const.

MoWf (m,ω) :=

4 (m + 1)2 γ2

 

 

 

 

 

(m + 1)2 γ2 + ω2

 

 

 

 

 

4

 

 

 

 

 

 

MoWf (0,ω)3.2

 

 

 

 

 

 

MoWf (1

,ω)

 

 

 

 

 

 

MoWf (2

,ω)2.4

 

 

 

 

 

 

MoWf (3

,ω)1.6

 

 

 

 

 

 

MoWf (4

,ω)

 

 

 

 

 

 

 

0.8

 

 

 

 

 

 

 

0

0

5

10

15

20

25

 

 

 

 

 

ω

 

 

3. Построить зависимость квадрата модуля частотной характеристики ортогонального фильтра от параметра γ; m = const.

m := 3

 

MoWf (γ,ω1) :=

4 (m + 1)2 γ2

(m + 1)2 γ2 + ω12

 

210

CreateMesh (MoWf )

4. Построить зависимость полосы пропускания ортогонального фильтра m - ого порядка от параметра γ; m = 0 - 4.

ω

c

(m,γ) := π γ (m + 1)

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

ωc(0,γ1 )160

 

 

 

 

 

 

ωc(1,γ1 )

 

 

 

 

 

 

ωc(2

120

 

 

 

 

 

 

,γ1 )

 

 

 

 

 

 

ωc(3

,γ1 ) 80

 

 

 

 

 

 

ωc(4

,γ1 )

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

 

0

0

5

10

15

20

25

 

 

 

 

 

 

 

γ1

 

 

211

Приложение 10

АНАЛИТИЧЕСКИЕ ВЫРАЖЕНИЯ ДЛЯ ВЫЧИСЛЕНИЯ КОЭФФИЦИЕНТОВ РАЗЛОЖЕНИЯ

Модель ρx (τ )= e−λ τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица П 10.1

ψk (τ ,γ / α)

 

 

 

 

 

 

 

 

 

 

βk

 

 

 

 

1

 

 

 

 

P(−1 2,0) (τ,γ)

γ (4k +1)∑k

Cks Cks+−s1−/

12/ 2

(−1)s

 

 

 

 

 

 

 

 

 

(λ +γ ((4s +1)/ 2))

k

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ (4k + 3)∑k

 

 

 

 

(−1)s

 

 

 

 

 

1

 

 

 

 

P(1 2,0 )(τ,γ)

Cks Cks++s1+/

12/ 2

 

 

 

 

 

 

 

(λ +γ ((4s + 3)/ 2))

k

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2γ (2k +1)∑k

Cks Cks+s (−1)s

 

 

 

 

1

 

 

 

 

 

P( 0,0 ) (τ,γ)

 

 

 

 

 

 

 

(λ +γ (2s +1))

k

 

 

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2γ (k +1)∑k

Cks Cks++s1+1 (−1)s

 

 

 

 

1

 

 

 

 

 

P( 1,0 ) (τ,γ)

 

 

 

 

 

 

 

 

 

(λ +γ (s +1))

k

 

 

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2γ (2k + 3)∑k

Cks Cks++s2+2 (−1)s

 

 

 

 

1

 

 

 

 

P( 2,0 ) (τ,γ)

 

 

 

 

 

 

 

 

 

(λ +γ (2s + 3))

k

 

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ ∑Cks

 

(−γ )

s

s+1

Lk (τ ,γ )

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s=0

 

 

(λ +γ / 2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

k −s (−γ )s (s +1)

(1)

(τ ,γ )

 

 

2

 

 

Lk

 

γ

 

∑Ck +1

 

 

 

 

 

 

 

 

 

 

 

 

(λ +γ / 2)

s+2

 

 

 

 

 

 

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ 3

k

 

 

k −s (−γ )s (s +1)(s + 2)

(2 )

(τ ,γ )

 

 

Lk

 

∑Ck +2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

(λ +γ / 2)

s+3

 

 

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

212

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]