No.2
- 回答日時:
この方程式は解析的には解けません。
この方程式の場合図的に解く方法が直感的で分かり易いですが、
プログラムに使うとなると図的解法は使えないですね。
数値解法を使うことになります。
プログラムに使うコンピュータがどんなものか分かりませんが、
大型機なら数値解法のモジュールが入っているはずです。
大して難しくありませんから、自分で組むことも十分可能です。
数値解法にはいくつか方法がありますので、問題に応じて適切なものを選択するとよいと思います。
No.3ベストアンサー
- 回答日時:
とてもスマートには行かないと思います。
本気で取り組むと、これだけで結構な仕事になっちゃう。この方程式をいろいろな係数について何度も解くのかどうか、係数の精度がどの位か、値の範囲がどうか、などによってかなり事情が変わりますが、とりあえず一般論。
A=0,w=0,B=0のどれかの場合は簡単ですから、これらの場合は除外します。
T=wt, b=-B/(wA), c=-C/Aとおいて
sin(T) + b T + c = 0
|b|≧1なら解は1個だけですが、|b|<1だと|b|が小さくなるほど解の数が多くなる。全部の解を求めるということだとなかなか大変です。
まず、|bT+c|≦1の範囲だけ調べれば良いことは自明です。(特定の範囲の解を求めるなら、そこだけ検討すればよい。)nπが求めたい解Tに近いような整数nを選んで、
nが偶数→T=x+nπ
nが奇数→T=nπ-x
と変数変換して代入し、
f(x)=sin(x)+p x + q
f'(x)=cos(x)+p
と表します。すると、x∈-π/2~π/2の間でf(x)=0の解を探せばよい。この範囲には最大3個の解があり得ます。そこで、近似的な出発値x[0]を見つけてこれをNewton法
x[n+1]=x[n]-f(x[n])/f'(x[n])
で改良するのが簡単です。ただし重解(y=0がy=f(x)の接線になっている)に近い場合は数値的に不安定になります。これをきちんと処理するのは難しい。(Tをπで割った余りを精密に計算する必要があり、数値の有効桁数が不足する現象(桁落ち)が起こるため、本当にy=f(x)がy=0に交差しているのか、接しているのか、微妙に離れているのか区別できなくなる恐れがあるからです。)
出発値x[0]を得るには、例えば適当な次数のマクローリン展開でも結構ですが、何度も解くのなら単にsineの数表を作っておいてlook upするのが簡単です。
この回答へのお礼
お礼日時:2001/02/23 10:28
丁寧な回答ありがとうございます.
プログラムで解く以外にはないのですね.
御指摘の通りの解の存在する|bT+c|≦1の範囲で,
sin関数と直線の交点はsin関数の1/4周期に最大1つあるようですので
これをもとに二分探索で探すことにいたします.
ありがとうございました.
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
Excelで合計値を基にデータを均...
-
a(x-α)(x-β)
-
16の4乗根は±2ではない!?
-
3次関数が極値をもつ必要十分条件
-
3次関数と直線が接する場合、...
-
tanX=Xの解
-
点P(x+y、xy)の軌跡を求めよ。...
-
2次不等式の解の答え方について
-
2次方程式X^2-3X-1=0の2つの...
-
解なし≠解はない
-
3x-7y=1の全ての整数解を求め...
-
答えを教えて
-
連立一次方程式における正のみの解
-
「この解は問題にあう」中学2...
-
三角関数
-
高校数学です。 m^2-11m-1が整...
-
数学の未解決問題が解けたら
-
つぎのステートメントはあって...
-
微分方程式の問題がわかりません
-
複数の品目での単価と全体の合...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
Excelで合計値を基にデータを均...
-
適正解と最適解
-
微分の重解条件は公式として使...
-
解なし≠解はない
-
微分方程式についてです。 xy'-...
-
tanX=Xの解
-
数学の質問です。 2つの2次方程...
-
微分方程式 定常解について・・・
-
3次関数と直線が接する場合、...
-
16の4乗根は±2ではない!?
-
複数の品目での単価と全体の合...
-
aの値に関係なくとよく問題で見...
-
解に3つ以上±や∓がある時複号...
-
必要条件、、十分条件、必要十...
-
数学についてです 「 aを定数と...
-
微分方程式で、分母=0の場合は...
-
一枚の板から何枚取れるか?
-
3次関数と1次関数が接するとき
-
何故グラフに接するとき重解に...
-
答えを教えて
おすすめ情報