2014 dxdy logo

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

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




 
 Помогите решить систему в Matlab или Mathematica
Сообщение27.01.2011, 16:48 
Здравствуйте,
помогите, пожалуйста, решить численно систему из 32 уравнений и 16 неравенств и нескольких условий.
Пытался решить с Математикой- не справляется, долго думает и не даёт результата..

Переменные a1,a2...a16 и b1,b2...b16 и q
все положительные, для q принято значение 0.3.
Мы ищем значения переменных a1,a2...a16 и b1,b2...b16.
Значок && это AND.

Если для этой системы есть бесконечное количество решений,
то можно принять a1=0.1 (если это поможет, конечно).

Буду очень благодарен любой помощи,
не могу решить это уже месяц.

=================================================
(a1 == 44/(b10 + b11 + b13 + b14 + b15 + b2 + b5 + b6 + b9 + b12 q +
b13 q + b16 q + b3 q + b4 q + b6 q + b7 q + b8 q +
b9 q)) &&

(b1 <
17/(a12 + a13 + a16 + a3 + a4 + a6 + a7 + a8 + a9 + a11 q + a13 q +
a14 q + a15 q + a2 q + a5 q + a6 q + a9 q)) &&

(a10 ==
1/(b1 q) (17 - a12 b1 - a13 b1 - a16 b1 - a3 b1 - a4 b1 - a6 b1 -
a7 b1 - a8 b1 - a9 b1 - a11 b1 q - a13 b1 q - a14 b1 q -
a15 b1 q - a2 b1 q - a5 b1 q - a6 b1 q - a9 b1 q)) &&

(a2 == 35/(b11 + b13 + b14 + b15 + b16 + b3 + b4 + b5 + b7 + b1 q +
b10 q + b11 q + b12 q + b3 q + b5 q + b6 q + b7 q + b8 q +
b9 q)) &&

(b2 <
15/(a10 + a11 + a12 + a3 + a5 + a6 + a7 + a8 + a9 + a11 q + a13 q +
a14 q + a15 q + a16 q + a3 q + a4 q + a5 q + a7 q)) &&

(a1 ==
1/b2 (15 - a10 b2 - a11 b2 - a12 b2 - a3 b2 - a5 b2 - a6 b2 -
a7 b2 - a8 b2 - a9 b2 - a11 b2 q - a13 b2 q - a14 b2 q -
a15 b2 q - a16 b2 q - a3 b2 q - a4 b2 q - a5 b2 q -
a7 b2 q)) &&

(a3 == 28/(b1 + b10 + b12 + b15 + b2 + b6 + b7 + b8 + b9 + b11 q +
b12 q + b13 q + b14 q + b15 q + b16 q + b2 q + b4 q + b5 q +
b8 q)) &&

(b3 <
14/(a11 + a12 + a13 + a14 + a15 + a16 + a2 + a4 + a5 + a8 + a10 q +
a12 q + a15 q + a2 q + a6 q + a7 q + a8 q + a9 q)) &&

(a1 ==
1/(b3 q) (14 - a11 b3 - a12 b3 - a13 b3 - a14 b3 - a15 b3 -
a16 b3 - a2 b3 - a4 b3 - a5 b3 - a8 b3 - a10 b3 q - a12 b3 q -
a15 b3 q - a2 b3 q - a6 b3 q - a7 b3 q - a8 b3 q -
a9 b3 q)) &&

(a4 == 33/(b1 + b12 + b14 + b15 + b3 + b6 + b7 + b8 + b9 + b10 q +
b11 q + b13 q + b14 q + b15 q + b16 q + b2 q + b5 q +
b8 q)) &&

(b4 <
20/(a10 + a11 + a13 + a14 + a15 + a16 + a2 + a5 + a8 + a12 q +
a14 q + a15 q + a3 q + a6 q + a7 q + a8 q + a9 q)) &&

(a1 ==
1/(b4 q) (20 - a10 b4 - a11 b4 - a13 b4 - a14 b4 - a15 b4 -
a16 b4 - a2 b4 - a5 b4 - a8 b4 - a12 b4 q - a14 b4 q -
a15 b4 q - a3 b4 q - a6 b4 q - a7 b4 q - a8 b4 q -
a9 b4 q)) &&

(a5 == 30/(b12 + b13 + b14 + b15 + b16 + b2 + b3 + b4 + b7 + b1 q +
b10 q + b11 q + b12 q + b2 q + b6 q + b7 q + b8 q +
b9 q)) &&

(b5 <
22/(a10 + a11 + a12 + a2 + a6 + a7 + a8 + a9 + a12 q + a13 q +
a14 q + a15 q + a16 q + a2 q + a3 q + a4 q + a7 q)) &&

(a1 ==
1/b5 (22 - a10 b5 - a11 b5 - a12 b5 - a2 b5 - a6 b5 - a7 b5 -
a8 b5 - a9 b5 - a12 b5 q - a13 b5 q - a14 b5 q - a15 b5 q -
a16 b5 q - a2 b5 q - a3 b5 q - a4 b5 q - a7 b5 q)) &&

(a6 == 34/(b1 + b10 + b13 + b14 + b15 + b16 + b2 + b5 + b8 + b1 q +
b11 q + b12 q + b14 q + b16 q + b3 q + b4 q + b7 q + b8 q +
b9 q)) &&

(b6 <
27/(a11 + a12 + a14 + a16 + a3 + a4 + a7 + a8 + a9 + a10 q +
a13 q + a14 q + a15 q + a16 q + a2 q + a5 q + a8 q)) &&

(a1 ==
1/(b6 + b6 q) (27 - a11 b6 - a12 b6 - a14 b6 - a16 b6 - a3 b6 -
a4 b6 - a7 b6 - a8 b6 - a9 b6 - a10 b6 q - a13 b6 q - a14 b6 q -
a15 b6 q - a16 b6 q - a2 b6 q - a5 b6 q - a8 b6 q)) &&

(a7 == 22/(b1 + b10 + b11 + b12 + b2 + b5 + b6 + b8 + b9 + b11 q +
b12 q + b13 q + b14 q + b15 q + b16 q + b2 q + b3 q + b4 q +
b5 q)) &&

(b7 <
21/(a11 + a12 + a13 + a14 + a15 + a16 + a2 + a3 + a4 + a5 + a10 q +
a11 q + a12 q + a2 q + a5 q + a6 q + a8 q + a9 q)) &&

(a1 ==
1/(b7 q) (21 - a11 b7 - a12 b7 - a13 b7 - a14 b7 - a15 b7 -
a16 b7 - a2 b7 - a3 b7 - a4 b7 - a5 b7 - a10 b7 q - a11 b7 q -
a12 b7 q - a2 b7 q - a5 b7 q - a6 b7 q - a8 b7 q -
a9 b7 q)) &&

(a8 == 26/(b1 + b10 + b11 + b15 + b2 + b3 + b4 + b5 + b6 + b9 +
b10 q + b12 q + b13 q + b14 q + b16 q + b3 q + b4 q + b6 q +
b7 q)) &&

(b8 <
29/(a10 + a12 + a13 + a14 + a16 + a3 + a4 + a6 + a7 + a10 q +
a11 q + a15 q + a2 q + a3 q + a4 q + a5 q + a6 q +
a9 q)) &&

(a1 ==
1/(b8 q) (29 - a10 b8 - a12 b8 - a13 b8 - a14 b8 - a16 b8 - a3 b8 -
a4 b8 - a6 b8 - a7 b8 - a10 b8 q - a11 b8 q - a15 b8 q -
a2 b8 q - a3 b8 q - a4 b8 q - a5 b8 q - a6 b8 q - a9 b8 q)) &&

(a9 == 22/(b1 + b10 + b11 + b13 + b14 + b15 + b16 + b2 + b5 + b6 +
b1 q + b11 q + b12 q + b13 q + b16 q + b3 q + b4 q + b7 q +
b8 q)) &&

(b9 <
26/(a11 + a12 + a13 + a16 + a3 + a4 + a7 + a8 + a10 q + a11 q +
a13 q + a14 q + a15 q + a16 q + a2 q + a5 q + a6 q)) &&

(a1 ==
1/(b9 + b9 q) (26 - a11 b9 - a12 b9 - a13 b9 - a16 b9 - a3 b9 -
a4 b9 - a7 b9 - a8 b9 - a10 b9 q - a11 b9 q - a13 b9 q -
a14 b9 q - a15 b9 q - a16 b9 q - a2 b9 q - a5 b9 q -
a6 b9 q)) &&

(a10 == 21/(b12 + b13 + b14 + b15 + b16 + b2 + b4 + b5 + b8 + b1 q +
b11 q + b12 q + b14 q + b15 q + b3 q + b6 q + b7 q + b8 q +
b9 q)) &&

(b10 <
24/(a11 + a12 + a14 + a15 + a3 + a6 + a7 + a8 + a9 + a12 q +
a13 q + a14 q + a15 q + a16 q + a2 q + a4 q + a5 q +
a8 q)) &&

(a1 ==
1/b10 (24 - a11 b10 - a12 b10 - a14 b10 - a15 b10 - a3 b10 -
a6 b10 - a7 b10 - a8 b10 - a9 b10 - a12 b10 q - a13 b10 q -
a14 b10 q - a15 b10 q - a16 b10 q - a2 b10 q - a4 b10 q -
a5 b10 q - a8 b10 q)) &&

(a11 == 16/(b10 + b13 + b16 + b2 + b3 + b4 + b5 + b6 + b7 + b9 +
b1 q + b12 q + b14 q + b15 q + b16 q + b2 q + b7 q + b8 q +
b9 q)) &&

(b11 <
21/(a12 + a14 + a15 + a16 + a2 + a7 + a8 + a9 + a10 q + a13 q +
a16 q + a2 q + a3 q + a4 q + a5 q + a6 q + a7 q +
a9 q)) &&

(a1 ==
1/b11 (21 - a12 b11 - a14 b11 - a15 b11 - a16 b11 - a2 b11 -
a7 b11 - a8 b11 - a9 b11 - a10 b11 q - a13 b11 q - a16 b11 q -
a2 b11 q - a3 b11 q - a4 b11 q - a5 b11 q - a6 b11 q -
a7 b11 q - a9 b11 q)) &&

(a12 == 23/(b1 + b10 + b11 + b2 + b3 + b5 + b6 + b7 + b8 + b9 +
b10 q + b13 q + b14 q + b15 q + b16 q + b3 q + b4 q + b5 q +
b7 q)) &&

(b12 <
27/(a10 + a13 + a14 + a15 + a16 + a3 + a4 + a5 + a7 + a10 q +
a11 q + a2 q + a3 q + a5 q + a6 q + a7 q + a8 q +
a9 q)) &&

(a1 ==
1/(b12 q) (27 - a10 b12 - a13 b12 - a14 b12 - a15 b12 - a16 b12 -
a3 b12 - a4 b12 - a5 b12 - a7 b12 - a10 b12 q - a11 b12 q -
a2 b12 q - a3 b12 q - a5 b12 q - a6 b12 q - a7 b12 q -
a8 b12 q - a9 b12 q)) &&

(a13 == 18/(b1 + b12 + b14 + b16 + b3 + b4 + b7 + b8 + b9 + b1 q +
b10 q + b11 q + b14 q + b15 q + b16 q + b2 q + b5 q + b6 q +
b9 q)) &&

(b13 <
30/(a10 + a11 + a14 + a15 + a16 + a2 + a5 + a6 + a9 + a12 q +
a14 q + a16 q + a3 q + a4 q + a7 q + a8 q + a9 q)) &&

(a1 ==
1/(b13 +
b13 q) (30 - a10 b13 - a11 b13 - a14 b13 - a15 b13 - a16 b13 -
a2 b13 - a5 b13 - a6 b13 - a9 b13 - a12 b13 q - a14 b13 q -
a16 b13 q - a3 b13 q - a4 b13 q - a7 b13 q - a8 b13 q -
a9 b13 q)) &&

(a14 == 19/(b10 + b11 + b12 + b13 + b16 + b3 + b4 + b6 + b7 + b8 +
b1 q + b10 q + b13 q + b15 q + b2 q + b4 q + b5 q + b6 q +
b9 q)) &&

(b14 <
37/(a10 + a13 + a15 + a2 + a4 + a5 + a6 + a9 + a10 q + a11 q +
a12 q + a13 q + a16 q + a3 q + a4 q + a6 q + a7 q +
a8 q)) &&

(a1 ==
1/b14 (37 - a10 b14 - a13 b14 - a15 b14 - a2 b14 - a4 b14 -
a5 b14 - a6 b14 - a9 b14 - a10 b14 q - a11 b14 q - a12 b14 q -
a13 b14 q - a16 b14 q - a3 b14 q - a4 b14 q - a6 b14 q -
a7 b14 q - a8 b14 q)) &&

(a15 == 11/(b10 + b11 + b12 + b13 + b14 + b16 + b3 + b4 + b7 + b1 q +
b10 q + b2 q + b3 q + b4 q + b5 q + b6 q + b8 q +
b9 q)) &&

(b15 <
30/(a10 + a2 + a3 + a4 + a5 + a6 + a8 + a9 + a10 q + a11 q +
a12 q + a13 q + a14 q + a16 q + a3 q + a4 q + a7 q)) &&

(a1 ==
1/b15 (30 - a10 b15 - a2 b15 - a3 b15 - a4 b15 - a5 b15 - a6 b15 -
a8 b15 - a9 b15 - a10 b15 q - a11 b15 q - a12 b15 q -
a13 b15 q - a14 b15 q - a16 b15 q - a3 b15 q - a4 b15 q -
a7 b15 q)) &&

(a16 == 13/(b1 + b11 + b12 + b13 + b3 + b4 + b6 + b7 + b8 + b9 +
b10 q + b11 q + b13 q + b14 q + b15 q + b2 q + b5 q + b6 q +
b9 q)) &&

(b16 <
35/(a10 + a11 + a13 + a14 + a15 + a2 + a5 + a6 + a9 + a11 q +
a12 q + a13 q + a3 q + a4 q + a6 q + a7 q + a8 q +
a9 q)) &&

(a1 ==
1/(b16 q) (35 - a10 b16 - a11 b16 - a13 b16 - a14 b16 - a15 b16 -
a2 b16 - a5 b16 - a6 b16 - a9 b16 - a11 b16 q - a12 b16 q -
a13 b16 q - a3 b16 q - a4 b16 q - a6 b16 q - a7 b16 q -
a8 b16 q - a9 b16 q)) &&

(a1 > 0 && a2 > 0 && a3 > 0 && a4 > 0 && a5 > 0 && a6 > 0 && a7 > 0 &&
a8 > 0 && a9 > 0 && a10 > 0 && a11 > 0 && a12 > 0 && a13 > 0 &&
a14 > 0 && a15 > 0 && a16 > 0 && b1 > 0 && b2 > 0 && b3 > 0 &&
b4 > 0 && b5 > 0 && b6 > 0 && b7 > 0 && b8 > 0 && b9 > 0 &&
b10 > 0 && b11 > 0 && b12 > 0 && b13 > 0 && b14 > 0 && b15 > 0 &&
b16 > 0 && q > 0) &&

(q == 0.3)

=================================================

 
 
 
 Re: Помогите решить систему в Matlab или Mathematica
Сообщение30.01.2011, 20:55 
Аватара пользователя
а свой код из mathematica выложи

 
 
 [ Сообщений: 2 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group