空气污染指数计算matlab

以下是一个简单的MATLAB代码来计算空气污染指数（AQI）：

% 输入污染物浓度和标准限值
pm25 = input('PM2.5浓度（μg/m³）：');
pm10 = input('PM10浓度（μg/m³）：');
co = input('CO浓度（mg/m³）：');
no2 = input('NO2浓度（μg/m³）：');
o3_8h = input('O3 8小时滑动平均值（μg/m³）：');
so2 = input('SO2浓度（μg/m³）：');

pm25_limit = 35;
pm10_limit = 50;
co_limit = 4;
no2_limit = 40;
o3_8h_limit = 160;
so2_limit = 50;

% 计算AQI
pm25_aqi = calc_aqi(pm25, pm25_limit, 0, 50, 0, 100);
pm10_aqi = calc_aqi(pm10, pm10_limit, 0, 50, 0, 100);
co_aqi = calc_aqi(co, co_limit, 0, 50, 0, 100);
no2_aqi = calc_aqi(no2, no2_limit, 0, 50, 0, 100);
o3_8h_aqi = calc_aqi(o3_8h, o3_8h_limit, 100, 150, 150, 200);
so2_aqi = calc_aqi(so2, so2_limit, 0, 50, 0, 100);

aqi_array = [pm25_aqi, pm10_aqi, co_aqi, no2_aqi, o3_8h_aqi, so2_aqi];
aqi = max(aqi_array);

% 输出AQI和空气质量等级
if aqi <= 50
    level = '优';
elseif aqi <= 100
    level = '良';
elseif aqi <= 150
    level = '轻度污染';
elseif aqi <= 200
    level = '中度污染';
elseif aqi <= 300
    level = '重度污染';
else
    level = '严重污染';
end

fprintf('AQI = %d, 空气质量等级为%s\n', aqi, level);

function aqi = calc_aqi(c, c_low, i_low, i_high, c_high, c_very_high)
% 计算单个污染物的AQI
if c <= c_low
    aqi = i_low;
elseif c <= c_high
    aqi = (i_high - i_low) / (c_high - c_low) * (c - c_low) + i_low;
elseif c <= c_very_high
    aqi = (100 - i_high) / (c_very_high - c_high) * (c - c_high) + i_high;
else
    aqi = 100;
end
end

使用时，输入各项污染物的浓度和标准限值即可计算出AQI和空气质量等级。其中，`calc_aqi`函数用于计算单个污染物的AQI值。