用MATLAB求極值

1、用MATLAB求極值靈活的運用MATLAB的計算功能,可以很容易地求得函數(shù)的極值._23x4x4.八例3.6.1求y3x24X4的極值xx1解首先建立函數(shù)關(guān)系：symssy=(3*xA2+4*x+4)/(xA A2+x+1);/然后求函數(shù)的駐點：dy=diff(y);/xz=solve(dy)/xz=0-2知道函數(shù)有兩個駐點x1=0和x2=-2,考察函數(shù)在駐點處二階導數(shù)的正負情況：d2y=diff(y,2);/z1=limit(d2y,x,0)/z1=-2z2=limit(d2y,x,-2)/z2=2/9于是知在x1=0處二階導數(shù)的值為z1=-2,小于0,函數(shù)有極大值；在x2=-2處二階導數(shù)的

2、值為Z2=2/9,大于0,函數(shù)有極小值.如果需要,可順便求出極值點處的函數(shù)值：y1=limit(y,x,0)/y=4y2=limit(y,x,-2)/y2=8/3事實上,如果知道了一個函數(shù)的圖形,那么它的極值情況和許多其它特性是一目了然的.而借助MATLAB的作圖功能,我們很容易做到這一點.例解symsx/y=(3*xA2+4*x+4)/(xA2+x+1);/得到如下列圖形ezplot(y)/如何用MATLAB求函數(shù)的極值點和最大值比方說y=x+xA2+1,怎樣用matlab來算它的極值和最大值？求極值：symsxyy=xA A3+xA A2+1diff(y)%求導ans=3*xA2+2*xs

3、olve(ans)%求導函數(shù)為零的點ans=-2/30極值有兩點.求最大值,既求-y的最小值：f=(x)(-xA3-xA2-1)f=(x)(-xA3-xA2-1)x=fminunc(f,-3,3)%在-3;-3范圍內(nèi)找Warning:Gradientmustbeprovidedfortrust-regionmethod;usingline-searchmethodinstead.Infminuncat354Optimizationterminated:relativeinfinity-normofgradientlessthanoptions.TolFun.x=-0.6667f(x)ans=-

4、1.1481在規(guī)定范圍內(nèi)白最大值是1.1481由于函數(shù)的局限性,求出的極值可能是局部最小(大)值.求全局最值要用遺傳算法.如何用matlab求多元函數(shù)的極值在0 xpi/2,0ydzy=diff(z,y)dzy=sin(x)*cos(y)*sin(x+y)+sin(x)*sin(y)*cos(x+y)xx,yy=solve(cos(x)*sin(y)*sin(x+y)+sin(x)*sin(y)*cos(x+y)=0,.sin(x)*cos(y)*sin(x+y)+sin(x)*sin(y)*cos(x+y)=0,x,y)%,y可以不用？A=diff(z,x,2);%對于定義域內(nèi)的駐點求解二階

5、偏導函數(shù)B=diff(diff(z,x),y);C=diff(z,y,2);D=A*C-BA2;wu別極大值點D1=subs(subs(D,x,xx(1),y,yy(1)%subs含義見補充局部A1=subs(subs(A,x,xx(1),y,yy(1)D2=subs(subs(D,x,xx(2),y,yy(2)A2=subs(subs(A,x,xx(2),y,yy(2)%極大值點和極大值xx(2)yy(2)JDZ=subs(subs(z,x,xx(2),y,yy(2)程序運行結(jié)果ans=1/3*pi,1/3*piJDZ=3/8*3A(1/2)補充matlab中的繪圖函數(shù)很多,三維的有好幾個呢

6、,下面我給出兩種繪制方法1.使用ezmesh或者ezsurf這兩個函數(shù)是簡易繪圖函數(shù),可以直接使用字符串繪圖ezmesh(x*y-(1/3)*(xA3)+1/4*(xA2)*(sqrt(2*(xA2-y*2)-1/2*(y*(sqrt(2*(xA2-y*2)2.先計算坐標點數(shù)據(jù),再繪圖,可以使用mesh和surf函數(shù)x,y=meshgrid(0:0.1:10);z=x.*y-(1/3)*(x.A A3)+1/4*(x.A A2).*(sqrt(2*(x.A A2-y*2)-1/2*(y.*(sqrt(2*(x.A A2-y*2);z(imag=0)=nan;mesh(x,y,z)subs是賦值

7、函數(shù),用數(shù)值替代符號變量替換函數(shù)例如：輸入subs(a+b,a,4)意思就是把a用4替換掉,返回4+b也可以替換多個變量,例如：subs(cos(a)+sin(b),a,b,sym(alpha),2)分別用字符alpha替換a和2替換b,返回cos(alpha)+sin(2)用法很靈活,仔細看幫助,會得到你想要的形式的matlabsubs如何用懸賞分：0-解決時間：2021-12-317:36在解一個模型,其中：Qd=(-2*x-2*wd+2*a-a*k+y*k+wf*k+td*k)/(4*b-2*b*kA2)現(xiàn)在求出來x=1/2*(48*a+32*k+16*wd-16*cdu-4*wf*k-

8、4*td*k-12*a*k-10*a*kA2+td*kA3-8*kA3+4*kA2*cdu-4*k*cfu+cfu*kA3+3*a*kA3-6*wd*kA2+wf*kA3)/(5*kA2-16)y=(-8*a+a*kA2-2*wd*k+10*a*k+8*wf-3*wf*kA2-3*td*kA2+2*cfu*kA2+4*kA2+8*td-8*cfu-2*k*cdu)/(5*kA2-16)要把x和y代回Qd,請教達人,用matlab的subs怎么寫??？萬分感謝！問題補充：Qd=(-2*x1-2*wd+2*a-a*k+y1*k+wf*k+td*k)/(4*b-2*b*kA2)x1=1/2*(48*a

9、+32*k+16*wd-16*cdu-4*wf*k-4*td*k-12*a*k-10*a*kA2+td*kA3-8*kA3+4*kA2*cdu-4*k*cfu+cfu*kA3+3*a*kA3-6*wd*kA2+wf*kA3)/(5*kA2-16)y1=(-8*a+a*kA2-2*wd*k+10*a*k+8*wf-3*wf*kA2-3*td*kA2+2*cfu*kA2+4*kA2+8*td-8*cfu-2*k*cdu)/(5*kA2-16)subs(Qd,x1,y1,x,y)以上輸入之后報錯為：Undefinedfunctionorvariablex1.請問是哪里除了問題？提問者：nxdviol

10、a一級最正確答案symsx1,y1QdQd=(-2*x1-2*wd+2*a-a*k+y1*k+wf*k+td*k)/(4*b-2*b*kA2)subs(Qd,x1,y1,x,y)為何Matlab用subs后給出的結(jié)果是個代數(shù)式而不是數(shù)字？懸賞分：0-解決時間：2021-10-822:16使用d=subs(f4,x,y,z,2,3,1)后,結(jié)果給出了一個很長的代數(shù)式,里面有pi和atan,但都沒有計算,為什么不給出一個最終的數(shù)字結(jié)果？提問者：再見長江-二級最正確答案subs函數(shù)就是替換符號表達式的函數(shù),要計算需要用eval函數(shù)急詢：Matlab中subs(S)的應(yīng)用疑問？懸賞分：5-解決時間：2

11、006-6-2912:57急詢：Matlab中subs(S)的應(yīng)用疑問？1.subs(S)到底是什么函數(shù)？其什么作用,謝謝！2.下段繪圖程序如下的話就畫出不來：symsxt;t=0:pi/60:2*pi;y1=int(sin(x),x,0,t);y2=int(cos(x),x,0,t);plot(y1,y2)而把plot(y1,y2)改成plot(subs(y1),subs(y2)就可以運行成功畫出來了.(1)為什么？(2)subs在此起到什么作用？(3)而不加subs為什么運行不了,畫不出來？急,多謝！(上次提問時候把y1寫成y1=sin(x)了,現(xiàn)在改為int(sin(x),x,0,t);

12、抱歉)問題補充：(上次提問時候把y1寫成y1=sin(x)了,現(xiàn)在改為int(sin(x),x,0,t);抱歉)請具體講一下subs(S)到底是什么意思？(我知道subs(f,x,t)是什么意思,但我不明白sub(S)在此到底具體是什么意思？請詳細一下,謝謝！)提問者：wanglicun-二級最正確答案subs(S)的意思是：S以前是sym變量,subs(S)為double變量,plot只能對double變量作圖.曲線極值的標注：用matlab中plot繪制的曲線時,把曲線上峰值的大小標注在上面你指的是離散點的最大和最小值吧？給你個簡單的例子看看：代碼：x=0:pi/50:2*pi;y=sin

13、(x);y1=max(y);x1=x(find(y=y1);y2=min(y);x2=x(find(y=y2);plot(x,y);holdonplot(x1,y1,r*);plot(x2,y2,g*);離散數(shù)據(jù)就極值點腦法一x=0:0,01:20;y=2*sin(x/2)+cos(2*x)/2;%indmax=find(diff(sign(diff(y)0)+1;%極小值點ind=find(diff(sign(diff(y)=0)+1;%極值點plot(x,y,x(ind),y(ind),ro)%r法二x=0:0.01:20;y=2*sin(x/2)+cos(2*x)/2;indmax=intersect(find(diff(y)0)+1,find(diff(y)0);%極大值點indmin=intersect(find(diff(y)0)