本文主要是介绍水环境学——利用最小二乘法对BOD降解好氧系数K1估算,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
由BOD的实验室化验资料估算讲解系数K1时,可以采用最小二乘法。
//"input.txt"BOD化验资料
t(d) y(mg/L)
1 55
2 81
3 102
4 119
5 131
6 140
7 147
8 153
9 159
10 162#include "stdafx.h"
#include <iostream>
#include <fstream>
#include <cmath>
void main()
{using namespace std;const int N = 10,//天数M = 100;//最大循环次数double y[N],t[N],f1[N],f2[N],Sum11 = 0, Sum22 = 0, Sum12 = 0, Sum1y = 0, Sum2y = 0,a, d, L0, h, K1;const double h0 = 0.0001;ifstream infile;infile.open("input.txt");for(int i = 0; i<N; i++)infile>>t[i]>>y[i];infile.close();K1 = 0.5;//假设K1的初始值ofstream outfile;outfile.open("output.txt");for(int k = 0; k < M; k++){for(int i = 0; i<N; i++){f1[i] = 1 - exp(-K1*t[i]);f2[i] = t[i]*exp(-K1*t[i]);Sum11 += pow(f1[i],2);Sum22+= pow(f2[i],2);Sum12
这篇关于水环境学——利用最小二乘法对BOD降解好氧系数K1估算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!