由三角的三个点返回三角形的角度计算参数

鬼子兵

;;由三角的三个点返回三角形的角度计算参数
(defun c:te()
        (setq fanhui '())
        ;三角
        (setq p_a (getpoint "\np1"))
        (setq p_b (getpoint "\np2"))
        (setq p_c (getpoint "\np3"))
        ;;三边
        (setq l_a (distance p_b p_c))
        (setq l_b (distance p_a p_c))
        (setq l_c (distance p_a p_b))
        ;计算P
        (setq js_p (* (+ l_a l_b l_c) 0.5))
        ;;三角面积
        (setq s_1 (sqrt (* js_p (- js_p l_a) (- js_p l_b) (- js_p l_c))))
        ;;sina\sinb\sinc
        (setq sin_a (/ s_1 (* l_b l_c 0.5)))
        (setq sin_b (/ s_1 (* l_a l_c 0.5)))
        (setq sin_c (/ s_1 (* l_a l_b 0.5)))
        ;;初算cosa\cosb\cosc
        (setq cos_a (sqrt (- 1 (* sin_a sin_a))))
        (setq cos_b (sqrt (- 1 (* sin_b sin_b))))
        (setq cos_c (sqrt (- 1 (* sin_c sin_c))))
        (setq j_a (/ (* 180.0 (atan (/ sin_a cos_a))) PI))
        (setq j_b (/ (* 180.0 (atan (/ sin_b cos_b))) PI))
        (setq j_c (/ (* 180.0 (atan (/ sin_c cos_c))) PI))
        ;;圆周率角度相等测试
        (if (not (equal (+ j_a j_b j_c) 180.0 0.000000001))
                (progn
                        (setq max_l (max l_a l_b l_c))
                        ;;l_a边大
                        (if (= max_l l_a)
                                (progn
                                        (setq cos_a (- 0.0 (sqrt (- 1 (* sin_a sin_a)))))
                                        (setq j_a  (/ (* 180.0 (atan (/ sin_a cos_a))) PI))
                                        (if (< j_a 0.0)
                                                (setq j_a (+ 180.0 j_a))
                                        )
                                )
                        )
                        ;;l_b边大               
                        (if (= max_l l_b)
                                (progn
                                        (setq cos_b (- 0.0 (sqrt (- 1 (* sin_b sin_b)))))
                                        (setq j_b  (/ (* 180.0 (atan (/ sin_b cos_b))) PI))
                                        (if (< j_b 0.0)
                                                (setq j_b (+ 180.0 j_b))
                                        )                                       
                                )
                        )
                        ;;l_c边大                       
                        (if (= max_l l_c)
                                (progn
                                        (setq cos_c (- 0.0 (sqrt (- 1 (* sin_c sin_c)))))
                                        (setq j_c  (/ (* 180.0 (atan (/ sin_c cos_c))) PI ))
                                        (if (< j_c 0.0)
                                                (setq j_c (+ 180.0 j_c))
                                        )                                       
                                )
                        )                       
                )
        )
        ;;a边对应的中线长
        (setq m_a_l (* 0.5 (sqrt (- (* 2.0 (+ (* l_b l_b) (* l_c l_c))) (* l_a l_a)))))
        ;;a边的中点
        (setq m_a_p (list
                                        (+ (/ (- (car p_b) (car p_c)) 2.0) (car p_c))                               
                                        (+ (/ (- (cadr p_b) (cadr p_c)) 2.0) (cadr p_c))
                                        (+ (/ (- (caddr p_b) (caddr p_c)) 2.0) (caddr p_c))                                       
                                )
        )                       
        (command "line" p_a p_b "")
        (command "line" p_b p_c "")       
        (command "line" p_a p_c "")       
        (command "line" p_a m_a_p "")
        (print j_a)
        (print j_b)       
        (print j_c)               
        (setq fanhui
                (list
                        ;;返回三点
                        (list "dian_a" p_a)
                        (list "dian_b" p_b)       
                        (list "dian_c" p_c)       
                        ;;返回三边
                        (list "bian_a" l_a)
                        (list "bian_b" l_b)
                        (list "bian_c" l_c)
                        ;;返回三个角的角度
                        (list "jiao_a" j_a)
                        (list "jiao_b" j_b)
                        (list "jiao_c" j_c)       
                        ;;返回三角形面积
                        (list "mianji" s_1)       
                )
        )
)