Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2008-12-06, 20:04   #1
killer166
初级会员
 
注册日期: 2008-12-06
年龄: 37
帖子: 2
声望力: 0
killer166 正向着好的方向发展
默认 [求助]求助,程序的错误

[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的最大值
望高手赐教
killer166 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码


相似的主题
主题 主题作者 版面 回复 最后发表
[分享]给大家传一个图像空间变换和直方图变换的数学推导 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


所有时间均为北京时间。现在的时间是 13:46


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.