shecheng.m

遗传算法得一种编译程序,不是很正统

a=30*pi/180;
b=45*pi/180;
r=30*pi/180;
w=0;
h=0.125;
g=9.8;
D=0.75*10^-2;     %喷嘴直径，m
p=400*10^3; %工作压力，
q=0.85*pi*0.375^2*10^-2*sqrt(2*g*40000);
n=3;
t1=20;
d=2.31*D^0.831*p^-0.257; %最大水滴直径，m
m=pi*d^3*10^3/6; %水滴质量
f=(1293.31-4.96*t1+0.02807*t1^2-0.000188*t1^3)/1000;   %空气密度，单位kg/m?其中f的大小随空气温度的升高而降低
%q=Q/n;   %Q为管道流量 q为每个喷头的流量，n为喷嘴数目
k=0.18*f*d^2;                 %水滴与空气的摩擦系数kg/m
u=q/(pi*0.375^2*10^-2);                   %水流初速度
t1=20;                   %空气温度
v1=u*sin(a)*cos(r)+w*sin(b);   %v1代表x轴方向的初速度,u为不考虑风速时喷头出口处水的速度w风速通常采用距地2m的速度
                               %a喷嘴处水流初速度与z轴的夹角
v2=u*sin(a)*sin(r)+w*sin(b);   %r喷嘴处水流的初始速度在水平面上的投影与y轴的夹角b水平风速与x轴（正北方向）的夹角                               
v3=u*cos(a);
t=sqrt(m)*(acos(cos(atan(sqrt(k/(m*g))*v3)/exp(h*k/m)))+atan(sqrt(k/(m*g))*v3))/sqrt(g*k);  %水滴在空气中运动的总时间
y=(m/k)*log(v2*k*t/m+1);    %x，y，z水滴在任一时刻的位移分别在三个方向上的投影
x=(m/k)*log(v1*k*t/m+1);    
s=sqrt(x^2+y^2);
disp('shecheng=')
disp(s)