在CAD内解方程，一个极为实用的程序

highflybir · 发表于 2006-9-8 08:22:00

本帖最后由作者于 2006-9-8 12:04:32 编辑

刚刚学了点LISP 程序，就要卖弄了:

这是一个解方程的LISP程序，很简短。

用法：１解压文件得到solve.lsp程序。

２，在ＣＡＤ中用AP命令加载此lisp程序。

３，在ＣＡＤ中输入solve按照提示输入表达式和初始值。（你可以先估计一个解作为初始值）。

４，最后出现的数字为解。

另外注明：

１，方程只接受x为一元变量，多元方程还在研究中。

２，方程接受ＣＡＤ中认可的表达式，如x^5+x^3-3*sin(x)+1；（如不清楚的请参考ＣＡＬ命令帮助中的语法和函数）注意了，方程只需写成多项式，无须写成f(x)＝0式。

３，方程解题精度可达１６位，只是表达实数根。

４，程序并没有检错判断，有可能无解现象出现，或其他错误出现。也没有用到对话框。

我把源代码发布上来，请各位大侠帮忙斧正，如要转载程序到他处，为尊重本人劳动成果见，请注明来源与原作者。

highflybir · 发表于 2006-9-8 12:44:00

这里是源代码：

(defun func (/ dx dy y1 y2 y3 y4)
   (setq y1 (cal p))
   (setq dx 0.00000000000001)
   (setq x (+ x dx))
   (setq y2 (cal p))
   (setq dy (- y2 y1))
   (setq y3 (/ dy dx))
   (setq x (- x dx))
   (setq y4 (- x (/ y1 y3)))
)
(defun C:solve (/ p x y)
   (arxload "geomcal.arx")
   (setq p (getstring "\n请输入函数表达式: "))
   (setq x (getreal "\n请输入初始值: "))
   (setq y (func))
(while (> (abs (- y x)) 0.000000000000001)
    (setq x y)
           (setq y (func))
   )
   (rtos y 2 20)
)

ahlzl · 发表于 2006-9-8 12:50:00

要力争实现解方程和画图的“无缝”连接。

Joseflin · 发表于 2006-9-9 08:20:00

希望能解联立方程式的LISP早日出现，如此几何解题法当会更加为大众接受。

quark · 发表于 2006-9-9 19:52:00

Paula · 发表于 2006-9-21 16:51:00

hehe 谢谢了!!

账号		自动登录	找回密码
密码			注册

在CAD内解方程，一个极为实用的程序

本帖子中包含更多资源