24小时热门版块排行榜    

查看: 1495  |  回复: 4

xujunp

金虫 (正式写手)

[求助] matlab含参数二重积分计算问题

%这是我的程序,结果是是带int的式子,求大神指点下
syms ro Z a c S0  n0 Aj Po2 lambda Sph
Z0=5; S0=0.025; n0=10^25; Aj=0.002; Po2=21278; lambda=0.053;
R=sqrt((a-ro)^2+(c-Z)^2);
ne=n0*exp((-(Z-Z0)^2-ro^2)/2*S0^2);
Sph=ne*Aj*Po2^2*exp(-lambda*Po2*R)/(4*pi*R);
A=int(int(Sph,Z,0,0.1),ro,0,4);
%积分不出结果,参数是a和c
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

长_青

木虫 (知名作家)

咱们好像是老乡?
2楼2017-11-01 16:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujunp

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 长_青 at 2017-11-01 16:14:13
咱们好像是老乡?

这都能知道,不过这个怎么处理呢
3楼2017-11-01 16:25:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

这种 一般很难有解析解的,可以考虑数值求解。

P.S.  你的参数a、c 是已知的 还是 未知的。 若已知可以代入直接计算,否则可以考虑用匿名函数来表述(参数传递)。
数值计算
4楼2017-11-06 09:55:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujunp

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 独孤神宇 at 2017-11-06 09:55:13
这种 一般很难有解析解的,可以考虑数值求解。

P.S.  你的参数a、c 是已知的 还是 未知的。 若已知可以代入直接计算,否则可以考虑用匿名函数来表述(参数传递)。

可以是已知的,如果是已知的,积分上限的函数可以用数值解吗

发自小木虫Android客户端
5楼2017-11-06 10:23:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xujunp 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 7/350 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +4 SCSIOyxguo 2026-09-03 7/350 2026-09-05 10:39 by Honglianlian
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见