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

Vychislitelny_praktikum

.pdf
Скачиваний:
12
Добавлен:
16.03.2015
Размер:
7.27 Mб
Скачать

2. Рассчитать значения ортогональных функций k-ого порядка в "нуле".

k := 0.. 5

P(k ,0)

1

-1

1

-1

1

-1

3. Определение нормы ортогональных функций:

3.1. Определить значение нормы ортогональных функций из выражения (2.3) (см. содержание лабораторной работы). Результат представить в виде матрицы значений с разрядностью (k,m), Привести графическую интерпретацию.

k := 6

i := 0.. k 1

m := k

μ(τ ,γ ) := 1

j := 0.. m 1

NormP(i,j ) := P(i,τ ,γ ) P(j ,τ ,γ ) μ(τ ,γ ) dτ

0

NormMP := matrix(k ,m,NormP)

 

0.5

0

0

0

0

0

 

 

 

0

0.25

0

0

0

0

 

 

 

0

0

0.167

0

0

0

 

NormMP =

 

 

 

0

0

0

0.125

0

0

 

 

 

 

 

 

 

0

0

0

0

0.1

0

 

 

0

0

0

0

0

0.083

173

NormMP

3.2. Определить значения нормы ортогональных функцийk-ого порядка, используя выражение, приведенное в таблице 2.1 (см. содержание лабораторной работы). Результат представить в виде вектора значений.

1 NormPTi := 2γ (i + 1)

NormPTi =

0.5

0.25

0.167

0.125

0.1

0.083

4. Рассчитать длительность, используя выражения (2.5) и (2.6) (см. содержание лабораторной работы). Построить графическую зависимость длительности ортогональной функции k-ого порядка от параметра масштаба.

k := 0.. 5

174

 

 

tau4(k )

P(k ,τ ,γ )2 dτ

tau2(k ) :=

P(k ,τ ,γ ) dτ

 

:=

 

0

 

 

0

 

tau2(k ) =

 

 

 

tau4(k ) =

 

 

1

 

 

 

0.5

 

 

0.5

 

 

 

0.25

 

 

0.333

 

 

 

0.167

 

 

0.25

 

 

 

0.125

 

 

0.2

 

 

 

0.1

 

 

0.167

 

 

 

0.083

 

 

γ1 := 0.1,0.1 + 0.1.. 2

 

 

 

 

i := 0.. 20

 

 

 

 

 

 

tau2V0i := tau2(0,i 0.1 + 0.1)

tau2V2i := tau2(2,i 0.1 + 0.1)

tau2V4i := tau2(4,i 0.1 + 0.1)

 

 

 

 

tau2V1i := tau2(1,i 0.1 + 0.1)

tau2V3i := tau2(3,i 0.1 + 0.1)

tau2V5i := tau2(5,i 0.1 + 0.1)

tau2M 0 := tau2V0

tau2M 2 := tau2V2

 

tau2M 4 := tau2V4

tau2M 1 := tau2V1

tau2M 3

:= tau2V3

 

tau2M 5 := tau2V5

(tau2M

0 )

 

 

 

 

 

 

i

 

 

 

 

 

(tau2M

1 )8

 

 

 

 

 

i

 

 

 

 

 

(tau2M

2 )

 

 

 

 

 

 

i6

 

 

 

 

(tau2M

3 )

 

 

 

 

 

 

i

 

 

 

 

 

(tau2M

4 )4

 

 

 

 

 

i

 

 

 

 

 

(tau2M

5 )

 

 

 

 

 

 

i2

 

 

 

 

 

 

0.5

 

1

1.5

2

 

 

 

 

i 0.1+0.1

 

 

175

 

 

0

1

2

3

4

5

 

0

10

5

3.333

2.5

2

1.667

 

1

5

2.5

1.667

1.25

1

0.833

 

2

3.333

1.667

1.111

0.833

0.667

0.556

 

3

2.5

1.25

0.833

0.625

0.5

0.417

 

4

2

1

0.667

0.5

0.4

0.333

 

5

1.667

0.833

0.556

0.417

0.333

0.278

 

6

1.429

0.714

0.476

0.357

0.286

0.238

tau2M =

7

1.25

0.625

0.417

0.313

0.25

0.208

 

8

1.111

0.556

0.37

0.278

0.222

0.185

 

9

1

0.5

0.333

0.25

0.2

0.167

 

10

0.909

0.455

0.303

0.227

0.182

0.152

 

11

0.833

0.417

0.278

0.208

0.167

0.139

 

12

0.769

0.385

0.256

0.192

0.154

0.128

 

13

0.714

0.357

0.238

0.179

0.143

0.119

 

14

0.667

0.333

0.222

0.167

0.133

0.111

 

15

0.625

0.312

0.208

0.156

0.125

...

 

 

 

 

 

 

 

 

tau2M

176

tau4V0i

:= tau4(0,i 0.1 + 0.1)

tau4V2i := tau4(2,i 0.1 + 0.1)

tau4V4i

:= tau4(4,i 0.1 + 0.1)

tau4V1i

:= tau4(1,i 0.1 + 0.1)

tau4V3i := tau4(3,i 0.1 + 0.1)

tau4V5i

:= tau4(5,i 0.1 + 0.1)

tau4M 0 := tau4V0

tau4M 2 := tau4V2

 

tau4M 4 := tau4V4

tau4M 1 := tau4V1

tau4M 3 := tau4V3

 

tau4M 5 := tau4V5

(tau4M

0 )

 

 

 

 

 

 

i

 

 

 

 

(tau4M

1 )4

 

 

 

 

 

 

i

 

 

 

 

(tau4M

2 )

 

 

 

 

 

 

i3

 

 

 

 

(tau4M

3 )

 

 

 

 

 

 

i

 

 

 

 

(tau4M

4 )2

 

 

 

 

 

 

i

 

 

 

 

(tau4M

5 )

 

 

 

 

 

 

i1

 

 

 

 

 

 

0.5

1

1.5

2

 

 

 

 

i 0.1+0.1

 

 

 

 

 

0

1

2

3

4

5

 

0

5

2.5

1.667

1.25

1

0.833

 

1

2.5

1.25

0.833

0.625

0.5

0.417

 

2

1.667

0.833

0.556

0.417

0.333

0.278

 

3

1.25

0.625

0.417

0.313

0.25

0.208

 

4

1

0.5

0.333

0.25

0.2

0.167

 

5

0.833

0.417

0.278

0.208

0.167

0.139

 

6

0.714

0.357

0.238

0.179

0.143

0.119

tau4M =

7

0.625

0.313

0.208

0.156

0.125

0.104

 

8

0.556

0.278

0.185

0.139

0.111

0.093

 

9

0.5

0.25

0.167

0.125

0.1

0.083

 

10

0.455

0.227

0.152

0.114

0.091

0.076

 

11

0.417

0.208

0.139

0.104

0.083

0.069

 

12

0.385

0.192

0.128

0.096

0.077

0.064

 

13

0.357

0.179

0.119

0.089

0.071

0.06

 

14

0.333

0.167

0.111

0.083

0.067

0.056

 

15

0.312

0.156

0.104

0.078

0.063

...

 

 

 

 

 

 

 

 

177

tau4M

178

Приложение 6

ВЫЧИСЛИТЕЛЬНЫЙ ПРАКТИКУМ 3. «ОПРЕДЕЛЕНИЕ МАКСИМАЛЬНОЙ ДЛИТЕЛЬНОСТИ И ИНТЕРВАЛА ДИСКРЕТИЗАЦИИ

ОРТОГОНАЛЬНЫХ ФУНКЦИЙ»

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

1. Задать вид ортогональных функций.

k

k!

 

 

(k + s + 1)!

 

s

 

(s+1) γ τ

P(k,τ,γ) :=

 

(1)

e

 

 

 

 

 

 

(k s)! s!

k! (s + 1)!

 

 

s = 0

 

 

 

 

 

 

 

 

 

2. Задать порядок ортогональных функций.

k:=

3.Задать значение параметра масштаба.

γ:= 1

4.Задать погрешность приближения.

δn :=

 

δ :=

0.02

if

δn

 

 

1

 

 

 

 

 

 

 

 

 

0.05

if

δn

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

0.1

if

δn

 

 

3

 

 

 

 

 

 

5. Построить график ортогональной функции.

0.5

 

 

 

 

0

 

 

 

δ

 

 

 

 

P(k,τ,γ)

 

 

 

 

 

P(k,τ,γ)0.5

 

 

 

 

 

P(k,τ,γ) −δ

 

 

 

 

 

1

 

 

 

 

 

1.5

0

2.5

5

7.5

10

 

 

 

τ

 

 

179

6. Определить максимальный интервал корреляции на основании графика.

τ1 := 10

Given

P(k1) δ

Find(τ1) = 5.641

7. Получить на основании общей формулы частные формулы ортогональной функции, ее первой, второй и третьей производной.

P(k,τ,γ) 6 exp(−τ) 105 exp(2 τ) + 560 exp(3 τ) 1260exp(4 τ) + 1260exp(5 τ) 462exp(6 τ)

d P(k,τ,γ) → −6 exp(−τ) + 210exp(2 τ) 1680exp(3 τ) + 5040exp(4 τ) 6300exp(5 τ) + 2772exp(6 τ)

dτ

d2 P(k,τ,γ) 6 exp(−τ) 420exp(2 τ) + 5040exp(3 τ) 20160exp(4 τ) + 31500exp(5 τ) 16632exp(6 τ) dτ2

d3 P(k,τ,γ) → −6 exp(−τ) + 840exp(2 τ) 15120exp(3 τ) + 80640exp(4 τ) 157500exp(5 τ) + 99792exp(6 τ) dτ3

8. Вычислить максимум второй производной.

500

375

d2 P(k,τ,γ) 250 dτ2

125

0

 

 

 

 

 

 

0

1

2

3

4

 

 

 

τ

 

 

 

180

2

fder (τ) := d P(k,τ,γ) dτ2

τ := 25

Maximize(fder ) = 0

fder (Maximize(fder )) = 666

9. Вычислить интервал дискретизации для ортогональной функции.

t :=

 

 

8 δ

 

fder (Maximize(fder ))

 

 

 

 

 

 

 

t= 0.015

10.Определить максимальную длительность ортогональной функции (см. п.6).

τmax := 5.64

11. Вычислить количество интервалов дискретизации.

τmax

 

n := trunc

 

 

+ 0.5

 

 

 

t

 

n= 364

12.Занести полученные значения в таблицу (Microsoft Word, Microsoft Excel).

13 - 15. Выполнить аналогично 1-12.

181

16. Построить зависимость максимального интервала корреляции от порядка при заданном параметре масштаба.

При параметре масштаба γ = 1, k = 0..5.

j := 0.. 5

3.9124.59

4.983

Vector_τmax := 5.259

5.47

5.641

6

 

 

 

 

 

 

5.5

 

 

 

 

 

 

5

 

 

 

 

 

 

Vector_τmax

 

 

 

 

 

 

j

 

 

 

 

 

 

4.5

 

 

 

 

 

 

4

 

 

 

 

 

 

3.5

0

1

2

3

4

5

 

 

 

 

j

 

 

17. Построить зависимость максимального интервала корреляции от

параметра масштаба при заданном порядке ортогональной функции.

При заданном порядке k = 0, γ = 0.5, 0.5+0.1..1.

ε := 0.5,0.5 + 0.1.. 1

7.8246.52

5.589 Vector_τmax := 4.894.347

3.912

182

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