matlab数学建模微分法与优化问题.pptxVIP

M文件;例：定义函数 f(x1,x2)=;关系与逻辑运算;2. 逻辑运算符;for循环 允许一组命令以固定的和预定的次数重复 for x=array commands； end 在for和end语句之间的命令串commands按数组(array)中的每一列执行一次. 在每一次迭代中，x被指定为数组的下一列，即在第n次循环中，x=array(：，n); while expression commands； end 只要在表达式(expression)里的所有元素为真，就执行while和end语句之间的命令串commands. ; if-else-end结构 （1）一个选择的一般形式是： if (expression) commands； end 如果在表达式(expression)里的所有元素为真，就执行if和end语句之间的命令串commands. ;（2）多个选择的一般形式是： if （expression1） commands1； else if （expression2） commands2； else if （expression3） commands3； else if …… ………………………………… else commands； end end end …… end ;continue命令的作用是结束本次循环，即跳过循环体中尚未执行的语句，接着进行下一次是否执行循环的判断。(当条件满足时，continue语句才被调用，跳出后面的语句，循环继续，条件不满足则继续执行后面的程序。);input命令用来提示用户从键盘输入数据、字符串或表达式，并接受输入值。 r=input(‘prompt’)在屏幕上显示提示信息prompt，等待用户输入，并将值赋给r，直接回车则输出空矩阵。 r=input(‘prompt’,’s’)返回的字符串作为文本变量而不作为变量名或者数值。;代数方程的求解;函数 fzero 非线性方程的数值解 格式 x = fzero (fun,x0) %用fun定义表达式f(x)，x0为初始解。 x = fzero (fun,x0,options) 说明 该函数采用数值解求方程f(x)=0的一个根。;非线性方程组数值解;微分方程的解析解; y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x);微分方程的数值解;使用泰勒公式，有龙格-库塔法、线性多步法等方法．;; 令 y1=x，y2=y1’; 1．建立M文件rigid．m如下： function dy=rigid(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-0.51*y(1)*y(2);;;由(1),(2)消去t, 整理得模型:;解法二(数值解法);解法三(建立参数方程求数值解);4． 解导弹运动轨迹的参数方程;用MATLAB优化工具箱解线性规划;3. 模型：min z=cX ;解 编写M文件xxgh.m如下： c=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6]; A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.