96年数模竞赛a题:最优捕鱼策略问题

来源:百度知道 编辑:UC知道 时间:2024/05/27 11:19:49
求这道题的详细答案,包括问题的假设,分析;模型的建立等

一定要有matlab的程序代码或者c语言的也可以
如果文件太大,请发到我的邮箱dinckywyle@163.com,不胜感激

%保存为M文件,文件名为"smax1.m",在主窗口中输入"[ss3,ss4,st]=smax1()"命令来运行!

%k为捕捞强度,通过改变k来改变结果。

function [ss3,ss4,st]=smax1()
s10=122e9;
s20=29.7e9;
s30=10.1e9;
s40=3.29e9;
s0=[s10,s20,s30,s40];
a1=0.8;
ss3=0;
ss4=0;
a=1.109e5;
t1=1;
syms t;
d=10;
k=10;
h=1e5;
while d>=0
sb1=s10*exp(-a1*(2/3)*t1);
sb2=s20*exp(-a1*(2/3)*t1);
sb3=s30*exp(-(a1+0.42*k)*(2/3)*t);
sb4=s40*exp(-(a1+k)*(2/3)*t);

%sb3q=s30*exp(-(a1+0.42*k)*(2/3)*t1);

s1=sb1*exp(-(1/3)*t1);
s2=sb2*exp(-(1/3)*t1);
%s3=sb3*exp(-(1/3)*t);
%s4=sb4*exp(-(1/3)*t);

s3q=s30*exp((-(a1+0.42*k)*(2/3)-(1/3))*t1);
s3=s30*exp((-(a1+0.42*k)*(2/3)-(1/3))*t);
s4=s40*exp((-(a1+k)*(2/3)-(1/3))*t);

sp3=quad