matlab频响优化,MATLAB中关于频响图函数最优化的程序问题
我是一名大四的學生,現在正在做畢業設計,因MATLAB從未學過,也是邊學邊做,我需要求出IGv函數的頻率響應圖,因IGv的表達式很復雜,這里我沒列出,在下面的程序中有的,以下是我的程序,請高手幫我看看程序有什么問題哈,萬分感謝!
%%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是變量,整個函數就是要在不同的k,c,omega,M取值時使IGv達到最小,也就是函數最優化的問題,我寫的程序不知道對不對,錯的地方我也不太會改,請各位高手幫幫忙哈.
總結
以上是生活随笔為你收集整理的matlab频响优化,MATLAB中关于频响图函数最优化的程序问题的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: mysql select null 0,
- 下一篇: matlab学习-线性规划