面向微分方程的数字仿真.docVIP

计算机仿真作业3 。 【作业1】 面向微分方程的数字仿真 1. 一个生长在罐中的细菌简单模型。假设细菌的出生率和当前细菌的总数成正比，死亡率和当前的总数的平方成正比。若 以代表当前细菌的总数，则细菌的出生率可以表示为： 细菌的死亡率可表示为： 细菌数量的变化率可以表示为出生率与死亡率之差。于是该系统可表示为如下微分方程 ，假设b=1，p=0.5，当细菌的总数为100，用simulink计算t的值。 解： 2. 求下列函数最小值（方法不限）。 s.t. 解：将目标函数写成二次函数的形式 其中： H=[1,0;0,1]; f=[-1;-2]; A=[2,3;1,4]; b=[6;5]; lb=zeros(2,1); [xopt,fopt]=quadprog(H,f,A,b,[],[],lb) xopt = 0.7647 1.0588 fopt = -2.0294 3. 用遗传算法工具箱求以下函数最大值（选作题） 解： function z=test(x) z=-x^3+60*x^2-900*x-100 当最大值x=10 函数最大值：4100 第 1 页 共 3 页