数学建模社区-数学中国
标题:
开放式计算软件OpenLu-工程计算助手
[打印本页]
作者:
forcal
时间:
2012-1-3 14:04
标题:
开放式计算软件OpenLu-工程计算助手
OpenLu是OpenFC的升级版本,OpenLu用Lu作为支持脚本。OpenLu力求消除冗余,简化计算。基本计算内容包括表达式计算、解非线性方程(组)、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
6 E$ H' C! y8 D# P1 [
OpenLu工作模式有三种,其即时编译计算模式特别适合于常用计算:在代码窗口写代码时,即时给出代码中的错误,若没有错误,将进行计算给出结果。
0 f* d$ h" I) y% k- N9 P. F' |
OpenLu不仅支持整数、实数、复数、三维向量等基本运算,还可以自定义常量、函数、命令菜单等。按需定制,简化计算是OpenLu的基本特点。
. F" `2 b8 _8 b1 F1 Y. s* A' V1 L
% [* q/ [2 M9 I/ M4 ~
本次升级,OpenLu简化了LuIMSL、LuOpt等扩展库,保留求方程(组)全部解、参数优化拟合等基本功能,函数用法也有诸多简化,以微分方程求解为例:
9 M5 {7 e8 W4 j' u4 M$ G
/ j! H' J/ i/ [
[例子] 设一阶微分方程组及初值为:
/ {( F0 P# u) F/ R+ @( J$ p6 ]
r'=2r-2rf, r(0)=1
$ _/ u9 ^, {! O
f'=-f+rf, f(0)=3
4 h% Q1 N7 \7 f/ k
计算t=1,2,...,10时的r、f的值。
- u N+ P# M3 I. p7 x: ~ k9 \" ? q
* B, U" J: \* s$ M
程序如下:
!!!using["IMSL","math"];
& s! P# f, |1 A d
f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
! ]" \3 i ?* H: b
ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];
复制代码
结果:
0. 1. 3.
1 U& n" d' Y% H8 S, o) E, `2 }
, M' m. J2 `5 g) c/ {* W8 X
1. 7.73453e-002 1.46445
; J5 B4 I6 g+ k2 _' _
, t4 G6 ^% j+ i4 C5 l. l
2. 8.49774e-002 0.577954
( ?- |- B- x% [# i1 @: W
# d* ~ W3 d x' s9 f# e4 [ O
3. 0.290891 0.249253
; ?5 H" ^. j/ W- o
0 g9 ?" r# n- E
4. 1.4466 0.187219
* L7 [- v; N W0 Q# H9 T- [
$ W$ Z# ^5 p) Z2 |/ N& K1 V9 z- D
5. 4.05146 1.43948
& g# o0 k7 f3 I: C8 Y
! t! `3 X5 w4 i( R
6. 0.175618 2.2586
, ]! D9 ~- M# q4 ^
/ D3 P; c2 Y2 O9 j; X6 @
7. 6.53112e-002 0.9088
5 a% b6 i& W3 V4 V4 v& U; e5 i
' v, R, `$ Z0 t
8. 0.147227 0.366718
6 s0 P- Q Y4 @ x+ e$ r w
% N9 \# z" ]2 V8 V- \+ L: w! U5 ?
9. 0.650596 0.187575
# W' P: C P, f3 b/ u: x7 ^# g
# k& H m( T* s
10. 3.14433 0.348821
复制代码
详细说明:
http://www.forcal.net/sysm/lu1/openlu_readme.htm
+ H0 C! |; _1 c/ P' n: F. \( L" {
下载试用:
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