|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
& r& V+ D8 U1 w: s g$ a ]- @
1 线性方程组求解
; U9 [* J/ a: i I左除法
7 L; _) A! E! s3 s9 l/ u, KAx=b====>x=A\b %
% T; F$ j( N4 w矩阵分解1 h- [9 U+ N* R. z v
8 C, t) Y3 m# ^5 s( D3 C3 e
2 非线性方程组求解5 }2 P8 l! h% L: m
fzero(‘fname’,x0,tol,trace) %x0为搜索起点,函数会求离x0最近那个根,tol控制相对精度,trace指代迭代信息,1时显示,0不显示
3 B) G- I! y0 e! K1 z( q/ k6 J g! rX=fsolve(‘fname’,x0,option) %option为最优化工具箱的选项设定。最优化工具箱提供了20多个选项,用户可以使用optimset命令将它们显示出来。如果想改变其中某个选项,则可以调用optimset()函数来完成。例如,Display选项决定函数调用时中间结果的显示方式,其中‘off’为不显示,‘iter’表示每步都显示,‘final’只显示最终结果。optimset(‘Display’,‘off’)将设定Display选项为‘off’. s1 E Y' [# L) I2 X
! k; N5 c8 |1 y7 X& U9 H' E* K3 常微分方程# s/ ?! Q7 r6 U# ~' C. S
[t,y]=ode23(‘fname’,tspan,y0)0 ?3 N( o% x! {2 }; {) |
[t,y]=ode45(‘fname’,tspan,y0) %其中fname是定义f(t,y)的函数文件名,该函数文件必须返回一个列向量。tspan形式为[t0,tf],表示求解区间。y0是初始状态列向量。t和y分别给出时间向量和相应的状态向量, \+ C9 K% [' d4 h7 `& @$ f7 H
; g3 V5 p; [( h5 V7 q) d) n
4 极值
B5 o: X; w# s3 f1 ?fmin(‘fname’,x1,x2) %x1,x2范围
& }% [% g S7 }4 \4 ?fmins(‘fname’,x0) %x0为初始值向量
6 B7 a9 H5 c2 k+ e4 z5 r
# | r! g _) Q* L+ C |
|