|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
/ L. U- J- f! u4 H$ O
1 线性方程组求解
. E; C: j! [) V3 b, ]$ W左除法
0 m5 b3 ` o: o' ZAx=b====>x=A\b %. A6 }$ O9 C, u: [( V2 }2 N
矩阵分解5 i2 Y. G9 [% V5 X
6 b; E C, T$ W. l2 非线性方程组求解
) @& V9 K. u0 X% j6 lfzero(‘fname’,x0,tol,trace) %x0为搜索起点,函数会求离x0最近那个根,tol控制相对精度,trace指代迭代信息,1时显示,0不显示
/ g2 y8 A( a% s1 mX=fsolve(‘fname’,x0,option) %option为最优化工具箱的选项设定。最优化工具箱提供了20多个选项,用户可以使用optimset命令将它们显示出来。如果想改变其中某个选项,则可以调用optimset()函数来完成。例如,Display选项决定函数调用时中间结果的显示方式,其中‘off’为不显示,‘iter’表示每步都显示,‘final’只显示最终结果。optimset(‘Display’,‘off’)将设定Display选项为‘off’+ n% a& u4 K4 k/ H2 c+ C; w
& i$ d; {: D. d+ [3 常微分方程
! A* z1 G3 o% s. o) j7 H[t,y]=ode23(‘fname’,tspan,y0)
1 j( ^7 j+ p3 c2 o' e; e/ P[t,y]=ode45(‘fname’,tspan,y0) %其中fname是定义f(t,y)的函数文件名,该函数文件必须返回一个列向量。tspan形式为[t0,tf],表示求解区间。y0是初始状态列向量。t和y分别给出时间向量和相应的状态向量" i% U$ W& g& L) o1 ]0 t
6 z2 w. v/ w' i
4 极值
- o3 f' ^3 S7 T& i, U, I$ ^# Tfmin(‘fname’,x1,x2) %x1,x2范围* T$ c* u2 e5 y- @* z
fmins(‘fname’,x0) %x0为初始值向量
# X# R5 X2 B1 f& S2 _( C4 ?, N+ n. E& f6 U" K1 L* \/ d9 U
|
|