Matematika 5.2. Не решается система ОДУ. ГЛЮК? как бороться?

[Konstantin_m25]

Ув. товарищи,
есть такая системка:

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`" - в хэлпе сказано - если эта ошибка возникает там где ее не должно быть - обратитесь к разработчику....



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

[Leierkastenmann]

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

[Konstantin_m25]

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

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 не работают почемуто. ХЕЛП...

[Gafield]

Все получилось при замене данных в десятичном формате на точные:

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;

[Konstantin_m25]

Gafield - спасибо огромное!!!!

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

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

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

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