当集成贝塞尔函数时发生错误

2视图(30天)
你好,
我试图解决贝塞尔函数的积分方程如附件所示的图像。
我得到一个错误”运算符“+”为操作数的类型不支持function_ha金宝appndle”。
任何人都可以指定如何纠正错误。
代码如下,
gammalv = 72 e - 3;
gammasv = 0.06;
h = 50 e-6;
ν= 0.5;
E = 3000;
R = 0.2 e - 3;
r = 0.0002;
= @ (s)(5 - 12 *ν+ 8 *怒。^ 2 + 2 * s。^ 2 * h。^ 2 +(3 - 4 *ν)* cosh (2 *。* h))。/((3 - 4 *ν)* sinh (2 *。* h) - 2 *。* h);
b = @ (s)(2 *(1 -ν^ 2)。*。* gammasv) / E;
Q = @ (s) 2 *(1 -(ν^ 2))/ (s。* E)。* (1。/ (a + b));
f = @(年代)。* (R . * besselj(0年代。* R) - ((2 * besselj (1, s。* R)) / s))。* Q (s)。* besselj(0年代。* R);
ans =积分(f, 0,正);
谢谢,
Alwar

接受的答案

明星黾
明星黾 2021年9月14日
现在的问题是:
Q = @ (s) 2 *(1 -(ν^ 2))/ (s。* E)。* (1. / ((s) + b (s)));
↑←在这里
这些都需要评估他们的论点(我提供)。
gammalv = 72 e - 3;
gammasv = 0.06;
h = 50 e-6;
ν= 0.5;
E = 3000;
R = 0.2 e - 3;
r = 0.0002;
= @ (s)(5 - 12 *ν+ 8 *怒。^ 2 + 2 * s。^ 2 * h。^ 2 +(3 - 4 *ν)* cosh (2 *。* h))。/((3 - 4 *ν)* sinh (2 *。* h) - 2 *。* h);
b = @ (s)(2 *(1 -ν^ 2)。*。* gammasv) / E;
Q = @ (s) 2 *(1 -(ν^ 2))/ (s。* E)。* (1. / ((s) + b (s)));
f = @(年代)。* (R . * besselj(0年代。* R) - ((2 * besselj (1, s。* R)) / s))。* Q (s)。* besselj(0年代。* R);
ans =积分(f, 0,正);
警告:正或NaN值。
无限的结果仍然有问题。第二项的 “f” (与 “年代” 分子和分母)可能是原因。我离开你的决心。
2的评论
明星黾
明星黾 2021年9月14日
我的荣幸!
返回的结果是 因为结果将 0/0 计算,而 结果的结果 x / 0 计算与 “x” 任何非零值。
如果我的回答帮助你解决你的问题,请 接受 它!

登录置评。

更多的答案(0)

类别

找到更多的在贝塞尔函数帮助中心文件交换

下载188bet金宝搏


释放

R2021a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的