【无标题】

 

%1lbf=4.45N
%ft是长度单位=30.48cm
%gal是体积单位
%圆桶在很短的时间内到达了终端速度，忽略不计加速路程
clear clc;
G=527.436;
Flo=470.327;
vg=0.08;
[e,a]=V(G,Flo,vg)

function [v,S] = V(G, Flo,vg)
%UNTITLED2 此处显示有关此函数的摘要
%   此处显示详细说明
v=0;
t=0;
Zu=0;
v2=(G-Flo)/vg
S=0;
t1=[0];
s1=[0];
v1=[0];
while  S <600 || v<v2
a=(G-Flo-Zu);
Zu=v*vg;
t=t+0.001;
S=S+double(a,v,v+a*t);
v=v+at;
t1=[t1 t];
s1=[s1 S];
v1=[v1 v];
end
plot(t1,s1,'b+',t1,v1,'r*')
end