数学建模社区-数学中国
标题:
开放式计算软件OpenLu-工程计算助手
[打印本页]
作者:
forcal
时间:
2012-1-3 14:04
标题:
开放式计算软件OpenLu-工程计算助手
OpenLu是OpenFC的升级版本,OpenLu用Lu作为支持脚本。OpenLu力求消除冗余,简化计算。基本计算内容包括表达式计算、解非线性方程(组)、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
) ~1 j! y" `& g" z% h4 \8 ?9 q H; I. [
OpenLu工作模式有三种,其即时编译计算模式特别适合于常用计算:在代码窗口写代码时,即时给出代码中的错误,若没有错误,将进行计算给出结果。
0 i- X: r. A: u7 z' U6 f, T" [
OpenLu不仅支持整数、实数、复数、三维向量等基本运算,还可以自定义常量、函数、命令菜单等。按需定制,简化计算是OpenLu的基本特点。
0 b0 F, I$ @' j
9 {& d& S1 r/ }0 M5 W6 X
本次升级,OpenLu简化了LuIMSL、LuOpt等扩展库,保留求方程(组)全部解、参数优化拟合等基本功能,函数用法也有诸多简化,以微分方程求解为例:
% t9 S! {% o1 m0 W) g: q
: e( k. I/ n: d8 S
[例子] 设一阶微分方程组及初值为:
& v# i; ^0 c# _4 o2 w
r'=2r-2rf, r(0)=1
! W% Y$ Z: D& M9 ]! N. E7 b
f'=-f+rf, f(0)=3
' B& @: v" K( b; g) M2 I5 p
计算t=1,2,...,10时的r、f的值。
; L: n8 N0 A, n. u6 L4 i4 q( B
3 }: t0 q( `7 C' T- M
程序如下:
!!!using["IMSL","math"];
) ?( O% F7 d! p! ^2 @
f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
0 F; M! h* s1 ?8 y; Z2 Y
ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];
复制代码
结果:
0. 1. 3.
! v( e1 t; i1 N3 | T/ q
/ C: M" y* n4 |! C
1. 7.73453e-002 1.46445
* y$ i& n) T1 }7 N. k9 G1 T7 g
' s2 }( E' v6 |$ O" _. ^
2. 8.49774e-002 0.577954
8 x/ ?3 ~3 o9 f$ p- n
9 S0 @' z6 U# E5 s. Y5 T d A
3. 0.290891 0.249253
8 V7 g6 L! A$ I i+ G, q( @
( T2 _1 ?. X( a
4. 1.4466 0.187219
3 f4 w. s: z" [; o; b- I' [ P/ R. T
3 @# z( p: F. Z7 Z- `5 w7 p3 r
5. 4.05146 1.43948
3 V4 ~9 e' g1 v S* V
9 T& M/ A+ t1 h' T) D" D3 m) {+ |
6. 0.175618 2.2586
5 e) n: s* w( a E X
# K. H: |* C9 _7 _( c
7. 6.53112e-002 0.9088
6 i- e; S( }% V# M2 o! J8 c
6 b7 Y2 f$ R* }. V( @8 ^9 ~
8. 0.147227 0.366718
9 J4 ]! l; E9 Y: z0 m+ w* p. P4 J- D
0 M( ~& s8 W6 W4 O7 u
9. 0.650596 0.187575
# L" ~9 c' f- v" i3 ~% k
' F- v2 V# B* D! V" m6 V6 C' n, ^
10. 3.14433 0.348821
复制代码
详细说明:
http://www.forcal.net/sysm/lu1/openlu_readme.htm
3 }0 a c( X% f. ^
下载试用:
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