| 查看: 1855 | 回复: 36 | ||
wgdd木虫 (正式写手)
|
[求助]
求助一个数值积分问题，用matlab的quadgk函数来计算，谢谢！
|
|
现求解一个半无限震荡积分问题： 被积函数为s*besselj(0,s*R)./(s.^2-k0^2), 积分区间为[0,inf], 其中R=6, k0=10, 均为已知常数。我打算用两种方法来计算这个问题。(解析解和数值解) 方法一，解析解 对于这样一个积分，可以查到数学手册中积分结果为pi*i/2*besselh(0,1,k0*R). 经过计算可以得到结果为-0.0744。 计算程序 % analytic method k0=10; R=6; p_inc=real(pi*1i/2*besselh(0,1,k0*R)); 方法二，采用matlab中的quadgk函数进行计算，由于k0为奇点，所以将积分区间分为两块。[0,k0]和[k0,inf]，计算结果为0.0390 计算程序 % numerical method k0=10; R=6; p_f=@(s)(s).*besselj(0,s*R)./(s.^2-k0^2); p1=quadgk(p_f,0,k0); p2=quadgk(p_f,k0,inf); p_incN=p1+p2; 求各位大牛支招，为啥两种方法差别这么大呀~正确的数值积分方法应该是怎么样的呢？十分感谢！@月只蓝@beefly |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有213人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有111人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有68人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有1人回复
 |
2楼2017-08-17 17:25:05
 |
3楼2017-08-17 19:47:06
|
|
4楼2017-08-17 22:59:40
|
|
5楼2017-08-18 01:07:01
|
|
6楼2017-08-18 02:06:45
 |
7楼2017-08-18 13:54:45
 |
8楼2017-08-18 15:56:36
|
|
9楼2017-08-18 16:43:44
|
|
10楼2017-08-19 00:30:27
10