新しく発言する  EXIT  インデックスへ
SIN,COS倍角式

  SIN,COS倍角式 しばっち 2008/02/08 20:45:02 
  EXTERNALSUBSINEX(F,COSA(),SINA(,),COSB()... しばっち 2008/02/08 20:45:37 
  │└!チェビシェフの多項式Tn(x)からn倍角の公... 山中和義 2008/02/09 15:09:21 
  │ └山中和義さんご返事ありがとうございます。 しばっち 2008/02/10 17:34:35  (修正1回)
  !tanのn倍角の公式 山中和義 2008/02/10 17:34:41  (修正1回)
   └つづき 山中和義 2008/02/10 17:35:22  (修正1回)

Re: !tanのn倍角の公式  返事を書く  ノートメニュー
山中和義 <drdlxujciw> 2008/02/10 17:35:22 ** この記事は1回修正されてます
つづき



!2倍角の公式
!●tan2θ=(tanθ+tanθ)/(1-tanθtanθ)

!●tanθ=(tan0+tanθ)/(1-tan0tanθ)=t
!●tan2θ=(tanθ+tanθ)/(1-tanθtanθ)=2*t/(1-t^2)
!●tan3θ=(tan2θ+tanθ)/(1-tan2θtanθ)=(3*t-t^3)/(1-3*t^2)
!●tan4θ=(tan3θ+tanθ)/(1-tan3θtanθ)=(4*t-4*t^3)/(1-6*t^2+t^4)
!  :
!  :

DIM t(2*N)
LET t(1)=1 !tanθの定義


DIM P00(2*N),Q00(2*N) !1つ前の分子、分母 tan(w-1)θ=P/Q
LET P00(0)=0 !tan0
LET Q00(0)=1

DIM P99(2*N),Q99(2*N) !分子、分母
FOR w=1 TO N
PRINT "tan";w;"θ="

!分子の部
DIM T1(2*N)
CALL poly_mul(N,Q00,t, T1) !tan(w-1)θ+tanθ=P/Q+t=(P+Qt)/Q
CALL poly_add(N,T1,P00, P99) !通分して分子のみ ※分母は同じため約せる
CALL disp_poly(n,P99,"tanθ")

PRINT REPEAT$("-",w*7)

!分母の部
CALL poly_mul(N,P00,t, T1) !1-tan(w-1)θtanθ=1-P/Q*t=(Q-P*t)/Q
CALL poly_sub(N,Q00,T1, Q99)
CALL disp_poly(n,Q99,"tanθ")
PRINT


MAT P00=P99 !次へ
MAT Q00=Q99
NEXT w


END



なお、sin、cosのn倍角の公式は加法定理を使って、
2変数(sinθ、cosθ)の多項式で導出できます。


 インデックスへ  EXIT
新規発言を反映させるにはブラウザの更新ボタンを押してください。