|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
) D% i1 y: I) q/ ]
1 线性方程组求解
( `) n O" B+ n% g% ?左除法& F1 c y. }6 c2 U i. A9 t
Ax=b====>x=A\b %5 e; v6 w/ a- r) o, \/ X: |
矩阵分解
6 h' s3 d$ Y$ l6 A9 K1 s- B
8 C/ z6 n4 d4 [* n( \0 }2 非线性方程组求解
2 l9 C6 P o7 B+ p1 R+ D6 \- ]fzero(‘fname’,x0,tol,trace) %x0为搜索起点,函数会求离x0最近那个根,tol控制相对精度,trace指代迭代信息,1时显示,0不显示+ f* |: e* M& M1 ?+ C B$ g* u; G
X=fsolve(‘fname’,x0,option) %option为最优化工具箱的选项设定。最优化工具箱提供了20多个选项,用户可以使用optimset命令将它们显示出来。如果想改变其中某个选项,则可以调用optimset()函数来完成。例如,Display选项决定函数调用时中间结果的显示方式,其中‘off’为不显示,‘iter’表示每步都显示,‘final’只显示最终结果。optimset(‘Display’,‘off’)将设定Display选项为‘off’0 P. e4 ]* q, H2 X1 q1 [* r
8 h* ~2 H2 C! q* K; d& y3 常微分方程3 K' I- @- Q z! L4 y4 M
[t,y]=ode23(‘fname’,tspan,y0)
7 t& d5 ~$ N$ U/ j S8 W[t,y]=ode45(‘fname’,tspan,y0) %其中fname是定义f(t,y)的函数文件名,该函数文件必须返回一个列向量。tspan形式为[t0,tf],表示求解区间。y0是初始状态列向量。t和y分别给出时间向量和相应的状态向量
* l( r. t5 h' X& {* [7 |( d, @2 Y+ ^
0 t1 g( k3 Z: \4 u' c7 u$ o4 极值
& p6 O: b( x+ D, B( G! w/ sfmin(‘fname’,x1,x2) %x1,x2范围
- F( A$ A G* {" X$ ]6 C; Efmins(‘fname’,x0) %x0为初始值向量/ c0 X' T0 @3 [% F2 w8 ]
4 F3 [, O0 M% J, q. X* B( H. D |
|