Mình có làm 1 bài mô phỏng đạn bay bằng matlab tương tự như bài này
http://vatlyvietnam.org/forum/showthread.php?p=10251#post10251nhưng nhiều chỗ chưa hiểu lắm vì cách kí hiệu đạo lượng hơi khác mình học, mong các bạn chỉ giúp:
ta giới hạn xét ở các vận tốc nhỏ , khi đó thì lực ma sát của không khí tỷ lệ bậc nhất với vận tốc : Fms=βr′
thời điểm ban đầu :
r(0)=0
V(0)=(V0cosα)+(V0sinα)
phương trình chuyển động :
mr′′=−βr′−mg
⇒ V′+βV/m=−g
nhân hai vế phương trình trên với e^(βt/m) :
V′.e^(βt/m)+βmV.e^(βt/m)=−g.e^(βt/m)
bây giờ thì ta thấy việc nhân một lượng như vậy rất là có lợi khi tích phân hai vế :
V.e^(βt/m)=−∫g.e^(βt/m)=−gm/β.e^(βt/m)+C1
xác định C1 dựa vào điều kiện ban đầu
⇒ C1=V(0)+gm/β
vậy ta được biểu thức của V(t) :
V=−gm/β.(1−e^(βt/m))+V(0).e^(βt/m)
tích phân lên ,được biểu thức của r(t):
r=−gm/β.(t+m/β.e^(βt/m))−V(0).m/β.e^(βt/m)+C2
bây giờ lại xác định C2 từ điều kiện ban đầu , được kết quả cuối ( tự làm ) [/b]
1. Nhân hai vế phương trình trên với e^(β*t/m) :
thì e^(β*t/m) là j vậy ?
2. Biểu thức V=−gm/β.(1−e^(βt/m))+V(0).e^(βt/m) mình không hiểu lắm :|
mong các bạn giúp mình, bỏ 2 ngày ra nháp rồi mà vẫn chưa hiểu dc :|
