Error message from integral calculation
显示 更早的评论
Hi, I'm calculating the survival function (lx) that is expressed with an integral such as
to calculate this I have created following program
% main.m
param_1938 = [2.06441912000572E-07/1000,0.197642212387667/100000,1.23947876070978/10];
mu_x=@(t) f_lx(param_1938);
l_x = exp(-integral(@(t) mu_x(t),0,106));
%f_lx.m
function res=f_lx(param)
x=(0:1:106)';
a=param(1);
b=param(2);
c=param(3);
res = zeros(size(x));
ind = x>100;
res(ind) = a+b*exp(c*100)+(x(ind)-100)*0.001;
res(~ind) =a+b*exp(c*x(~ind));
end
This result in the error message
Error using integralCalc/finalInputChecks (line 515)
Output of the function must be the same size as the input. If FUN is an array-valued integrand, set the 'ArrayValued' option to true.
Error in integralCalc/iterateScalarValued (line 315)
finalInputChecks(x,fx);
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);
I'm not seeing where this goes wrong. Please can anyone shed some light on this??
回答（1 个）
Steven Lord
2018-4-7
0 个投票
You want to integrate with respect to x, not with respect to your parameters. See the "Parameterized Function" example on the documentation page for the integral function and the "Parameterizing Functions" page linked in the Topics section on that page to see how to integrate a parameterized function like yours.
2 个评论
Orongo
2018-4-9
Steven Lord
2018-4-9
And you changed f_lx to accept the first input argument as x and does not overwrite the array integral passed into your integrand function with a vector of a fixed length?
类别
在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息
2018-4-9
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
