| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) | 
![]()  | 
	
| 		
			
			 | 
		#1 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2010-05-31 
				
				年龄: 34 
				
					帖子: 3
				 
				
				
				声望力: 0 ![]()  | 
	
	
	
		
		
			
			 
			
			我做一个人口增长的拟合,在确定参数初值时遇到了严重的麻烦,求高手给与帮助!!!万分感谢!!! 
		
		
		
		
		
		
		
	以下是我写的程序: >> t=[1790:10:2000]; >> population=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4]; >> ft=inline('a(1).*a(2)./(a(1)+exp(-a(3).*t).*a(2)-exp(-a(3).*t).*a(1))','a','t'); >> a=nlinfit(t,population,ft,[1 1 1]) Warning: Rank deficient, rank = 1, tol = 2.6167e-014. > In nlinfit>LMfit at 296 In nlinfit at 168 Warning: Rank deficient, rank = 1, tol = 2.6050e-014. > In nlinfit>LMfit at 296 In nlinfit at 168 Warning: Rank deficient, rank = 1, tol = 2.6038e-014. > In nlinfit>LMfit at 296 In nlinfit at 168 Warning: The Jacobian at the solution is ill-conditioned, and some model parameters may not be estimated well (they are not identifiable). Use caution in making predictions. > In nlinfit at 225 a = 1.0000 94.7727 1.0000 当然[1 1 1]是我瞎写的,显然最后的结果不对啊!但到底应该是什么呢?怎么确定?:handshake  | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#2 | 
| 
			
			 高级会员 
			
			
			
			注册日期: 2007-06-24 
				
				年龄: 70 
				
					帖子: 188
				 
				
				
				声望力: 22 ![]()  | 
	
	
	
		
		
			
			 
			
			a=[1.391057891*10^(-11), 44657285.16, 0.02154729647]
		 
		
		
		
		
		
		
		
	 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#3 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2010-05-31 
				
				年龄: 34 
				
					帖子: 3
				 
				
				
				声望力: 0 ![]()  | 
	
	
	
		
		
			
			 
			
			哥,请问你是如何算出的?能给与解释吗?非常感谢啊!!!
		 
		
		
		
		
		
		
		
	 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#4 | 
| 
			
			 高级会员 
			
			
			
			注册日期: 2007-06-24 
				
				年龄: 70 
				
					帖子: 188
				 
				
				
				声望力: 22 ![]()  | 
	
	
	
		
		
			
			 
			
			本例的非线性拟合有较大难度,原因在于对初值的要求太高。随机初值与最优值之间相差较大,用matlab自带的拟合命令(lsqcurvefit, nlinfit) 很难实现全局最优拟合。本人自编了一个小程序,能在任意初值条件下稳定实现一般曲线曲面的拟合问题。
		 
		
		
		
		
		
		
		
	 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#5 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2010-05-31 
				
				年龄: 34 
				
					帖子: 3
				 
				
				
				声望力: 0 ![]()  | 
	
	
	
		
		
			
			 
			
			功力深厚啊!能不能将您那个程序发给我用用呢?我邮箱[email protected]  QQ是349041471,如果可以的话非常感谢!! 
		
		
		
		
		
		
		
	 
		 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#6 | |
| 
			
			 初级会员 
			
			
			
			注册日期: 2011-04-23 
				
				年龄: 40 
				
					帖子: 4
				 
				
				
				声望力: 0 ![]()  | 
	
	
	
		
		
			
			 引用: 
	
 slgu,您好。我也遇到一个对初值非常敏感的非线性拟合问题,麻烦您也传我一份您的程序吧!还请您出手相救啊!! [email protected] 或者 [email protected] 叩谢!  | 
|
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#7 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2013-01-04 
				
				
				
					帖子: 1
				 
				
				
				声望力: 0 ![]()  | 
	
	
	
		
		
			
			 
			
			您好!  刚进实验室,导师让我编个程序对数据进行高斯拟合,也遇到了这种问题。能把您的程序发给我试试么?导师交的第一个任务就不完成以后再实验室都混不下去了。。。急求,谢谢啦!!我的邮箱是[email protected]
		 
		
		
		
		
		
		
		
	 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#8 | 
| 
			
			 高级会员 
			
			
			
			注册日期: 2007-06-24 
				
				年龄: 70 
				
					帖子: 188
				 
				
				
				声望力: 22 ![]()  | 
	
	
	
		
		
			
			 
			
			可将数据给出,或以文件的形式置上,并将可能的模型给出,否则我代为寻找模型需花费较长时间。我将尽力给出合适的结果。很抱歉,源程序尚不想公示。
		 
		
		
		
		
		
		
		
	 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 |