24小时热门版块排行榜    

查看: 2219  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 pang_2020 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 337一志愿华南理工0805材料求调剂 +4 mysdl 2026-03-11 6/300 2026-03-11 23:42 by BruceLiu320
[考研] 一志愿华科071000生物学 338求调剂 +4 九月九里 2026-03-05 4/200 2026-03-11 20:22 by 求调剂zz
[考研] 311求调剂 +6 zchqwer 2026-03-10 6/300 2026-03-11 15:54 by 30660438
[考研] 调剂 +5 调剂的考研学生 2026-03-09 5/250 2026-03-10 22:13 by peike
[考研] 086000生物与医药319分求调剂 +4 Tolkien 2026-03-07 8/400 2026-03-10 21:34 by Tolkien
[考研] 2026考研求调剂-材料类-本科211一志愿985-初试301分 +10 虫友233 2026-03-07 10/500 2026-03-10 17:10 by Demonsssss
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 材料与化工,291,求调剂 +12 咕噜咕噜123123 2026-03-05 13/650 2026-03-10 13:50 by allen-yin
[考研] 材料与化工 一志愿山大 321分 求调剂 +6 每天散步 2026-03-09 7/350 2026-03-10 10:52 by 星空星月
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
[基金申请] 面上项目还需要AI说明吗? +3 liyundong 2026-03-08 3/150 2026-03-09 22:30 by kingkocxr
[考研] 310 070300化学求调剂 +4 扑风铃的猫 2026-03-08 5/250 2026-03-09 19:46 by 30660438
[考研] 301求调剂 +11 朝天椒的双马尾 2026-03-05 11/550 2026-03-08 21:50 by 知足常乐的乐
[考研] 0701-322 求调剂 +3 jiliuxian 2026-03-06 8/400 2026-03-08 19:31 by jiliuxian
[考研] 求调剂 +3 拾柒12。 2026-03-08 3/150 2026-03-08 17:33 by bingxueer79
[考研] 求调剂,不管什么专业,我是可塑造的人才一枚,希望遇到知己老师捞捞我 +4 13102137290 2026-03-06 5/250 2026-03-07 21:21 by lissomchan
[考研] 08工科求调剂 +3 隆LLL 2026-03-06 4/200 2026-03-06 15:26 by 隆LLL
[考研] 268求调剂 +4 刘合华 2026-03-05 4/200 2026-03-06 11:18 by lature00
[考研] 材料与化工354调剂 +4 Lucy-xiao 2026-03-06 7/350 2026-03-06 09:57 by hxuq8429
[考研] 085602 293分求调剂 +3 SivanNano. 2026-03-05 3/150 2026-03-05 19:41 by caszguilin
信息提示
请填处理意见