Vychislitelny_praktikum

m

ρ(τ ,m,χ) :=

k = 0

Δτ := 0.081649

Nx := 37

τ := 0,0 + Δτ

⌠

∞

− λ τ

β5Mχ1m_1k

:= 2 χ1c(m_1) (k + 1)

P7(k ,τ ,χ1c(m_1)) e

cos (ω τ ) dτ

⌡0

k

m

s

2 (−1)

(k +

∑

1) 1 −

(−1)

β5Mχ1m_1s

b5Mχ1m_1k := β5Mχ1m_1k +

s = 0

m_1

ρ5χ1(τ ) := ∑

k = 0

⌠∞

− λ τ

β5Mχ2m_1k := 2 χ2c(m_1) (k + 1)

P7(k ,τ ,χ2c(m_1)) e

cos (ω τ ) dτ

⌡0

k

m

s

2 (−1)

− ∑

(k + 1) 1

(−1)

β5Mχ2m_1s

s = 0

b5Mχ2m_1k := β5Mχ2m_1k +

(m + 1) (m + 2)

m_1

ρ5χ2(τ ) :=

∑

k = 0

⌠∞

β5Mχm_1k

:= 2 χminMm_1−2 (k + 1)

⌡0

k

m

s

2 (−1)

(k

− ∑

+ 1) 1

(−1)

β5Mχm_1s

b5Mχm_1k := β5Mχm_1k +

s = 0

(m + 1) (m + 2)

m_1

ρ5χ(τ ) := ∑

k = 0

ρ5χ1 (τ)

ρ5χ2 (τ) 0.5

ρ5χ (τ)

0

0

1

2

ρt(τ)

− 0.5

− 1

τ

δ3V8i := δ1(i 0.1 + 0.01,10,5)

δ3M 8 := δ3V8

(

8 )

= 0.143

min δ3M

(

(

δ3M

8 )

,δ3M

8 )

0.1 + 0.01 = ( 1.11)

match

min

χminM

:= match

(

(

8 )

,δ3M

8 )

0 0.1 + 0.01

8

min δ3M

⌠

∞

− λ τ

β5Mχ1m_2k

:= 2 χ1c(m_2) (k + 1)

P7(k ,τ ,χ1c(m_2)) e

cos (ω τ ) dτ

⌡0

k

m

s

2 (−1)

∑

(k + 1) 1 −

(−1)

β5Mχ1m_2s

b5Mχ1m_2k := β5Mχ1m_2k +

s = 0

m_2

ρ10χ1(τ ) :=

∑

k = 0

⌠∞

− λ τ

β5Mχ2m_2k := 2 χ2c(m_2) (k + 1)

P7(k ,τ ,χ2c(m_2)) e

cos (ω τ ) dτ

⌡0

k

m

s

2 (−1)

(k

− ∑

+ 1) 1

(−1)

β5Mχ2m_2s

s = 0

b5Mχ2m_2k := β5Mχ2m_2k +

m_2

ρ10χ2(τ ) :=

∑

k = 0

⌠∞

P7(k ,τ ,χminMm_2−2) e− λ τ cos (ω τ ) dτ

β5Mχm_2k

:= 2 χminMm_2−2 (k + 1)

⌡0

k

m

s

2 (−1)

(k

+

∑

1) 1 −

(−1)

β5Mχm_2s

s = 0

b5Mχm_2k := β5Mχm_2k +

1

ρ10χ1 (τ)

ρ10χ2 (τ) 0.5

ρ10χ (τ)

0

0

1

ρt(τ)

2

− 0.5

− 1

τ

k

m

s

2 (−1)

(k + 1) 1 − ∑

(−1)

βM s

b5(k ,γ ) := βM k +

s = 0

(m + 1) (m + 2)

∑

k + 1

γ π

k = 0

k −1

otherwise

2

∑

φ (ω,s)

s = 0

λ

1

1

St(ω) := 2π

2

+ (ω − ω0)

2

+ (ω + ω0)

2

λ

λ

ω := 0,0 + 0.1.. 10

0.2

0.16

Sx(ω)

0.12

Snx(ω)

St(ω) 0.08

0.04

0

2

4

6

8

10

ω

Проверка условия нормировки модели нормированной корреляционной

функции при ограничении на модель.

мощности. Проверить условие нормировки.

⌠

∞

β5s (k ,γ ) := 2γ (k + 1)

P7(k ,ω,γ ) St(ω) dω

⌡0

βsM k := β5s (k ,γ )

m

βsM

γ2 − ∑ s + 1s

c5(k ,γ ) := βsM k

+

s = 0

m

∑ s +1

1

s = 0

m

Sa(ω) := ∑ (c5(k ,γ ) P7(k ,ω,γ ))

k = 0

0.2

0.16

Sa(ω)

0.12

Snx(ω)

St(ω) 0.08

0.04

0

2

4

6

8

10

ω

m

c5(k ,γ )

2 ∑

= 0.5

2γ (k + 1)

k = 0

1

m

k

s

ω

Fx(ω) :=

π ∑

b5(k) ∑

combin(k ,s) combin(k + s + 1,s + 1)

(−1) atan

γ (s + 1)

k = 0

s = 0

Ft(ω) :=

1

(φ1 (ω) + φ2 (ω))

2π

ω := 0,0 + 0.1.. 15

0.5

0.5

0.4

Fx(ω)0.3

Ft(ω)

0.2

0.1

0

3

6

9

12

15

ω

ортогональной моделью. Построить график.

γ пр(ω) := Fx(ω) − Ft(ω)

0.5

0.01

5×10−

3

γ пр(ω)

0

5

10

15

− 5×10−

3

− 0.01

− 0.015

ω