- ベストアンサー
高精度乱数関数
Fookyの回答
- Fooky
- ベストアンサー率71% (59/82)
> 0.<1回目のrand()のbitパターン><2回目のrand()のbitパターン> > としたいという意味です。（よね？） その解釈は成り立たないでしょう。 その理由は、 １．除算前にrand()をlong doubleにキャストしている したがって、除算は４倍精度の浮動小数点形式で行われます。 ２の乗数による除算が単純な右シフトで表現されるデータ形式は、 符合無し整数型だけです。 また、もし整数での除算であったとしても、 ２．RAND_MAX+1での除算に相当する右シフトの結果はかならず０になる 非除数は最大でもRAND_MAXなので、８ビット目より上位にビットが 立つことはないんですから > 0.<1回目のrand()のbitパターン><2回目のrand()のbitパターン> このような意図であれば、bitパターンの生成は整数型で行ない、 また、除算ではなく、１回目のrand()にRAND_MAX+1を掛けて 左シフトし、2回目を足し合わせる。ということをした後に、 RAND_MAXの２乗で割る、という手順を踏まないと行けないでしょう。 long double x; int i; unsigned long k; for( i = 0, k = 0 ; i < 4 ; i++ ){ k = k << 15; k += (unsigned long)random() & 0x7fff; } x = (long double)k; for( i = 0 ; i < 4 ; i++ ) x /= (long double)(RAND_MAX + 1); random()の出力が７ビット以上、longが８バイト と仮定しています。
関連するQ&A
- 標準正規分布の乱数
RAND()関数は ((double)rand() / (1.0 + RAND_MAX))と定義します。 中心極限定理により、一様乱数を足し合わせると正規分布に近づくことから、 x = 分散 * (Σ[1～12]RAND() - 6) + 平均 で正規乱数が作れる。標準正規分布は分散1、平均0なのでその乱数は x = Σ[1～12]RAND() - 6 ですよね。この乱数を例えば100個羅列するにはどうしたらいいのでしょうか？ もし間違ってたら指摘してください。 参考文献「Cによるシミュレーションプログラム 石川宏」 #include <stdio.h> #include <stdlib.h> #define RAND() ((double)rand() / (1.0 + RAND_MAX)) #define NUMBER 10000 /* 発生させる乱数の数 */ main(void) { int j; double u, x; srand(5); for (j = 0; j <= 11; j++) { u = u + RAND(); } x = u - 6.0; }
- ベストアンサー
- C・C++・C#
- 大きな数の乱数を作るには
C 初心者です。 表題のように、unsigned longのスケールの乱数をつくりたいんですが、以下のように記述すると値がいつも同じになります。この理由と、正しく動作するにはどう直したらいいのか教えてください。 unsigned long ul; ul = 4294967295UL * rand() / (RAND_MAX + 1); 値は常に131071でした。
- ベストアンサー
- C・C++・C#
- 乱数について
乱数の分布を見るために以下のようなプログラムを書きました。 #include <stdio.h> #include <stdlib.h> #include <math.h> int main() { int i,imax, S[RAND_MAX], r; double x,y; FILE *output1; output1=fopen("random2.data","w"); imax=100000; for(i=0;i<=imax;i++){ r = rand(); S[r] += 1; } for(i=0;i<=RAND_MAX;i++){ fprintf(output1,"%d %d \n",i,S[i]); } return 0; } するとコンパイルできて実行もできるのですが、なぜか乱数が30000を 超えるくらいのところでおかしな値になりました。 原因がわからないのでどなたか教えてください。
- ベストアンサー
- C・C++・C#
- 乱数について(Ｃ言語)
C言語において，乱数の範囲を 0 ≦ r < 1 とする場合には double r=(double)rand()/(RAND_MAX+1); とするのは知っているのですが0 < r ≦ 1にする場合の方法がわからず困っています． アドバイスいただきたいです．
- ベストアンサー
- その他（インターネット・Webサービス）
- 乱数のdouble型について
JSPから下記のクラスファイルを呼び出し、戻り値を返すように作りたいのですがうまくいきません。 1.Math.floor(Math.random()*100)-50 上記で実行しても小数点以下がでてしまいます。 出ないようにするのは無理なのでしょうか？ 2.乱数はdouble型以外だめなのでしょうか？ コンパイルするとdouble型なので間接参照できません。というようなエラーが出てしまいます。 これは結果を文字型に変更して戻り値として返したいのですができません。 どうすれば理想どおりにできるようになりますか？ public String getR(){ double rnd = Math.floor(Math.random()*100)-50; if(rnd.length = 4){ rnd = rnd.substring(0,2) ; } else{ rnd = rnd.substring(0,3) ; } _R = rmd; return _R; }
- ベストアンサー
- Java
- 正規分布の乱数生成
C言語で正規分布の乱数を発生させたいのですがどうすればいいのでしょうか？ 自分なりにネットで検索して調べたのですが void gaussrand() { static double V1, V2, S; static int phase = 0; double X; if(phase == 0) { do { double U1 = (double)rand() / RAND_MAX; double U2 = (double)rand() / RAND_MAX; V1 = 2 * U1 - 1; V2 = 2 * U2 - 1; S = V1 * V1 + V2 * V2; } while(S >= 1 || S == 0); X = V1 * sqrt(-2 * log(S) / S); } else X = V2 * sqrt(-2 * log(S) / S); phase = 1 - phase; } こうありました。 例えば平均50の分散9の正規分布の乱数を1000個発生させて、配列seiki[1000]に代入したいときは、このプログラミングをどのようにすればいいのでしょうか？ もちろん、このソースではなく、他のもので説明していただけても全然構いません。 また、もしよろしければ、正規分布の他に、二項分布など他の分布でのデータの生成方法もお教えいただけたら幸いです。 よろしくお願いいたします。
- 締切済み
- C・C++・C#
- 乱数を発生させるプログラムを教えてください。
タイトルのままなのですが、１から１００までの乱数を発生させるプログラムを知りたいです。 乱数をｘとおくと、ｘの値は、０＜ｘ＜１の範囲内でお願いします。 C言語で、rand関数を用いて、どうかお願いします。
- ベストアンサー
- C・C++・C#
- 精度を上げたいのですが…
#include <stdio.h> #include <time.h> #include <stdlib.h> #define MAX 1000 main(void) { int i; float x1, x2, en, sum=0.0, s; srand( (unsigned)time( NULL ) ); for(i=0;i<MAX;i++) { x1=((float)rand()/(float)RAND_MAX); x2=((float)rand()/(float)RAND_MAX); if(en=(x1-0.5)*(x1-0.5)+(x2-0.5)*(x2-0.5)<=(0.5)*(0.5)) { sum++; } } s=sum/MAX; printf("円の面積:%15.6e\n",s); } この方法で円の面積を求めたんですが、もう少し精度を上げたいと思います。ただそのプログラムをどうやって書けばいいのかさえわからずとまどっています。円全体でなくその一部を考え、またその部分を簡単に面積が求められるようにわけるプログラムを組みたいのですがどのようにすればいいのか教えてください。
- ベストアンサー
- C・C++・C#
お礼
右シフトと言うよりは、指数部分の減算を期待しているのですが、これは間違いなのでしょうか。 確かに、ビット合成すれば正確な値は得られそうですが、それにはlong doubleの内部表現を完全に理解し、ビットフィールドを用いなければ、ならないような気がするのですが。