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θ") MAT P00=P99 !次へ MAT Q00=Q99 NEXT w END なお、sin、cosのn倍角の公式は加法定理を使って、 2変数(sinθ、cosθ)の多項式で導出できます。 |