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


所有时间均为北京时间。现在的时间是 15:41


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