投稿

N88-BASICで懸垂線 (1回目)

イメージ
2024/7/3(水) N88-BASICで懸垂線 (1回目)   (catenary) 懸垂線(カテナリー、紐を垂らしたときの曲線)   ■ 前提 ▼ 参照 https://ulprojectmail.blogspot.com/2024/07/catenary-1.html 懸垂線 (1回目) より   懸垂線はy = cosh(x)のグラフになる   ■ 解説 y = sinh(x) y = cosh(x) y = tanh(x) のグラフを描画しました     VL,NL,XL-BASICとdlg~.zip(cate00 1 .bas)は このブログ(以下のリンク)から ダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい  

懸垂線 (1回目)

イメージ
2024/ 7 /1( 月 ) 懸垂線 (1回目)   ( c atenary) 懸垂線 ( カテナリー、紐を垂らしたときの曲線 )   ■ 導出 ▼   定義   T : 紐の張力 θ:水平方向からの張力方向への角度 ds : x~x +dx 間 (点 A,B 間 )の紐の長さ ρ:紐の線密度 g : 重力加速度 f (x): 紐の高さ   ▼   微分方程式の導出 y  = y (x)   …   紐の形のグラフ y '( x) = dy/dx = tan θ (x)   …   ①   点 A : θ( x), T( x) 点 B : θ( x+dx) ,   T( x+dx) 点 A B 間の紐の長さ ds ds  = √ (dx 2  + dy 2 )   = dx √ {1 + (dy/dx) 2 }   水平方向 T (x)cos θ (x) = T(x+dx)cos θ( x+dx) どの xでも同じ値なので T (x)cos θ (x) = const. H  = T (x)cos θ (x)   …   ② と置く(水平張力)   垂直方向 d y/dx = y '( x) T (x)sin θ (x) +   dsρ g = T(x+dx)sin θ( x+dx) ρ gdx √ {1 + (dy/dx) 2 } = T(x+dx)sin θ( x+dx) - T (x)sin θ (x) {T(x+dx)sin θ( x+dx)- T (x)sin θ (x)}/dx = ρ g √ {1 + (dy/dx) 2 }   (d/dx){ T (x)sin θ (x)} = ρ g √ {1 + (dy/dx) 2 }   …   ③   ①,②より T (x)sin θ (x) = T(x)cos θ (x)sin θ (x)/cos θ (x) = T(x)cos θ (x)tan θ (x) = T(x)cos θ (x)(dy/dx) =   H (dy/dx) を ③に代入 (d/dx){ H (dy/dx)} = ρ g √ {1 + (dy/dx) 2 }   (d/dx) 2 y = ( ρ g/ H ) √ {1 + (dy/dx) 2 }   λ  = ρ g/ H , f(x) = dy/dx と置く (

三角関数 (5回目)

2024/6/28(金) 三角関数 (5回目)   ■ 定義 ▼ 双曲線関数 y = sinh(t) = {exp(t)-exp(-t)}/2 y = cosh(t) = {exp(t)+exp(-t)}/2   ▼ 備考 tanh -1 (x)については(3回目)を参照   ■ 結果 ▼ 逆双曲線関数 t = sinh -1 (y) =  log(y + √{y 2  + 1}) t = cosh -1 (y) = ±log(y + √{y 2  + 1})   ▼ 双曲線関数(加法定理) sinh(a±b) = sinh(a)cosh(b)±sinh(b)cosh(a) cosh(a±b) = cosh(a)cosh(b)±sinh(a)sinh(b) sinh(2a) = 2sinh(a)cosh(a) cosh(2a) = cosh 2 (a) + sinh 2 (a)   ▼ 双曲線関数(その他) cosh(a)+1 = 2cosh 2 (a/2) sinh(a) = 2sinh(a/2)cosh(a/2) {cosh(a)+1}/sinh(a) = tanh(a/2)     ■ 導出 ▼ 逆双曲線関数 y = sinh(t)  , y = cosh(t) t = sinh -1 (y) , t = cosh -1 (y) をlogで表す (log = log e  とする)   y = exp(x) log(y) = x   y = sinh(t) = {exp(t)-exp(-t)}/2 2y = exp(t)-exp(-t) 2yexp(t) = exp(2t) - 1 exp(2t) - 2yexp(t) - 1 = 0 exp(t) = y±√{y 2  + 1} ≧ 0 = y + √{y 2  + 1} ≧ 0 t = log(y + √{y 2  + 1})   y = sinh(t) t = sinh -1 (y) = log(y + √{y 2  + 1})     y = cosh(t) = {exp(t)+exp(-t)}/2 2y = exp(t)+exp(-t) 2yexp(t) = exp(2t) + 1 exp(2t) - 2yexp(t) + 1 = 0 exp(t) = y±√{y 2  - 1} t = log(y±

N88-BASIC,Pythonで円周率 (3回目)

イメージ
2024/6/24(月) N88-BASIC,Pythonで円周率 (3回目)   ■ N88-BASIC https://ulprojectmail.blogspot.com/2022/07/n88-basicpi-3.html N88-BASICで円周率 (3回目)   NL-BASIC(N88-BASIC互換?)で走らせ 1040桁まで正しく表示されている事を確認しました それ以降の桁の確認はしていません     ■ Python 桁数(5の倍数)を入力し円周率を表示する N88-BASICをPythonに書き換えました   表示は桁数が分かり易いように 3.14…を 314…と表示しています   Pythonの整数型は無限桁なので もっと効率の良い書き方ができますが あえてほぼそのまま書き換えています   式やプログラムのミスなどにより 正しく計算出来ていない可能性はありますので 使用には十分注意して下さい     VL,NL,XL-BASICと d lg~.zip ( pi 00 3 . py )は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい  

N88-BASIC,Pythonで円周率 (2回目)

2024/6/16(日) N88-BASIC,Pythonで円周率 (2回目)   ■ N88-BASIC https://ulprojectmail.blogspot.com/2022/07/n88-basicpi-2.html N88-BASICで円周率 (2回目)   円に内接(外接)する正多角形を利用して 円周率の範囲を求めます     ■ Python nを入力し 円に内接(外接)する正n角形を利用して 円周率の範囲を表示する N88-BASICをPythonに書き換えました     VL,NL,XL-BASICと d lg~.zip ( pi 00 2 . py )は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい  

N88-BASIC,Pythonで円周率 (1回目)

イメージ
2024/6/12(水)   N88-BASIC,Pythonで円周率 (1回目)   ■ N88-BASIC https://ulprojectmail.blogspot.com/2022/04/n88-basicpi.html N88-BASICで円周率     ■ Python 乱数を使用して円周率を求める N88-BASICをPythonに書き換えました     VL,NL,XL-BASICと d lg~.zip ( pi 001. py )は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい  

量子力学2 (2回目)

2024/6/9(日) 量子力学2 (2回目)   (Quantum mechanics)   ディラック方程式     ■ 前提 https://ulprojectmail.blogspot.com/2023/05/quantum-1.html 量子力学 (1回目) https://ulprojectmail.blogspot.com/2023/05/relativity-2.html 相対性理論 (2回目) https://ulprojectmail.blogspot.com/2023/09/electro-2.html 電磁気学 (2回目) https://ulprojectmail.blogspot.com/2023/09/electro-4.html 電磁気学 (4回目)   ▼ 定義 A :振幅(m) λ:波長(m)(m/回) k :角波数(rad/m) [k = 2π/λ] ω:各振動数(rad/s) [ω = 2πf] t :時間(s) ν:物質波などの振動数(Hz) c :光速(m/s) x :位置ベクトル(m) m :質量(kg) E :エネルギー(J) p :運動量(kg・m/s) [p = mv] h :プランク定数(6.62607015×10 -34 J・s) ℏ :ディラック定数 [ℏ = h / (2π)]   ナブラ∇ = (∂/∂ x ) = (∂/∂x,∂/∂y,∂/∂z) grad f = ∇f , div E  = ∇・ E  , rot E  = ∇× E   ラプラシアン∇・∇ = ∇ 2  = Δ ダランべルシアン □ = ∂ μ ∂ μ  = -(1/c 2 )(∂ 2 /∂t 2 ) + Δ   素粒子 p = h/λ = ℏk k = p/ℏ p = E/c、λ=c/ν E = pc = pνλ = hν = ℏω ω = E/ℏ   ▼ クライン・ゴルドン方程式 E = iℏ(∂/∂t) p  = -iℏ(∂/∂ x ) E 2  = (mc 2 ) 2  + | p c| 2     {iℏ(∂/∂t)} 2  = (mc 2 ) 2  + {-iℏc(∂/∂ x )} 2   -ℏ 2 (∂ 2 /∂t 2 ) = -ℏ 2 c 2 (∂ 2 /∂ x 2 ) + m 2 c 4     -ℏ 2 (∂ 2