[求助]：谁能给一个求点到直线的垂点的代码

aeo

要 代数 求的那种.

libo2000

2. 
   
3.   [FONT=courier new]
   
4. ;垂足：点到线
   
5. (defun chengs_getpoint_per (p0 p1 p2 / r12 p3)
   
6.   (setq
   
7.            r12 (angle p1 p2)
   
8.            r0 (+ r12 (/ pi 2))
   
9.            p3 (polar p0 r0 1.0)
   
10.   )
    
11.   (inters p1 p2 p0 p3 nil)
    
12. )
    
13. ;垂足2：线上点p0,线点p1 p2；w垂足距线长度
    
14. (defun chengs_getpoint_per2 (p0 p1 p2 w / r12 p3)
    
15.   (setq
    
16.            r12 (angle p1 p2)
    
17.            r0 (+ r12 (/ pi 2))
    
18.            p3 (polar p0 r0 w)
    
19.   )
    
20. )
    
21. ;垂足3：点－>线
    
22. (defun chengs_getpoint_per3 (p0 ln1 / p1 p2  ln1_list r12 p3)
    
23.   (if (cdr ln1) (setq ln1 (car ln1)))
    
24.   (setq ln1_list (entget ln1))
    
25.   (setq p1 (cdr (assoc 10 ln1_list)))
    
26.   (setq p2 (cdr (assoc 11 ln1_list)))
    
27.   (setq
    
28.            r12 (angle p1 p2)
    
29.            r0 (+ r12 (/ pi 2))
    
30.            p3 (polar p0 r0 1.0)
    
31.   )
    
32.   (inters p1 p2 p0 p3 nil)
    
33. )
    
34. [/FONT]
    

<br />   [FONT=courier new]<br /> ;垂足：点到线<br /> (defun chengs_getpoint_per (p0 p1 p2 / r12 p3)<br />   (setq<br />            r12 (angle p1 p2)<br />            r0 (+ r12 (/ pi 2))<br />            p3 (polar p0 r0 1.0)<br />   )<br />   (inters p1 p2 p0 p3 nil)<br /> )<br /> ;垂足2：线上点p0,线点p1 p2；w垂足距线长度<br /> (defun chengs_getpoint_per2 (p0 p1 p2 w / r12 p3)<br />   (setq<br />            r12 (angle p1 p2)<br />            r0 (+ r12 (/ pi 2))<br />            p3 (polar p0 r0 w)<br />   )<br /> )<br /> ;垂足3：点－>线<br /> (defun chengs_getpoint_per3 (p0 ln1 / p1 p2  ln1_list r12 p3)<br />   (if (cdr ln1) (setq ln1 (car ln1)))<br />   (setq ln1_list (entget ln1))<br />   (setq p1 (cdr (assoc 10 ln1_list)))<br />   (setq p2 (cdr (assoc 11 ln1_list)))<br />   (setq<br />            r12 (angle p1 p2)<br />            r0 (+ r12 (/ pi 2))<br />            p3 (polar p0 r0 1.0)<br />   )<br />   (inters p1 p2 p0 p3 nil)<br /> )<br /> [/FONT]<br />

aeo

最初由 libo2000 发布
[B]
  
;垂足：点到线
(defun chengs_getpoint_per (p0 p1 p2 / r12 p3)
  (setq
           r12 (angle p1 p2)
           r0 (+ r12 (/ pi 2))
           p3 (polar p0 r0 1.0)
  )
  (inters p1... [/B]

inters,谁都会想到的.

我要的:
已知 p0 p1 p2坐标.
用* / + - sin cos 等 数学公式求.

panzhongw

看看我的程序

nuance

(setq a (distance p1 p2))
(setq b (distance p0 p1))
(setq c (distance p0 p2))
(setq cosang (/ (- (+ (* a a)(* b b))(* c c))(* 2 a b)))
(setq l (* b cosang))
(setq pt-per (polar p1 (angle p1 p2) l))

余弦定理是a*a+b*b-2abcosang=c*c吧？我忘了
不知道你做什么用？

mmmm

不用那么复杂，根据三条边长计算面积（由数学公式）
在用标准三角形面积计算方法算：高

aeo

好象谁也逃不掉 inters distance polar ,都不能用,就数学公式可用

比如求p1(x1 y1) p2(x2 y2)的距离: (sqrt (+ (* (- x2 x1)(- x2 x1))(*(- y2 y1)(- y2 y1))))
就要这种方法.
举个例子:
晓东贴过的,求三角形面积:

2. 
   
3. (defun tarea (p1 p2 p3)
   
4.   (- (+ (* (car p1) (cadr p2)) (* (car p2) (cadr p3)) (* (car p3)
   
5.                              (cadr p1)
   
6.                               )
   
7.     ) (* (car p1) (cadr p3)) (* (car p2) (cadr p1)) (* (car p3)
   
8.                             (cadr p2)
   
9.                             )
   
10.   )
    
11. )
    

<br /> (defun tarea (p1 p2 p3) <br />   (- (+ (* (car p1) (cadr p2)) (* (car p2) (cadr p3)) (* (car p3) <br />                              (cadr p1) <br />                               ) <br />     ) (* (car p1) (cadr p3)) (* (car p2) (cadr p1)) (* (car p3) <br />                             (cadr p2) <br />                             ) <br />   ) <br /> ) <br />

(可能要除于2 .)

eachy

2. 
   
3. ;;; 返回点到直线距离
   
4. ;;;返回值：点到直线的距离，点在线段的左(上)侧为负值, 在右(下)侧为正值。;
   
5. ;;;方向根据给定的起始点和终止点定。
   
6. ;;; Request distance from a point(p0) to a line(p1 p2).
   
7. (defun pntoln (p0 p1 p2 / x1 y1 x2 y2 c1 c2 c3)
   
8.   (setq        x1 (car p1)
   
9.         y1 (cadr p1)
   
10.         x2 (car p2)
    
11.         y2 (cadr p2)
    
12.         c1 (- y2 y1)
    
13.         c2 (- x1 x2)
    
14.         c3 (- (* x2 y1) (* x1 y2))
    
15.   )
    
16.   (/ (+ (* c1 (car p0)) (* c2 (cadr p0)) c3)
    
17.      (sqrt (+ (* c1 c1) (* c2 c2)))
    
18.   )
    
19. )

<br /> ;;; 返回点到直线距离<br /> ;;;返回值：点到直线的距离，点在线段的左(上)侧为负值, 在右(下)侧为正值。;<br /> ;;;方向根据给定的起始点和终止点定。<br /> ;;; Request distance from a point(p0) to a line(p1 p2).<br /> (defun pntoln (p0 p1 p2 / x1 y1 x2 y2 c1 c2 c3)<br />   (setq        x1 (car p1)<br />         y1 (cadr p1)<br />         x2 (car p2)<br />         y2 (cadr p2)<br />         c1 (- y2 y1)<br />         c2 (- x1 x2)<br />         c3 (- (* x2 y1) (* x1 y2))<br />   )<br />   (/ (+ (* c1 (car p0)) (* c2 (cadr p0)) c3)<br />      (sqrt (+ (* c1 c1) (* c2 c2)))<br />   )<br /> )

aeo

最初由 eachy 发布
[B]
;;; 返回点到直线距离
;;;返回值：点到直线的距离，点在线段的左(上)侧为负值, 在右(下)侧为正值。;
;;;方向根据给定的起始点和终止点定。
;;; Request distance from a point(p0) to a line(p1 p2).
(d... [/B]

以前见过了,我就是靠它已经解决问题了.

panzhongw

能说一下为什么不用INTERS等函数，有什么特别的意义吗？

netbee

;;;===============================
;;;(ppt ptt1 ptt2 ptt3)返回点pt1至直线pt2 pt3的垂点坐标
(defun ppt (ptt1 ptt2 ptt3 / pt osmode_old angtmp)
  (if (/= (getvar "osmode") 0)
    (progn (setq osmode_old (getvar "osmode"))
           (setvar "osmode" 0)
    )
  )
  (setq angtmp (-(angle ptt2 ptt1)(angle ptt2 ptt3)))
  (if (> angtmp pi)
    (setq angtmp (- (* pi 2) angtmp))
  )

  (if (> angtmp (/ pi 2))
    (setq pt (polar ptt2
                    (+ (angle ptt2 ptt3) pi)
                    (sqrt (- (^ (distance ptt1 ptt2) 2)
                             (^ (sort ptt1 ptt2 ptt3) 2)
                          )
                    )
             )
    )
    (setq pt (polar ptt2
                    (angle ptt2 ptt3)
                    (sqrt (- (^ (distance ptt1 ptt2) 2)
                             (^ (sort ptt1 ptt2 ptt3) 2)
                          )
                    )
             )
    )
  )
  (if osmode_old
    (setvar "osmode" osmode_old)
  )
  pt
)