Boltzmann 函数曲线拟合的 Lisp 程序.doc

Boltzmann 函数曲线拟合的 Lisp 程序1 原程序 在计算机中，将以下原代码写入记事本 中并保 存文件名为“bzlm.lsp” (setq smx (lambda ( / k wi a1 a2 b1 b2 c1 c2 sx) (setq wi (mapcar (lambda ( x / ) (expt 2.718282 (/ (- x m3) m4) xi) k 0 a1 (apply + (mapcar (lambda ( y / w) (setq w (nth k wi) k (1+ k) (/ y (+ 1 w) yi) k 0 a2 (apply + (mapcar (lambda ( y / w) (setq w (nth k wi) k (1+ k) (/ (* y w) (+ 1 w) yi)b1 (apply + (mapcar (lambda ( w / ) (/ 1 (expt (+ 1 w) 2) wi) b2 (apply + (mapcar (lambda ( w / ) (/ w (expt (+ 1 w) 2) wi)c1 b2 c2 (apply + (mapcar (lambda ( w / ) (expt (/ w (+ 1 w) 2) wi) m1 (/ (- (* a1 c2) (* a2 c1) (- (* b1 c2) (* b2 c1) m2 (/ (- (* b1 a2) (* b2 a1) (- (* b1 c2) (* b2 c1) k 0 sx (apply + (mapcar (lambda ( x / y w) (setq w (nth k wi) y (nth k yi) k (1+ k) (expt (- y (+ (/ (- m1 m2) (+ 1 w) m2) 2) xi) ) (if (car s_min) (if ( sx (car s_min) (setq s_min (list sx m1 m2 m3 m4) nil) (setq s_min (list sx m1 m2 m3 m4) ) ) (setq mmc (lambda (range / m1 m2 s_min m3 m4 rm3 rm4 tm4 q3 q4) (setq rm3 (abs range) rm4 rm3 q3 rm3 q4 rm4 m3 0 m4 0) (repeat 5 (setq rm3 (+ m3 q3) m3 (- m3 q3) rm4 (+ m4 q4) m4 (- m4 q4) tm4 m4 q3 (* q3 0.1) q4 (* q4 0.1) (while (= m3 rm3) (while (= m4 1) (smx) (setq m4 (+ m4 q4) (setq m3 (+ m3 q3) m4 tm4) ) (setq m3 (nth 3 s_min) m4 (last s_min) ) s_min ) ) (setq cy (lambda ( / m1 m2 m3 m4) (if sc (progn (setq m1 (nth 1 sc) m2 (nth 2 sc) m3 (nth 3 sc) m4 (nth 4 sc) yc (+ (/ (- m1 m2) (+ 1 (expt 2.718282 (/ (- xc m3) m4) m2) )(set_tile cy (vl-princ-to-string yc) ) (setq cx (lambda ( / m1 m2 m3 m4 tm) (if sc (progn (setq m1 (nth 1 sc) m2 (nth 2 sc) m3 (nth 3 sc) m4 (nth 4 sc) (if (and ( (setq tm (- (/ (- m1 m2) (- yc m2) 1) 0) (progn (setq xc (+ m3 (* m4 (log tm) (set_tile cx (vl-princ-to-string xc) ) (alert 超出取值范围!) ) ) ) ) ) (setq draw (lambda ( / lx ly x1 x2 tem m1 m2 m3 m4 i) (setq lx (max (abs (setq x2 (apply max xi) (abs (setq x1 (apply min xi) ly (max (abs (apply max yi) (abs (apply min yi) m1 (nth 1 sc) m2 (nth 2 sc) m3 (nth 3 sc) m4 (nth 4 sc) ) (entmake (list (0 . LINE) (100 . AcDbEntity) (cons 10 (0 -10) (cons 11 (list 0 ly) (62 . 4) (entmake (list (0 . LINE) (100 . AcDbEntity) (cons 10 (list (* -1 lx) 0) (cons 11 (list lx 0) (62 . 4) (setq tem (0 . SPLINE) (100 . AcDbEntity) (100 . AcDbSpline) (70 . 8) (71 . 3) (72 . 13) (73 . 9) (74 . 7) (62 . 1) i 0 lx nil ) (while (= (length xi) 3) (= (length xi) (length yi) (progn (setq sc (mmc 100) i 0 pc (apply max (mapcar (lambda ( x / y m1 m2 m3 m4) (setq y (nth i yi) i (1+ i) m1 (nth 1 sc) m2 (nth 2 sc) m3 (nth 3 sc) m4 (nth 4 sc) (abs (- y (+ (/ (- m1 m2) (+ 1 (expt 2.718282 (/ (- x m3) m4) m2) ) xi) ) (set_tile ls (vl-princ-to-string pc) (set_tile lcs (vl-princ-to-string (nth 0 sc) (set_tile m1 (vl-princ-to-string (nth 1 sc) (set_tile m2 (vl-princ-to-string (nth 2 sc) (set_tile m3 (vl-princ-to-string (nth 3 sc) (set_tile m4 (vl-princ-to-string (nth 4 sc) ) (setq xi nil yi nil sc nil) ) ) (action_tile cx (progn (setq xc (read $value) yc nil) (acti