C言語の課題について

質問

何度も申し訳ありません。
http://okwave.jp/qa/q7877142.html
や
http://okwave.jp/qa/q7878947.html
より、回答を受けてプログラムを作成しているのですが、
衝突時の時間と座標をどのように作成すれば良いか分からなく、
一秒後の点Pの位置を算出したのですが、合っているかどうかが分かりません。
非常に難しくて困っています。お助け下さい。

ソースコード
#include <stdio.h>
#include<math.h>

int main(void)
{
float p,v,vx,vy,x,y,x1,y1,katamuki,gyaku,seiki_x,seiki_y,n,touei_x,touei_y,idoux,idouy,px,py;

x=1.0;
y=1.0;
vx=1.0; //x方向の速さ
vy=-2.0; //y方向の速さ
y1=(1/y); //逆数
x1=(1/x); //逆数
seiki_x=x1/(sqrt(x1*x1+y1*y1));
seiki_y=y1/(sqrt(x1*x1+y1*y1));
n=-(vx*seiki_x)+(vy*seiki_y); //nは内積
touei_x=n*seiki_x; //touei_xは投影x軸
touei_y=n*seiki_y; //touei_yは投影y軸
idoux=vx+2*touei_x; //衝突後、点Pのxの移動方向
idouy=vy+2*touei_y; //衝突後、点Pのyの移動方向
px=-idoux; //点Pのx座標
py=-idouy; //点Pのy座標

printf("1秒後の位置は点P(x,y)=(%f,%f)",px,py);

return 0;
}

出力結果
1秒後の位置は点P(x,y)=(2.000000,5.000000)

何かおかしなところがあれば追加や修正等をしてくれれば幸いです。

回答No.1

f272 · ベストアンサー

ちゃんとやる気はあるようですね。
これでも見て理解してみてね。

#include<stdio.h>
#include<math.h>

int main(void)
{
double x=1.0,y=1.0;//初期位置
double vx=1.0,vy=-2.0;//初期速度
double ax=1.0/3.0,ay=-1.0;//壁は（1/3）x-y=0
double nx,ny,px,py,dot,ss,tc,t1=1.0;

//t秒後の位置(x+vx*t,y+vy*t)だからax*(x+vx*t)+ay*(y+vy*t)=0のとき衝突する
//つまりt=-(ax*x+ay*y)/(ax*vx+ay*vy)のとき
tc=-(ax*x+ay*y)/(ax*vx+ay*vy);
if (tc<1.0) {//衝突は1秒以内
x+=vx*tc;//衝突位置
y+=vy*tc;
ss=sqrt(ax*ax+ay*ay);
nx=ax/ss;//壁の正規化した法線
ny=ay/ss;
dot=(vx*nx+vy*ny);//内積
px=-dot*nx;//投影ベクトル
py=-dot*ny;
vx+=2.0*px;//衝突後速度
vy+=2.0*py;
x+=vx*(t1-tc);//t1秒後位置
y+=vy*(t1-tc);
} else {//衝突は1秒よりあと
x+=vx*t1;//t1秒後位置
y+=vy*t1;
}

printf("1秒後の位置は点P(x,y)=(%f,%f)\n",x,y);

return 0;
}

お礼

すごく丁寧な回答ありがとうございます。
おかげさまでかなり理解を深めることができました。

また、何度も質問してしまって、申し訳ありませんでした。