matlab的问题.之前使用这个公式,我采用[x,y] =meshgrid(0:1:2000,-200:1:200);by0 =0.32*x.*(1+0.0004.*x).^(-1/2);bz0 =0.24*x.*(1+0.0001.*x).^(-1/2);tempy1 =(-1/2).*((y.*y./by0./by0)+(H^2.by0.bz0));tempy2 =2.718282.^(tempy1);c =Q/pi/u.*((by0.
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/25 12:24:59
matlab的问题.之前使用这个公式,我采用[x,y] =meshgrid(0:1:2000,-200:1:200);by0 =0.32*x.*(1+0.0004.*x).^(-1/2);bz0 =0.24*x.*(1+0.0001.*x).^(-1/2);tempy1 =(-1/2).*((y.*y./by0./by0)+(H^2.by0.bz0));tempy2 =2.718282.^(tempy1);c =Q/pi/u.*((by0.
matlab的问题.
之前使用这个公式,我采用
[x,y] =meshgrid(0:1:2000,-200:1:200);
by0 =0.32*x.*(1+0.0004.*x).^(-1/2);
bz0 =0.24*x.*(1+0.0001.*x).^(-1/2);
tempy1 =(-1/2).*((y.*y./by0./by0)+(H^2.by0.bz0));
tempy2 =2.718282.^(tempy1);
c =Q/pi/u.*((by0.*bz0).^(-1)).*tempy2;
meshc(x,y,c);
xlabel('x轴向距离(m)');
ylabel('y轴向距离(m)');
title('铅')
这个编程.数据代入啥的我自己输入.
但是建模的时候,我对这个式子进行了改动
变成了这个样子,就是说好几个一样的式子,但是得乘以个频率,然后加一块,水平不够,因为要画图,上面X.Y已经用了矩阵了,再算的时候就用不了矩阵了,能帮我解决下这个问题吗,
式中X(x,y,0,H)为下风向x米、横向y米、地面处的扩散的气体浓度,单位为kg/m;
Q为源强(即源释放速率),单位为kg/s;
u为平均风速,单位为m/s;
σy为水平扩散系数,计算公式为0.16*x.*(1+0.0001.*x).^(-1/2)
σz为垂直扩散系数,计算公式为0.12x
H为泄漏源有效高度,单位为m;y为横向距离,单位为m;z为垂直方向距离,单位为m.
因为风速的关系,我们按频率来对式子进行修正,就变成了
其中Hi和fi分别是风速为ui是对应的高度和频率
数据:
Q=0.015079644737231kg/s
风速的倒数1/u=[0.763762 1.013766 0.805282 1.981871 0.435155 0.502842 0.861758 0.431514]
频率f=[0.067692 0.033846 0.070769 0.043076 0.080000 0.233846
0.083076 0.387692]
matlab的问题.之前使用这个公式,我采用[x,y] =meshgrid(0:1:2000,-200:1:200);by0 =0.32*x.*(1+0.0004.*x).^(-1/2);bz0 =0.24*x.*(1+0.0001.*x).^(-1/2);tempy1 =(-1/2).*((y.*y./by0./by0)+(H^2.by0.bz0));tempy2 =2.718282.^(tempy1);c =Q/pi/u.*((by0.
请说清你的问题,你是数学有问题还是语文有问题?你这个问题的表述谁看谁不明白.先把你的已知参数列出,并最好具体赋值,再把已知的变量和要计算的变量列出,最后你的求和公式是从多少求到多少个啊,是i在循环吗,有递推吗?要画图最多可有4个变量,也就是说X,x,y和o或H,o和H最好为定值,不是变量,或起码一个不是变量.这样连X,x,y和另一个变量o或H,总共4个,多一个都没法画图了.fi是变量还是参数?怎么其它数学建模的都做得挺顺利,这里的问题就层出不穷呢?挺简单一问题搞这么复杂呀!