Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2008-12-06
年龄: 37
帖子: 2
声望力: 0 ![]() |
![]()
[catchfish.m]
function y=catchfish(x) syms k f10 f20 f30 f40 total; %f10,f20,f30,f40表示4个鱼龄的鱼群的初始量,k表示捕捞强度。 global f10 f20 f30 f40 total k; total表示捕捞总量。 %求解方程组 x1=dsolve('Dx1=-0.7*x1','x1(0)=f10'); t=1;f20=subs(x1); x2=dsolve('Dx2=-0.7*x2','x2(0)=f20'); t=1;f30=subs(x2); x31=dsolve('Dx31=-(0.7+k)*x31','x31(0)=f30'); t=2/3;f31=subs(x31); x32=dsolve('Dx32=-0.7*x32','x32(2/3)=f31'); t=1;f40=subs(x32); x41=dsolve('Dx41=-(0.7+3*k)*x41','x41(0)=f40'); t=2/3,f41=subs(x41); x42=dsolve('Dx42=-0.7*x42','x42(2/3)=f41'); nn=1.1*10^5*(0.5*f31+f41); eq=f10-nn*1.2*10^11/(1,2*10^11+nn);%求一龄鱼的存活数量 s=solve(eq,f10);f10=s(2); sym t; t3=subs(subs(int(k*x31,t,0,2/3))); t4=subs(subs(int(3*k*x41,t,0,2/3))); total=15.6*t3+20.8*t4; k=x; y=subs(-total); [bestcatchfish.m] global a10 a20 a30 a40 total; [k,mtotal]=fminbnd('catchfish',0,50); ezplot(total,0,50); xlabel('捕捞强度系数k'); ylabel('总捕捞量'); title('捕捞强度-总捕捞量关系图'); format long; k total=-total; a10=eval(a10); a20=eval(a20); a30=eval(a30); a40=eval(a40); format short; clear 运行显示 Undefined function or variable 'k' 或Undefined function or variable ‘X3’ 我想求K的最大值 望高手赐教 |
![]() |
![]() |
![]() |
|
|
![]() |
||||
主题 | 主题作者 | 版面 | 回复 | 最后发表 |
[分享]给大家传一个图像空间变换和直方图变换的数学推导 | johnny8088 | MATLAB论坛 | 0 | 2009-01-31 03:26 |
有哪位做过BEESSEL法数字下变频? | yqgame365 | MATLAB论坛 | 0 | 2008-12-07 23:53 |
[求助]在SIMULINK里,如何将inegrator上方的箭头移到下面去 | yijianmei | MATLAB论坛 | 1 | 2008-07-24 10:15 |
求助~!matlab 初级问题~~急要~~谢谢啦 | embyg | MATLAB论坛 | 0 | 2008-07-08 15:35 |
Simplify Instrument Communication: Using Instrument Drivers with the Instrument Contr | TechnicalArticles | MATLAB技术文章 | 0 | 2008-01-06 16:32 |