版块导航: 正在加载中...

登录注册

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

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

hzlhm

至尊木虫 (著名写手)

应助: 387 (硕士)
金币: 18110.5
红花: 53
帖子: 2879
在线: 603.7小时
虫号: 1517335
注册: 2011-11-30
性别: GG
专业: 常微分方程与动力系统

【答案】应助回帖

感谢参与，应助指数 +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

123a.jpg

赞一下

回复此楼

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

caoao22

QQ：2120156492

2楼2021-07-23 21:45:00

已阅回复此楼关注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的回帖