Question about converting integral formula to matlab code

1 次查看（过去 30 天）
형준 이
형준 이 2022-9-23
评论： Star Strider 2022-9-23
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更早的评论
형준 이
형준 이 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 个）

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

产品


版本

R2022a

Community Treasure Hunt

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

Start Hunting!

Translated by