请指教三元非线性函数拟合 - MATLAB爱好者论坛-LabFans.com

lgdca · 2009-07-27, 22:13

请求高手指教！
有如下式子：v=k*P^a*(1-(pe/p)^a*(n/ne)^b)
其中待求参数为k,a,b,其他p,n,v(都为22行的列向量)；pe=1106.1, ne=2.4607
p=[2686
2675
2670
2659
2653
2648
2642
2637
2631
2626
2620
2615
2604
2598
2593
2587
2582
2576
2571
2565
2560
2554];
n=[1.52157
1.72677
1.82006
2.0253
2.13727
2.23058
2.34256
2.43588
2.54788
2.64121
2.75322
2.84657
3.05195
3.16399
3.25736
3.36942
3.4628
3.57487
3.66827
3.78036
3.87377
3.98587];
v=[0.10259
0.1026
0.046645
0.10262
0.055985
0.046655
0.05599
0.04666
0.056
0.046665
0.056005
0.046675
0.051345
0.05602
0.0233425
0.028015
0.015563333
0.056035
0.02335
0.018681667
0.015568333
0.018683333];
我试过fminunc，设置了option.TolFun=1e-20，老是从中间跳出，不知还有什么函数或方法，恳求高手指教！

slgu · 2009-07-29, 06:57

[k,a,b]=[-12.5751052758, -0.6558990751, -0.600379582476]
RSS =0.004825553373
MSe = 0.00026809
R^2 = 0.67342
v与p，n的关系不太密切，所用模型并非非常恰当。

2009-07-27, 22:13	#1
lgdca 初级会员注册日期: 2009-06-03 年龄: 43 帖子: 1 声望力: 0	请指教三元非线性函数拟合请求高手指教！有如下式子：v=kP^a(1-(pe/p)^a(n/ne)^b) 其中待求参数为k,a,b,其他p,n,v(都为22行的列向量)；pe=1106.1, ne=2.4607 p=[2686 2675 2670 2659 2653 2648 2642 2637 2631 2626 2620 2615 2604 2598 2593 2587 2582 2576 2571 2565 2560 2554]; n=[1.52157 1.72677 1.82006 2.0253 2.13727 2.23058 2.34256 2.43588 2.54788 2.64121 2.75322 2.84657 3.05195 3.16399 3.25736 3.36942 3.4628 3.57487 3.66827 3.78036 3.87377 3.98587]; v=[0.10259 0.1026 0.046645 0.10262 0.055985 0.046655 0.05599 0.04666 0.056 0.046665 0.056005 0.046675 0.051345 0.05602 0.0233425 0.028015 0.015563333 0.056035 0.02335 0.018681667 0.015568333 0.018683333]; 我试过fminunc，设置了option.TolFun=1e-20，老是从中间跳出，不知还有什么函数或方法，恳求高手指教！此帖于 2009-07-29 19:50 被 lgdca 编辑。原因: k,a,b都为正数，范围应当在（0，10）*

2009-07-29, 06:57	#2
slgu 高级会员注册日期: 2007-06-24 年龄: 70 帖子: 188 声望力: 21	回复: 请指教三元非线性函数拟合 [k,a,b]=[-12.5751052758, -0.6558990751, -0.600379582476] RSS =0.004825553373 MSe = 0.00026809 R^2 = 0.67342 v与p，n的关系不太密切，所用模型并非非常恰当。