我正在尝试在MATLAB 2017a中集成一个常量函数,但我被卡住了。首先,当我使用以下脚本进行集成时,我得到了正确的输出。因此,该脚本适用于x0
,取决于t
。
function E=sol(n,k)
x0 = @(t) t^(2);
j = 0;
E = zeros(n,1);
while j < n+1 ;
K = matlabFunction(subs(po(j,k))) ;
eval(sprintf('x%d = integral(K,0,1);',j+1)) ;
E(j+1,1) = subs(sprintf('x%d',j+1))
j = j+1;
end
end
函数po(j,k)
如下所示,
function A_j = po(j,k) % Adomian polynomials
if j >0
x = sym('x',[1 j]);
syms p; % Assinging a symbolic variable for p
syms x0;
S = x0+ sum(p.^(1:j) .* x) ; % Sum of p*x up to order j
Q =f(S,k); % Taking the k-th power of S, i.e.
A_nc = diff(Q,p,j)/factorial(j); % Taking the j-th order derivative
A_j = subs(A_nc,p,0) ; % Filling in p=0
else
syms x0;
S = x0;
A_j =f(S,k); % Taking the k-th power of S,
end
end
f(x,k)
的位置,
function F = f(x,k) % Nonlinear function of k power
F = x^k ;
end
现在,当我{cal sol(n,k)
时,它确实有效。但是,当我尝试在x0
中使用常量函数更改sol(n,k)
函数时,
function E=solcon(n,k)
x0 = @(t) 2.*ones(size(t));
j = 0;
E = zeros(n,1);
while j < n+1 ;
K = matlabFunction(subs(po(j,k))) ;
eval(sprintf('x%d = integral(K,0,1);',j+1)) ;
E(j+1,1) = subs(sprintf('x%d',j+1))
j = j+1;
end
end
它不起作用,正如您所看到的,我添加*ones(size(t));
只是为了使其成为t
的函数。但不幸的是,当我打电话时,它仍然不起作用,
K = matlabFunction(subs(po(j,k))) ;
我明白了,
@()4.0
然后我打电话时出错,
eval(sprintf('x%d = integral(K,0,1);',j+1))
任何人都可以帮助我尝试整合常量吗?
我致电solcon(10,2)
时收到的错误是
Error using symengine>@()4.0
Too many input arguments.
Error in integralCalc/iterateScalarValued (line 314)
fx = FUN(t);
Error in integralCalc/vadapt (line 132)
[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);
Error in integralCalc (line 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);
Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);
Error in solcon1 (line 7)
eval(sprintf('x%d = integral(K,0,1);',j+1)) ;
编辑2 我使用了以下脚本,
function E=solcon(n,k)
x0 = @(t) 2.*ones(size(t));
j = 0;
E = zeros(n,1);
while j < n+1 ;
K = matlabFunction(subs(po(j,k))) ;
fstr= func2str(K)
if fstr(3) == ')';
x{j+1} = K*(1-0)
else x{j+1} = integral(K,0,1)
end
E(j+1,1) = subs(x{j+1},1);
j = j+1
end
end
但发生以下错误,
Undefined operator '*' for input arguments of type 'function_handle'.
Error in solcone1 (line 9)
x{j+1} = K*(1-0);
答案 0 :(得分:11)
我将忽略使用eval
的可怕选择,尤其是当你可以
x{j+1} = integral(K,0,1);
Read more on why dynamic variables and eval are terrible
您的问题是matlabFunction
是智能手机。当它检测到您的函数与x
没有任何依赖关系时,它会为您提供一个空输入参数@()4.0
的函数。如您所见,integral
不喜欢这样。
解决问题的一种方法是在调用integral
之前检测到这个问题。您可以检查它是否有输入参数,如果没有,则可以“手动”评估积分
...
while j < n+1 ;
K = matlabFunction(subs(po(j,k))) ;
fstr=func2str(K);
if fstr(3)==')'
x{j+1}=K()*(1-0); % evaluate the integral yourself
else
x{j+1} = integral(K,0,1);
end
E(j+1,1) = subs(x{j+1});
j = j+1;
end
...
虽然我是这样,但问题要困难得多。要么重写整个事物,要么使用eval:
...
while j < n+1 ;
K = matlabFunction(subs(po(j,k))) ;
fstr=func2str(K);
if j==0
eval(sprintf('x%d = K()*(1-0);;',j+1)) ;
else
eval(sprintf('x%d = integral(K,0,1);',j+1)) ;
end
E(j+1,1) = subs(sprintf('x%d',j+1));
j = j+1;
end
...