24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1528  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

astringent

铜虫 (著名写手)


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

有两个点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
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

匿名

用户注销 (知名作家)


★ ★
astringent(金币+10): 2011-02-17 02:57:30
余泽成(金币+2): 谢谢参与应助,欢迎常来程序语言版! 2011-02-20 18:35:06
本帖仅楼主可见
4楼2011-02-16 17:10:15
已阅   申请程序强帖   回复此楼   编辑   查看我的主页
查看全部 4 个回答
★ ★
余泽成(金币+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基准轴方向的位移。我写了一个程序, ...

程序本身语法似乎没有问题,但能不能达到你的目的就不清楚了,这是算法的问题,呵呵……
2楼2011-01-21 09:12:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)


★ ★
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方向。
3楼2011-02-16 16:25:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +7 璃茉一定上岸 2026-04-10 7/350 2026-04-10 15:25 by jiajinhpu
[考研] 298求调剂 +8 钉叮咚冬瓜 2026-04-07 10/500 2026-04-10 11:25 by 高维春
[考研] 材料专硕283求调剂 +18 试试看呗 2026-04-04 19/950 2026-04-10 10:05 by may_新宇
[考研] 一志愿华中农微生物,288分,三年实验经历 +10 代fish 2026-04-09 10/500 2026-04-10 09:49 by potato妹
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 085600,321分求调剂 +19 大馋小子 2026-04-04 20/1000 2026-04-09 14:12 by Delta2012
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 一志愿211,0703化学305分求调剂 +20 严西西戏 2026-04-06 26/1300 2026-04-08 22:12 by 凯凯要变帅
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 338求调剂 +5 小猪红色 678 2026-04-06 6/300 2026-04-07 21:18 by 乔哒哒哒
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 296求调剂 +3 汪!?! 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 292求调剂 +11 2022080213 2026-04-04 13/650 2026-04-04 18:38 by macy2011
[考研] 22408,264求调剂 +3 ywh729 2026-04-03 4/200 2026-04-04 11:04 by ywh729
信息提示
请填处理意见