| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) | 
![]()  | 
	
| 
	 | 
| 		
			
			 | 
		#1 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2008-03-20 
				
				年龄: 38 
				
					帖子: 1
				 
				
				
				声望力: 0 ![]()  | 
	
	
	
		
		
			
			 
			
			各位大侠帮我看看我的程序哪有问题! 
		
		
		
		
		
		
		
	clc; clear; syms Kd Ka n1 n2 y1 y2; X=[Kd Ka]; N=[n1 n2]; y=[y1 y2]; Z=[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*8/4)-exp(-X(2)*8/4))-8.5]^2+[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*28/4)-exp(-X(2)*28/4))-7]^2+[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*36/4)-exp(-X(2)*36/4))-6.1]^2+[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*56/4)-exp(-X(2)*56/4))-7.2]^2; x=[1 2]; old=subs(Z,[X(1) X(2)],[x(1) x(2)]); Z1=[diff(Z,X(1));diff(Z,X(2))]; H=[diff(Z,X(1),2) diff(Z1(1),X(2));diff(Z1(2),X(1)) diff(Z,X(2),2)]; step=Z1.'*Z1\(Z1.'*H*Z1); F=subs(step,[X(1) X(2)],[x(1) x(2)]); N(1)=x(1)-F*subs(Z1(1),[X(1) X(2)],[x(1) x(2)]); N(2)=x(2)-F*subs(Z1(2),[X(1) X(2)],[x(1) x(2)]); new=subs(Z,[X(1) X(2)],[N(1) N(2)]); while (abs(new-old)/old)>exp(-6), x(1)=N(1); x(2)=N(2); old=new; F=subs(step,[X(1) X(2)],[x(1) x(2)]); N(1)=x(1)-F*subs(Z1(1),[X(1) X(2)],[x(1) x(2)]); N(2)=x(2)-F*subs(Z1(2),[X(1) X(2)],[x(1) x(2)]); new=subs(Z,[X(1) X(2)],[N(1) N(2)]); end X(1)=N(1); X(2)=N(2); 为什么总说??? Error using ==> >= Function '>=' is not defined for values of class 'sym'. 要怎么改改啊?谢谢各位!  | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
![]()  | 
	
	
| 主题工具 | |
| 显示模式 | |
		
  | 
	
		 | 
			 
			相似的主题
		 | 
	||||
| 主题 | 主题作者 | 版面 | 回复 | 最后发表 | 
| BMP图像转换成时间序列 | zlnhappy | MATLAB论坛 | 0 | 2008-04-28 17:53 | 
| matlab编写板卡驱动 | emirlee | MATLAB论坛 | 2 | 2008-04-18 09:16 | 
| 【求助】如何做FFT变换? | ceeshy | MATLAB论坛 | 1 | 2007-12-21 11:27 | 
| 求救-----基于小波域的数字水印 | tcgqb111 | MATLAB论坛 | 0 | 2007-05-24 07:27 |