集成在一个恒定的函数上

时间:2018-06-08 08:19:39

标签: matlab symbolic-math numerical-integration

我正在尝试在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);

1 个答案:

答案 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
...