问题描述:
请问我这个matlab程序哪里错了啊,
Q=input('请输人泄漏源强(mg/s):Q=');
H=input('输人计算风速(m/s):H=');
d=input('请输人计算精度(m):d=');
Z0=input('请输人地面粗糙长度(m):Z0:');
Yj=meshgrid(50:d:1000,-100:d:100);%定义解空间和计算精度
byO=0.08*(1+0.0001)*(-1/2);%计算Y轴向的基本扩散参数
bzO=0.06*(1+0.O015)*(-1/2);%计算z轴向的扩散参数
by:byo.(1+0.38 Zo);%对Y轴向的扩散参数按地面粗糙长度进行修正.
= (2.53-0.13 log( )).(0.55+0.042 log
( )).(-1)..(0.35-0.03 log( ));%按地面粗糙长度计算z轴向的扩散参数修正系数
= bzO."fz"%对 轴向的扩散参数按地面粗糙长度进行修正.
%高斯扩散浓度模拟计算%
tempyl= -Y.Y./by./by./2;
tempy2=2.718282.^(tempy1);
C:Q/pi/u ((69,.).^(-1)).tempy2
Cs=input('请输人所有求解浓度(mg/m3):');
%所输人数以”[”和”]”作为开头和结束.
contour( ,Y,C,Cs);%以 ,Y,C变量分别作为的轴,绘制扩散浓度分布图
shading interp;
colorbar;
d;
)dabe1(' 轴向距离(m)');
ylabel('Y轴向距离(m)');
title('气体扩散下风向浓度分布图')
Q=input('请输人泄漏源强(mg/s):Q=');
H=input('输人计算风速(m/s):H=');
d=input('请输人计算精度(m):d=');
Z0=input('请输人地面粗糙长度(m):Z0:');
Yj=meshgrid(50:d:1000,-100:d:100);%定义解空间和计算精度
byO=0.08*(1+0.0001)*(-1/2);%计算Y轴向的基本扩散参数
bzO=0.06*(1+0.O015)*(-1/2);%计算z轴向的扩散参数
by:byo.(1+0.38 Zo);%对Y轴向的扩散参数按地面粗糙长度进行修正.
= (2.53-0.13 log( )).(0.55+0.042 log
( )).(-1)..(0.35-0.03 log( ));%按地面粗糙长度计算z轴向的扩散参数修正系数
= bzO."fz"%对 轴向的扩散参数按地面粗糙长度进行修正.
%高斯扩散浓度模拟计算%
tempyl= -Y.Y./by./by./2;
tempy2=2.718282.^(tempy1);
C:Q/pi/u ((69,.).^(-1)).tempy2
Cs=input('请输人所有求解浓度(mg/m3):');
%所输人数以”[”和”]”作为开头和结束.
contour( ,Y,C,Cs);%以 ,Y,C变量分别作为的轴,绘制扩散浓度分布图
shading interp;
colorbar;
d;
)dabe1(' 轴向距离(m)');
ylabel('Y轴向距离(m)');
title('气体扩散下风向浓度分布图')
问题解答:
我来补答展开全文阅读