10.01.2015 Views

Київський національний університет - Геологічний факультет

Київський національний університет - Геологічний факультет

Київський національний університет - Геологічний факультет

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Тема_14_MathCad_7.doc 65<br />

Sumkv ( C)<br />

:=<br />

n−1<br />

∑<br />

i<br />

=<br />

0<br />

( ( )) 2<br />

yy i − Signal xx i , C<br />

- сума квадратів відхилень<br />

Функція-підпрограма знаходженя методом випадкових випробувань того значення<br />

XminVek, при якому у m-вимірному паралелепіпеді, при якому приймає найменше значення<br />

задана функція Sumkv. Протилежні вершини паралелепіпеда задані векторами Cmin та<br />

Cmax. Nprob - кількість статистичних випробувань, m - розмірність простору, у якому<br />

відшукується точка мінімуму .<br />

XminVek( Sumkv , Cmin, Cmax, Nprob , m) := Xmin ← Cmin<br />

Ymin ←<br />

for<br />

return<br />

Sumkv ( Cmin)<br />

i ∈ 1..<br />

Nprob<br />

X ← Rndvekt ( Cmin, Cmax,<br />

m)<br />

p ←<br />

Sumkv ( X)<br />

Xmin ← X if p < Ymin<br />

Ymin ← p if p < Ymin<br />

Xmin<br />

Обчислення оцінок параметрів С 0<br />

, С 1<br />

, С 2 , С 3<br />

:<br />

Cmin :=<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

−1<br />

2<br />

70<br />

3<br />

⎞ ⎟<br />

⎟<br />

⎟<br />

⎟<br />

⎠<br />

- нижні границі<br />

параметрів<br />

сигналу<br />

Cmax :=<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

2<br />

10<br />

90<br />

9<br />

⎞<br />

⎟<br />

⎟<br />

⎟<br />

⎟<br />

⎠<br />

- верхні границі<br />

параметрів<br />

сигналу<br />

m := 4<br />

- кількість параметрів сигналу<br />

Nprob := 4000<br />

- кількість випобувань для визначення оцінок параметрів спостережень при<br />

мінімізації суми квадратів відхилень<br />

Cpr := XminVek( Sumkv , Cmin, Cmax, Nprob , m)<br />

Cpr =<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

0.406<br />

6.107<br />

79.056<br />

7.603<br />

⎞<br />

⎟<br />

⎟<br />

⎟<br />

⎟<br />

⎠<br />

x := 0, 2..<br />

110<br />

- фон<br />

- амплитуда<br />

- положение<br />

- размах<br />

- оцінка параметрів сигналу методом<br />

найменших квадратів<br />

Sumkv ( Cpr) = 2.235568<br />

Рис. 77. Приклад використання підпрограми-функції для обчислення оцінки відразу кількох параметрів<br />

корисного сигналу методом найменших квадратів. Частина 2.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!