3、lve('Dy=y+2*x','y(0)=1','x');ezplot(y,[01])Result:1.2Code:clear[ty]=ode23('wf',[010],[01]);plot(t,y(:,1))M文件functionyd=wf(t,y)yd=[y(2);-y(1)*cos(t)]Result:以下为部分数值解和图像3.Rossler微分方程组:当固定参数b=2,c=4时,试讨论随参数a由小到大变化(如a∈(0,0.65))而方程解的变化情况,并且画出空间曲线图形,观察空间曲线是否形成混沌状?Code:clearglo
4、balafora=0:0.05:0.65;[tx]=ode45('P3',[010],[111]);plot(t,x)holdonendM文件:functionequ=fun(t,x)globalab=2;c=4;equ=[0-1-1;1a0;x(3)0-c]*x+[0;0;b];endResult:可见,图像呈混沌状。4.Apollo卫星的运动轨迹的绘制Code:clearglobaluu1u=1/82.45;u1=1-u;[tx]=ode45('P4',[010],[1.200-1.04935751]);plot(t,x)M文件
5、:functioneqt=P4(t,x)globaluu1r1=sqrt(((x(1)+u))^2+x(3)^2);r2=sqrt(((x(1)-u1)^2)+x(3)^2);eqt(1,1)=x(2);eqt(2,1)=2*x(4)+x(1)-u1*(x(1)+u)/r1^3-u*(x(1)-u1)/r2^3;eqt(3,1)=x(4);eqt(4,1)=-2*x(2)+x(3)-u1*x(3)/r1^3-u*x(3)/r2^3;Result:5.盐水的混合问题一个圆柱形的容器,内装350升的均匀混合的盐水溶液。如果纯水以每秒14
6、升的速度从容器顶部流入,同时,容器内的混合的盐水以每秒10.5升的速度从容器底部流出。开始时,容器内盐的含量为7千克。求经过时间t后容器内盐的含量。Code:clear[tx]=ode23('P5',0:0.1:10,[7350]);count=1;fort=0:0.1:10plot(t,x(count,1)/x(count,2),'--');xlabel('ʱ¼ä');ylabel('º¬ÑÎÁ¿£¨kg/L£©');holdoncount=count+1;endM文件:functionequ=P5(t,x)equ=[-10.5
7、*x(1)/x(2);3.5];result:9.River-bay系统水污染问题一条河流和河湾与大湖相连,位于湾上游的小河是造成湾污染的主要因素,另有一座铝厂恰好建在湾旁,也造成污染。当湾中污染物平均浓度达到1.6mg/l时,铝厂将被迫暂时关闭。假使该湾的容量为4,000,000公升,流入和流出河湾的水流速度均为40,000公升/天,若当前该河湾水中的污染物浓度为0.8mg/l,河水中污染物的浓度为0.5mg/l。要求1)建立湾中水污染状况的模型;2)计算30天后该河湾水的污染物浓度;3)该河湾水的污染物浓度是否能达到一个稳定值?
8、4)如将4,000,000mg污染物瞬间排入河水中,求铝厂必须关闭多长时间?5)列出并讨论影响河湾污染的模型中未考虑到的因素至少四种。(1)假设t时刻,河湾中污染物浓度为W,河水中污染物浓度为W1则浓度变化量等于流入的污染物减去流出的