24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1265  |  回复: 2

莺呼起2020

铜虫 (小有名气)

[求助] 利用matlab对贝塞尔函数进行积分并求解未知变量已有1人参与

利用matlab对贝塞尔函数进行积分,使用离散相加的方法,程序能够跑通,但是与正确结果有较大误差,烦请大神帮忙看程序是否存在问题,指点迷津。感谢感谢!!
所计算的公式如图所示。

程序如下所示:

syms e0

xl = linspace(1,1000,10001);
xl(end) = [];
dxl = xl(2)-xl(1);

f7  = sum(e0.* ((((besselj(0,xl*a)) - (besselj(0,xl*b))).^2)./xl) )*dxl;
f8 = (t0/(1i*k0*sqrt(t0)))*((1-t1)/(1+t1));

f9 = f7 - f8;
e0 = vpasolve(f9 == 0);

q = vpa(e0);

其中e0为符号数,a、b、k0均为实数,t0、t1为复数,i=sqrt(-1)。求解e0。
感谢感谢!

利用matlab对贝塞尔函数进行积分并求解未知变量


发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

johnliu1983

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
莺呼起2020: 金币+20 2020-10-06 09:28:18
你用的离散积分取的数值点不够多,需要再多一点儿。  另外离散积分你需要用离散积分公式计算,不能这样直接求和。最简单的你可以用trapz这个函数试一试。如果你要求的精度更高,你需要找一下更好的离散积分公式。
2楼2020-09-14 09:10:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

莺呼起2020

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by johnliu1983 at 2020-09-14 09:10:14
你用的离散积分取的数值点不够多,需要再多一点儿。  另外离散积分你需要用离散积分公式计算,不能这样直接求和。最简单的你可以用trapz这个函数试一试。如果你要求的精度更高,你需要找一下更好的离散积分公式。

好的,我试一下,谢谢!!
3楼2020-10-06 09:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 莺呼起2020 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[留学DIY] 10043以后咋出去啊 +3 RuiLove 2024-04-22 3/150 2024-04-27 12:16 by icm639
[教师之家] 博士论文被抄袭 +26 和尚敲小木鱼 2024-04-22 43/2150 2024-04-27 11:14 by liuyuzhou
[论文投稿] AC审稿 +3 马winner 2024-04-26 6/300 2024-04-27 09:08 by 马winner
[硕博家园] 聊天 +12 暮色恋伊人 2024-04-22 13/650 2024-04-27 08:40 by WASM
[硕博家园] 博士白读了 +47 Da_Meng_Zi 2024-04-21 52/2600 2024-04-27 08:25 by shl2112501
[论文投稿] LWT投 +3 AChen92 2024-04-26 3/150 2024-04-26 22:16 by hizifu
[考研] 没学上 +6 季向阳 2024-04-26 12/600 2024-04-26 21:06 by 季向阳
[考博] 申博求助 +4 dskabdh 2024-04-24 11/550 2024-04-26 15:54 by dskabdh
[考研] 381求调剂 +4 小刺猬987654321 2024-04-25 6/300 2024-04-26 10:57 by czl12138
[考博] 取博导收留 5+4 zzb777888 2024-04-20 10/500 2024-04-26 08:52 by polymerfriend
[教师之家] 某种做法不行。说过几遍了。同学还那样做。再那样做就给低分 +4 河西夜郎 2024-04-24 4/200 2024-04-26 08:51 by Quakerbird
[论文投稿] Nature一直在编辑手里,考虑好几天了,是悬了吗 +12 彩虹初见 2024-04-24 12/600 2024-04-25 19:21 by 雪径踏青
[考博] 求博导 +6 好okjh 2024-04-21 10/500 2024-04-25 14:04 by 好okjh
[论文投稿] 一直找不到审稿人 +5 lizhengke06 2024-04-21 6/300 2024-04-25 14:01 by chongdong
[基金申请] 国社科项目,你们学校都限额申报吗? +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[考博] 24年 申博 化学/材料 一作6篇sci +9 wangyp123 2024-04-23 11/550 2024-04-24 19:01 by bangbangbiu
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
[教师之家] 大家访学都是怎么找的啊? +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[教师之家] 找工作过程中碰到好多人问南京某新换校长211的政策,规定如下 +15 ZHONGWU_U 2024-04-21 17/850 2024-04-24 05:12 by juta256
[论文投稿] 研究光催化的,好中的三四区 20+3 sl.0117 2024-04-20 3/150 2024-04-22 09:53 by bobvan
信息提示
请填处理意见