Question

matlab画自定义函数图象

myfunc=inline('beta(1).*(1+beta(2).*exp(beta(3).*(beta(4).-x)))^(-1./beta(2))','beta','x');
上面这个公式是Richards模型函数
哪位大神能告诉我错在哪儿了？

Answer 1

错在多写了“.”和")"
myfunc=inline('beta(1).*(1+beta(2).*exp(beta(3).*(beta(4)-x))^(-1/beta(2))','beta','x');

追问

但是还是用不起来诶
x=[1:1:45];
y=[3 7 9 14 16 18 20 24 28 33 38 43 49 60 63 77 82 87 91 96 102 104 108 108 109 112 115 124 126 126 127 127 127 128 128 129 130 130 130 131 131 131 131 131 131];
myfunc=inline('a(1).*(1+a(2).*exp(a(3).*(a(4)-x))^(-1/a(2))','a','x');
a=nlinfit(x,y,myfunc,[0.1 0.1 0.1 0.1]);
plot(x,y,'*',a(1).*(1+a(2).*exp(a(3).*(a(4)-x))^(-1/a(2)))

能帮看一下么

追答

我试一下，有几个问题。

1、你给出的拟合函数不对，从散点图可以看到，所拟合的函数应该是S型曲线函数。

S型曲线函数 y=1/(a+b*exp(-cx))

2、函数使用不当。应该使用lsqcurvefit（）而不是nlinfit（）。

3、绘图函数使用不正确。

x=1:1:45;
y=[3 7 9 14 16 18 20 24 28 33 38 43 49 60 63 77 82 87 91 96 102 104 108 108 109 112 115 124 126 126 127 127 127 128 128 129 130 130 130 131 131 131 131 131 131];
myfunc=inline('1./(a(1)+a(2).*exp(-a(3).*x))','a','x');
b=[rand(1)/1000    rand(1)/1000    1]
a=lsqcurvefit(myfunc,b,x,y);
vpa(a)
xi=1:1:45;
yi=1./(a(1)+a(2).*exp(-a(3).*xi));
plot(x,y,'*',xi,yi,'r-')

Answer 2

是不是括号太多了