python实现阿列金分类法
输入河流主要离子含量(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实现阿列金分类法相关推荐
- python:克里金插值
python:克里金插值 最近写代码遇到了使用样本数据做克里金插值的事情.于是将Excel保存的[x坐标,y坐标,样本值]数据结合tif数据做了克里金插值,并将代码记录下来. 克里金插值结果: 输入数 ...
- python pandas dataframe 列 转换为离散值
python pandas dataframe 列 转换为离散值 import pandas as pd import numpy as np df1 = pd.DataFrame({'日期': [' ...
- [转载] python处理数据列_Python中基于跨列的数据处理
参考链接: Python 克里斯蒂安Cristian算法 python处理数据列 Let's suppose your manager gives you a random dataset and t ...
- 使用Python通过拉马努金公式快速求π
使用Python通过拉马努金公式快速求π 一.前言 π是一个数学常数,定义为:圆的周长与直径的比值. π是一个无理数,也是一个超越数,它的小数部分无限不循环. π可以用来精确计算圆周长.圆面积.球体积 ...
- python excel增加一列_(用Python修改excel中一列数据)python新增一列
python怎么从excel中读取数据? ⒈ #import παέτο import xlrd #Ρύθμι διαδρής path='C:\\Users\\jyjh\\Desktop\\data ...
- python 写入第二列_python读写Excel表格的实例代码(简单实用)
这篇文章主要介绍了python读写Excel表格的方法,本文通过实例代码给大家介绍的非常详细,具有一定的参考借鉴价值,需要的朋友可以参考下.需要先安装两个库:pip install xlrd.pip ...
- python查看各列数据类型_pandas中查看数据类型的几种方式
在获得数据之后.分析数据之前,我们一般需要对数据总体进行一个概览,如有哪些字段,每个字段的类型,值是否缺失等,以下列出了几种方法,供我们方便快捷的查看dataframe的数据类型. 1.维度查看:df ...
- python dataframe 新列_Python第二十四课:Pandas库(四)
Python第二十四课:Pandas库(四)点击上方"蓝字",关注我们. 不知不觉,我们已经跨越了千难万险,从零开始,一步步揭开了Python神秘的面纱.学到至今,回过头,才晓得自 ...
- python增加一列数据_Python编程给numpy矩阵添加一列方法示例
首先我们有一个数据是一个mn的numpy矩阵现在我们希望能够进行给他加上一列变成一个m(n+1)的矩阵 import numpy as np a = np.array([[1,2,3],[4,5,6] ...
最新文章
- Html Agility Pack基础类介绍及运用
- (更新中)篮球相关英语积累与2020年NBA东部决赛:迈阿密热火与波士顿凯尔特人
- java综合项目推荐,推荐20个5月最热门的Java开源项目
- mysql --force -f_Mysql_mysql force Index 强制索引
- python分析工具有哪些_常用Python数据分析工具汇总
- Elasticsearch查询性能优化
- OpenCV形态学操作
- Unity热更新之资源的上传与下载
- linux 查看磁盘分区,文件系统,使用情况的命令和相关工具介绍,新型linux 查看磁盘分区、文件系统、使用情况的命令和相关工具介绍.doc...
- Myeclipse破解后报错解决
- 常见的工业通讯协议有哪些
- 大数据物流项目:概述及Docker入门(一)
- 《分布式微服务电商源码》-项目简介
- html是l面包屑效果,CSS制作面包屑
- 高手速成android开源项目[View篇]
- Win10下怎么查看WIFI密码
- Python提取Word文件中的目录标题保存为Excel文件
- 【热搜】想卷深度学习必会的10题【最全AI面经】
- Lucas定理和拓展Lucas定理
- Linux下载hfs文件,linux读写mac HFS+
热门文章
- 广告加载失败,请稍后重试 ks: 310004,gdt:5006
- 机房环境监控系统的作用,机房监控的主要目的是什么
- Unity XCode自动发布ipa
- 故宫景点功课17:内廷西路(中)
- [知识梳理]虚拟语气
- 路由器vueRouter
- 【2022HVV系列】蓝队手册更新版(建议收藏)
- php开发勋章插件,勋章中心 2.0 For Discuz!7.0
- 戴尔DELL SCV/SC系列存储故障 Storage Center停机错误的解决方案
- 算法导论 练习5.4-4