24小时热门版块排行榜

返回列表

astringent

铜虫 (著名写手)

应助: 16 (小学生)
金币: 2210.7
帖子: 1180
在线: 354.8小时
虫号: 900839

[交流] 【求助】关于向量投影的问题

有两个点A（x1,y1,z1）,B（x2,y2,z2）,现在以由点A和B确定的直线为基准轴，我想求的是当以物体由点C（x3,y3,z3）运动到点D（x4,y4,z4）时，沿AB基准轴方向的位移以及沿垂直于AB基准轴方向的位移。我写了一个程序，劳烦大家帮忙看看，这个程序有没有问题以及能不能实现我的目的?请大家不吝赐教。
program get_projection
implicit none
real x1,y1,z1
real x2,y2,z2
real x3,y3,z3
real x4,y4,z4
real theta
real dx1,dy1,dz1,dx2,dy2,dz2
real px1,py1,pz1,px2,py2,pz2
real r1,r2

! 求基准轴AB
  dx1=x2-x1
  dy1=y2-y1
  dz1=z2-z1
!求CD的位移（前保存的结果）
dx2=x4-x3
  dy2=y4-y3
  dz2=z4-z3
!求AB和CD的模
  r1=dx1**2+dy1**2+dz1**2
r2=dx2**2+dy2**2+dz2**2
！由点乘公式计算AB与CD所形成的角theta
theta=acos ((dx1*dx2+dy1*dy2+dz1*dz2)/(sqrt(r1)*sqrt(r2)))
!计算沿AB基准轴方向的位移
  px1=dx2*cos(theta)
  py1=dy2*cos(theta)
  pz1=dz2*cos(theta)
!沿垂直于AB基准轴方向的位移
  px2=dx2*sin(theta)
  py2=dy2*sin(theta)
  pz2=dz2*sin(theta)
end

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有54人回复
垃圾破二本职称评审标准已经有19人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

高维数组投影到子空间已经有4人回复
向量的投影已经有4人回复
【求助】关于经典物理中时间各向同性的理解已经有13人回复
【求助】关于向量投影的问题已经有3人回复
【求助】国外的一道作业题，不会做！已经有10人回复
【讨论】ARCGIS中如何实现同一区域不同投影和比例尺的两张矢量图的要素叠合？已经有4人回复
【求助】transiesta中关于轨道的问题已经有18人回复
【转帖】理解矩阵已经有51人回复
【讨论】NEB寻找过渡态专题讨论已经有49人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-01-21 04:03:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
余泽成(金币+2): 谢谢参与应助！ 2011-01-21 13:28:24
astringent(金币+5): 2011-01-22 05:27:21

引用回帖:

Originally posted by astringent at 2011-01-21 04:03:30:
有两个点A（x1,y1,z1）,B（x2,y2,z2）,现在以由点A和B确定的直线为基准轴，我想求的是当以物体由点C（x3,y3,z3）运动到点D（x4,y4,z4）时，沿AB基准轴方向的位移以及沿垂直于AB基准轴方向的位移。我写了一个程序， ...

程序本身语法似乎没有问题，但能不能达到你的目的就不清楚了，这是算法的问题，呵呵……

赞一下(1人)

回复此楼

2楼2011-01-21 09:12:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

leedobb

金虫 (正式写手)

程序强帖: 1
应助: 77 (初中生)
金币: 152.4
帖子: 856
在线: 270.4小时
虫号: 1199535

★ ★
astringent(金币+5): 2011-02-17 02:57:21
astringent(金币+10): 2011-02-17 03:11:16
余泽成(金币+2): 谢谢参与应助，欢迎常来程序语言版！ 2011-02-20 18:34:56

发现一个小问题，
把下面一段改一下

" px1=dx2*cos(theta)
  py1=dy2*cos(theta)
  pz1=dz2*cos(theta)
!沿垂直于AB基准轴方向的位移
  px2=dx2*sin(theta)
  py2=dy2*sin(theta)
  pz2=dz2*sin(theta)
“

==》
”
px1= dx1/sqrt(r1) * sqrt(r2) *cos(theta)
py1= dy1/sqrt(r1) * sqrt(r2) *cos(theta)
pz1= dz1/sqrt(r1) * sqrt(r2) *cos(theta)
!沿垂直于AB基准轴方向的位移
  px2=dx2-px1
  py2=dy2-py1
  pz2=dz2-pz1
“

原程序里p1求的向量是沿cd方向的，所以不对
p2也沿cd方向。

赞一下(1人)

回复此楼

3楼2011-02-16 16:25:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖