通过程序计算出来的数值和手动计算的数值不一样，是不是我的程序写错了？

gprspda

Andyhon 发表于 2011-10-22 22:04
(setq α (getangle "\n输入转弯坡度角: "))
==>
  (setq α (getReal "\n输入转弯坡度角: "))

打扰你这么长时间真是不好意思，非常感谢你的帮助

Andyhon

比较版

2. 
   
3. 
   
4. (defun c:flx (/ p p0 p1 p2 p3 p4 p5 p6 v c h w r1 r2 b1 b2 b3 α e t1)
   
5.   (setq v (getReal "\n输入指示空速: "))
   
6.   (setq c (getReal "\n输入与ISA的差值: "))
   
7.   (setq h (getReal "\n输入高度（米）: "))
   
8.   (setq w (getReal "\n输入风速: "))
   
9.   (setq p (getpoint "\n确定起始点: "))
   
10.   (setq α (getAngle "\n输入转弯坡度角: "))
    
11.   
    
12.   (setq t1 (/ (* v 171233 (expt (- (+ 288 c) (* 0.006496 h)) 0.5))
    
13.             (expt (- 288 (* 0.006496 h)) 2.628)
    
14.          )
    
15.   )
    
16. 
    
17.   ;;; 公式是r2=3431*tan(25°）/(π t1).
    
18.   ;;; (setq r2 (/ (* 3431 (tan α))  pi t1))
    
19.   (setq r2 (/ (* 3431 (/ (sin α) (cos  α))) Pi t1))
    
20.   
    
21.   (if (> r2 3)
    
22.     (setq r2 3)
    
23.   )
    
24.   
    
25.   (setq r1 (/ t1 (* 20 pi r2)))
    
26.   (setq p0 (polar p 0 r1))
    
27.   (setq e (/ w (* 40 r2)))
    
28.   (setq p1 (polar p0 (/ pi 2) e))
    
29.   (setq p2 (polar p0 0 e))
    
30.   (setq p3 (polar p0 (/ (- pi) 2) e))
    
31.   (setq b1 (expt (+ (* r1 r1) (* e e)) 0.5))
    
32.   (setq b2 (+ r1 e))
    
33.   (setq b3 (+ r1 (* 2 e)))
    
34.   (setq p4 (polar p1 (/ pi 2) b1))
    
35.   (setq p5 (polar p2 0 b2))
    
36.   (setq p6 (polar p3 (/ (- pi) 2) b3))
    
37.   
    
38.   (setVar "OsMode" 0)
    
39. 
    
40. 
    
41.   ;|
    
42.   (command "arc" "c" p1 p4 p)
    
43.   (command "arc" "c" p2 p5 p4)
    
44.   (command "arc" "c" p3 p6 p5)
    
45.   |;
    
46. 
    
47.   (command "arc" "c" p1 p4 p)
    
48.   (setq e1 (entlast))
    
49.   (command "arc" "c" p2 p5 p4)
    
50.   (setq e2 (entlast))
    
51.   (command "arc" "c" p3 p6 p5)
    
52.   (setq e3 (entlast))
    
53.   
    
54.   (command "Pedit" "m" e1 e2 e3 "" "j" "" "")
    
55. 
    
56.   (setq fn (open (getfiled "将文件保存为：" "gc.txt" "txt" 1) "w"))
    
57.   (write-line (strcat "t1=" (rtos t1)) fn)
    
58.   (write-line (strcat "r2=" (rtos r2)) fn)
    
59.   (write-line (strcat "r1=" (rtos r1)) fn)
    
60.   (write-line (strcat "b1=" (rtos b1)) fn)
    
61.   (write-line (strcat "b2=" (rtos b2)) fn)
    
62.   (write-line (strcat "b3=" (rtos b3)) fn)
    
63.   (close fn)
    
64.   (command "_zoom" "0.5X")
    
65.   (princ)
    
66. )
    

<br /> <br /> (defun c:flx (/ p p0 p1 p2 p3 p4 p5 p6 v c h w r1 r2 b1 b2 b3 α e t1)<br />   (setq v (getReal "\n输入指示空速: "))<br />   (setq c (getReal "\n输入与ISA的差值: "))<br />   (setq h (getReal "\n输入高度（米）: "))<br />   (setq w (getReal "\n输入风速: "))<br />   (setq p (getpoint "\n确定起始点: "))<br />   (setq α (getAngle "\n输入转弯坡度角: "))<br />   <br />   (setq t1 (/ (* v 171233 (expt (- (+ 288 c) (* 0.006496 h)) 0.5))<br />             (expt (- 288 (* 0.006496 h)) 2.628)<br />          )<br />   )<br /> <br />   ;;; 公式是r2=3431*tan(25°）/(π t1).<br />   ;;; (setq r2 (/ (* 3431 (tan α))  pi t1))<br />   (setq r2 (/ (* 3431 (/ (sin α) (cos  α))) Pi t1))<br />   <br />   (if (> r2 3)<br />     (setq r2 3)<br />   )<br />   <br />   (setq r1 (/ t1 (* 20 pi r2)))<br />   (setq p0 (polar p 0 r1))<br />   (setq e (/ w (* 40 r2)))<br />   (setq p1 (polar p0 (/ pi 2) e))<br />   (setq p2 (polar p0 0 e))<br />   (setq p3 (polar p0 (/ (- pi) 2) e))<br />   (setq b1 (expt (+ (* r1 r1) (* e e)) 0.5))<br />   (setq b2 (+ r1 e))<br />   (setq b3 (+ r1 (* 2 e)))<br />   (setq p4 (polar p1 (/ pi 2) b1))<br />   (setq p5 (polar p2 0 b2))<br />   (setq p6 (polar p3 (/ (- pi) 2) b3))<br />   <br />   (setVar "OsMode" 0)<br /> <br /> <br />   ;|<br />   (command "arc" "c" p1 p4 p)<br />   (command "arc" "c" p2 p5 p4)<br />   (command "arc" "c" p3 p6 p5)<br />   |;<br /> <br />   (command "arc" "c" p1 p4 p)<br />   (setq e1 (entlast))<br />   (command "arc" "c" p2 p5 p4)<br />   (setq e2 (entlast))<br />   (command "arc" "c" p3 p6 p5)<br />   (setq e3 (entlast))<br />   <br />   (command "Pedit" "m" e1 e2 e3 "" "j" "" "")<br /> <br />   (setq fn (open (getfiled "将文件保存为：" "gc.txt" "txt" 1) "w"))<br />   (write-line (strcat "t1=" (rtos t1)) fn)<br />   (write-line (strcat "r2=" (rtos r2)) fn)<br />   (write-line (strcat "r1=" (rtos r1)) fn)<br />   (write-line (strcat "b1=" (rtos b1)) fn)<br />   (write-line (strcat "b2=" (rtos b2)) fn)<br />   (write-line (strcat "b3=" (rtos b3)) fn)<br />   (close fn)<br />   (command "_zoom" "0.5X")<br />   (princ)<br /> )<br />

gprspda

Andyhon 发表于 2011-10-22 22:22
比较版

谢谢，两个版本运行的结果都是一样的...
有没有一个程序，可以让我画出来的弧连成一个整体？现在画出的弧是分三段的

Andyhon

  (command "arc" "c" p1 p4 p)
  (command "arc" "c" p2 p5 p4)
  (command "arc" "c" p3 p6 p5)


===>
  (command "arc" "c" p1 p4 p)
  (setq e1 (entlast))
  (command "arc" "c" p2 p5 p4)
  (setq e2 (entlast))
  (command "arc" "c" p3 p6 p5)
  (setq e3 (entlast))
  
  (command "Pedit" "m" e1 e2 e3 "" "j" "" "")

gprspda

Andyhon 发表于 2011-10-23 09:57
(command "arc" "c" p1 p4 p)
  (command "arc" "c" p2 p5 p4)
  (command "arc" "c" p3 p6 p5)

程序运行之后出来的选项应该选哪个？我挨个试了一下，好像都不可以啊

Andyhon

修订了 32#

gprspda

Andyhon 发表于 2011-10-23 10:17
修订了 32#

太感谢你了，谢谢，终于可以交差了...

gprspda

Andyhon 发表于 2011-10-23 10:17
修订了 32#

不好意思，还要打扰一下，程序运行之后为什么不能把结果显示在gc.txt这个文件里了？以前还是可以的

Andyhon

这会儿重测正常,
建议您重新启动后再调试

gprspda

Andyhon 发表于 2011-10-23 15:24
这会儿重测正常,
建议您重新启动后再调试

我重新启动了，按照正常步骤，加载，运行，图形结果正确，但是参数的计算结果没有显示..是不是有什么地方冲突？