应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » MATLAB解积分方程组
51/1返回列表
查看: 3304  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者caoao22将赠送您 200 个金币

caoao22

木虫 (正式写手)

[求助] MATLAB解积分方程组 已有1人参与

原程序：

r_lower=0.99；
r_point=1；
alpha=0.05;

myfun=@(a,b,r) 1/beta(a,b)*r^(a-1)*(1-r)^(b-1);
myfun1=@(a,b,r) myfun(a,b,r)*r;
myfun2=@(a,b) integral(@(r)myfun(a,b,r),r_lower,1)-alpha;
myfun3=@(a,b) integral(@(r)myfun1(a,b,r),0,1)-r_point;
myfun4=@(r)[myfun2(a,b);myfun3(a,b)];
options=optimset('TolFun',1e-10,'TolX',1e-10);
[x,fval]=fsolve(myfun4,[200;200],options);

程序报错，始终找不到原因，求大神指教

MATLAB解积分方程组
图1.jpg


MATLAB解积分方程组-1
图2.jpg
回复此楼

» 猜你喜欢
1楼 2021-07-22 22:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
在这里用句柄函数来自定义积分方程是不合适的，应创建以函数体的结构形式来自定义积分函数。
x0=[200;200];
options=optimset('TolFun',1e-10,'TolX',1e-10);
[x,fval]=fsolve(@fun,x0,options);
a = x(1)
b = x(2）

function y=fun(x)
a=x(1);b=x(2);
。。。。。。
end
MATLAB解积分方程组-2
123a.jpg
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

caoao22
QQ：2120156492
2楼2021-07-23 21:45:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caoao22

木虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by hzlhm at 2021-07-23 21:45:00
在这里用句柄函数来自定义积分方程是不合适的，应创建以函数体的结构形式来自定义积分函数。
x0=;
options=optimset('TolFun',1e-10,'TolX',1e-10);
=fsolve(@fun,x0,options);
a = x(1)
b = x(2）

functio ...

能不能附上整个程序过程呢，谢谢
赞一下 回复此楼
3楼2021-07-25 00:56:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by caoao22 at 2021-07-25 00:56:42
能不能附上整个程序过程呢，谢谢...

求解方法已上传，希望题主想一想，动动手，不要什么都想要。想要的话，可以商榷。
赞一下 回复此楼
QQ：2120156492
4楼2021-07-25 17:55:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江山印999
牛

发自小木虫IOS客户端
5楼2022-07-28 15:15:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 caoao22 的主题更新
51/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸！ 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
信息提示
关闭
请填处理意见
关闭