matlab函数求极值

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,*,函数的极值,1、一元函数的极值,函数命令：fminbnd,调用格式：x,feval,exitflag,output=fminbnd(fun,x1,x2,options),%求fun在区间(x1,x2)上的极值.,返回值：,x:函数fun在(x1,x2)内的极值点,feval:求得函数的极值,exitflag:,exitflag0,函数收敛于解x处,exitflag=0,已达最大迭代次数,exiflag,A=fsolve(myfsolve2,ones(3)%,ones（3）初始矩阵，直全为1,myfsolve2(A),4、数值积分,函数命令：quad quadl odbquad,调用格式:,q1=quad(fun,a,b,tol),%采用自适应Simpson算法计算积分,q1=quadl(fun,a,b,tol),%采用自适应Lobatto算法计算积分,q2=dblquad(fun,xmin,xmax,ymin,ymax,tol),%二重积分,q3=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol),%三重积分,例1：求,I,=,.,(1)利用quad,quadl计算,format long,I=quad(exp(-x.2),0,1,1e-8),Il=quadl(exp(-x.2),0,1,1e-8),I=,0.74682413285445,Il=,0.74682413398845,(2)采用编程计算:,function I=myquad1(a,b,n),x=linspace(a,b,n);%,把ab区间平均分成n等份,y=exp(-x.2)*(b-a)/n;%,高底=每个取边梯形的面积,I=sum(y);,I1=myquad1(0,1,10000),I2=myquad1(0,1,100000),I1=,0.74681784375801,I2=,0.74682350396218,例,2,：计算,利用,dblquad,s1=dblquad(x.y,0,1,1,2),%,被积函数写成数组的形式,db_(double),二重积分，即两次积分,s1=,0.40546626724351,s2=dblquad(x,y)x.y,0,1,1,2),%,采用匿名函数表示被积函数,x y,的区间,s2=,0.40546626724351,s3=dblquad(inline(x.y),0,1,1,2),s3=,0.40546626724351,求,s1=,s2=,s3=,例,3,：,syms x,x,符号变量,f1=1/(1+x2);,f2=1/(x2+2*x+3);,f3=1/(x2+2*x-3);,s1=int(f1,1,inf),1,到正无穷,s2=int(f2,-inf,inf),int符号积分,s3=int(f3,-inf,inf),s1=,1/4*pi,s2=,1/2*pi*2(1/2),s3=,NaN,不确定的结果,s4=quad(1./(1+x.2),1,inf),%quad,只能计算闭型积分,Warning:Infinite or Not-a-Number function value encountered.,In quad at 95,s4=,NaN,s5=quad(1./(1+x.2),1,10),s5=,0.6857,s6=int(f1,1,10),s6=,atan(10)-1/4*pi,vpa(s6,10),ans=,.6857295105,