2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 Matematika 5.2. Не решается система ОДУ. ГЛЮК? как бороться?
Сообщение05.03.2008, 16:41 


05/03/08
4
Ув. товарищи,
есть такая системка:

Vb==(L3 + L2 + 2*M23)*i32'[t] + (R3 + R2)*i32[t] + (emf3[t] - emf2[t]) + (L2 +
M12 + M23 - M13)*i12'[t] + R2*i12[t],
Vb==(L1 + L2 + 2*M12)*i12'[t] + (R1 + R2)*i12[t] + (emf1[t] - emf2[t]) + (L2 +
M23 + M12 - M13)*i32'[t] + R2*i32[t],i12[to] == i12o, i32[to] == i32o

Предварительно определены функции:
emf1[t_] = Kt11*w*Sin[1*(w*t + offset1)] + Kt13*w*Sin[3*(w*t + offset1)],
emf2[t_] и emf3[t_] в таком же духе.
Всем параметрам заданы числовые значения.

Если в DSolve подставить дифуры как они записаны, то получаем "Программа выполнила недопустимую операцию...". Если составить уравнения отдельно с подставленными параметрами, то получаем сообщение: "Unable to find principal vectors for eigenvalue -104.16666666666667`" - в хэлпе сказано - если эта ошибка возникает там где ее не должно быть - обратитесь к разработчику....

:( :(

помогите плиз.
П.С. Мэпл такое оч. быстренько решает. но нужно в математике.

 Профиль  
                  
 
 
Сообщение13.03.2008, 20:10 
Аватара пользователя


15/01/06
200
Вы бы приводили точно код, который пытаетесь выполнить в Математике, а то сидеть и руками определять все коэффициенты и функции, пока вобьешь все и желание искать решение проблемы отпадет :wink:

 Профиль  
                  
 
 вот пожалуйста код
Сообщение24.03.2008, 11:32 


05/03/08
4
Извиняюсь за задержку с кодом. Вот:

Код:
Kt11 = 0.001102658; Kt13 = 0.0001102658;
Kt21 = Kt11; Kt23 = Kt13;
Kt31 = Kt11; Kt33 = Kt13;
L1 = 0.00005; L2 = L1; L3 = L1;
R1 = 0.005; R2 = R1; R3 = R1;
M12 = (L1 + L2)/n; M13 = (L1 + L3)/n; M23 = (L2 + L3)/n; n = 10;
offset1 = Pi/6; offset2 = offset1 + 2*Pi/3; offset3 = offset1 + 4*Pi/3;
RPM = 60; w = 2*Pi*RPM/60; Vb = 10; T = 2*Pi/w;
emf1[t_] = Kt11*w*Sin[1*(w*t + offset1)] + Kt13*w*Sin[3*(w*t + offset1)];
emf2[t_] = Kt21*w*Sin[1*(w*t + offset2)] + Kt23*w*Sin[3*(w*t + offset2)];
emf3[t_] = Kt31*w*Sin[1*(w*t + offset3)] + Kt33*w*Sin[3*(w*t + offset3)];
tokon20 = DSolve[{Vb == (L3 + L2 + 2*M23)*i32'[
        t] + (R3 + R2)*i32[t] + (emf3[t] - emf2[t]) + (L2 + M12 +
           M23 - M13)*i12'[t] + R2*i12[t],
          Vb == (L1 + L2 + 2*M12)*i12'[t] + (R1 + R2)*i12[
          t] + (emf1[t] - emf2[t]) + (L2 + M23 + M12 - M13)*i32'[t] + R2*i32[
          t]}, {i32[t], i12[t]}, t]
tokon21 = DSolve[{Vb == (L3 + L2 + 2*M23)*i32'[
        t] + (R3 + R2)*i32[t] + (emf3[t] - emf2[t]) + (L2 + M12 +
           M23 - M13)*i12'[t] + R2*i12[t],
          Vb == (L1 + L2 + 2*M12)*i12'[t] + (R1 + R2)*i12[
          t] + (emf1[t] - emf2[t]) + (L2 + M23 + M12 - M13)*i32'[t] + R2*i32[
          t], i12[to] == i12o, i32[to] == i32o}, {i32[t], i12[t]}, t]


Вот 2 последних DSolve не работают почемуто. ХЕЛП...

 Профиль  
                  
 
 
Сообщение24.03.2008, 15:59 
Заслуженный участник


22/01/07
605
Все получилось при замене данных в десятичном формате на точные:
Код:
Kt11 = 1102658 10^-9; Kt13 = 1102658 10^-10;
Kt21 = Kt11; Kt23 = Kt13;
Kt31 = Kt11; Kt33 = Kt13;
L1 = 5 10^-5; L2 = L1; L3 = L1;
R1 = 5 10^-3; R2 = R1; R3 = R1;

 Профиль  
                  
 
 С П А С И Б О !!!!!!!!!!!!!!!!!!!!!
Сообщение25.03.2008, 10:45 


05/03/08
4
Gafield - спасибо огромное!!!! :D :D

я и не знал о такой фиче как возможность представления в виде дроби.

и если позволите тогда еще вопрос - есть ли функция встроенная для перевода числа с точкой в натуральную дробь.?

Еще раз спасибо.

П.С. Кажется нашел - Rationalize[...]

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 5 ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: dgwuqtj


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group