vba

###前提・実現したいこと
現在Visual basicのプログラムをC言語に書き換えています
そこでRandomizeの関数がでてくるのですが、これをどのように書き換えればいいのかわかりません
自分なりにプログラムを作りましたが、同じような値が出力されません
特にわからない箇所が
Randomize4の意味は毎回同じシートを使って乱数を指定しているので乱数の出方のパターンは同じということですか？
また、Rnd1とRnd2に代入される乱数は同じ値ですか？

最終的にはこの部分と同じプログラムをC言語で書きたいです
なので書き換えの部分で変な箇所があればどのように直すのか教えていただきたいです

Visual basic
以下はプログラムの途中部分です
Randomize 4
PI = 4 * Atn(1)
For i = 1 To N / 2 - 1
For j = 1 To N / 2 - 1

        If i <> 0 Or j <> 0 Then
        
        kx = i / (region)
        ky = j / (region)
        K = Sqr(kx ^ 2 + ky ^ 2)
        Rnd1 = Rnd
        Rnd2 = Rnd
        
        gauss = Abs(Sqr(-2 * Log(Rnd1)) * Cos(2 * PI * Rnd2))
                    
         '****************************
         'gaussの上限，下限 (04.6.20)
            
            If gauss < 0.001 Then
                gauss = 0.001
            ElseIf gauss > 10 Then
                gauss = 10
            End If
         '****************************
                    
        Spec = ParaA1 * (K ^ (-ParaB1))
        Spec2 = Spec * (gauss ^ flag)
        rad = Sqr(Spec2 * 0.25)
             
       
           If Sqr(kx ^ 2 + ky ^ 2) >= lim Then
              phase = PI / 2
           Else
              phase = 2 * PI * Rnd
           End If
                

        A(i, j) = rad * Cos(phase)
        B(i, j) = rad * Sin(phase)
        A(N - i - 1, N - j - 1) = rad * Cos(phase)
        B(N - i - 1, N - j - 1) = rad * Sin(phase)
        A(i, N - j - 1) = rad * Cos(phase)
        B(i, N - j - 1) = rad * Sin(phase)
        A(N - i - 1, j) = rad * Cos(phase)
        B(N - i - 1, j) = rad * Sin(phase)
        End If

以下C言語の書き換えです
visual basicでRnd1,Rnd2に代入される乱数が同じとして作ってみました
srand(3);
pi=4atan(1)
rnd = (double)rand()/RAND_MAX;
for(i = 0;i <= n/2-1;i++){
for(j = 0;j <= n/2-1;j++){
kx = i/region;
ky = j/region;
k = sqrt(kxkx + ky*ky);
gauss = fabs( sqrt(-2 * log(rnd)) * cos(2 * pi * rnd) );
/gaussの上限，下限 (04.6.20)/
if(gauss < 0.001){
gauss = 0.001;
}else{
if(gauss > 10){
gauss = 10;
}
}

  if(k == 0){
    spec = 1000;
  }else{
    spec = paraa1 * pow(k,-parab1);
   }

  spec2 = spec * pow(gauss,flag);
  rad = sqrt(spec2 * 0.25);
  if(sqrt(kx*kx + ky*ky) >= lim){
    phase = pi / 2;
  }else{
     phase = 2 * pi * rnd;
   }

  a[i][j] = rad * cos(phase);
  b[i][j] = rad * sin(phase);
  a[n-i-1][n-j-1] = rad * cos(phase);
  b[n-i-1][n-j-1] = rad * sin(phase);
  a[i][n-j-1] = rad * cos(phase);
  b[i][n-j-1] = rad * sin(phase);
  a[n-i-1][j] = rad * cos(phase);
  b[n-i-1][j] = rad * sin(phase);
}

}

行動規範の内容に同意します

回答2件

ベストアンサー

Randomize4の意味は毎回同じシートを使って乱数を指定しているので乱数の出方のパターンは同じということですか？

そうです。

また、Rnd1とRnd2に代入される乱数は同じ値ですか？

違います。Rndは関数なので、使う度に異なる値が返ってきます。

VBのRndは変数のように見えますが関数なので、当然C言語側もRndが書かれていた場所で乱数を取得しないといけません。

rnd()関数を用意します。

C
1double rnd(void)
2{
3	return (double)rand()/RAND_MAX;
4}

VBのRnd関数を使っている箇所をrnd()関数に置き換えます。

C
1// gaussの計算のところ
2gauss = fabs( sqrt(-2 * log(rnd())) * cos(2 * pi * rnd()) );
3
4// 中略
5
6// phaseの計算のところ
7phase = 2 * pi * rnd();

ループ手前のrnd = (double)rand()/RAND_MAX;は不要なので削除してください。変数宣言も削除してください。

投稿2016/05/27 11:11

編集2016/05/27 11:15

catsforepaw

総合スコア5944

退会済みユーザー

2016/05/27 12:52

visual basicのrndの戻り値が0以上1未満になっているのですが C言語で書き換えてみたら戻り値が整数になっているようで計算ができません乱数の戻り値を1より小さくするにはどうしたらいいですか？

catsforepaw

2016/05/27 12:55

上のrnd関数のコード例のようにやれば、0以上1未満の値が得られます。

退会済みユーザー

2016/05/27 13:10 編集

すいません勘違いしました気にしないでください

catsforepaw

2016/05/27 13:15

多分そうだろうなとは思いました。ご自身のソースでもやってましたし。気にしてませんよ。

退会済みユーザー

2016/05/27 13:41

ありがとうございます値を確かめるために printf("kx=%d ",kx); printf("gauss=%d ",gauss);のプログラムを追加したのですが表示されません変数iを含んだ計算式はこの形では表示されないんですか？

catsforepaw

2016/05/27 13:48

printfの書式は変数の型に合わせないといけません。kxもgaussもおそらくdouble型だと思うのですが、「%d」だとint型を表示するためのものなので正しく表示してくれません。double型を表示する場合は「%f」か「%g」を使ってください。

退会済みユーザー

2016/05/27 14:02

for(i = 0;i <= n/2-1;i++){ j = 0; if(i != 0){ printf("i=%d ", i); printf("region=%f ", region); kx = i / (region); ky = j / (region); k = sqrt(kx*kx + ky*ky); printf("kx=%f ", kx); printf("k=%f ", k); printf("parab1=%lf ", parab1); gauss = fabs( sqrt(-2 * log(rnd)) * cos(2 * pi * rnd) ); /*gaussの上限，下限　(04.6.20)*/ if(gauss < 0.001){ gauss = 0.001; }else{ if(gauss > 10){ gauss = 10; } } spec = 5.14 * paraa1 * pow(k,-parab1); spec2 = spec * pow(gauss,flag); rad = sqrt(spec2 * 0.25); if(sqrt(kx*kx) >= lim){ phase = pi / 2; }else{ phase = 2 * pi * rnd; } a[i][0] = rad * cos(phase); b[i][0] = rad * sin(phase); a[n-i-1][0] = rad * cos(phase); b[n-i-1][0] = rad * sin(phase); } ｝上記のように%fを使っても表示されませんでした書き方に何か問題があるのでしょうか？

catsforepaw

2016/05/27 14:15

n には 4 以上の値が入っていますか？　n が3以下だと n/2-1 は0になって`if(i != 0)`の条件に合致することなくfor文が終了してしまいます。

退会済みユーザー

2016/05/27 14:24

nには512が入るので4以上にはなります他に何か理由がみあたりますか？

catsforepaw

2016/05/27 14:31

そうなると問題は根深いですね。`printf("i=%d ", i);`をfor文の直後に書いてみてください。そもそもfor文自体が実行されているかどうか疑わしいです。

退会済みユーザー

2016/05/27 14:54

printf("i=%d ", i);をfor文の直後においたらなぜかkxが出ました