请教各位大侠，这个用MATLAB怎样编程？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

现在有个定积分的隐函数求解，r为积分变量，a=0.99,b=5,pe=10,pw=4,rw=0.1,re=100,要求pe*，请教各位高手matlab 的fsolve 可以直接求解吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB matla atlab 怎样编程 各位大侠 怎样编程

为什么要用fsolve那么复杂？这个积分可以自己手算的诶……

把所有p和r带下标的玩意统统视为常数，也就是被积分函数是
\[\frac{2b - ar + r + \sqrt{ar - r - C + 4abrD}}{2br}\]
然后，将分子拆开，从而变成
\[\frac{1}{r} + \frac{1 - a}{2b} + \frac{\sqrt{(a - 1 + 4abD) r - C}}{2br}\]
前两项都可以求出积分来，下面考虑第三项

那么记常数\[\\k = a - 1 + 4 abD\\\]，并做变量代换得\[\\t = \sqrt{(a - 1 + 4abD)r - C}\\\]
因此
\[r = \frac{t^2 + C}{k}\]
\[dr = 2ktdt\]
上面被积函数第三项变成\[\int 2k^2 \frac{t^3}{t^2 + C} dt\]
这个积分你要是算不出我也没办法了……

BTW，有时候多算一点，matlab毕竟是数值运算，直接解误差一般比较大~

非常感谢！你这个处理方法对我的文章非常有借鉴之处！我准备使用！！
程序昨晚做出来了
function ff=djf(p1)
pe=10;pw=4;la=0.001;lc=0.2;a=(la-lc)/la;b=1/la;re=100;rw=0.1;

  g=@(r,p1)((2*b-a*r+r+sqrt((a*r-r-b*(p1-pw)/log(re/rw)).^2+4*a*b*r*(p1-pw)/log(re/rw)))./(2*b*r));
   ff=pe-pw-quadl(@(r)g(r,p1),0.1,100);
end

   
   

%主程序调用文件
% [x,feval]=fsolve(@djf,6,optimset('Display','off','TolX',1e-8,'TolFun',1e-8))