当集成贝塞尔函数时发生错误
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
0评论
接受的答案
明星黾
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,正);
无限的结果仍然有问题。第二项的
“f”
(与
“年代”
分子和分母)可能是原因。我离开你的决心。
。
2的评论
明星黾
2021年9月14日
我的荣幸!
返回的结果是
南
不
正
因为结果将
0/0
计算,而
正
结果的结果
x / 0
计算与
“x”
任何非零值。
如果我的回答帮助你解决你的问题,请
接受
它!
。