第五章Matlab在高等数学中的应用讲解.pptVIP

MATLAB中的函数lsqnonlin用来求解如下形式的非线性最小二乘优化问题: 调用格式有： (1)x=lsqnonlin(fun,x0); (2)x=lsqnonlin(fun,x0,lb,ub); (3)x=lsqnonlin(fun,x0,lb,ub,options); (4)x=lsqnonlin(problem); (5)[x,resnorm]=lsqnonlin(...)返回值resnorm为残差的平方，即最优值 (6) [x,resnorm,residual]=lsqnonlin(...)返回值residual为残差 (7) [x,resnorm,residual,exitflag]=lsqnonlin(...) (8) [x,resnorm,residual,exitflag,output]=lsqnonlin(...) (9) [x,resnorm,residual,exitflag,output,lambda]=lsqnonlin(...) (10) [x,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(...)jacobian是最优值x处的雅克比矩阵 例题 建立方程的m文件 function y=lsqfun(x) y=[sin(x(1)+x(2)-2); 1/(2-(x(1)-3)^2); exp(2*x(1))+exp(2*x(2)); x(1)^2+x(2)^2-x(1)*x(2) ... +x(1)+1]; 在MATLAB R2009a中用于求解线性规划的函数式linprog，能求解下面形式线性规划问题： 调用形式如下： (1)x=linprog(f,A,b) 求解 (2)x=linprog(f,A,b,Aeq,beq) (3)x=linprog(f,A,b,Aeq,beq,lb,ub) (4) x=linprog(f,A,b,Aeq,beq,lb,ub,x0) (5) x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options) (6)x=linporog(problem) 结构problem用来指定线性规划问题 (7)[x,fval]=linprog(...) (8)[x,fval,exitflag]=linprog(...) (9) [x,fval,exitflag,output]=linprog(...) (10)[x,fval,exitflag,output,lambda]=linprog(...)lambda是包含最优解处拉格朗日乘子的结构体 例题： 在MATLAB中，提供了函数bintprog来求解0-1整数规划，其优化模型为： 它的调用形式有： x=bintprog(f) x=bintprog(f,A,b) x=bintprog(f,A,b,Aeq,beq) x=bintprog(f,A,b,Aeq,beq,x0) x=bintprog(f,A,b,Aeq,beq,x0,options) x=bintprog(problem) [x,fval]=bintprog(...) [x,fval,exitflag]=bintprog(...) [x,fval,exitflag,output]=bintprog(...) 例题： MATLAB中的函数quadprog可用来解决如下二次规划模型： 调用形式为： (1)x=quadprog(H,f,A,b); (2)x=quadprog(H,f,A,b,Aeq,beq); (3)x=quadprog(H,f,A,b,Aeq,beq,lb,ub); (4) x=quadprog(H,f,A,b,Aeq,beq,lb,ub,x0); (5) x=quadprog(H,f,A,b,Aeq,beq,lb,ub,x0,options); (6) x=quadprog(problem); (7)[x,fval]=quadprog(...) (8)[x,fval,exitflag]=quadprog(...) (9) [x,fval,exitflag,output]=quadprog(...) (10) [x,fval,exitflag,output,lambda]=quadprog(...) 例题： 计算结果见下页 注：在默认情况下，如果指定了如下两种约束中的任意一种，quadprog采用大型算法(large-scale)。 (1)自变量的上界和下界约束； (2)线性方程约束，约束矩阵Aeq的列向量线性无关，且Aeq是稀疏矩阵。 5.5. 积分的计算 5.5.1 定积分和不定积分 符号积分用函数int来实现，其调用格式为：