integral(exp(double integral)) ?===>integral of riemann

2 次查看(过去 30 天)
Hi,
i have an equation to 3 unknown (x,y theta)
i try this code
a =2.2;
b =0:10:100;
r = 50;
for i = 1:length(b)
bin=100;
resultat=0;
for k=1:bin
y=k*r/bin;
eq1= @(x,theta) exp(-b(i).*(sqrt(x.^2.*(cos(theta)).^2+(x.*sin(theta)-y).^2).*x)) ;
In2=integral2(eq1,0,r,0,2*pi);
resultat=resultat+ exp(-a* (1-In2))*2*y/(r^2);
end
end
==== >result does not work correctly
  2 个评论
John D'Errico
John D'Errico 2019-7-20
Without even looking more carefully at the numbers, as soon as I see this start:
exp(-3.7154e+05*d(i).*(
I will predict your problem is in the form of exponential underflows. The result will be numerical garbage.
Marwen Tarhouni
Marwen Tarhouni 2019-7-20
编辑:Marwen Tarhouni 2019-7-20
Without this number ,i removed it. It's important to me to solve this equation ,i try this code as solution of my equation but i dont know if it's correctly.

请先登录,再进行评论。

采纳的回答

Sulaymon Eshkabilov
Hi,
In your code, the variable d is assigned instead of b on line 2: d= 0:10:100; When this is fixed then everything works ok.
You can also start your simulation with the command: clearvars
Good luck

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Programming 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by