Ошибка в функции

@Mutlu · Регистрация: 26.09.2011

Здравствуйте,
Проблема в том, что не могу понять где ошибка в самой функции.
Вот задание:
Система ОДУ (Модель Ресслера)
x'=-y-z
y'=x+Ay
z'=B+z(x-C)
Начальные условия (x(0) = -0.7, y(0) = -0.7, z(0) = 1)
Интервал времени(t = [0 50])
Параметры(A = 0.2, B = 0.2, C = 5.7)
Вот *.m файл – l.m

Matlab M

1 2	function y=l(t,x) y = [(-x(1)-x(2)); (x(3)+0.2.x(1)); (0.2+x(2)(x(3)-5.7))];

В командной строке ввожу:

Matlab M

x0=[-0.7 -0.7 1];
ts=[0 50];
[t,x]=ode45('l',ts,x0);
figure(1), plot(t,x),grid on
figure(2), comet3(x(:,1),x(:,2),x(:,3))
figure(2), comet3(x(:,1),x(:,2),x(:,3))

Всё вроде выводит, но препод говорит, что в функции ошибка, не могли бы проверит, где имеено там ошибка?
Спасибо!

Зосима · 18.01.2013, 12:06

Чтоб не путаться рекоммендую в комментариях писать где какая переменая что значит:

Matlab M

y=l(t,x)
% x(1) - это y
% x(2) - это dy/dt
% x(3) - это dy/dt
A = 0.2;
B = 0.2;
C = 5.7;
 
y= zeros(3,1);
y(1) = -x(2)-x(3);  
y(2) = x(1) - A*x(2);
y(3) = B + x(3)*(x(1)-C);
end

Как я вижу у тебя x(1) это y, x(2) - z, x(3) - x - может это запутало препода?

@Mutlu · 18.01.2013, 12:14 **[ТС]**

Сообщение от Зосима

Чтоб не путаться рекоммендую в комментариях писать где какая переменая что значит:

Matlab M

y=l(t,x)
% x(1) - это y
% x(2) - это dy/dt
% x(3) - это dy/dt
A = 0.2;
B = 0.2;
C = 5.7;
 
y= zeros(3,1);
y(1) = -x(2)-x(3);  
y(2) = x(1) - A*x(2);
y(3) = B + x(3)*(x(1)-C);
end

Как я вижу у тебя x(1) это y, x(2) - z, x(3) - x - может это запутало препода?

Так получается у меня всё правильно?

Зосима · 18.01.2013, 13:08

Mutlu, вроде да

только обозначено чуток по другому и нет пояснений

@Mutlu · 18.01.2013, 13:15 **[ТС]**

Сообщение от Зосима

Mutlu, вроде да

только обозначено чуток по другому и нет пояснений

Спасибо, напишу тогда пояснение.

Зосима 5210 / 3542 / 370 Регистрация: 02.04.2012 Сообщений: 6,455 Записей в блоге: 17
	18.01.2013, 13:08	4
	Mutlu, вроде да только обозначено чуток по другому и нет пояснений 1

@Mutlu 11 / 11 / 2 Регистрация: 26.09.2011 Сообщений: 185
	18.01.2013, 13:15 [ТС]	5
	Сообщение от Зосима Mutlu, вроде да только обозначено чуток по другому и нет пояснений Спасибо, напишу тогда пояснение. 0

Опции темы