matlab - 使用 for 循环计算累积和
问题描述
我试图在 27 层的模型字段中得出变量的总和。大多数变量适用于每一层,但对于其中一个变量,我正在测量高度的变化,因此从给定层的总高度中减去前一层的“顶部高度”。
基本上,我只是不确定如何使用 for 循环在任何时候表示累积和。
根据我的代码,我目前正在尝试使用两个 for 循环来计算累积和并获取错误。
"Attempted to access flheight(299,162,0,12); index must be a positive integer or logical."
我知道这是因为“flheight(299,162,0,12)”不存在,因为第三维没有“layer = 0”。
no2molcm2 = 0;
dh = 0;
patm = 0;
no2ppm = 0;
for n=0:26
for i=1:27
T = Temp(299,162,i,12); % K
dh = (flheight(299,162,i,12)*100) -flheight(299,162,n,12)*100;
patm = sum(Pres(299,162,i,12))*(1/101325); %atm
R = 82.06; % cm3*atm/(k*mol)
av = 6.022140857747*(10^23); % 1/mol
no2ppm = sum(no2(299,162,i,12));
no2molcm2 = cumsum(((no2ppm*av*patm)/(R*T))*dh);
end
end
我的问题是我到底怎么能告诉matlab当它看到这个输入(或错误输出)时将它设置为零?
解决方案
因此,根据评论,我意识到我应该采取简单的方法并手动指定第一层,然后将其余层相加。这是修复它的代码:
latco = 1;
lonco = 200;
layer1 = (no2(latco,lonco,1,12)*av*(Pres(latco,lonco,1,12)*
(1/101325)))/(R*Temp(latco,lonco,1,12)*h1);
for n=1:26
for i=2:27
T = Temp(latco,lonco,i,12); % K
dh1 = flheight(latco,lonco,i,12)*100;
dh2 = flheight(latco,lonco,n,12)*100;
dh = dh1 - dh2;
patm = sum(Pres(latco,lonco,i,12))*(1/101325); %atm
no2ppm = no2(latco,lonco,i,12);
no2molcm2_26 = sum(((no2ppm*av*patm)/(R*T))*dh);
end
end
no2final = layer1 + no2molcm2_26
推荐阅读
- javascript - Google 广告 dfp js 代码未呈现所有广告位
- java - 如何通过从第一类调用第二类的方法来调用第三类的方法
- c# - 获取 ListView 内开关的选中状态
- r - 如何更换
特定列的值? - postgresql - 当名称实际上是大写时,如何使用确保表名是小写的查询?
- javascript - TypeError:无法读取未定义 bcryptjs 会话的属性“_id”
- c# - 如何使用反射以优雅的方式访问嵌套对象?
- angular - 如何在不覆盖的情况下加载数据?
- python - 如何从气流连接中保存为字符串的字节中获取字节值?
- android - 如何创建一些图像较小而一些图像较大的 Instagram 搜索布局