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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » MATLAB解积分方程组
51/1返回列表
查看: 3343  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

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

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的回帖
查看全部 5 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学070300 求调剂 +29 哈哈哈^_^ 2026-04-12 29/1450 2026-04-18 15:56 by Equinoxhua
[考研] 一志愿华中农业071010，320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考博] 求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师 +3 猪猪人Zzz 2026-04-12 3/150 2026-04-17 19:10 by 阳阳阳^_^
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 294求调剂 +7 淡然654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[考研] 0854求调剂 +21 门路摸摸 2026-04-15 25/1250 2026-04-17 15:45 by qzxyhcsy
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 5/250 2026-04-17 10:02 by bobvan
[考研] 291求调剂 +9 关忆北. 2026-04-14 9/450 2026-04-16 22:49 by cfdbai
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-16 15:18 by jiahl2024
[考研] 322求调剂 +8 123安康 2026-04-12 15/750 2026-04-16 11:07 by Espannnnnol
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 通信工程求调剂！！！ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 105500药学求调剂 +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[考研] B区0809 ，数一英一，290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
信息提示
关闭
请填处理意见
关闭