Warning: Reached the limit on the maximum number of intervals in use. Approximate bound on error is 2.4e+02. The integral may not exist, or it may be difficult to approximate numerically to the requested accuracy.
16 ビュー (過去 30 日間)
古いコメントを表示
Hi all, I'm try to solve the integral in the code and I have the error. Any suggestions? The function seems to not have singularities.
Roi=916;
Row=1024;
r=50;
d=50;
H0=50;
va=9.46
vi=30.99
b=(H0-vi*vi/(2*9.806))*2
a=(vi/va-1)*(3/2*Roi/Row*b)
N=1000
t=linspace(0,30,N)
for i=1:length(t)
E1=@(x)besselj(0,r/d.*x).*sqrt(x.*tanh(x)).*x.*sin(sqrt(9.81*t(i)^2/d.*x.*tanh(x))).*(sin(a/d.*x)-a./d.*x.*cos(a/d.*x))./(a/d.*x).^3
E(i)= integral(E1,0,inf)
end
E2=-2*(a^3)/(d^3)*sqrt(2*H0/d-a/d*b/a)/(1+2/3*Row/Roi*a/b)*E*d;
max(E2)
plot(t,E2)
0 件のコメント
回答 (1 件)
Torsten
2019 年 2 月 22 日
E1=@(x)besselj(0,r/d.*x).*sqrt(x.*tanh(x)).*x.*sin(sqrt(9.81*t(i)^2/d.*x.*tanh(x))).*(sin(a/d.*x)-a./d.*x.*cos(a/d.*x))./(a/d.*x).^3
At x=0, you divide by 0 here.
参考
カテゴリ
Help Center および File Exchange で Bessel functions についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!