【问题标题】:Is it possible to define an integral function and integrate it using `integral`?是否可以定义一个积分函数并使用“integral”对其进行积分?
【发布时间】:2017-05-19 04:57:37
【问题描述】:

我想做以下整合:

为了在matlab中使用函数integral进行积分,我定义了函数

先申请integral。以下是我的功能:

OrderEpsilon2IntegrandIntegralIntegrand=@(T,s)lambda.*(T-s).*exp(-1.*lambda.*s);
OrderEpsilon2IntegrandIntegral=...
@(T,u)integral(@(t)OrderEpsilon2IntegrandIntegralIntegrand(T,t),u,T);

OrderEpsilon2Integrand=@(T,u)...
(T-u).*sigma(u).*v_0-...
OrderEpsilon2IntegrandIntegral(T,u).*exp(lambda.*u).*sigma(u).*v_0;

OrderEpsilon2Integral=...
@(t)arrayfun(@(T)integral(OrderEpsilon2Integrand(T,u),InitialTime,T),t);

但是,我收到以下错误消息:

Error using integral (line 85)
A and B must be floating-point scalars.

Error in (T,u)integral(@(t)OrderEpsilon2IntegrandIntegralIntegrand(T,t),u,T)

Error in @(T,u)(T-u).*sigma(u).*v_0-OrderEpsilon2IntegrandIntegral(T,u).*exp(lambda.*u).*sigma(u).*v_0

Error in @(u)OrderEpsilon2Integrand(2.3,u)

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);

出现这种错误的问题是什么?我对消息没有任何想法。另外,我想知道如何更正代码。但是,纠正代码有一些限制:

首先,不要写

明确。

其次,不要使用integral2。如果可能的话,我想定义具有积分的被积函数。

【问题讨论】:

    标签: matlab integral


    【解决方案1】:

    为什么不呢?:

    syms s u T t sg lm
    %%sg=u^2+u+1; %%Define sigma expression in here
    f=int((T-u)*sg-int(lm*(T-s)*exp(-lm*s),s,u,T)*exp(lm*u),u,t,T) 
    

    导致:

    T*t + T/lm + (T^2*sg)/2 - t/lm + (sg*t^2)/2 - T^2/2 - 1/lm^2 - t^2/2 + (exp(-T*lm)*exp(lm*t))/lm^2 - T*sg*t
    

    【讨论】:

    • 这是因为最终我想计算,例如 OrderEpsilon2Integral([1,2,3,4])。定义符号不能帮助我计算 OrderEpsilon2Integral([1,2,3,4]),对吗?
    • 不,你最终可以创建一个for 循环,如果没有矢量化,它可以通过任何标准处理器毫无问题地解决?
    猜你喜欢
    • 2016-08-11
    • 2016-12-13
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2013-05-05
    • 1970-01-01
    • 1970-01-01
    • 2017-07-05
    相关资源
    最近更新 更多