Question

matlab 定义常量

编写程序总是出问题，将问题简单化后如下，总是报错。请高人看下格式哪里出问题了。要求是Z是一个可变的常数，可以随时更改。
x=[1,2,3,4,5];
y=[7,16,31,0,79];
z=2;
myfunc=inline('beta(2)+beta(1).*x.*z','beta','z','x');
beta=nlinfit(x,y,myfunc,[10,0.1]);
a=beta(1),b=beta(2);

Caused by:
Error using inline/subsref (line 13)
Not enough inputs to inline function.

Answer 1

用nlinfit(x,y,@fun,beta) 只能让函数有两个输入fun（beta，x）

其中beta是参数，x是自变量
如果fun有多于两个的变量，第三个变量就不能通过调用nlinfit传递给fun
结果报错你的inline function输入参数数量不足。

而z也不能作为beta（3）传进去，因为拟合过程会改变他的值，最纵不能起到固定常数的作用

不过我们还是有办法的，把需要的额外不参与拟合的参数通过x传给fun
因为nlinfit不理x和y具体有多少个值，只要把x传递给fun，能够得到长度和y相同的数列就可以
所以，你要实现的功能可以通过以下程序实现

x=[1,2,3,4,5];
y=[7,16,31,0,79];
z=2;
x=[x z];
myfunc=inline('beta(2)+beta(1).*x(1:end-1).*x(end)','beta','x');
beta=nlinfit(x,y,myfunc,[10,0.1]);
a=beta(1);b=beta(2);

将z塞到x屁股后面，从x传递给myfunc，x变为[1,2,3,4,5,2];
但是要注意在函数内部x的最后一个元素作为z，前面的5个元素参与运算
那样myfunc的计算结果仍然是和y序列等长的5个数

顺便讲句，其实可以通过这样的办法实验二维曲面的拟合
设现有三组等长列向量对应二维曲面的x，y，z坐标
x=[x1 x2 ...... xn]';
y=[y1 y2 ...... yn]';
z=[z1 z2 ...... zn]';
要拟合为曲线方程z=f(beta,x,y),其中beta是参数，x，y是两个输入
以二维旋转对称高斯曲面为例子z=A*exp((x^2+y^2)/B^2)+C

可以这样操作
X=[x y];
Y=z;
f=inline('beta(1)+exp((X(:,1).^2+X(:,2).^2)./beta(2)^2)+beta(3)','beta','X');
beta=nlinfit(X,Y,f,beta0); %参数的初始值beta0
A=beta(1);B=beta(2);C=beta(3);