| 查看: 2290 | 回复: 8 | |||
| 本帖产生 1 个 程序强帖 ，点击这里进行查看 | |||
[交流]
【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
|
|||
|
已经知道直线向量a上面的2个点（x1,y1,z1), (x2,y2,z2)和直线向量b上面的2个点（x3,y3,z3), (x4,y4,z4),求直线a，b的夹角。（0~180degree） 夹角公式是cos 其实我最大困难是在fortran中怎么定义这个新的函数-直线向量(恩，这个向量类似一个2组数组)。本人是 对fortran用的不熟。 请大家帮忙写写这段程序，重点就是这个向量的定义，最好能注明清楚~~谢谢。 |
回复此楼» 猜你喜欢
材料科学与工程320求调剂，080500
已经有9人回复
-
0703化学调剂 348分
已经有10人回复
-
307求调剂
已经有12人回复
-
求调剂
已经有13人回复
-
一志愿211，0703化学305分求调剂
已经有15人回复
-
284求调剂
已经有14人回复
-
材料工程调剂
已经有5人回复
-
287求调剂
已经有12人回复
-
362求调剂
已经有12人回复
-
273求调剂
已经有42人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 物理模拟是用c语言好还是fortran好？
已经有30人回复
- 如何用fortran编写高斯白噪声程序
已经有6人回复
- 求初学者能用的fortran编译器
已经有6人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】fortran中怎么存储不确定长度的数据【已完结】
已经有6人回复
- 【求助】这样的公式，在fortran里该怎么写了，谢谢【求助成功】
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
汕头大学海洋科学、生物学、生物与医药接受调剂
+1/1494
-
长江大学化工学院2026年硕士研究生调剂公告（理、工均收）
+1/91
-
天津商业大学接收2026生物、食品类硕士研究生调剂
+2/56
-
东南大学—有机多孔功能材料团队（杰青团队）招收 AI+COFs 方向科研助理
+1/43
-
华侨大学-发光材料与信息显示研究院-招研究生
+1/40
-
【通知】北京信息科技大学仪器科学与光电工程学院招收调剂生（2026）
+1/40
-
安徽工业大学有调剂名额
+1/38
-
福建医科大学生信系饶世涛老师2026年招收调剂优秀硕士生1名-生物医学信息学
+1/38
-
海南医科大学（海南省医学科学院）李东伟课题组招聘研究助理
+2/36
-
盐城工学院化工学院招收07及08开头生化环材食品等相关专业调剂生
+2/32
-
渤海大学储能材料课题组招收调剂研究生
+1/7
-
中国地震局地质研究所+地学/数学相关专业+接收调剂
+1/7
-
福建理工大学材料学院招收专业代码08开头学术型硕士研究生
+1/7
-
找好工作来我这-广东唯一石化院校——资源与环境专硕招生
+1/5
-
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”学硕与专硕有缺额
+1/5
-
北大口腔招人工智能、生物再造科研助理（优秀者可留校或推荐读博）
+1/5
-
【上海调剂】985联合培养！别错过能送你去丹麦/交大/同济/中科院的神仙导师！
+1/4
-
海水鱼类遗传育种课题组招收硕士调剂生
+1/4
-
总分338求调剂，本科大连理工一志愿北航生医工083100
+1/4
-
北京高校副校长团队招收机械类，环境类学硕和专硕
+1/3
★ ★ ★
tpp001(金币+5):很感谢您的帮忙。不过我这里的数据有点复杂，我根据您的这个程序转换不过来，不知道能不能在麻烦您一下，可以吗？输入文件详细的情况在下面回复中，谢谢 2010-12-18 15:38:27
余泽成(金币+3):辛苦了！ 2010-12-19 19:27:55
tpp001(金币+5):很感谢您的帮忙。不过我这里的数据有点复杂，我根据您的这个程序转换不过来，不知道能不能在麻烦您一下，可以吗？输入文件详细的情况在下面回复中，谢谢 2010-12-18 15:38:27
余泽成(金币+3):辛苦了！ 2010-12-19 19:27:55
|
写了一个简单的小程序，可以满足你使用的要求：   输入文件的含义如下：  下面是三个输入文件以及测试结果  测试输入对应的两个矢量如下图所示：   测试输入对应的两个矢量如下图所示：   测试输入对应的两个矢量如下图所示：  [ Last edited by maomao1210 on 2010-12-18 at 14:59 ] |
2楼2010-12-18 14:56:39
|
真的不好意思，我的看懂了您上面的代码，可是我还是不会用来写我这个复杂点的情况，可能我有点笨了，不知道能不能麻烦您看看我这个情况，谢谢。。。。 总共有29条直线向量，每条直线向量上都知道2个点的坐标，A（a1,a2,a3），B（b1,b2,b3）,【所以输入文件中有58个点的坐标】而我要得出他们任意2条直线向量的夹角。。。 输入文件的格式是这样的 a1 a2 a3 b1 b2 b3 #这是第一条直线向量上面的2个点坐标 c1 c2 c3 d1 d2 d3 #这是第2条直线向量上面的2个点坐标 ...... ..... (n-1)1 (n-1)2 (n-1)3 n1 n2 n3 #这是第n条直线向量上面的2个点的坐标 一共58行数据 |
3楼2010-12-18 15:51:00
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zyj8119(金币+10, 程序强帖+1):希望maomao1210继续努力。。。 2010-12-18 16:52:40
小木虫(金币+0.5):给个红包，谢谢回帖交流
zyj8119(金币+10, 程序强帖+1):希望maomao1210继续努力。。。 2010-12-18 16:52:40
|
按照你的要求，写的如下程序，可以符合你的要求：    下面是我用四条线的简单测试结果：  那个黑色的方阵 是计算结果，两两LINE之间的，构成了一个对称矩阵。 C C ==================================================== C www.muchong.com C ==================================================== C Implicit None Integer Idicp,Jdicp,Kdicp,Ldicp Integer Num_Max,Num_Line Real*8 X1,Y1,Z1,X2,Y2,Z2 Parameter(Num_Max=1000) Real*8 A(3,Num_Max) Real*8 Theta(Num_MAx,Num_MAx),Output Open(10,File='input',Status='Old') Open(20,File='Output',Status='Unknown') C C Read necessary data C Read (10,*) Num_Line Do Idicp =1, Num_Line Read(10,*) X1,Y1,Z1 Read(10,*) X2,Y2,Z2 A(1,Idicp)= X2-X1 A(2,Idicp)= Y2-Y1 A(3,Idicp)= Z2-Z1 EndDo C C Find the theta between every possible combin two lines C Do Idicp=1, Num_Line Do Jdicp=1,Idicp Theta(Idicp,Jdicp)= 0.0D0 Call Sub(A(1,Idicp),A(1,Jdicp),Output) Theta(Idicp,Jdicp)= Output Theta(Jdicp,Idicp)= Theta(Idicp,Jdicp) EndDo Enddo C C Output the matrice C Do Idicp=1, Num_Line Write(*,'(30f12.3)') (Theta(Idicp,Jdicp),Jdicp=1,Num_Line) Write(20,'(30f12.3)') (Theta(Idicp,Jdicp),Jdicp=1,Num_Line) Enddo End Subroutine Sub(A,B,Output) Implicit None Real*8 A(3),B(3) Real*8 Dicp,AA,BB,AB,Output AB=A(1)*B(1)+A(2)*B(2)+A(3)*B(3) AA=A(1)*A(1)+A(2)*A(2)+A(3)*A(3) BB=B(1)*B(1)+B(2)*B(2)+B(3)*B(3) Dicp = AB / ( Dsqrt(AA) * Dsqrt(BB) ) Output= Acos(Dicp)*180.0d0/3.1415926D0 Return End Enjoy It !!!! [ Last edited by maomao1210 on 2010-12-18 at 16:37 ] |
4楼2010-12-18 16:34:04
|
非常感谢版主您的帮忙，回复的也非常及时。我有一点疑问向您请教。 在调用子程序的时候，在主程序中式这样的写的： Call Sub(A(1,Idicp),A(1,Jdicp),Output) 而在子程序中： Subroutine Sub(A,B,Output) Implicit None Real*8 A(3),B(3) Real*8 Dicp,AA,BB,AB,Output AB=A(1)*B(1)+A(2)*B(2)+A(3)*B(3) AA=A(1)*A(1)+A(2)*A(2)+A(3)*A(3) BB=B(1)*B(1)+B(2)*B(2)+B(3)*B(3) 在子程序中需要用到的自变量A(3)是一个一维数组，大小是3。 而您在主程序中设定的自变量A(1,Idicp)是一个二位数组。我就是这一点想不通，因为：这个A(1,Idicp)在主程序中我认为它对应了一个具体的数据：A(1,Idicp)= X2-X1。而在子程序中会用到A(1),A(2),A(3),这3个数据的，所以我老是觉得它们对应不上去。  还有一个问题是：在输出结果的时候，我希望输出的时候，没有重复的夹角出现比如line 1与 line2 之间的夹角，和line 2 与line 1之间的夹角是一样，只要输出一个就行。还有就是 同一直线的夹角也不用输出来。因为我后面还个程序输出的数据进行一些分析的。谢谢 |
5楼2010-12-19 10:34:05
6楼2010-12-19 11:02:45
|
您好，首先谢谢您上面的程序，最近我才用这个程序来处理我的数据。出现了一些问题，我不太清楚是我的数据的问题还是在数学模型中出现了问题。当然程序肯定是没有问题的。 我把我的数据（输入文件和程序）传上来。请您帮我分析一下是什么问题。 我的输入文件中有27条直线，因为每条直线有2个点，所以我的输入文件中有54个点。当我设定只是计算开头的2条直线的夹角时，程序可以正常运行。可是当我计算开头的4条直线中任意2条直线的夹角时，程序出现错误。。。会是什么问题呢，我搞不懂？我在下面还把我的出错信息给贴出来  [ Last edited by tpp001 on 2011-1-6 at 20:20 ] |
7楼2011-01-06 20:17:44
8楼2011-01-07 21:01:50
9楼2011-01-07 22:16:52
