No.1
- 回答日時:
No.2
- 回答日時:
コンパイラ・・・プログラムですか?
^が使えなくても、Cのpow関数みたいに、べき剰を計算してくれる関数はないんでしょうか?
pow(x, 0.5)みたいにかけるような・・・・
No.3
- 回答日時:
sight さんが仰るように何か関数があるような気はしますが、
どうしても近似が必要であれば次のテイラー展開の結果を使ってください。
(1+x)^α = 1 + (α,1)*x + (α,2)*x^2 + … + (α,n)*x^n + o(x^n)
ここで、(α,n) は2項係数でαは自然数でなくてもよいですが n は自然数です。
(通常カッコの中のαと n は縦に並べて書きますが今は横に並べています。)
つまり、
(α,n) = α(α-1)…(α-n+1)/n!
です。ただし、(α,0)=1 です。
今は、α=1/2 ですから
(1/2,1) = (1/2)/1! = 1/2
(1/2,2) = (1/2)(1/2-1)/2! = -1/8
(1/2,3) = (1/2)(1/2-1)(1/2-2)/3! = 1/16
などになります。
したがって、
(1+x)^0.5 = 1 + x/2 - x^2 /8 + x^3 /16 …
と展開できます。ご質問の x^0.5 であれば、x+1 ⇒ x とずらして
x^0.5 = 1 + (x-1)/2 - (x-1)^2 /8 + (x-1)^3 /16 …
ですね。
あとは近似の精度に気をつけて下さい。
No.4ベストアンサー
- 回答日時:
√x = x^0.5 は言語によってはsqrt(x), sqr(x), power(x,0.5) などと書かれます。
そういうのがない場合には、数値計算をする関数(functionあるいはprocedure, subroutine, method)を作ってやれば良いですね。 要するに、大抵の言語でコンパイラに作りつけの関数sqrtがやっていることを自分で書けばよい。
~~~~~~~~~~~~~~~~~~~~~
√xを高い精度で計算するには、漸化式を使うのが速いと思います。
もちろんx<0のときはエラー、x=0なら√x=0です。
x>0のとき、y=√x に近づく近似値の列y[0], y[1], ...., y[k], .... を次々と作り、誤差が小さくなったら終わります。
●漸化式
適当な出発値y[0]から始めて
y[n+1]=y[n]/2+x/(2y[n])
という漸化式でy[1], y[2], ... を作ります。(Newton法)
この式は、繰り返しの度に有効桁数が倍になる、「2次収束」という性質を持っているので、出発値y[0]が正解に近ければ、ほんの数回の繰り返しで十分な精度(10進法で10桁以上)に達します。
●繰り返しを何処で打ち切るかの判断
筋から言うと、正解√xに対する相対誤差 δ
δ=|(y[n]^2)/x - 1|/2
が十分小さくなったら打ち切る、ということで良いでしょう。
このδがなぜ相対誤差なのか。近似値y[n]が絶対誤差εを含んでいて
y[n]=sqrt(x)+ε
であるとすると、
δ=|(y[n]^2)/x - 1|/2
= |(x+2y[n]ε+(ε^2))/x - 1|/2
= |(2ε√x+ε^2)/x |/2
≒|(2ε√x)/x |/2 (|ε|は小さいのでε^2 は無視できる。)
= |(ε/√x |
だからです。
幾らなら十分小さくなったと言えるか。これは要求される精度に依存しますが、一般に計算機内部での浮動小数点数値の有効桁数を超えたら、それ以上計算したってしょうがないですね。
これを逆手に取って、
y[n]とy[n+1]が同じになったら打ち切る
というやり方もあります。
●出発値y[0]
出発値y[0]があまりにもでたらめだと、上記の漸化式が旨く働かないことがあります。
y[0]は、数値xの計算機内部での表現を利用すると簡単に決められます。
浮動小数点の数値xは符号部S、仮数部D、指数部Eで表され、たいてい16進法で
x = S×D×(16^E)
となっている。Eは整数です。
ここで、x≠0の場合には
1>D≧1/16
です。
この問題ではx≠0だし、符号部SはS=1に決まっています。従って
√x = √D ×(4^E)
ですね。
√D=√(1+(D-1))≒1+(D-1)/2 ≒1
と荒っぽい近似ができます。また
Eが偶数の場合には
√x = √D ×16^(E/2)
Eが奇数の場合には
√x = (√D )/4 ×16^((E+1)/2)
です。
だから
Eが偶数のとき y[0] = 符号部1 仮数部(1+(D-1)/2) 指数部(E/2)
Eが奇数のとき y[0] = 符号部1 仮数部(1+(D-1)/2)/4 指数部((E+1)/2)
とでもすれば良いでしょう。もっと手抜きして
Eが偶数のとき y[0] = 符号部1 仮数部1 指数部(E/2)
Eが奇数のとき y[0] = 符号部1 仮数部1 指数部((E+1)/2)
でも大丈夫です。(漸化式の繰り返しが少し増えるだけです。)
内部表現の調べかたが分からない場合には、計算でDとEを出すこともできます。
16^(E-1) < x < 16^(E)
となるEを見つければ良いのです。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
百分率の統計解析について
-
最小二乗法は、なぜ「二乗」な...
-
有効数字について教えて下さい...
-
信号対雑音比について
-
有効数字が整数部分の一桁で表...
-
測定誤差や有効数字は厳しく指...
-
データの相対誤差について
-
集積公差について教えて下さい。
-
誤差率 理論値が0の時
-
最小二乗法の重みづけについて...
-
20人対象のアンケートに統計...
-
プラスマイナス1.5の範囲に...
-
近似値、誤差の限界について
-
重回帰分析で手持ちデータがま...
-
有効数字について 6×3.0だった...
-
アップ率の求め方について
-
パーセントの合計と平均について
-
パーセントの平均の計算式。 42...
-
日本人男性の平均身長172cmと言...
-
小数点以下
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報