N88-BASICで振り子 (3回目)

2024/1/9(火)
N88-BASICで振り子 (3回目)
 
二重振り子(Double pendulum)
 
■ 前提
https://ulprojectmail.blogspot.com/2024/01/pendulum-5.html
振り子 (5回目)
より
 
▼ 定義
ℓ_i :紐の長さ(m)
m_i :質点の質量(kg)
θ_i:鉛直下方向からの角度(rad)
g  :重力加速度
x_i :支点からの右方向の変位(m)
y_i :支点からの下方向の変位(m)
 
▼ 微分方程式
{m₁ + m₂ - m₂cos²(θ₁ - θ₂)}ℓ ₁θ(･･)₁ 
= {ℓ₁θ(・)₁²cos(θ₁ - θ₂) - ℓ₂θ(・)₂²}m₂sin(θ₁ - θ₂)
+ {m₂sinθ₂cos(θ₁ - θ₂) - (m₁ + m₂)sinθ₁}g
 
{m₁ + m₂ - m₂cos²(θ₁ - θ₂)}ℓ₂θ(･･)₂ 
= {m₂ℓ₂θ(・)₂²cos(θ₁ - θ₂) + (m₁ + m₂)ℓ₁θ(・)₁²}sin(θ₁ - θ₂)
+ {sinθ₁cos(θ₁ - θ₂) - sinθ₂}(m₁ + m₂)g
 
■ 動作
数値計算の2重振り子を表示します
(積分θ(・)=∫θ(･･)dt等は長方形近似で計算しています)
 
■ PC-8801
1070 VSYNC = &HA0 '--- &HA0:PC-98 , &H40:PC-88
を
1070 VSYNC = &H40 '--- &HA0:PC-98 , &H40:PC-88
に変更するとPC-8801用になります
 
 
VL,NL,XL-BASICとblg～.zip(pend003.bas)は
以下のリンクからダウンロードできます
https://ulprojectmail.blogspot.com
Readme.txtを読んで遊んで下さい