Vychislitelny_praktikum
.pdf2. Рассчитать значения ортогональных функций 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. Определить максимальный интервал корреляции на основании графика.
τ1 := 10
Given
P(k,τ1,γ) δ
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