我正在尝试运行一个包含函数句柄的循环。如果我对 n 使用任何特定值,则代码运行良好。但是,当我尝试在循环中运行代码时,会出现错误:
Nonscalar arrays of function handles are not allowed; use cell
arrays instead.
我对 MATLAB 比较陌生,所以在尝试了几个小时后,我的想法已经用完了......任何帮助都将不胜感激。
代码如下
for m=1:N
for n=1:N
int(n) = @(w,s0,v,r,t) chfun_norm(s0,v,r,t,w(n)-(alpha+1)*1i)/(alpha^2 + alpha - w(n)^2 + 1i*(2*alpha+1)*w(n));
int(n) = @(w)int(w,s0,v,r,t);
g(n) = exp(1i*(b-log(s0))*w(n)-alpha*k(m)-r*t)*int(n)*pond(n);
g(n) = exp(-1i*2*pi/N*(n-1)*(m-1))*g(n);
end
y(m) = real(sum(g))*stepw;
end
编辑
我已经使用单元格数组重新制定了循环,但我现在收到以下错误
Undefined function 'mtimes' for input arguments of type 'function_handle'
以下是更新后的代码:
y = zeros(N,1);
intCell = {zeros(N,1)};
gCell = {zeros(N,1)};
for m=1:N
for n=1:N
intCell{n} = @(w,s0,v,r,t) chfun_norm(s0,v,r,t,w(n)-(alpha+1)*1i)/(alpha^2 + alpha - w(n)^2 + 1i*(2*alpha+1)*w(n));
intCell{n} = @(w)intCell{n}(w,s0,v,r,t);
gCell{n} = exp(1i*(b-log(s0))*w(n)-alpha*k(m)-r*t)*intCell{n}*pond(n);
gCell{n} = exp(-1i*2*pi/N*(n-1)*(m-1))*g(n);
end
y(m) = real(sum(g))*stepw;
end