#include #include using namespace std; class index { private: double te;//乾球 public: int s , rh;//s:選択,rh:相対湿度 double *pte , t , w , h , WB , SR , WS ,THI;//*pte:乾球ポインタ,t:乾球飽和水蒸気圧,w:湿球,h:湿球飽和水蒸気圧WB:WBGT,SR:夏=0.2,冬0.1,WS:平均風速,THI:不快指数 int info();//タイトル表示 double temp(double);//温度を飽和水蒸気圧に変換 double wet(double , int );//乾球の飽和水蒸気圧と相対湿度から湿球を求める int judge(double , double);//飽和水蒸気密度からインフルエンザ環境判定 int wbgt(double, int );//乾球と相対湿度などからwbgtや不快指数を判断する int select(); //相対湿度か湿球のどちらが分かるか選択 }; int index::info() {//タイトル表示 cout << "\033[36m ###################################### \033[m" << endl; cout << "\033[32m Environmental temperature and humidity \033[m" << endl; cout << "\033[36m ###################################### \033[m" << endl; cout << endl; cout << "Select1:Wet bulb temparature measurement" << endl; cout << "Select2:Humidity measurement" << endl; return 0; } double index::temp(double T) {//温度Tを飽和水蒸気圧Pに変換 double P =1-(273.15+T)/647.3 ; P = 221200 * exp((-7.76451 * P + 1.45838 * pow( P , 1.5) - 2.7758 * pow(P , 3) -1.23303 * pow(P , 6)) / (1 - P)); cout << "Satulation vapor pressure Pws = " << P << "hPa" << endl; return P; } double index::wet(double P , int RH) {//乾球の飽和水蒸気圧Pと相対湿度RHから湿球Tを求める double T =P * RH / 100; T = log10(T / 6.1078); T = (237.3 * T) / (7.5 - T); cout << "Wet bulb temperature (℃) =" << T <<"℃" <> s; if (s == 1) { SR = 0; WS = 0; } else { cout << "When season summer ?(1:Yes , Other:No)"; cin >> s; cout << "How wind speed ?(m/s)"; cin >> WS; if (s == 1){ SR = 0.2; } else { SR = 0.1; } } WB = 0.735 * T + 0.0374 * RH + 0.00292 * T * RH + 7.619 * SR - 4.557 * SR * SR - 0.572 * WS - 4.064; cout << "WBGT=" << WB << "℃"; if (WB>=31) { cout << "\033[5;31m Danger !! \033[m" <=28) { cout << "\033[31m Strict caution !! \033[m" <=25) { cout << "\033[5;33m Warning! \033[m" <=21) { cout << "\033[33m Caution \033[m" <= 85) { cout << "THI:" << "\033[5;31m Extreame heat !!! \033[m"; } else if (THI >= 80) { cout << "THI:" << "\033[31m Very hot !! \033[m"; } else if (THI >= 75) { cout << "THI:" << "\033[33m Hot ! \033[m"; } else if (THI >= 70) { cout << "THI:" << "\033[32m Not hot \033[m"; } else if (THI >= 65) { cout << "THI:" << "\033[32m Confortable \033[m"; } else if (THI >= 60) { cout << "THI:" << "\033[32m Not feel \033[m"; } else if (THI >= 55) { cout << "THI:" << "\033[36m Chilly ! \033[m"; } else { cout << "THI:" << "\033[34m Cold !! \033[m"; } return 0; } int index::select(){//相対湿度か湿球のどちらが分かるか選択 pte = &te; cout << "Choose select 1 or 2 ?" << endl; cin >> s; switch (s) { case 1: cout << "(Dry bulb)Temparature ? (-20~50 ℃) " << endl; cin >> *pte; t = temp(*pte); cout << "Humidity ? (%)" << endl; cin >> rh; w = wet(t , rh); h = temp(w); judge(h , *pte); wbgt(*pte , rh); break; case 2: cout << "(Dry bulb)Temparature ? (-20~50 ℃) " << endl; cin >> *pte; t = temp(*pte); cout << "Wet bulb temperature? (℃)" << endl; cin >> w; h = temp(w); rh =int(round(( h / t ) * 100)); cout << "Humidity =" << rh << "(%)" << endl; judge(h , *pte); wbgt(*pte , rh); break; default: cout << "Sorry try again!" << endl; select(); } return 0; } int main() { class index myobj; myobj.info(); myobj.select(); }

environmental temperature humidity