我是一名大四的學(xué)生,現(xiàn)在正在做畢業(yè)設(shè)計(jì),因MATLAB從未學(xué)過,也是邊學(xué)邊做,我需要求出IGv函數(shù)的頻率響應(yīng)圖,因IGv的表達(dá)式很復(fù)雜,這里我沒列出,在下面的程序中有的,以下是我的程序,請(qǐng)高手幫我看看程序有什么問題哈,萬分感謝!

%%igmax is global maximum

Igmin=1000000000;

%%exmperimantal values

M=0:1:100;

C=0:1:100000;

For omiga=0:1:100;

K=0:1:1000;

%%igvlmax is local maximum;

Igvlmax=0;

%%igv’s calculational equation

Igv=sqrt(16*c.^4*omiga.^4+16*k.^4+4*M.^2*k.^2*omiga.^4+32*c.^2*omiga.^2*k.^2-16*M*k.^3*omiga.^2+4*M.^2*c.^2*omiga.^6)/((-M*omiga.^2+2*k).^2+4*c.^2*omiga.^2);

If igv>=igvlmax;

Igvlmax=igv;

End;

End;

If igvlmax<=igvmin;

Igvmin=igvlmax;

Kmin=k; cmin=c; mmin=m;

End

End

End

End

If IGvlmin>=IGvmax

K=Kmin

c=Cmin

m=Mmin

For omiga=0:1:100;

Igv(i)=sqrt(16*c.^4*omiga.^4+16*k.^4+4*M.^2*k.^2*omiga.^4+32*c.^2*omiga.^2*k.^2-16*M*k.^3*omiga.^2+4*M.^2*c.^2*omiga.^6)/((-M*omiga.^2+2*k).^2+4*c.^2*omiga.^2);

End;

For i=0:1:100

Plot (I, igv(i))

其中k,c,omega,M是變量,整個(gè)函數(shù)就是要在不同的k,c,omega,M取值時(shí)使IGv達(dá)到最小,也就是函數(shù)最優(yōu)化的問題,我寫的程序不知道對(duì)不對(duì),錯(cuò)的地方我也不太會(huì)改,請(qǐng)各位高手幫幫忙哈.

總結(jié)

以上是生活随笔為你收集整理的matlab频响优化,MATLAB中关于频响图函数最优化的程序问题的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。