| Цитата(W4FhLF @ 27.3.2010, 19:00) | ты за 3 дня в гугл своё вопрос вбить не пытался?
http://drobotenko.com/code_rus.html
|
Ну как же не пытался, я все перепробовал. Дело в том что тут как показано в примере програмирование без main, так что тот код в ссылке не работает. А во вторых у меня все время выкидывает ошибку: template <int N> error C2143: syntax error : missing ',' before 'constant' и я не пойму в чем дело.
Я пытался сделать свое собственный код, но он работает только на небольших матрицах:
| Код |
MyMatrix InvertMatrix(MyMatrix& CovMatrix)
{
int N=CovMatrix.rows();
int j, i, r;
double factor;
MyMatrix Inverted(N,N);
Inverted=CreateUnitMatrix(N);
for (i=0; i < N; i++)
{
factor=1.0/CovMatrix(i,i); //Making sure that the Diagonal is 1.
for (j=0; j < N; j++)
{ CovMatrix(i,j)=CovMatrix(i,j)*factor;
Inverted(i,j)=Inverted(i,j)*factor;
}
//Making zeros in column i
for (r=0; r < N; r++)
{
double Numb=CovMatrix(r,i);
for (size_t j=0; j < CovMatrix.rows(); j++)
{
if (r==i) {}//Do nothing (this is the row used for creating the other zeros
else
{
CovMatrix(r,j)=CovMatrix(r,j)-Numb*CovMatrix(i,j);
Inverted(r,j)=Inverted(r,j)-Numb*Inverted(i,j);
}
}
}
}
RoundingZeros(Inverted);
return Inverted;
}
MyMatrix CreateUnitMatrix(int Dimension)
{
MyMatrix UnitMatrix(Dimension,Dimension);
for(int j=0; j<Dimension; j++)
{
for(int i=0; i<Dimension; i++)
{
if (j==i)
{UnitMatrix(i,j)=1;}
else
{UnitMatrix(i,j)=0;}
}
}
return UnitMatrix;
}
MyMatrix RoundingZeros(MyMatrix& Matrix)
{
for (size_t i=0; i < Matrix.rows(); i++)
{
for (size_t j=0; j < Matrix.columns(); j++)
{
if(fabs(Matrix(i,j))<0.000000000001)
{Matrix(i,j)=0.0;}
}
}
return Matrix;
|
Тут похоже не хватает пивотирования, то бишь если в диагонале стоит 0 то просто надо перепрыгивать этот ряд. У меня уже просто нет фантазии на это дело )) Если кто может помоч советом буду очень признателен !!! |
