| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) | 
![]()  | 
	
| 		
			
			 | 
		#2 | 
| 
			
			 高级会员 
			
			
			
			注册日期: 2008-11-07 
				住址: 湖南长沙 
				
				
					帖子: 233
				 
				
				
				声望力: 22 ![]()  | 
	
	
	
		
		
			
			 
			
			这样就可以了,其实程序没有什么大问题,就是少了几个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 
				__________________ 
		
		
		
		
	坚持就是胜利,努力就有奇迹。  | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 |