java马斯京根法计算汇流系数P

package game;import java.util.*;
public class hydrological_forecast2 {public static void main(String[] args) {Scanner sc=new Scanner(System.in);// 参数Ke、Xe和n的确定
//      System.out.print("请输入河长L值：");
//      int L = sc.nextInt();
//      System.out.print("请输入△t值：");
//      int bt = sc.nextInt();
//      System.out.print("请输入x值：");
//      double x = sc.nextDouble();
//      System.out.print("请输入K值：");
//      int K = sc.nextInt();int L = 512;int bt = 3;double x = 0.34;int K = 12;int Ke = bt;int n = K/bt;double Le = L/n;double Xe = 0.5-n*(1-2*x)/2;System.out.println("Ke值为"+Ke);System.out.println("Xe值为"+Xe);System.out.println("n值为"+n);
//       C0 C1 C2 的确定double C0 = (0.5*bt-Ke*Xe)/(0.5*bt+Ke-Ke*Xe);double C1 = (0.5*bt+Ke*Xe)/(0.5*bt+Ke-Ke*Xe);double C2 = (-0.5*bt+Ke-Ke*Xe)/(0.5*bt+Ke-Ke*Xe);double CC = C0+C1+C2;System.out.println("C0值=："+C0);System.out.println("C1值=："+C1);System.out.println("C2值=："+C2);System.out.println("C0+C1+C2值=："+CC);// 创建数组列表sites储存PArrayList<Double> sites = new ArrayList<Double>();double A = C1+C0*C2;// 将0段汇流系数P加入数组列表sites.add(Math.pow(C0, 4));// 其余段汇流系数pmdouble sum;for (int m = 1; m <20; m++) {sum = 0;for (int i = 1; i <= n ; i++) {if (m-i >= 0) {double A1 = Math.pow(A, i);double B = Factorial(n)*Factorial(m-1)/(Factorial(i)*Factorial(i-1)*Factorial(n-i)*Factorial(m-i));double C01 = Math.pow(C0, n-i);double C21 = Math.pow(C2, m-i);sum += A1*B*C01*C21;}}sites.add(sum);}// 输出System.out.println("汇流系数保留六位小数：");double falg=0;int p = 0;for (int i = 0; i < sites.size(); i++) {System.out.println(String.format("%.6f", sites.get(i)));falg += sites.get(i);// 定义条件，跳出循环if (falg >= 1) {p = i;break;}}outQ(sites,p);}// 阶乘计算public static double Factorial(int n) {double sum=1;for (int i = 1; i <= n ; i++) {sum*=i;}return sum;}// 输入输出Q  注意：此方法未完善public static void  outQ(ArrayList<Double> sites,int p) {Scanner sc = new Scanner(System.in);System.out.println("请输入上游入流量个数：");int x = sc.nextInt();System.out.println("请输入上游入流量：");double[] I = new double[x];for (int i = 0; i < x; i++) {I[i] = sc.nextDouble();}// 计算入流乘于汇流系数double[][] qp = new double[I.length][p];for (int j = 0; j < I.length; j++) {for (int k = 0; k < p; k++) {qp[j][k] = I[j]*sites.get(k);}}// 输出下游流量double sum;for (int i = 0; i < p; i++) {sum = 0;for (int j = qp.length-1; j >= 0 ; j--) {sum += qp[i][j];sum += qp[i+1][j-1];System.out.println(sum);}}}
}

java马斯京根法计算汇流系数P相关推荐

马斯京根法matlab程序,（一）马斯京根流量演算法
(1)x值试算表: 经过对取不同的x值作图可知,当x＝0.1时,直线较为接近作直线求的K＝△s/△Q'=7.5,其中x是流量比重因子,k是河段发生稳定流的传播时间;根据求的x和k求的流量系数Co,C1 ...
java之利用行列式展开法，计算行列式的值
一.提出问题如何利用行列式展开法计算行列式的值二.基本理论与操作说明三.举例说明 1.代码 public class first {static int[][] geta(int m[][],i ...
电动力学每日一题 2021/10/15 Fourier变换法计算均匀电流密度产生的磁场
电动力学每日一题 2021/10/15 Fourier变换法计算均匀电流密度产生的磁场无限长均匀电流无限长圆柱面均匀电流密度无限长均匀电流假设z轴上有一根非常细的电线,携带均匀电流I0I_0I ...
nssl1163-小x游世界树【树形dp,二次扫描和换根法】
正题题目大意一棵树,一条边的权是原本的权值减去出发点的加速. 求一个点使得这个点到所有点路径边权和最小. 解题思路我们先求出以1为根时的答案然后用换根法我们从1转移到2,我们会发现红色的部 ...
Java实现标题相似度计算，文本内容相似度匹配，Java通过SimHash计算标题文本内容相似度
目录一.前言二.关于SimHash 补充知识一).什么是海明距离二).海明距离的应用三).什么是编辑距离三.SimHash算法的几何意义和原理一).SimHash算法的几何意义二). ...
eaxsinbx_用“三行列表法”计算两次分部积分题
用"三行列表法"计算两次分部积分题徐群飞Ξ(绍兴文理学院上虞分院数学系 ,浙江上虞 312300)摘要 :两次分部积分题的计算比较复杂 ,若用"三行列表法&quo ...
matlab用jc法计算可靠度,用MATLAB实现JC法计算结构可靠度程序
软件设计开发本栏目责任编辑 :谢媛媛 ComputerKnowledgeand Technology 电脑知识与技术第 5 卷第 29 期 (2009 年 10 月) 用 MATLAB 实现 JC 法 ...
【软件测试】使用边界值分析法和等价类划分法计算佣金
[软件测试]使用边界值分析法和等价类划分法计算佣金前言 1.边界值分析法 1.1 边界分析 1.1.1 设计测试用例 1.2 程序源码 2.等价类划分法 2.1 划分等价类 2.2 为有效等价类设计 ...
Java五子棋Ai-权值法
Java五子棋Ai-权值法一. 什么是权值法? 二.实现权值法Ai会面临的问题 1.如何表达棋盘上未下处的横竖斜的棋子情况?(类似01120) 2.如何将未下处的横竖斜情况和权值联系起来?(Hash ...

java马斯京根法计算汇流系数P

java马斯京根法计算汇流系数P相关推荐

最新文章

热门文章