Добро пожаловать, гость
:: алгоритмы  и методы :: :: олимпиадные задачи :: :: связь :: :: о сайте :: :: форум ::

Форум работает в режиме архива, только для чтения и поиска.
Архив 2004 Архив 2007 Архив 2013

 
 
Опции темы Поиск в этой теме Опции просмотра
  #1  
Старый 08.05.2009, 22:15
Новичок

Отправить личное сообщение для ciaonataha Посмотреть профиль Найти все сообщения от ciaonataha
 
Регистрация: 09.02.2009
Сообщений: 15

Помогите Матрица Жакобиана
Матрица Жакобиана
Помогите записать алгоритм на C++ - "матрицу Жакобиана для нелинеиной системы уравнении":
//Y - неизвестные переменные
//F - вектор уравнении
//equilibre_Fe_Cr_C(n,Y,F) - функция с уравнениями
У меня не получается. И самое страшное, что не могу найти ошибки.
Код:
Код:
void  Jacobian(std::vector <double> &Y,  std::vector <double> &F, std::vector < std::vector<double> > &Jacobian2D)
{

 int i,j;
 int n;
 int i1 = Y.size();
 int i2 = i1;
 n = i1;
 std::vector <double> Ynew(i2);
 std::vector < std::vector<double> > Fnew;

 Fnew.resize(i1);
 Jacobian2D.resize(i1);

 for (i = 0; i < Jacobian2D.size(); i++)
   {
     Fnew[i].resize(i2);
     Jacobian2D[i].resize(i2);
   }
 
 for( i=0 ; i < Jacobian2D.size() ; i++)
   for(j=0 ; j < Jacobian2D[i].size() ; j++) 
     {
       Ynew[j] = Y[j] + kdelta;
       equilibre_Fe_Cr_C(n,Ynew,F);
       Fnew[i][j] = F[i];
       equilibre_Fe_Cr_C(n,Y,F); 
       Jacobian2D[i][j] = (Fnew[i][j] - F[i])/kdelta; 
     }
}
  #2  
Старый 08.05.2009, 22:16
Новичок

Отправить личное сообщение для ciaonataha Посмотреть профиль Найти все сообщения от ciaonataha
 
Регистрация: 09.02.2009
Сообщений: 15

Код:
kdelta = 0.001
  #3  
Старый 24.12.2009, 15:59
гость

 
Сообщений: n/a

Жакобиан в штанах у тебя. Якобиан считай мб?
 


Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск
Опции просмотра


Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
мкэ: матрица жесткости стержня (растяжение + изгиб) BreakPoint Математические алгоритмы 2 21.01.2010 02:26
Матрица Жакобиана ciaonataha Задачи 0 08.05.2009 18:36
Матрица достижимости из матрицы инцидентности Лимонадный Джо Графы 5 02.09.2008 03:52
динамические массивы. матрица гость Математические алгоритмы (другое) 4 27.11.2007 11:38