mathlab решение системы ду 2го порядка.помогите найти ошибку

Apathy · 4 May 2012

Apathy said:
где v1,v2-предельные скорости тела

мю здесь динамическая проницаемость среды
r-радиус тела

здесь с-безразмерный коэффициент формы,для шара =0.4
p0-плотность среды
p-плотность тела
Помогите пожалуйста составить и решить систему средствами матлаб
сам пытался делать так
Code:
global V0 alpha
% формирование вектора начальных условий
Y0 = [0; V0*cos(alpha);0;V0*sin(alpha)];
% вызов солвера от файл-функции oscil, начального и конечного
% момента времени и вектора начальных условий
[T, Y] = ode45('oscil', [0 5], Y0);
Code:
function F = oscil(t, y)
global m alpha 
v1=m*9.8/(6*pi*0.0182*0.1);
v2=sqrt(2*m*9.8/(7.8*1000));
F = [y(2);...
-9.8*(abs(sqrt(y(2)^2+y(3)^2))./v1+(sqrt(y(2)^2+y(3)^2)./v2)^2).*cos(alpha);...
y(3);...
-9.8-9.8*(abs(sqrt(y(2)^2+y(3)^2))./v1+(sqrt(y(2)^2+y(3)^2)./v2)^2).*sign(y(3)*sin(alpha))*sin(alpha)];
спасибо за помощь
ps мой вариант работает не верно, тело улетает в бесконечность. [/img]
Click to expand...

Useful Searches

mathlab решение системы ду 2го порядка.помогите найти ошибку

Apathy New Member