MATLAB的数值计算.ppt

建立m文件 function xdot=wf(t,x) xdot=zeros(2,1) xdot(1)=x(2) xdot(2)=x(2)*(1-x(1)^2)-x(1) 给定区间、初始值;求解微分方程 t0=0; tf=20; x0=[0 0.25]; [t,x]=ode23(wf, t0, tf, x0) plot(t,x), figure(2),plot(x(:,1),x(:,2)) 命令格式: [T,Y] = ODE23(ODEFUN,TSPAN,Y0) 建立m文件 function dxdt=wf(t,x) dxdt=[x(2);x(2)*(1-x(1)^2)-x(1)]; 求解微分方程 [t,x]=ode23(@wf,[0 30],[0 0.25]); plot(t,x); figure(2) plot(x(:,1),x(:,2)) 七、函数优化 寻优函数： fmin —— 单变量函数 fmins —— 多变量函数 constr —— 有约束条件 无约束条件 例1：f(x)=‘x2+3x+2’在[-5 5]区间的最小值 f=fmin(x^2+3*x+2,-5,5) 例2：f(x)=100(x2-x12)2+(a-x1)2在x1=a, x2=a2处有最小值 function f=xun(x,a) f=100*(x(2)-x(1).^2).^2+(a-x(1