Имеет ли решение эта система уравнений?

@Halford · Регистрация: 19.11.2014

Уравнений, как и переменных - 19. Вопрос - можно ли найти эти переменные?

Matlab M

>> syms t11 t12 t13 t14 t22 t23 t31 t33 t34 t42 x1 y1 x2 y2 x3 y3 x4 y4 z real
Solution = solve (-log(t11)-1-x1-y1-5*z,-log(t12)-1-x1-y2-5*z,-log(t13)-1-x1-y3-5*z,-log(t14)-1-x1-y4-5*z,-log(t22)-1-x2-y2-3*z,-log(t23)-1-x2-y3-3*z,-log(t31)-1-x3-y1-5*z,-log(t33)-1-x3-y3-5*z,-log(t34)-1-x3-y4-5*z,-log(t42)-1-x4-y2-5*z,40-t11-t12-t13-t14,20-t11-t31,10-t22-t23,25-t12-t22-t42,30-t31-t33-t34,25-t13-t23-t33,10-t42,20-t14-t34,430-5*t11-5*t12-5*t13-5*t14-3*t22-3*t23-5*t31-5*t33-5*t34-5*t42);
t11 = eval(Solution.t11);
t12 = eval(Solution.t12);
t13 = eval(Solution.t13);
t14 = eval(Solution.t14);
t22 = eval(Solution.t22);
t23 = eval(Solution.t23);
t31 = eval(Solution.t31);
t33 = eval(Solution.t33);
t34 = eval(Solution.t34);
t42 = eval(Solution.t42);
x1 = eval(Solution.x1);
y1 = eval(Solution.y1);
x2 = eval(Solution.x2);
y2 = eval(Solution.y2);
x3 = eval(Solution.x3);
y3 = eval(Solution.y3);
x4 = eval(Solution.x4);
y4 = eval(Solution.y4);
z = eval(Solution.z);

Warning: Explicit solution could not be found.
> In solve at 81
??? Comma separated list expansion has cell syntax for an array that
is not a cell.

Error in ==> sym.sym>sym.subsref at 1364
[inds{k},refs{k}] = privformat(inds{k});

@Nata1979 · 24.11.2014, 20:13

Если не решается в символьной форме, то попробуй решить в числовой форме с использованием
функции fsolve с начальной точкой (пробуй задавать).
И смотри выходной аргумент exitflag.

@Halford · 25.11.2014, 19:31 **[ТС]**

Я пробовал в числовой форме. Так как я новичок, то в приведенном ниже примере (несколько модифицированном в переменных) есть ошибки. Начальной точки нет по условию. Можно ли задавать самому?

>> function OdeMy
%
[t,y]=ode15s(odefun1, [], []);
%
function dy=odefun1(t,y)
%
t11=y(1);
t12=y(2);
t13=y(3);
t14=y(4);
t22=y(5);
t23=y(6);
t31=y(7);
t33=y(8);
t34=y(9);
t42=y(10);
lambda11=y(11);
lambda12=y(12);
lambda21=y(13);
lambda22=y(14);
lambda31=y(15);
lambda32=y(16);
lambda41=y(17);
lambda42=y(18);
beta=y(19);
%
dy=zeros(19,1);
%
dy(1)=(-log(t11)-1-lambda11-lambda12-5*beta);
dy(2)=(-log(t12)-1-lambda11-lambda22-5*beta);
dy(3)=(-log(t13)-1-lambda11-lambda32-5*beta);
dy(4)=(-log(t14)-1-lambda11-lambda42-5*beta);
dy(5)=(-log(t22)-1-lambda21-lambda22-3*beta);
dy(6)=(-log(t23)-1-lambda21-lambda32-3*beta);
dy(7)=(-log(t31)-1-lambda31-lambda12-5*beta);
dy(8)=(-log(t33)-1-lambda31-lambda32-5*beta);
dy(9)=(-log(t34)-1-lambda31-lambda42-5*beta);
dy(10)=(-log(t42)-1-lambda41-lambda22-5*beta);
dy(11)=(40-t11-t12-t13-t14);
dy(12)=(20-t11-t31);
dy(13)=(10-t22-t23);
dy(14)=(25-t12-t22-t42);
dy(15)=(30-t31-t33-t34);
dy(16)=(25-t13-t23-t33);
dy(17)=(10-t42);
dy(18)=(20-t14-t34);
dy(19)=(430-5*t11-5*t12-5*t13-5*t14-3*t22-3*t23-5*t31-5*t33-5*t34-5*t42)
??? function OdeMy
|
Error: Function definitions are not permitted at the prompt or in scripts.

@Nata1979 0 / 0 / 1 Регистрация: 24.11.2014 Сообщений: 8
	24.11.2014, 20:13	2
	Сообщение было отмечено Зосима как решение Решение Если не решается в символьной форме, то попробуй решить в числовой форме с использованием функции fsolve с начальной точкой (пробуй задавать). И смотри выходной аргумент exitflag. 0

@Halford 0 / 0 / 0 Регистрация: 19.11.2014 Сообщений: 3
	25.11.2014, 19:31 [ТС]	3
	Я пробовал в числовой форме. Так как я новичок, то в приведенном ниже примере (несколько модифицированном в переменных) есть ошибки. Начальной точки нет по условию. Можно ли задавать самому? >> function OdeMy % [t,y]=ode15s(odefun1, [], []); % function dy=odefun1(t,y) % t11=y(1); t12=y(2); t13=y(3); t14=y(4); t22=y(5); t23=y(6); t31=y(7); t33=y(8); t34=y(9); t42=y(10); lambda11=y(11); lambda12=y(12); lambda21=y(13); lambda22=y(14); lambda31=y(15); lambda32=y(16); lambda41=y(17); lambda42=y(18); beta=y(19); % dy=zeros(19,1); % dy(1)=(-log(t11)-1-lambda11-lambda12-5beta); dy(2)=(-log(t12)-1-lambda11-lambda22-5beta); dy(3)=(-log(t13)-1-lambda11-lambda32-5beta); dy(4)=(-log(t14)-1-lambda11-lambda42-5beta); dy(5)=(-log(t22)-1-lambda21-lambda22-3beta); dy(6)=(-log(t23)-1-lambda21-lambda32-3beta); dy(7)=(-log(t31)-1-lambda31-lambda12-5beta); dy(8)=(-log(t33)-1-lambda31-lambda32-5beta); dy(9)=(-log(t34)-1-lambda31-lambda42-5beta); dy(10)=(-log(t42)-1-lambda41-lambda22-5beta); dy(11)=(40-t11-t12-t13-t14); dy(12)=(20-t11-t31); dy(13)=(10-t22-t23); dy(14)=(25-t12-t22-t42); dy(15)=(30-t31-t33-t34); dy(16)=(25-t13-t23-t33); dy(17)=(10-t42); dy(18)=(20-t14-t34); dy(19)=(430-5t11-5t12-5t13-5t14-3t22-3t23-5t31-5t33-5t34-5t42) ??? function OdeMy \| Error: Function definitions are not permitted at the prompt or in scripts. 0

Опции темы

Имеет ли решение эта система уравнений?

Решение