24小时热门版块排行榜    

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 329  |  回复: 0
【悬赏金币】回答本帖问题,作者engine9817将赠送您 5 个金币

engine9817

铜虫 (小有名气)

[求助] 求设计一个被积分函数中需要解方程的积分程序

如果有一个被积方程为G(x,y),积分区间假设为(1,2)
其中y的值与x是由另一个函数决定如f(x,y)=0; 这个式子不能直接得到y的解析式,复杂点的话,可能一个x有多个y值,但只有一个合适,如正的或实的。
这样在求积G(x,y)时,对应于每个x都要求解一个y值,还可能要对y值进行筛选。
我一开始在function文件里面设计了G(x,y)的程序.方便起见,我们假设
f(x,y)=x-y^2=0; 假设这里把y解析表示,要求用fsolve解出y

function integrandfun=G(x)

y=fsolve(@(t) x-t.^2, 1);  

integrandfun=x.^2+y;

主函数中调用 这个方程积分如

In=quadgk(@(x) integrandfun, 0,1);

结果是: 3.558078204809361与直接积分x.^2+sqrt(x)的结果3.552284749830793有差别
提示:
警告: Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm
instead.
> In fsolve at 285
  In fkk at 3
  In @(x)fkk(x)
  In quadgk>evalFun at 330
  In quadgk>f1 at 348
  In quadgk>vadapt at 249
  In quadgk at 188
  In figures at 12

No solution found.

fsolve stopped because the last step was ineffective. However, the vector of function
values is not near zero, as measured by the default value of the function tolerance.

<stopping criteria details>

请问有进一步的优化吗?事实上如果把f(x,y)变得复杂点,积分存在更多问题。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 engine9817 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程与技术(学硕)297求调剂 5+5 @z张同学 2024-03-24 13/650 2024-03-29 15:12 by fyyybin
[考研] 290工学材料类求调剂 +4 L需要调剂 2024-03-29 4/200 2024-03-29 15:03 by 凭楼听雨
[考研] 330材料与化工求调剂 +5 哈—hhh 2024-03-28 7/350 2024-03-29 12:56 by 伊人葭儿
[教师之家] 毕业6年,感觉科研的路越走越窄了,不知道怎么办 +17 余姚和谐 2024-03-27 17/850 2024-03-29 11:01 by zju2000
[考研] 304材料专硕求调剂 +3 烟云墨雨飞啊 2024-03-28 3/150 2024-03-29 10:45 by feizhang1122
[考研] 335化工求调剂合适院校 +3 给你补个蛋 2024-03-28 3/150 2024-03-29 02:39 by 储能小站cc
[基金申请] MDPI的期刊 +5 zju2000 2024-03-26 5/250 2024-03-28 23:00 by zzahkj
[考研] 282求调剂 +8 朱朱女朱 2024-03-24 14/700 2024-03-28 20:06 by 麦田里的守望者
[精细化工] 有机硅改性环氧树脂 5+3 广州小苏 2024-03-22 3/150 2024-03-28 16:03 by CDMACD
[考研] 315求调剂 +4 无昵称- 2024-03-28 4/200 2024-03-28 15:20 by vivian_sfw
[考研] 309求调剂 +3 化工qxam 2024-03-23 3/150 2024-03-28 10:35 by hangyuan13
[有机交流] 溴素溴化反应生成溴苄 100+4 tanjishuang 2024-03-25 6/300 2024-03-28 09:37 by 高树桩
[考研] 287求调剂 +9 王杨具牢 2024-03-24 11/550 2024-03-27 20:55 by siwuhen
[有机交流] 胺类化合物萃取 30+3 Zephyr@ 2024-03-25 3/150 2024-03-27 18:15 by organic2020
[考研] 284求调剂 +5 她与残局皆遗憾 2024-03-26 5/250 2024-03-27 13:34 by zhyzzh
[考研] 求调剂 +7 黎明~太阳 2024-03-26 8/400 2024-03-27 09:53 by 黎明~太阳
[论文投稿] 研0,大佬们SCI在哪看呀 5+3 王二狗demo 2024-03-26 5/250 2024-03-27 09:40 by topedit
[考研] 化学工程与技术320求调剂 +7 Wjrjxmt 2024-03-26 8/400 2024-03-27 08:50 by 实验室@
[考研] 317求调剂 +4 (?????)????? 2024-03-25 5/250 2024-03-25 23:58 by DJ35
[考研] 293食品科学与工程求调剂 +3 懒懒er.。 2024-03-22 5/250 2024-03-23 12:06 by 打个小呼噜
信息提示
请填处理意见