Математическая модель метода главных компонент

| Категория реферата: Рефераты по статистике
| Теги реферата: ответы по контрольной, курсовые работы бесплатно
| Добавил(а) на сайт: Serdjukov.

_b_[i,j]:=-2; for q:=2 to m do
{----вычисляем p[q] и определитель-----} begin

{----вычисляем A[q]----} for i:=1 to m do for j:=1 to m do begin s:=0; for k:= 1 to m do s:=s+r[i,k]*_b_[k,j]; a_[i,j]:=s end;

{------вычисляем p[q]-------} s:=0; for i:=1 to m do s:=s+a_[i,i]; p[q]:=s/q;

{----вычисляем B[q]-----} for i:=1 to m do for j:=1 to m do if ij then b_[i,j]:=a_[i,j] else b_[i,j]:=a_[i,j]-p[q];

{----присваиваем предыдущим переменным значения текущих-----} for i:= 1 to m do for j:=1 to m do begin

_a_[i,j]:=a_[i,j];

_b_[i,j]:=b_[i,j] end end;
{---------===решаем характеристическое уравнение===----------} p[0]:=1; for i:=1 to m do p[i]:=-p[i]; for i:=1 to m do for j:=1 to m do l[i,j]:=0;
{------задаем начальные приближения------} for i:=1 to m do l[i,i]:=-p[i]/p[i-1];
{------выполняем итерационный процесс по методу Ньютона--------} repeat w:=0; for i:=1 to m do begin b:=0;_b:=0;

{-----вычисляем значение полинома в i-й точке-------} for j:=0 to m do begin s:=1; for k:=0 to m-j-1 do s:=s*l[i,i]; b:=b+p[j]*s end;

{------находим максимальную невязку-------} if b>w then w:=b;

{------вычисляем значение производной в i-й точке------} for j:=0 to m-1 do begin s:=1; for k:=0 to m-j-2 do s:=s*l[i,i];

_b:=_b+(m-j)*p[j]*s end;

{------вносим поправку для i-й точки-------} l[i,i]:=l[i,i]-(b/_b) end

{----выходим из процесса при достижении требуемой точности----} until w