不确定度分析

 则

 所以

#include<iostream>
#include <fstream>
using namespace std;
int main()
{
	double Qmax = 200;
	double delT1 = 0.5;
	double delT2 = 0.1;
	double Q;
	double temp1,temp2,min,max;
	double T1, T2;
	double Un;
	double Qreal;
	double H = 0.15;
	double lam = 0.18;
	double pi = 3.1415926;
	std::ofstream out("uncertainty_Un0.1.txt");
	if (!out)
	{
		// 打开文件失败
		std::cerr << "Error opening file!" << std::endl;
		return 1;
	}
	out << "uncertainty"<<"    "<<"Q"<<"    "<<"Qreal" <<"    "<<"ratioQloss" << "    " << "h" << "    " << "T1" << "    " << "T2" << "    " << "delT" << std::endl;

	for (Q = 10; Q <= Qmax; Q = Q + 10)
	{
		double delQ = Q * 0.05;
		for (double h=30;h<=300;h=h+10)
		{
			for (T1=313;T1<=333;T1++) 
			{
				for (T2=283;T2<=313;T2++) 
				{
					temp2 = 2 * pi * H * (T1 - T2) * lam;
					double Qloss = temp2 / log(0.006);
					Qreal = Q + Qloss;
					double ratioQloss = 0.05 - Qloss / Q;
					double delT = T1 - T2;
					temp1 = (delQ * delQ) + ((Qreal * delT1) / (delT)) * ((Qreal * delT1) / (delT)) + ((Qreal * delT2) / (delT) * (Qreal * delT2) / (delT));
					Un = sqrt(temp1)/ Qreal;
					if (Un <= 0.1) {
						out << Un << "    " << Q<<"    " << Qreal<<"    "<<ratioQloss << "    " << h << "    " << T1 - 273 << "    " << T2 - 273 << "    " << delT << std::endl;
					};
					delT = 0;
					Qloss = 0;
				};
			};
		};
				temp1 = 0;
				Un = 0;
				delQ = 0;
				Qreal = 0;

	};
	out.close();
	return 0;
}