Диф. уравнение MatLab

@vikuk · Регистрация: 01.03.2011

Author24 — интернет-сервис помощи студентам

Нужно написать м-функцию для решения примерно такого уравнения
du/dt=(a/b)*((u+b+x)^2-(u-b+x)^2)) (1)

//я только начинаю работать с матлабом поетому извените за простые вопросы, хочу использовать матлабовские dll-ки.

в это уравнение передаются значения a,b,x.

я пробовал следующи вариант

Matlab M

function [t,u] ode45_1(f,tspan(масив 2х елементов Т1,Т2), u0,x,a,b)
[t,u]=ode45(@f,tspan,u0,x,a,b)
end
 
function f (u,x,a,b)
f=(a/b)*((u+b+x)^2-(u-b+x)^2)) 
end

ну это какой то схематический вариант вышел, надеюсь вы поможете.
спасибо

@Галина Борисовн · 01.03.2011, 17:50

Matlab M

clear
clc
T0=0;
global a b x
a=0.2;b=3;x=4;
interval=[0:0.05:1];
[t,u]=ode45(@FF1,interval,T0);
plot(t,u,' r  ');
legend('u(t)');
grid on
 
function F=FF1(t,u)
global a b x 
F=(a./b).*((u+b+x).^2-(u-b+x).^2);

Оба файла сохранить в одной папке

@vikuk · 01.03.2011, 19:37 **[ТС]**

спасибо,
мне нужно сделать в одном м-файле, в который я передаю T0, a, b, x, interval

Matlab M

function DU (a,b,x,interval,T0)
[t,u]=ode45(@FF1,interval,T0)
end
function F=FF1(t,u)
global a b x 
F=(a./b).*((u+b+x).^2-(u-b+x).^2);
end

вот так как нибуть можна сделать???

Добавлено через 6 минут
или так?

Matlab M

function [u,t]=DU(a,b,x,interval,T0)
FF1=(a./b).*((u+b+x).^2-(u-b+x).^2);
[t,u]=ode45(@FF1,interval,T0)
end

@Галина Борисовн · 01.03.2011, 20:24

В одном файле подпрограмма-функция и вызывающая программа не получится.
Можно решить в символьном виде так:
dsolve('Du=(a/b)*((u+b+x)^2-(u-b+x)^2)')
ans =
-x+exp(4*a*t)*C1

@vikuk · 01.03.2011, 21:09 **[ТС]**

что значит
ans =
-x+exp(4*a*t)*C1
???

в символьном виде так

Matlab M

1
2
3

function [u,t]=DU(a,b,x,interval,T0)
[t,u]=ode45(@'Du=(a/b)*((u+b+x)^2-(u-b+x)^2)',interval,T0)
end

???

@Галина Борисовн · 01.03.2011, 21:15

Сообщение от vikuk

что значит
ans =
-x+exp(4*a*t)*C1

Это результат интегрирования диф.уравнения u(t)=-x+exp(4*a*t)*C1, C1-константа, значение которой определяется из начального условия.

@vikuk · 02.03.2011, 12:26 **[ТС]**

есть такая м-функция, для решения диф. ур. 1 порядка:

Matlab M

function [t,u] ode45_1(strfunk, tspan, u0)
[t,u]=ode45(@odefunc, tspan,u0,[],strfunk)
end
 
function dudt (t,u,strfunc)
strfunction=strcat(strfunk,'+0*t+*u')
F=inline(strfunction);
dudt=feval(F,t,u)
end

-----------
мне нужно передать в функцию еще значения трех параметров.
есть вариант при формирувание strfunk конвертировать флоат то стрінг, таким образом формируя то что надо. но я думал может есть какой нибуть прямой метод.
вообще то я хочу сделать матлабовскую dll, которую подключаю в VS2008

@Галина Борисовн 2834 / 2131 / 86 Регистрация: 02.05.2010 Сообщений: 3,195
	01.03.2011, 20:24	4
	В одном файле подпрограмма-функция и вызывающая программа не получится. Можно решить в символьном виде так: dsolve('Du=(a/b)((u+b+x)^2-(u-b+x)^2)') ans = -x+exp(4at)C1 0

@Галина Борисовн 2834 / 2131 / 86 Регистрация: 02.05.2010 Сообщений: 3,195
	01.03.2011, 21:15	6
	Сообщение от vikuk что значит ans = -x+exp(4at)C1 Это результат интегрирования диф.уравнения u(t)=-x+exp(4at)C1, C1-константа, значение которой определяется из начального условия. 0

	02.03.2011, 12:26