Вектор параметров для ode45, MATLAB |
Здравствуйте, гость ( Вход | Регистрация )
Вектор параметров для ode45, MATLAB |
![]()
Сообщение
#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.05.2025, 23:41 |