24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3338  |  回复: 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的回帖

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朵)

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的回帖
5楼2022-07-28 15:15:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 caoao22 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22专硕求调剂 +10 haoyun上岸 2026-04-11 12/600 2026-04-16 22:21 by 猪会飞
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 310求调剂 +16 666真好 2026-04-11 18/900 2026-04-15 13:28 by 黑科技矿业
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 材料085601调剂 +32 何润采123 2026-04-10 34/1700 2026-04-14 08:47 by 木木mumu~
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
信息提示
请填处理意见