computing finite time grammian

Question

Srighakollapu Vali M 2018-11-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/432969-computing-finite-time-grammian

回答： madhan ravi 2018-11-30

Hi When I am trying to find finite grammian I getting the following error:

-----------------------------------------------------
>>integral(f,0,1)
Matrix dimensions must agree.
Error in @(tau)expm(A.*tau)*(b*b')*expm(A'.*tau)
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);
 
________________________________________
My A=[-1/2 0;0 -1]
b=[1;1]
f = @(tau) expm(A.*tau)*(b*b')*expm(A'.*tau)