数学建模社区-数学中国

标题: 开放式计算软件OpenLu-工程计算助手 [打印本页]

作者: forcal    时间: 2012-1-3 14:04
标题: 开放式计算软件OpenLu-工程计算助手
    OpenLu是OpenFC的升级版本,OpenLu用Lu作为支持脚本。OpenLu力求消除冗余,简化计算。基本计算内容包括表达式计算、解非线性方程(组)、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。; M# O$ n; L% j. Z( Z7 D; Q3 N
    OpenLu工作模式有三种,其即时编译计算模式特别适合于常用计算:在代码窗口写代码时,即时给出代码中的错误,若没有错误,将进行计算给出结果。( w) h4 ]. h: [# h
    OpenLu不仅支持整数、实数、复数、三维向量等基本运算,还可以自定义常量、函数、命令菜单等。按需定制,简化计算是OpenLu的基本特点。
7 m; F, ]8 {# x. V/ }3 ^- j6 |* C; ^8 z  C3 g" ?# O/ G
    本次升级,OpenLu简化了LuIMSL、LuOpt等扩展库,保留求方程(组)全部解、参数优化拟合等基本功能,函数用法也有诸多简化,以微分方程求解为例:' l: `* Q, @7 {/ m9 j3 X+ }

4 Y7 g+ q% R. o' U/ L2 f    [例子] 设一阶微分方程组及初值为:
- I- H% F" X# g& d        r'=2r-2rf, r(0)=1. ?# B* p; [3 E, Q; X* k
        f'=-f+rf,  f(0)=3
7 Z& u, x- R% |% w* W! J/ b    计算t=1,2,...,10时的r、f的值。# e% N1 q5 y. e' N7 g
$ W/ n1 {9 u0 e2 L6 z6 @7 a
    程序如下:
  1. !!!using["IMSL","math"];  W9 n  d. H" j+ W' U
  2. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义* a1 I! W4 E7 T- o% a# d
  3. ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];
复制代码
结果:
  1. 0.  1.            3.
    ! @% b7 p& c: O5 n
  2. * b, D$ l4 E* `$ K  e3 j
  3. 1.  7.73453e-002  1.46445' Z5 B" e2 }/ N/ `4 ~
  4. & b' ]$ C/ A, G" R
  5. 2.  8.49774e-002  0.577954/ T" R! V3 }% X4 v9 }! e# x" {
  6. # \7 a) R' x  X0 s2 Q& N0 d
  7. 3.  0.290891      0.249253+ J; D. Z, u) A3 ^2 g
  8. . L) y% [+ S+ I! p
  9. 4.  1.4466        0.187219  U: i/ e' ?5 I3 b, x
  10. $ c8 X5 `$ M  c6 U9 O
  11. 5.  4.05146       1.439481 L) Y- C  B0 L1 J3 `- ~

  12. + j. h* f7 c; m3 l, |4 {8 M3 r7 F
  13. 6.  0.175618      2.25866 x3 f0 B4 S$ R

  14. 6 y: p; `  |1 U9 F" U& Y
  15. 7.  6.53112e-002  0.9088) h" g5 F4 A2 K

  16. + z4 r; l$ Q- n
  17. 8.  0.147227      0.3667180 v( y6 V* V' Q

  18. ' j1 g* T* d) q: n% V
  19. 9.  0.650596      0.1875754 D4 K) r+ P) k- I; t

  20. 8 H5 H+ J4 l( r9 l) Y6 D
  21. 10. 3.14433       0.348821
复制代码
详细说明:http://www.forcal.net/sysm/lu1/openlu_readme.htm) s: u4 d/ a0 ]% w+ S3 ~! x. I
    下载试用:http://www.forcal.net/xiazai/lu1/openlu32.rar
作者: 孤寂冷逍遥    时间: 2012-1-5 16:05

作者: liupeng723911    时间: 2012-11-29 09:43
说的不错!




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5