MATLAB解积分方程组
原程序：
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
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
在这里用句柄函数来自定义积分方程是不合适的，应创建以函数体的结构形式来自定义积分函数。
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
，
能不能附上整个程序过程呢，谢谢
求解方法已上传，希望题主想一想，动动手，不要什么都想要。想要的话，可以商榷。
牛