四个圆包围矩形且与矩形边相切

荒野孤行

想起刚开始学AutoCAD的时候，都是画的一些有数学关系的图，那时有很多图形都没画出来。现在用lisp写出来，感觉so easy！

源码如下：

1. ;;;*****矩形4圆相切 程序开始*****
   
2. (defun c:t1 ()
   
3.   (setvar "cmdecho" 0)
   
4.   (setvar "osmode" 15359)
   
5.   (princ "\n功能：通过矩形绘制出与其四边相切，且四圆两两相切的圆.")
   
6.   (initget 7)
   
7.   (setq l_len (getdist "\n请输入矩形的长度："))
   
8.   (initget 7)
   
9.   (setq l_wid (getdist "\n请输入矩形的宽度："))
   
10.   (initget 3)
    
11.   (setq pt1 (getpoint "\n请指定插入点："))
    
12.   (setvar "osmode" 0)
    
13.   (setq L1 (* l_wid 0.5))
    
14.   (setq L2 (* l_len 0.5))
    
15.   (setq a (* 2.0 L1))
    
16.   (setq b (+ (* L2 L2) (* L1 L1 -1.0) (* 4.0 L1 L2)))
    
17.   (setq c (* 2.0 L1 L1 L2))
    
18.   (setq r2 (/ (+ b (sqrt (+ (* b b) (* 4.0 a c)))) (* 2.0 a) 1.0))
    
19.   (setq r1 (+ (/ (* L1 L1) r2 1.0) (* 2.0 L1)))
    
20.   (setq pt2 (polar pt1 0 l_len))
    
21.   (setq pt12 (polar pt1 0 (* l_len 0.5)))
    
22.   (setq pt_r2 (polar pt12 (angtof "90") r2))
    
23.   (setq pt3 (polar pt2 (angtof "-90") l_wid))
    
24.   (setq pt4 (polar pt1 (angtof "-90") l_wid))
    
25.   (setq pt34 (polar pt1 (angtof "-90") (* l_wid 0.5)))
    
26.   (setq pt_r1 (polar pt34 (angtof "180") r1))
    
27.   (command "color" 1)
    
28.   (command "PLINE" pt1 pt2 pt3 pt4 "c")
    
29.   (command "color" 3)
    
30.   (command "CIRCLE" pt_r2 r2)
    
31.   (command "MIRROR" (entlast) "" pt34 pt_r1 "")
    
32.   (command "color" 4)
    
33.   (command "CIRCLE" pt_r1 r1)
    
34.   (command "MIRROR" (entlast) "" pt12 pt_r2 "")
    
35.   (command "color" "Bylayer")
    
36.   (setvar "osmode" 15359)
    
37.   (princ)
    
38. )
    
39. ;;;*****矩形4圆相切 程序结束*****
    

;;;*****矩形4圆相切 程序开始*****<br /> (defun c:t1 ()<br />   (setvar "cmdecho" 0)<br />   (setvar "osmode" 15359)<br />   (princ "\n功能：通过矩形绘制出与其四边相切，且四圆两两相切的圆.")<br />   (initget 7)<br />   (setq l_len (getdist "\n请输入矩形的长度："))<br />   (initget 7)<br />   (setq l_wid (getdist "\n请输入矩形的宽度："))<br />   (initget 3)<br />   (setq pt1 (getpoint "\n请指定插入点："))<br />   (setvar "osmode" 0)<br />   (setq L1 (* l_wid 0.5))<br />   (setq L2 (* l_len 0.5))<br />   (setq a (* 2.0 L1))<br />   (setq b (+ (* L2 L2) (* L1 L1 -1.0) (* 4.0 L1 L2)))<br />   (setq c (* 2.0 L1 L1 L2))<br />   (setq r2 (/ (+ b (sqrt (+ (* b b) (* 4.0 a c)))) (* 2.0 a) 1.0))<br />   (setq r1 (+ (/ (* L1 L1) r2 1.0) (* 2.0 L1)))<br />   (setq pt2 (polar pt1 0 l_len))<br />   (setq pt12 (polar pt1 0 (* l_len 0.5)))<br />   (setq pt_r2 (polar pt12 (angtof "90") r2))<br />   (setq pt3 (polar pt2 (angtof "-90") l_wid))<br />   (setq pt4 (polar pt1 (angtof "-90") l_wid))<br />   (setq pt34 (polar pt1 (angtof "-90") (* l_wid 0.5)))<br />   (setq pt_r1 (polar pt34 (angtof "180") r1))<br />   (command "color" 1)<br />   (command "PLINE" pt1 pt2 pt3 pt4 "c")<br />   (command "color" 3)<br />   (command "CIRCLE" pt_r2 r2)<br />   (command "MIRROR" (entlast) "" pt34 pt_r1 "")<br />   (command "color" 4)<br />   (command "CIRCLE" pt_r1 r1)<br />   (command "MIRROR" (entlast) "" pt12 pt_r2 "")<br />   (command "color" "Bylayer")<br />   (setvar "osmode" 15359)<br />   (princ)<br /> )<br /> ;;;*****矩形4圆相切 程序结束*****<br />

有空我就来

学习学习！谢谢

sfjlx

学习学习！谢谢

水仙的错

支持下,希望有更多好的程序

hbgsw

看看源码情况。

czcxxx

谢谢分享！！！

并州驴友

楼主的好代码，受教了，回去好好研究研究！

lrd1861

看看藏着啥  楼主是不是藏错了

winerfjy

看隐藏着什么

東洺Nile

看看代碼 學習 研究