画出如图的蛇形LINE图形

aimisiyou

newer 发表于 2016-10-23 14:10

此题只是求相关点坐标问题就能解决，用递归方法反而复杂了。递归的思路如下：1、初始值f(pt,0)=(pt3 pt2 pt1)
2、f(pt ,k)=F(f(polar pt (/ pi 2),(- k 1)  )),其中F运算为根据上一次运算的结果进行对称、平移运算得出此次三点（pt3'   pt2'   pt1'）
3、(setq i 0 ptlst nil )(while (< i n )  (setq ptlst (cons (f pt i) ptlst))   (setq i (+ i 1)))
4、对点集（reverse   ptlst）画多段线   

aimisiyou

(defun c:tt ()
  (defun disfun (ptc r n k)
    (setq h (/ (* 2.0 r) n))
    (setq dd (- (sqrt (- (* r r) (* (- r (* (+ k 1) h)) (- r (* (+ k 1) h)))))
                (sqrt (- (* r r) (* (- r (* k h)) (- r (* k h)))))
             )
    )
  )
  (defun f (ptc r n d k)
    (if (= k 0)
      (setq ptlst (list (list (+ (car ptc) (disfun ptc r n 0)) (+
                                                                  (cadr ptc)
                                                                  r
                                                                  (/ r -0.5
                                                                     n
                                                                  )
                                                               )
                        ) (list (+ (car ptc) (disfun ptc r n 0) d)
                                (+ (cadr ptc) r (/ r -1.0 n))
                          ) (list (+ (car ptc) (disfun ptc r n 0))
                                  (+ (cadr ptc) r)
                            )
                  )
      )
      (progn
        (setq ptlst (mapcar
                      '(lambda (x)
                         (list (+ (car x) (* (if (= (rem k 2) 0)
                                               1
                                               -1
                                             ) (if (and
                                                     (= (rem n 2) 0)
                                                     (= k (/ n 2))
                                                   )
                                                 0
                                                 (disfun ptc r n
                                                         (if (<= k
                                                                 (/ n 2)
                                                             )
                                                           k
                                                           (- k 1)
                                                         )
                                                 )
                                               )
                                          )
                               ) (- (cadr x) (/ r 0.5 n))
                         )
                       )
                      (mapcar
                        '(lambda (y)
                           (list (- (* 2 (car ptc)) (car y)) (cadr y))
                         )
                        (f ptc r n d (- k 1))
                      )
                    )
        )
      )
    )
  )
  (setq ptc (getpoint "\n圓心：")
        r (getreal "\n半徑：")
        n (getint "\n等分數:")
        d (getreal "\n尖角長度:")
  )
  (command "circle" ptc r)
  (setq i 0
        ptlist (list (polar ptc (* pi 0.75) (* r (sqrt 2))))
  )
  (while (< i n)
    (setq ptlist (append
                   (f ptc r n d i)
                   ptlist
                 )
    )
    (setq i (+ i 1))
  )
  (if (= (rem n 2) 0)
    (setq ptlist (cons (polar ptc (* pi -0.25) (* r (sqrt 2))) ptlist))
    (setq ptlist (cons (polar ptc (* pi 1.25) (* r (sqrt 2))) ptlist))
  )
  (entmake (append
             (list '(0 . "LWPOLYLINE") '(100 . "AcDbEntity") '(100 . "AcDbPolyline")
                   (cons 90 (length ptlist))
             )
             (mapcar
               '(lambda (pt)
                  (cons 10 pt)
                )
               ptlist
             )
           )
  )
)

雪山飞狐(lzh)

总觉得lisp括号看着头疼...

1.         [CommandMethod("TTTT")]
   
2.         public static void Test()
   
3.         {
   
4. 
   
5.             var doc = AcadApp.DocumentManager.MdiActiveDocument;
   
6.             var ed = doc.Editor;
   
7.             var db = doc.Database;
   
8.             var resCenter = ed.GetPoint("\n请输入圆心:");
   
9.             if (resCenter.Status != PromptStatus.OK) return;
   
10.             var resRadius = ed.GetDouble("\n请输入半径:");
    
11.             if (resRadius.Status != PromptStatus.OK) return;
    
12.             var resNum = ed.GetInteger("\n请输入等分数:");
    
13.             if (resNum.Status != PromptStatus.OK) return;
    
14.             var resLen = ed.GetDouble("\n请输入尖角水平距离:");
    
15.             if (resLen.Status != PromptStatus.OK) return;
    
16. 
    
17.             var radius = resRadius.Value;
    
18.             var diameter = radius * 2;
    
19.             var num = resNum.Value;
    
20.             var len = resLen.Value;
    
21. 
    
22.             Polyline pline = new Polyline();
    
23.             var x0 = -radius;
    
24.             var y0 = radius;
    
25.             var i = 0;
    
26.             pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);
    
27.             var setp = diameter / num;
    
28.             var y = radius - setp;
    
29.             var sign = 1;
    
30.             while (y >= -radius)
    
31.             {
    
32.                 var x = sign * Math.Sqrt(radius * radius - y * y);
    
33.                 x0 = y >= 0 ? x : -x0;
    
34.                 pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);
    
35.                 pline.AddVertexAt(i++, new Point2d(x0 + sign * len, y0 - setp / 2), 0, 0, 0);
    
36.                 y0 -= setp;
    
37.                 pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);
    
38.                 if (y < 0) x0 = x;
    
39.                 sign = -sign;
    
40.                 y -= setp;
    
41.             }
    
42.             x0 = -radius * sign;
    
43.             pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);
    
44. 
    
45.             using (var tr = db.TransactionManager.StartTransaction())
    
46.             {
    
47.                 var center = resCenter.Value;
    
48.                 var btr = (BlockTableRecord)tr.GetObject(db.CurrentSpaceId, OpenMode.ForWrite);
    
49.                 var cir = new Circle(center, Vector3d.ZAxis, radius);
    
50.                 pline.TransformBy(Matrix3d.Displacement(center - Point3d.Origin));
    
51.                 btr.AppendEntity(cir);
    
52.                 tr.AddNewlyCreatedDBObject(cir, true);
    
53.                 btr.AppendEntity(pline);
    
54.                 tr.AddNewlyCreatedDBObject(pline, true);
    
55.                 tr.Commit();
    
56.             }
    
57.         }

        [CommandMethod("TTTT")]<br />         public static void Test()<br />         {<br /> <br />             var doc = AcadApp.DocumentManager.MdiActiveDocument;<br />             var ed = doc.Editor;<br />             var db = doc.Database;<br />             var resCenter = ed.GetPoint("\n请输入圆心:");<br />             if (resCenter.Status != PromptStatus.OK) return;<br />             var resRadius = ed.GetDouble("\n请输入半径:");<br />             if (resRadius.Status != PromptStatus.OK) return;<br />             var resNum = ed.GetInteger("\n请输入等分数:");<br />             if (resNum.Status != PromptStatus.OK) return;<br />             var resLen = ed.GetDouble("\n请输入尖角水平距离:");<br />             if (resLen.Status != PromptStatus.OK) return;<br /> <br />             var radius = resRadius.Value;<br />             var diameter = radius * 2;<br />             var num = resNum.Value;<br />             var len = resLen.Value;<br /> <br />             Polyline pline = new Polyline();<br />             var x0 = -radius;<br />             var y0 = radius;<br />             var i = 0;<br />             pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);<br />             var setp = diameter / num;<br />             var y = radius - setp;<br />             var sign = 1;<br />             while (y >= -radius)<br />             {<br />                 var x = sign * Math.Sqrt(radius * radius - y * y);<br />                 x0 = y >= 0 ? x : -x0;<br />                 pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);<br />                 pline.AddVertexAt(i++, new Point2d(x0 + sign * len, y0 - setp / 2), 0, 0, 0);<br />                 y0 -= setp;<br />                 pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);<br />                 if (y < 0) x0 = x;<br />                 sign = -sign;<br />                 y -= setp;<br />             }<br />             x0 = -radius * sign;<br />             pline.AddVertexAt(i++, new Point2d(x0, y0), 0, 0, 0);<br /> <br />             using (var tr = db.TransactionManager.StartTransaction())<br />             {<br />                 var center = resCenter.Value;<br />                 var btr = (BlockTableRecord)tr.GetObject(db.CurrentSpaceId, OpenMode.ForWrite);<br />                 var cir = new Circle(center, Vector3d.ZAxis, radius);<br />                 pline.TransformBy(Matrix3d.Displacement(center - Point3d.Origin));<br />                 btr.AppendEntity(cir);<br />                 tr.AddNewlyCreatedDBObject(cir, true);<br />                 btr.AppendEntity(pline);<br />                 tr.AddNewlyCreatedDBObject(pline, true);<br />                 tr.Commit();<br />             }<br />         }

aimisiyou

雪山飞狐(lzh) 发表于 2016-10-31 00:10

武学的最高境界是如何将不起眼的落叶变成夺命利器！

雪山飞狐(lzh)

更高境界是手中无剑 呵呵 开个玩笑 不用为语言的好坏争论 只有算法是精髓

aimisiyou

雪山飞狐(lzh) 发表于 2016-10-31 09:12

赞同！万变不离其宗！