non linear differential equation

3 ビュー (過去 30 日間)
mohan
mohan 2012 年 11 月 30 日
please help me how to solve
dq/dt + q = a1 + b1*u(t) + 1/u(t)^2

回答 (1 件)

Azzi Abdelmalek
Azzi Abdelmalek 2012 年 11 月 30 日
編集済み: Azzi Abdelmalek 2012 年 11 月 30 日
function dq=myeq(t,q,t1,u1,a1,b1)
u=interp1(t1,u1,t)
dq=-q+a1+b1+1/u^2;
%Then call myeq
a1=1;
b1=2;
q0=0;
t1=1:0.1:10; % vector time
u1=t1.^2 % your vector u
[t q]=ode45(@(t,q) myeq(t,q,t1,u1,a1,b1),t1,q0)
plot(t,q)

カテゴリ

Help Center および File ExchangeProgramming についてさらに検索

製品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by