[原创]：在CAD内解方程，一个极为实用的程序

highflybird

刚刚学了点LISP 程序，就要卖弄了:
这是一个解方程的LISP程序，很简短。

用法：１解压文件得到solve.lsp程序。
２，在ＣＡＤ中用AP命令加载此lisp程序。
３，在ＣＡＤ中输入solve按照提示输入表达式和初始值。（你可以先估计一个解作为初始值）。
４，最后出现的数字为解。

另外注明：
１，方程只接受x为一元变量，多元方程还在研究中。
２，方程接受ＣＡＤ中认可的表达式，如x^5+x^3-3*sin(x)+1；（如不清楚的请参考ＣＡＬ命令帮助中的语法和函数）注意了，方程只需写成多项式，无须写成f(x)＝0式。
３，方程解题精度可达１６位，只是表达实数根。
４，程序并没有检错判断，有可能无解现象出现，或其他错误出现。也没有用到对话框。
我把源代码发布上来，请各位大侠帮忙斧正，如要转载程序到他处，为尊重本人劳动成果见，请注明来源与原作者。

huyuan1984

呵呵
我也是刚学的

花清扬

请把你的LISP文件源文件直接贴出来吧，这样的小文件用这种方法更方便，我下载你的rar文件试了好几次都没有成功。

highflybird

(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)
  )

花清扬

我测试过了，就以你第一贴中的函数式x^5+x^3-3*sin(x)+1测的，得到的解为“－0.8471650345690628”，但有以下几个问题：

1、要正确运行这个LSP程序，还需要在运行前手工加载一次CAL命令，否则提示：“; 错误: no function definition: CAL”。因此需要在程序里添加语句以自动调用CAL功能。
2、就这个函数，输入的初始值为0.6时，出现提示：“; 错误: 除数为零”。这种错误信息因该很容易被程序捕捉到，从而采用出错处理方法自动处理，防止程序异常终止。
3、程序的求解精度如能自定义则更好。
4、有的函数式有多个解，希望能自定义求解的区间，即由外部输入初始值的最小值和最大值；如果给定区间内无解，希望能给出提示，并由用户选择结束命令或重新输入。
5、如果程序能实现以上的基本功能了，建议再增加并完善对话框功能；至于多元方程求解功能，它属于比较麻烦的数值计算问题，等“1.0”版出来后再升级吧。

highflybird

谢谢楼上的指点：
现在已经更改一个问题，就是不需要手动加载CAL命令了，另外减少了“除数为零”的错误情况，但没有完全避免，我想可能是这个函数在零点导数为零。

水炮眼大侠

做这个有用吗?要做就做点与工作相关的.