Форум программистов, компьютерный форум, киберфорум
Наши страницы
Matlab
Войти
Регистрация
Восстановить пароль
 
klaxor2014
13 / 13 / 1
Регистрация: 13.02.2011
Сообщений: 101
1

Решение системы диф.уравнений (Движение снаряда)

13.11.2014, 21:07. Просмотров 479. Ответов 4
Метки нет (Все метки)

Здравствуйте, уважаемые форумчане, помогите пожалуйста разобраться в следующей проблеме.Решаю задачу о движении снаряда.Имеется следующая система: http://www.cyberforum.ru/cgi-bin/latex.cgi?\left\{\begin{matrix} &\frac{dx}{dt}={v}_{x} \\  &\frac{d{v}_{x}}{dt}=-\frac{k}{m}*{v}_{x} \\  &\frac{dy}{dt}={v}_{y} \\  &\frac{d{v}_{y}}{dt}=-\frac{k}{m}*{v}_{y}-g \end{matrix}\right.

Matlab M
1
2
3
4
5
6
7
8
function ds = syst(t,MP)
ds = zeros(4,1);
k=1;m=1;g=1;
ds(1)=MP(3)
ds(2)=MP(4);
ds(3)=-k/m*MP(3);
ds(4)=-k/m*MP(4)-g;
return;
Matlab M
1
2
3
clc;clear all;close all;
[T,Y] = ode45(@syst,[0 5],[1 1 1 1]);
figure;plot(T(:,1),Y(:,3),'x',T(:,1),Y(:,4),'y');grid on;

Не понятно, с каким индексами скармливать массив MP и какие массивы после вычислений ode выводить (ведь там выводятся графики производных также).
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
13.11.2014, 21:07
Ответы с готовыми решениями:

Решение системы диф.уравнений
Всем привет! Столкнулась с проблемой при решении системы обыкновенных...

Решение задачи Коши для системы 2 диф. уравнений методом Рунге-Кутты 4 порядка
Здравствуйте, помогите с написанием это программой, вообще не знаю как писать.

Решение диф. Уравнений
Друзья, хелп ми плиз... Найти аналитическое решение уравнения : y′ + y2 =...

Обычные системы диф. уравнений
Доброе время суток. Подскажите пожалуйста как нужно правильно расписать систему...

Решение диф.уравнений Mathlab
Помогите, пожалуйста, написать программу в Mathlab

4
nuHrBuH
253 / 253 / 150
Регистрация: 04.03.2011
Сообщений: 598
13.11.2014, 22:44 2
Matlab M
1
2
3
4
5
6
7
8
9
x = MP(1);
y = MP(2);
Vx = MP(3);
Vy = MP(4);
 
ds(1) = Vx;
ds(2) = Vy;
ds(3) = -k/m*Vx;
ds(4) = -k/m*Vy-g;
1
klaxor2014
13 / 13 / 1
Регистрация: 13.02.2011
Сообщений: 101
13.11.2014, 22:49  [ТС] 3
Спасибо, насколько понял в выходном массиве первые два столбца и есть искомые функции ускорения по осям x и y?
Не подскажете как построить график?Дело в том, что имею аналогичную систему в Simulink и там в качесве ответа выступают незатухающие гармонические функции, что я вообще не могу понять, как объяснить...
0
Зосима
4937 / 3310 / 313
Регистрация: 02.04.2012
Сообщений: 6,207
Записей в блоге: 15
Завершенные тесты: 1
14.11.2014, 11:34 4
Цитата Сообщение от klaxor2014 Посмотреть сообщение
Не подскажете как построить график?
plot(Y(:,1), Y(:,2))
1
klaxor2014
13 / 13 / 1
Регистрация: 13.02.2011
Сообщений: 101
14.11.2014, 13:55  [ТС] 5
Спасибо за подсказки, проблема решена, по крайней мере с кодом матлаб, проверил с помощью dsolve.Графики сходятся, буду решать проблему с решением аналогичной задачи в Simulink.
0
14.11.2014, 13:55
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
14.11.2014, 13:55

Построить график по решению системы диф уравнений
Есть решение системы диф уравнений: syms x(t) y(t) z(t) t Dx=diff(x);...

Решение системы уравнений (5 уравнений, 3 неизвестные)
мучаюсь и не могу решить систему: b0 = (A*( (A+1) +...

Решение системы уравнений
Добрый вечер! Помогите разобраться с решением системы двух уравнений с тремя...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
5
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru