请指教三元非线性函数拟合

[lgdca]

请求高手指教！
有如下式子：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，老是从中间跳出，不知还有什么函数或方法，恳求高手指教！