IPB                

Здравствуйте, гость ( Вход | Регистрация )


ФорУм - для ума ©
БСЭ; DJVU Библиотека - Основное книгохранилище
 
Ответить в данную темуНачать новую тему
Вектор параметров для ode45, MATLAB
korson
сообщение 26.03.2011, 14:45
Сообщение #1


Бакалавр
***

Группа: Member
Сообщений: 101
Регистрация: 16.8.2009
Из: Москва
Пользователь №: 36894
Поблагодарили: 28 раз(а)




Друзья, кто нибудь передавал ode45 вектор параметров вот так:
Код
----------------------------------
clear all
close all
clc

param = [1 2];
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options,param);

% Plotting the columns of the returned array Y versus T shows the solution

figure('color','w')
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')

----------------------------------

function dy = rigid(t, y, param)

% dy = zeros(3,1);    % a column vector

dy1 = param' * y(2) * y(3);
dy2 = -y(1) .* y(3);
dy3 = -0.51 * y(1) .* y(2);

dy = [dy1'; dy2'; dy3'];




Мне нужно получить решения для вектора параметров [1 2].
Перейти в начало страницы
+Цитировать сообщение

Ответить в данную темуНачать новую тему
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 



Текстовая версия Сейчас: 3.07.2020, 1:51


Rambler's Top100