Solving a Differential equation plus plotting results: Fluid Mechanics
3 views (last 30 days)
Show older comments
I am trying to create a code to: 1) Solve this differential Equation for h(t) (or h2) in terms of t. 2) Solve for t at a certain h2. 3) Plot h2 versus t
the equation is (dh/dt)=(-36290.323)*[(((2*mw*g*t)/(Gama*Atank))+(2*g*h))^(1/2)]
Thank you for your help.
0 Comments
Answers (1)
Youssef Khmou
on 12 Apr 2013
hi Kyle,
You can try as the following :
1) Create a function in M-file inwhich you define the equation :
function dh=MYFunc(t,h)
%(dh/dt)=(-36290.323)*[(((2*mw*g*t)/(Gama*Atank))+(2*g*h))^(1/2)]
A=-36290.323;
mw=2.33;
g=9.81;
Gamma=1.23;
Atank=4;
dh(1)=A*(((2*mw*g*t)/(Gamma*Atank))+(2*g*h(1))).^(1/2);
2) In the workspace you call the function "ode23" or "ode45" as :
t0=0; % staring time
tf=10; % t final
Initial_value=2; % the initial value for h
[t,h]=ode45('MYFunc',t0,tf,[Initial_value]);
figure, plot(t,abs(h));
Adjust the parameters, and try ...because in this case h is complex .
2 Comments
Youssef Khmou
on 12 Apr 2013
hi Kyle,
no i can not tell, that may be possible if we already know the Sampling frequency of t .
See Also
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!