Cプログラミングの質問なのですが，

質問

Cプログラミングの質問なのですが，
以下のプログラムで正規乱数を発生させたいのですが，どこがおかしいのでしょうか？
fp1のransuuはきちんとtxtで作成されています。

至急お助けください。

#include <stdio.h>
#include<stdlib.h>
#include<math.h>
#define PI 3.141592653589793238

int main (void)
{
FILE *fp1,*fp2;
int i,n;
unsigned int x1,x2;
double y1,y2;
fp1=fopen("ransu.txt","r");
fp2=fopen("seikiransu.txt","w");

for(i=0;i<n;i++)
{
fscanf(fp1,"%lf",&x1);
fscanf(fp1,"%lf",&x2);

y1=sqrt(2)*sqrt(-2*log(x1))*cos(2*PI*x2);
fprintf(fp2,"%lf\n",y1);

}
fclose(fp1);
fclose(fp2);
return 0;
}

回答No.5

TreatMeGently · ベストアンサー

＞for(i=0;i<n;i++)
　まず、前回答者さんも御指摘ですが、変数ｎが初期化されていませので無限ループかメモリ破壊
　を起こす可能性が有ります。
＞fscanf(fp1,"%lf",&x1);
　ｆですのでDoubleに変換しますので、x1もx2も破壊されます。
　倍精度浮動小数点ですのでスタック上の変数が破壊されます。
　CPUによって異なりますが、８バイト以上だと思いますので
　スタックを壊して変換します。
　従いましてDouble　x1、x2宣言が良いのでは？
＞y1=sqrt(2)*sqrt(-2*log(x1))*cos(2*PI*x2);
　ここでもcastしないと桁落ちするかも。
＞fprintf(fp2,"%lf\n",y1);
　このｆの前は”いち”ですか？”エル”ですか？
　”エル”なら”％ｌｆ”でDoubleでokですが、”いち”なら
　１桁だけの表示になります。

回答No.6

Tacosan

なんで x1 や x2 を unsigned int にしたの?
Box-Muller 法は完全に理解できている?

回答No.4

D-Matsu

「エラーがある」だけではなく「どんなエラーが」「どこで」出るかも書きましょう。
コンパイルエラーならエラーメッセージ本体と出ている行、実行時エラーならデバッガなどで出る箇所を特定しておくと回答側も追跡がかなり楽になります。

今のコードですぐに気付く分としては、nが不定値（初期化されていない）なので何回forが回るかは全く未定義です。
nをx1, x2の数で初期化すれば動作すると思いますが、そもそも読み終わり判定ならfeof()とwhileの組み合わせの方が確実です。

回答No.3

Tacosan

なんというか, Box-Muller 法をいろいろと勘違いしているように見える.
求める正規乱数の平均と分散はそれぞれどのような値にしたい?

お礼

平均0、分散2にしたいです

回答No.2

D-Matsu

#1でも言われてますが、scanf系の%lfはdoubleに対応しますがx1, x2の型はunsigned intになっていて不整合が起きています。
また、sizeof(double) > sizeof(unsigned int)の環境下では正しく値が取れないのみならず、メモリ破壊により予期せぬバグの原因となるかsegmentation faultの原因になります。
x1, x2の型が正しいのならunsigned intに対応する%uで取らなければならないし、型が間違っているのなら正しくdoubleで宣言しなければなりません。

お礼

ありがとうございます。
#include <stdio.h>
#include<stdlib.h>
#include<math.h>
#define PI 3.141592653589793238

int main (void)
{
FILE *fp1,*fp2;
int i,n;
unsigned int x1,x2;
double y1,y2;
fp1=fopen("ransu.txt","r");
fp2=fopen("seikiransu.txt","w");

for(i=0;i<n;i++)
{
fscanf(fp1,"%u",&x1);
fscanf(fp1,"%u",&x2);

y1=sqrt(2)*sqrt(-2*log(x1))*cos(2*PI*x2);
fprintf(fp2,"%lf\n",y1);

}
fclose(fp1);
fclose(fp2);
return 0;
}

上記のように変更したところまだエラーがあるようです。
ほかになにが問題点はあるのでしょうか？

回答No.1

Wr5

>unsigned int x1,x2;
>fscanf(fp1,"%lf",&x1);
>fscanf(fp1,"%lf",&x2);

型が合っていないように見えますが？