本文运用伯努力方程求解喷头处水喷洒速度,在解决喷洒半径时,考虑到了摩擦力对于快速运动的小水滴的作用。 喷洒水滴的半径随机性很大,不可能为一固定的常数,但在模型中认为是一个常数,是为了使问题得到简化。单个水滴的半径与喷洒距离之间存在着一定的关系,当喷洒的距离越远,水滴的半径越小。水滴半径越小,在喷洒过程中的损失越大。根据农田水利学中管道灌溉系统知识,在喷灌半径50-60%的范围内,即使各喷头水量不重叠,灌水量也能充分满足植株生长。而在60%以外,即喷头射程的后40%部分,随着距离的增大,水量越来越小,便不能满足植物的生长需要,而我们设计的喷灌方案中那些水量小的部分大部分分布在田的外部,具有一定的合理性。
附录: 解微分方程模型的源程序: function dx=rigit(t,x) dx=zeros(2,1); dx(1)=x(2); dx(2)=-3*0.6*1.25*x(2)^2/8/0.0025/1000;
function dy=rigit(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=9.8-3*0.6*1.25*y(2)^2/8/0.0025/1000;
[t,x]=ode45(@rigit,[0,4],[0,25.2607*sqrt(2)/2]) plot(t,x(:,1),'*')
喷头个数为n=2时的源程序: >> ezplot('(x-15)^2+(y-30)^2-400',[0,30,0,80]) >> ezplot('(x-15)^2+(y-50)^2-400',[0,30,0,80]) >> hold on >> ezplot('(x-15)^2+(y-30)^2-400',[0,30,0,80]) >> plot(15,30,'*') >> plot(15,50,'*')
喷头个数为n=1时的源程序: >> ezplot('(x-15)^2+(y-13.2288)^2-400',[0,30,0,80]) >> hold on >> ezplot('(x-15)^2+(y-13.2288*3)^2-400',[0,30,0,80]) >> plot(15,13.2288,'*') >> plot(15,13.2288*3,'*') 喷头个数为n=3时源程序: >> ezplot('(x-15)^2+(y-50)^2-400',[0,30,0,80]) >> hold on >> ezplot('(x-15)^2+(y-30)^2-400',[0,30,0,80]) >> ezplot('(x-15)^2+(y-40)^2-400',[0,30,0,80]) >> plot(15,40,'*') >> plot(15,30,'*') >> plot(15,50,'*')
首页 上一页 1 2 下一页 尾页 2/2/2