Double integral with complex numbers in a loop

Question

Maria445 2017-5-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/338395-double-integral-with-complex-numbers-in-a-loop

评论： Maria445 2017-5-6

Hi everyone, I'm trying to run this code (probably there are better ways to write it, but I'm not a pro at Matlab yet).

params is a vector of 4 parameters, Nstep =365*25, and p_sims is a matrix with Nstep+1 columns and as many raws as the number of simulations (1000 in my case).

I need to integrate a function in two variables, u and s and then use the result to create another matrix.

The problem is that I need to insert in the integrand a vector (the jjj column of p_sims), but this gives me an error (Inner matrix dimensions must agree), whether I use .* to multiply by something containing p_sims, or * How can I solve this?

(I tried with a loop to compute the jjj column of CPrice, but it took ages to run the code, so I was wondering if there was a more efficient way)

for jjj=Nstep+1:-1:1
        CFun = @(s,u) exp(-r*(s-jjj)).*p_sims(:,jjj).^gamma.*floor^(1-gamma)/pi.* ...
            exp(-1i*u.*log(p_sims(:,jjj)/floor)).*exp(1i*params(2).*s.*(-u-1i*gamma));   
CPrice(:,jjj) = X.*real(integral2(CFun,jjj,Nstep+1,0,inf)); 
phi(:,jjj) = p_sims(:,jjj).*(exp((params(2)-r)*(Nstep+1-jjj))*(1-(params(4)+1/2*params(3)^2)/params(1))^ ...
    (-params(1)*(Nstep+1-jjj))-1)  + CPrice(:,jjj);
end

3 个评论
显示 1更早的评论隐藏 1更早的评论

Maria445 2017-5-2

More precisely, with .* the error is "Matrix dimensions must agree", while using * the error is "Inner matrix dimensions must agree"

Maria445 2017-5-2

Cattura.PNG

This is the screenshot

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2017-5-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/338395-double-integral-with-complex-numbers-in-a-loop#answer_265398

The only * operation in your line of code is -r*(s-jjj) . We can see from your code that jjj will be a scalar.

Your s argument is input from the integral2() call. We know from https://www.mathworks.com/help/matlab/ref/integral2.html#inputarg_fun that

"The function fun must accept two arrays of the same size and return an array of corresponding values. It must perform element-wise operations."

so we know that s is a 2D array, of unspecified and variable size.

We are not given any information about the size of r. But clearly since the size of s is variable, if r is not a scalar, there will come a time where size(r,2) ~= size(s,1) . That would lead to the error about inner dimensions not matching. If r were a scalar, there would be no problem at that point, and since there is no other * operation in the line of code, we deduce that whatever r is, it is not a scalar.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Walter Roberson 2017-5-5

One option that could be experimented with (it might not be any faster at all) would be to convert your 2D integral into a pair of 1D integral() in which the outer integral calls upon the inner one. The inner one could be coded with 'ArrayValued', 'on'. To take care of the fact that the outer integral() will be passed a vector, you might need to arrayfun() the inner integral over its inputs.

... Or possibly you could just calculate your function over a mesh of s and u values, getting out a 3D array (s, u, p_sims column) and then trapz() along the appropriate dimensions. That would not be adaptive, but it might possibly be accurate enough for your purposes.

Maria445 2017-5-6