ÐиÑвÑÑкий наÑÑоналÑний ÑнÑвеÑÑиÑÐµÑ - ÐеологÑÑний ÑакÑлÑÑеÑ
ÐиÑвÑÑкий наÑÑоналÑний ÑнÑвеÑÑиÑÐµÑ - ÐеологÑÑний ÑакÑлÑÑеÑ
ÐиÑвÑÑкий наÑÑоналÑний ÑнÑвеÑÑиÑÐµÑ - ÐеологÑÑний ÑакÑлÑÑеÑ
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.