投稿

NEWS

2025/4/8(火 ) 波の合成 (1回目)  を公開しました    2025/4/4(金) VL,NL,XL-BASIC ver~28x1 を公開しました ダウンロードファイルはウイルスチェックしていますが ウイルスチェックソフトに引っかかる場合 誤認識かどうかは自己責任でご判断して下さい    2025/3/24(月 ) N88-BASICでフーリエ変換 (2回目)(離散フーリエ変換) を公開しました

波の合成 (1回目)

202 5 / 4 /8 (火 ) 波の合成   ( 1 回目 )   ( wave)   ■   2波の合成 ▼ 振幅が同じ場合     sin(α+β) = sinαcosβ + cosαsinβ -) sin(α-β) = sinαcosβ - cosαsinβ     sin(α+β) + sin(α-β) = 2sinαcosβ    Asin(α+β) + Asin(α-β) = 2Asinαcosβ … ① a = α+β , b = α-β と置くと α = ( a+b )/2   , β = ( a-b )/2   を ①式に代 入 Asina + Asinb = 2Acos ( (a-b) /2) sin( ( a+b) /2) 包絡線 2Acos ( (a-b) /2)   ▼ 振動数が同じ場合 底辺 X、高さY、斜辺R、角θ ( XR間 ) の 直角三角形 θ = Tan -1 (Y/X)   , R = √(X 2 +Y 2 )   , X = Rcosθ   , Y = Rsinθ  Rcosθsina   +   Rsinθcosa = Rsin(a+θ)   加法定理 に代入 Xsina + Ycosa = Rsin(a+θ) Xsina + Ycosa   = √(X 2 +Y 2 )sin{a + Tan -1 (Y/X)}   ▼ 振幅 も 振動数 も 異なる 場合 を含む Asin(α+β) = Asinαcosβ+Acosαsinβ   加法定理 Bsin(α-β) = Bsinαcosβ-Bcosαsinβ   加法定理 を加えて X=(A+B)cosβ、Y= ( A-B ) sinβと置くと Asin(α+β) + Bsin(α-β)   =   (A+B)sinαcosβ+(A-B)cosαsinβ = Xsinα + Ycosα   = √(X 2 +Y 2 )sin{α + Tan -1 (Y/X)} X 2 +Y 2  ...

四平方の定理 (2回目)

イメージ
2025/3/31(月) 四平方の定理 (2回目)   (De Gua's theorem)   ■ 任意の形の四平方の定理(デカルト・グアの定理、ド・グアの定理) A 2   + B 2   + C 2   = D 2   この式は任意の形の DをA,B,C面に投影しても成り立つ   図の三角錐の面の直角三角形の面積をそれぞれ A,B,C、 手前の面の三角形の面積を Dとする (x, y, zは座標軸とする)   ■ 任意の形の四平方の定理の導出 A,B,C,Dの法線ベクトル(単位ベクトル)をそれぞれ a,b,c,dとするとa,b,cは a = (1, 0, 0), b = (0, 1, 0), c = (0, 0, 1) となり dは d = (x, y, z)  … (x 2   + y 2   + x 2   = d・d = |d| 2   = 1) と置く   x軸の真上からDを見てAの面に投影すると 直角三角形 Aに見えるので A = D(a・d)となる … (a・dは内積でa//dなら1, a⊥dなら0) B,Cも同様にy軸,z軸の真上から見てDを投影すると B = D(b・d) C = D(c・d) (a・d) = (1, 0, 0)・(x, y, z) = 1・x + 0・y + 0・z = x (b・d) = (0, 1, 0)・(x, y, z) = 0・x + 1・y + 0・z = y (c・d) = (0, 0, 1)・(x, y, z) = 0・x + 0・y + 1・z = z より A = Dx, B = Dy, C = Dx A 2   + B 2   + C 2   = D 2 x 2   + D 2 y 2   + D 2 x 2   = D 2 (x 2   + y 2   + x 2 ) = D 2    … (上記よりx 2   + y 2   + x 2   = 1)   A 2   + B 2   + C 2  ...

四平方の定理 (1回目)

イメージ
2025/3/28(金) 四平方の定理 (1回目)   (De Gua's theorem)   ■ 四平方の定理(デカルト・グアの定理、ド・グアの定理) A 2   + B 2   + C 2   = D 2   A = ab/2  … △abdの面積 B = bc/2  … △bceの面積 C = ca/2  … △cafの面積 D         … △defの面積 とする   ■ 四平方の定理の導出 d 2   = a 2   + b 2   e 2   = b 2   + c 2   f 2   = c 2   + a 2   θ=∠de  … 二辺間の角度とする 余弦定理より f 2   = d 2   + e 2   - 2decosθ cosθ = (d 2   + e 2   - f 2 )/(2de) sinθ = √{1 - (d 2   + e 2   - f 2 ) 2 /(2de) 2 } = √{(2de) 2   - (d 2   + e 2   - f 2 ) 2 }/(2de)   D = (1/2)desinθ = (1/2)de√{(2de) 2   - (d 2   + e 2   - f 2 ) 2 }/(2de) = √{(2de) 2   - (d 2   + e 2   - f 2 ) 2 }/4   d 2   + e 2   - f 2   = a 2   + b 2   + b 2   + c 2   - c 2   - a 2   = 2b 2   d 2 e 2   = (a 2   + b 2 )(b 2   + c 2 ) ...

N88-BASICでフーリエ変換 (2回目)

イメージ
2025/3/24(月) N88-BASICで フーリエ変換 ( 2 回目 ) (fourier)   ■ 離散フーリエ変換(DFT)( Discrete Fourier Transform ) ▼ 前提 https://ulprojectmail.blogspot.com/2025/03/fourier-6.html フーリエ変換 (6回目) より   T:サンプル計測時間 N:サンプル総数 t k :サンプルの計測時刻( k   = 0, 1, 2, …, N-1) , t k   =   k( T/N ) y k :各サンプルの計測値( k   = 0, 1, 2, …, N-1)   周波数 f n   = n /T ( n   = 1 , …, N) y(t k ) = 2 Σ n = 0 N {A ( f n )exp(i 2πf n t k ) } A(f n ) = (2 /N ) Σ k = 0 N-1 { y k e xp(-i 2πf n t k ) } A = |A(f n )|  … |exp(-ix)| =  √{cos 2 x + sin 2 x)}   ▼ 動作 y(t)をcos,sinの合成波として定義し A = |A(f k )|   を求めて 周波数 f n   と、振幅 A を表示しています   VL,NL,XL-BASICとdlg~.zip( fou 00 2 .bas)は このブログ (以下のリンク)から ダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい  

フーリエ変換 (6回目)

2025/3/14(金) フーリエ変換 ( 6 回目 )   (fourier   t ransform )   ■ 離散フーリエ変換(DFT)( Discrete Fourier Transform ) ▼ 結果 T:サンプル計測時間 N:サンプル総数 t k :サンプルの計測時刻( k   = 0, 1, 2, …, N-1) , t k   =   k( T/N ) y k :各サンプルの計測値( k   = 0, 1, 2, …, N-1)   周波数 f n   = n /T ( n   = 1 , …, N)   y(t k ) = 2 Σ n = 0 N {A ( f n )exp(i 2πf n t k ) } A(f n ) = (2 /N ) Σ k = 0 N-1 { y k e xp(-i 2πf n t k ) }   A = |A(f n )|  … |exp(-ix)| =  √{cos 2 x + sin 2 x)}   ▼ 式 フーリエ変換 (FT) y:変位, t:時間, A:振幅, f:周波数 y(t) = ∫ -∞ ∞ {A ( f )e i2πft } d f A(f) = ∫ -∞ ∞ {y(t)e -i2πft } dt   T:サンプル計測時間 N:サンプル総数 t k :各サンプルの計測時刻( k   = 0, 1, 2, …, N-1) y k :各サンプルの計測値 ( k   = 0, 1, 2, …, N-1)   y(t k ) = y k   , Δt = T/N   , t k   = k Δt = k( T/N )   Δ f = 1/T 周波数 f n   = n Δ f = n/T ( n   = 1, …, N)   A(f) = ∫ -∞ ∞ {y(t)e -i2πft } dt = 2 ∫ 0 ∞ {y(t)e -i2πft } dt   + 2・0  … y(t)偶関数+奇関数 A(f n ) = 2 Σ k = 0...