24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3039  |  回复: 4
【悬赏金币】回答本帖问题,作者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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
5楼2022-07-28 15:15:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 caoao22 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见