Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#2 |
高级会员
注册日期: 2008-11-07
住址: 湖南长沙
帖子: 233
声望力: 21 ![]() |
![]()
这样就可以了,其实程序没有什么大问题,就是少了几个end和多了几处分号
%%IGvmin is global minimum IGvmin=1000000000000000; %%exmperimantal values M=1000000; for c=0:1:100000 for omiga=0:1:100 for k=100000:100:100000000 %%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 IGvmin>=IGvlmax IGvlmax=IGvmin; end end if IGvlmax<=IGvmin IGvmin=IGvlmax; Kmin=k; cmin=c; end end end if IGvlmin>=IGvmax k=kmin c=Cmin end 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)) hold on; end
__________________
坚持就是胜利,努力就有奇迹。 |
![]() |
![]() |