输入河流主要离子含量(mg/L)

一共是HCO3 CO3 SO4 Cl Ca Mg Na7个离子的含量数据,用空格隔开

思路:先把单位mg/L换成mmol/L

然后根据阿列金分类法的步骤去搞,花了一个小时

测试输入的是黄河中游的离子含量:

212 3.2 78.1 49.5 49.3 23.4 55.7
mg_L = input().split(' ')
mmol = []
mmol.append("{:.2f}".format(eval(mg_L[0]) / 61))
mmol.append("{:.2f}".format(eval(mg_L[1]) / 30))
mmol.append("{:.2f}".format(eval(mg_L[2]) / 48))
mmol.append("{:.2f}".format(eval(mg_L[3]) / 35.5))
mmol.append("{:.2f}".format(eval(mg_L[4]) / 20))
mmol.append("{:.2f}".format(eval(mg_L[5]) / 12))
mmol.append("{:.2f}".format(eval(mg_L[6]) / 25))
print(mmol)
# 212 3.2 78.1 49.5 49.3 23.4 55.7
HCO3 = eval(mmol[0])
CO3_12 = eval(mmol[1])
SO4_12 = eval(mmol[2])
Cl = eval(mmol[3])
Ca_12 = eval(mmol[4])
Mg_12 = eval(mmol[5])
Na = eval(mmol[6])aa = HCO3 + CO3_12list1 = []
list1.append(aa)
list1.append(SO4_12)
list1.append(Cl)dd = max(list1)
if dd == aa:print("碳酸盐类")
elif dd == SO4_12:print("硫酸盐类")
elif dd == Cl:print("氯化物类")list12 = []
list12.append(Ca_12)
list12.append(Mg_12)
list12.append(Na)dd = max(list12)
if dd == Ca_12:print("钙组")
elif dd == Mg_12:print("镁组")
elif dd == Na:print("钠组")bb = HCO3 + CO3_12 + SO4_12
cc = Ca_12 + Mg_12  # 1/2Ca + 1/2Mgif aa>cc:print("I型")
elif aa<cc<bb:print("II型")
elif cc>bb:print("III型")
elif aa == 0:print("IV型")
else:print("有错")

python实现阿列金分类法相关推荐

  1. python:克里金插值

    python:克里金插值 最近写代码遇到了使用样本数据做克里金插值的事情.于是将Excel保存的[x坐标,y坐标,样本值]数据结合tif数据做了克里金插值,并将代码记录下来. 克里金插值结果: 输入数 ...

  2. python pandas dataframe 列 转换为离散值

    python pandas dataframe 列 转换为离散值 import pandas as pd import numpy as np df1 = pd.DataFrame({'日期': [' ...

  3. [转载] python处理数据列_Python中基于跨列的数据处理

    参考链接: Python 克里斯蒂安Cristian算法 python处理数据列 Let's suppose your manager gives you a random dataset and t ...

  4. 使用Python通过拉马努金公式快速求π

    使用Python通过拉马努金公式快速求π 一.前言 π是一个数学常数,定义为:圆的周长与直径的比值. π是一个无理数,也是一个超越数,它的小数部分无限不循环. π可以用来精确计算圆周长.圆面积.球体积 ...

  5. python excel增加一列_(用Python修改excel中一列数据)python新增一列

    python怎么从excel中读取数据? ⒈ #import παέτο import xlrd #Ρύθμι διαδρής path='C:\\Users\\jyjh\\Desktop\\data ...

  6. python 写入第二列_python读写Excel表格的实例代码(简单实用)

    这篇文章主要介绍了python读写Excel表格的方法,本文通过实例代码给大家介绍的非常详细,具有一定的参考借鉴价值,需要的朋友可以参考下.需要先安装两个库:pip install xlrd.pip ...

  7. python查看各列数据类型_pandas中查看数据类型的几种方式

    在获得数据之后.分析数据之前,我们一般需要对数据总体进行一个概览,如有哪些字段,每个字段的类型,值是否缺失等,以下列出了几种方法,供我们方便快捷的查看dataframe的数据类型. 1.维度查看:df ...

  8. python dataframe 新列_Python第二十四课:Pandas库(四)

    Python第二十四课:Pandas库(四)点击上方"蓝字",关注我们. 不知不觉,我们已经跨越了千难万险,从零开始,一步步揭开了Python神秘的面纱.学到至今,回过头,才晓得自 ...

  9. python增加一列数据_Python编程给numpy矩阵添加一列方法示例

    首先我们有一个数据是一个mn的numpy矩阵现在我们希望能够进行给他加上一列变成一个m(n+1)的矩阵 import numpy as np a = np.array([[1,2,3],[4,5,6] ...

最新文章

  1. Html Agility Pack基础类介绍及运用
  2. (更新中)篮球相关英语积累与2020年NBA东部决赛:迈阿密热火与波士顿凯尔特人
  3. java综合项目推荐,推荐20个5月最热门的Java开源项目
  4. mysql --force -f_Mysql_mysql force Index 强制索引
  5. python分析工具有哪些_常用Python数据分析工具汇总
  6. Elasticsearch查询性能优化
  7. OpenCV形态学操作
  8. Unity热更新之资源的上传与下载
  9. linux 查看磁盘分区,文件系统,使用情况的命令和相关工具介绍,新型linux 查看磁盘分区、文件系统、使用情况的命令和相关工具介绍.doc...
  10. Myeclipse破解后报错解决
  11. 常见的工业通讯协议有哪些
  12. 大数据物流项目:概述及Docker入门(一)
  13. 《分布式微服务电商源码》-项目简介
  14. html是l面包屑效果,CSS制作面包屑
  15. 高手速成android开源项目[View篇]
  16. Win10下怎么查看WIFI密码
  17. Python提取Word文件中的目录标题保存为Excel文件
  18. 【热搜】想卷深度学习必会的10题【最全AI面经】
  19. Lucas定理和拓展Lucas定理
  20. Linux下载hfs文件,linux读写mac HFS+

热门文章

  1. 广告加载失败,请稍后重试 ks: 310004,gdt:5006
  2. 机房环境监控系统的作用,机房监控的主要目的是什么
  3. Unity XCode自动发布ipa
  4. 故宫景点功课17:内廷西路(中)
  5. [知识梳理]虚拟语气
  6. 路由器vueRouter
  7. 【2022HVV系列】蓝队手册更新版(建议收藏)
  8. php开发勋章插件,勋章中心 2.0 For Discuz!7.0
  9. 戴尔DELL SCV/SC系列存储故障 Storage Center停机错误的解决方案
  10. 算法导论 练习5.4-4