24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者caoao22将赠送您 200 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

caoao22

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2064
帖子: 337
在线: 35.5小时
虫号: 4098136
注册: 2015-09-24
性别: GG
专业: 概率论与随机分析

[求助] 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);

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

图1.jpg

图2.jpg

回复此楼

» 猜你喜欢

1楼2021-07-22 22:55:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

caoao22

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2064
帖子: 337
在线: 35.5小时
虫号: 4098136
注册: 2015-09-24
性别: GG
专业: 概率论与随机分析

送红花一朵

引用回帖:

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的回帖

相关版块跳转我要订阅楼主 caoao22 的主题更新

返回列表