Question about converting integral formula to matlab code

2022 9 23
0 个回答
更新时间：2022 9 23
8 次查看（30 天）

0 个投票

I'm trying to code an integral formula as a sum of discrete data.
That seems correct, but the result is different.
Please let me know if there are any changes to my code.
Please let me know if there is a new way.
The code I made is below.
~~~~~~~~~~~~~~~~~~~~
for r = 1:M
for p = 1:N
for f0 = 1:M0
for p0 = 1:N0-1
PSI(r,p) = PSI(r,p)+...
(1/(pi^2*Rho0^2))*...
K0(f0)*...
d_K0*d_Phi0_rad*...
U(f0,p0)*...
(Rho0*(Rho0-Rho(r)*cos(Phi_rad(p)-Phi0_rad(p0)))/(Rho0^2-2*Rho0*Rho(r)*cos(Phi_rad(p)-Phi0_rad(p0))+Rho(r)^2))*...
(norm([Rho(r)*cos(Phi_rad(p)) Rho(r)*sin(Phi_rad(p))]-[Rho0*cos(Phi0_rad(p0)) Rho0*sin(Phi0_rad(p0))])^2)*...
(exp(2*1i*K0(f0)*((norm([Rho(r)*cos(Phi_rad(p)) Rho(r)*sin(Phi_rad(p))]-[Rho0*cos(Phi0_rad(p0)) Rho0*sin(Phi0_rad(p0))]))-Rho0)));
end
end
end
end

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

Torsten
Torsten 2022-9-23
Any reason why you don't use "integral2" ?
Torsten
Torsten 2022-9-23
I think you should set up your problem for the use with "integral2".
형준 이
형준 이 2022-9-23
I thought about it, but it seems wrong to use that code.
Because U(f,pi) in term 4 is discrete matrix data of M0 by N0.
Therefore, U(f,pi) cannot be substituted into the corresponding code.
"integral2(fun,xmin,xmax,ymin,ymax)"
형준 이
형준 이 2022-9-23
And, I missed that the result of the integral is not a specific value, but the N by M matrix data. * PIS(rho,phi)
Star Strider
Star Strider 2022-9-23
It may be necessary to use integral iteratively with the 'ArrayValued' flag set if you are integrating an array. Check the inner integral before calculating the outer integral to be certain that it is doing what you want.

请先登录，再进行评论。

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

回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

产品

版本

R2022a

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by