式をC言語で立てれません
式をC言語で立てれません
式をC 言語で立てれません
こんにちは。私は今実験で円板振動子の中心軸上の音圧分布を求めています。
そこで写真のような図をプロットするために与えられた式をC言語で書いています。
式の詳細は
P/2ρcV ＝ |sinπ(√(r/λ)^2 + (a/λ)^2)(√はここまで)-(r/λ))|
です。
しかしうまくプロットできません。
プログラムは
-------------------------------------
#include <stdio.h>
#include <math.h>
#define F 120 //周波数(Hz)
#define C 1500 //音速(m/s)
#define A 0.2 // 半径(m)
int main()
{
double x, m, s, p, o, w, l;
int i;
for (i=0; i<100; i++)
{ x= (double)i; //観測点までの距離を表しました。距離が変化することで音圧が変わるためです。
l = C/F; //λをlとして音速÷周波数で定義できます。
m = x/l; //x軸です。
o = m*m + (A/l)*(A/l); //式のsinの中身です。
s = sin(M_PI*(sqrt (o) - m)); //sin全体を定義します。
w = fabs(s); //sin全体を絶対値に置き換えます。
p = w; //新しい変数に入れます。
printf("%9.9f %9.9f＼n", m, p);
}
return 0;
}
------------------------------
といった感じです。先生からはx軸のプログラムだけいじれば勝手に
y軸も出てくるとヒントを頂いているのでx軸についてのプログラムを
書いています。コンパイルはできても写真のようにいきません。
半径やλや周波数はa/λ=2.5と書いてあったので値は推測です。
C言語はあまり得意ではないので困っています。
お分かりになる方、ご教授お願い致します。
※式が間違っていたので新しく質問します。
