前端技术
HTML
CSS
Javascript
前端框架和UI库
VUE
ReactJS
AngularJS
JQuery
NodeJS
JSON
Element-UI
Bootstrap
Material UI
服务端和客户端
Java
Python
PHP
Golang
Scala
Kotlin
Groovy
Ruby
Lua
.net
c#
c++
后端WEB和工程框架
SpringBoot
SpringCloud
Struts2
MyBatis
Hibernate
Tornado
Beego
Go-Spring
Go Gin
Go Iris
Dubbo
HessianRPC
Maven
Gradle
数据库
MySQL
Oracle
Mongo
中间件与web容器
Redis
MemCache
Etcd
Cassandra
Kafka
RabbitMQ
RocketMQ
ActiveMQ
Nacos
Consul
Tomcat
Nginx
Netty
大数据技术
Hive
Impala
ClickHouse
DorisDB
Greenplum
PostgreSQL
HBase
Kylin
Hadoop
Apache Pig
ZooKeeper
SeaTunnel
Sqoop
Datax
Flink
Spark
Mahout
数据搜索与日志
ElasticSearch
Apache Lucene
Apache Solr
Kibana
Logstash
数据可视化与OLAP
Apache Atlas
Superset
Saiku
Tesseract
系统与容器
Linux
Shell
Docker
Kubernetes
[IEEE 754标准与JSON数值表示问...]的搜索结果
这里是文章列表。热门标签的颜色随机变换,标签颜色没有特殊含义。
点击某个标签可搜索标签相关的文章。
点击某个标签可搜索标签相关的文章。
Shell
...本,快速定位服务异常问题,极大地提升了运维响应速度和系统稳定性。同时,数据分析师也在利用awk处理CSV、JSON等多种格式的数据源,结合Python或R等高级编程语言进行深度分析和可视化呈现,为业务决策提供强有力的支持。 此外, awk不仅仅局限于处理结构化文本,它还可以结合正则表达式实现复杂模式匹配,这在网络安全领域同样大有可为,比如用于恶意流量的日志识别和追踪。 总的来说,awk作为一款经典且功能强大的文本处理工具,其价值在当今时代并未因新型技术的崛起而减弱,反而在与各类现代技术和场景的融合中焕发新生,持续为数据处理与分析工作带来便利与高效。因此,掌握awk并深入了解其在不同领域的实践案例,对于提升个人技能和工作效率具有显著的意义。
2023-05-17 10:03:22
67
追梦人-t
Mahout
...,继承自 Java 标准库中的 IllegalArgumentException。当调用 Mahout 库的方法或构造函数时,如果传入的参数不符合预期条件或者违反了方法执行的前提约束(例如矩阵维度不匹配或索引超出范围),该异常就会被抛出,用于提示开发者检查并修正错误的输入参数。 RandomAccessSparseVector , 在 Apache Mahout 中,RandomAccessSparseVector 是一种稀疏向量的实现类,特别适用于大部分元素为零的大维度向量场景。这种数据结构仅存储非零元素及其对应的索引,从而极大地节省了内存空间。相较于密集向量(如 DenseVector),稀疏向量在进行数值计算和存储时更加高效,尤其适合于大规模机器学习和数据挖掘任务中的特征向量表示。
2023-10-16 18:27:51
115
山涧溪流
Beego
...ego中如何解决这些问题。 一、UUID生成 在分布式系统中,我们常常需要生成全局唯一的标识符,也就是我们常说的UUID。UUID是一个128位的数字,可以用来表示一个特定的对象。在Go语言中,我们可以使用标准库中的math/rand包和time包来生成UUID。 go import ( "crypto/rand" "encoding/hex" "math/big" "time" ) func NewUUID() string { var b [16]byte _, err := rand.Read(b[:]) if err != nil { panic(err) } now := time.Now().UnixNano() b[6] = byte((now >> 40) & 0xf) b[7] = byte(now >> 32) b[8] = byte(now >> 24) b[9] = byte(now >> 16) b[10] = byte(now >> 8) b[11] = byte(now) return hex.EncodeToString(b[:]) } 二、自增ID生成 自增ID是一种常见的数据库主键生成方式,它通过不断增加一个整数值来保证数据的唯一性。在Beego这个框架里头,如果你想实现自动增长ID的功能,完全可以这样做:先定义一个模型,然后在这个模型里头添加一个类型为uint的ID字段,这就搞定了自增ID的需求。就像是给每一条记录分配一个独一无二的数字身份证一样,每次新增记录时,这个ID会自动加一,省去了手动指定ID的麻烦。 go type User struct { ID uint orm:"column(id);auto" Name string Email string Phone string Address string } 以上代码中,我们在User模型中定义了一个名为ID的字段,并设置了它的类型为uint和auto。这样,每次插入一条新的用户记录时,ID字段都会自动递增。 三、UUID和自增ID的选择 在实际开发中,我们常常需要根据具体的需求来选择生成哪种类型的ID。如果我们正在捣鼓一个分布式系统,那么选用UUID绝对是个更酷的选择。为啥呢?因为它可以在全球这个大舞台上保证每个ID都是独一无二的,就像每个人都有自己的指纹一样独特。假如我们正在捣鼓一个单机应用,那么选择自增ID可能是个更省心省力的办法。为啥呢?因为它生成的速度贼快,而且出岔子的概率也低得多,这样一来,我们就不用在这方面费太多心思啦! 四、总结 总的来说,生成UUID或自增ID是我们在开发Web应用时经常会遇到的问题。在Beego中,我们可以通过简单的代码就能实现这两种ID的生成。不过呢,具体要用哪种类型的ID,咱们还得根据实际需求来掂量决定。无论我们挑哪一个,只要能把数据的唯一性和安全性稳稳地守住,那就都是个没毛病的选择。
2023-11-17 22:27:26
589
翡翠梦境-t
JSON
JSON查询第二条记录:深入探索与实践 1. 引言 --- 在日常的Web开发和数据交互中,JSON(JavaScript Object Notation)扮演着至关重要的角色。这玩意儿就是个轻巧便捷的数据交换格式,瞅着贼容易让人理解,写起来也倍儿顺手;对机器来说,解析和生成它更是小菜一碟,轻松加愉快。本文将围绕“如何在JSON数据中查询第二条记录”这一主题进行探讨,通过实例代码演示,带您逐步揭开这个看似简单实则富含技巧的问题。 2. JSON基础认知 --- 首先,让我们温习一下JSON的基础知识。JSON数据呢,平常就像个小管家,喜欢把信息一对对地配好放在一起,这一对就叫键值对。这些“小对对”聚在一起,就成了一个“大对象”。而当很多个这样的“大对象”手牵手串成一串的时候,我们就称它为数组啦。例如: json { "employees": [ { "id": 1, "name": "John Doe", "position": "Manager" }, { "id": 2, "name": "Jane Smith", "position": "Developer" }, // 更多员工记录... ] } 在这个例子中,employees 是一个包含多个员工对象的数组,我们想要的目标是获取并查询数组中的第二条员工记录。 3. 查询JSON中的第二条记录 --- 那么,如何从上述JSON数据中提取出第二条记录呢?这就需要借助编程语言提供的JSON解析功能,这里我们以JavaScript为例,因为JSON的设计灵感就来源于JavaScript的对象表示法。 javascript let jsonData = { "employees": [ // 员工记录... ] }; // 获取第二条记录 let secondEmployee = jsonData.employees[1]; console.log(secondEmployee); 在这段代码中,jsonData.employees[1]就是我们获取到的第二条员工记录。注意,数组索引是从0开始的,所以索引1对应的是数组中的第二个元素。 4. 深入理解与思考 --- 细心的你可能已经注意到,这里的“第二条记录”实际上是基于数组索引的概念。要是有一天,JSON结构突然变了样儿,比如员工们不再像以前那样排着整齐的数组队列,而是藏在了其他对象的小屋里,那咱们查询的方法肯定也得跟着变一变啦。 json { "employeeRecords": { "record1": { "id": 1, "name": "John Doe", "position": "Manager" }, "record2": { "id": 2, "name": "Jane Smith", "position": "Developer" }, // 更多记录... } } 对于这种情况,由于不再是有序数组,查找“第二条记录”的概念变得模糊。我们无法直接通过索引定位,除非我们知道特定键名,如"record2"。不过,在现实操作里,咱们经常会根据业务的具体需求和数据的组织架构,设计出更接地气、更符合场景的查询方法。比如,先按照ID从小到大排个序,再捞出第二个记录;或者给每一条记录都标上一个独一无二的顺序标签,让它们在队列里乖乖站好。 5. 结论与探讨 --- 总的来说,查询JSON中的第二条记录主要取决于数据的具体结构。在处理JSON数据时,理解其内在结构和关系至关重要。不同的数据组织方式会带来不同的查询策略。在实际动手操作的时候,我们得把编程语言处理JSON的那些技巧玩得溜溜的,同时还要瞅准实际情况,琢磨出最接地气、最优解决方案。 最后,我鼓励大家在面对类似问题时,不妨像侦探破案一样去剖析JSON数据的构造,揣摩其中的规律和逻辑,这不仅能帮助我们更好地解决问题,更能锻炼我们在复杂数据环境中抽丝剥茧、寻找关键信息的能力。
2023-04-13 20:41:35
459
烟雨江南
SeaTunnel
...!本文将深入探讨这类问题的成因,并通过丰富的代码实例演示如何在SeaTunnel中妥善解决这些问题。 2. Parquet/CSV文件解析常见问题及其原因 2.1 数据类型不匹配 Parquet和CSV两种格式对于数据类型的定义和处理方式有所不同。比如,你可能会遇到这么个情况,在CSV文件里,某个字段可能被不小心认作是文本串了,但是当你瞅到Parquet文件的时候,嘿,这个同样的字段却是个整数类型。这种类型不匹配可能导致解析错误。 python 假设在CSV文件中有如下数据 id,name "1", "John" 而在Parquet文件结构中,id字段是int类型 (id:int, name:string) 2.2 文件格式规范不一致 Parquet和CSV对空值、日期时间格式等有着各自的约定。如CSV中可能用“null”、“N/A”表示空值,而Parquet则以二进制标记。若未正确配置解析规则,就会出现错误。 3. 利用SeaTunnel解决文件格式解析错误 3.1 配置数据源与转换规则 在SeaTunnel中,我们可以精细地配置数据源和转换规则以适应各种场景。下面是一个示例,展示如何在读取CSV数据时指定字段类型: yaml source: type: csv path: 'path/to/csv' schema: - name: id type: integer - name: name type: string transform: - type: convert fields: - name: id type: int 对于Parquet文件,SeaTunnel会自动根据Parquet文件的元数据信息解析字段类型,无需额外配置。 3.2 自定义转换逻辑处理特殊格式 当遇到非标准格式的数据时,我们可以使用自定义转换插件来处理。例如,处理CSV中特殊的空值表示: yaml transform: - type: script lang: python script: | if record['name'] == 'N/A': record['name'] = None 4. 深度思考与讨论 处理Parquet和CSV文件解析错误的过程其实也是理解并尊重每种数据格式特性的过程。SeaTunnel以其灵活且强大的数据处理能力,帮助我们在面对这些挑战时游刃有余。但是同时呢,我们也要时刻保持清醒的头脑,像侦探一样敏锐地洞察可能出现的问题。针对这些问题,咱们得接地气儿,结合实际业务的具体需求,灵活定制出解决问题的方案来。 5. 结语 总之,SeaTunnel在应对Parquet/CSV文件格式解析错误上,凭借其强大的数据源适配能力和丰富的转换插件库,为我们提供了切实可行的解决方案。经过实战演练和持续打磨,我们能够更溜地玩转各种数据格式,确保数据整合和ETL过程一路绿灯,畅通无阻。所以,下次你再遇到类似的问题时,不妨试试看借助SeaTunnel这个好帮手,让数据处理这件事儿变得轻轻松松,更加贴近咱们日常的使用习惯,更有人情味儿。
2023-08-08 09:26:13
76
心灵驿站
转载文章
... 1486:分数这一问题后,我们可以延伸至教育评估领域中关于考试设计与数据分析的最新研究进展。近日,美国教育考试服务中心(ETS)发布了一项关于利用大数据优化试题难度与区分度的研究报告。该研究表明,在大规模标准化测试中,运用机器学习算法和统计模型能够有效分析考生答题数据,精确调整题目难度和区分度,从而提高考试结果的信度和效度。 具体而言,研究人员借鉴了单峰函数优化方法,并创新性地结合三分法策略来动态调整试题参数,以实现得分分布的最佳匹配。这种方法不仅适用于编程竞赛的评分系统优化,更在各类资格认证、入学选拔等高风险考试设计中展现出了巨大潜力。同时,报告强调了保留有效数字的重要性,确保成绩计算和排名的公平性和准确性。 此外,随着我国新高考改革的深入推进,考试评价体系也在不断升级和完善。例如,部分地区引入智能化考试系统,通过实时监测和分析学生作答数据,动态生成适合不同层次学生的考题,实现了对考试难度和区分度的精细化管理,有力推动了教育公平与质量提升。 总之,从DTOJ 1486:分数这一具体的编程问题出发,我们看到了现代科技如何赋能传统考试评价方式,使其在保持公正严谨的同时,更加科学高效。未来,随着人工智能和大数据技术的持续发展,考试设计与数据分析将深度融合,进一步推动教育评价体系的现代化进程。
2023-08-30 11:55:56
154
转载
AngularJS
...置过滤器如date、json等的使用方法。 2. Angular 9/10过滤器新特性解读:尽管AngularJS已进入长期支持阶段,但其后续版本Angular仍保留了对数据处理的强大支持。在Angular 9/10中,管道(Pipe)作为过滤器的进化形态,提供了更丰富的功能和更高的性能。例如,通过自定义管道实现复杂的数据格式化需求,以及利用pure和impure管道优化性能表现。 3. 实战教程:构建响应式表单结合自定义过滤器:一篇近期的技术博客详细介绍了如何在Angular应用中结合自定义过滤器与响应式表单,实现实时数据验证和格式化显示,这为开发者解决实际项目中的具体问题提供了极具时效性的解决方案。 4. 案例分享:电商网站商品筛选功能实现:参考某知名电商平台近期公开的技术文章,其中详述了如何运用AngularJS(或Angular)过滤器进行多条件商品列表筛选,展示了过滤器在大规模数据处理场景下的高效应用。 5. 社区讨论:过滤器在状态管理库NGXS中的创新实践:随着状态管理库NGXS在Angular社区的广泛应用,有开发者提出并分享了如何将过滤逻辑融入到状态管理中,从而简化视图层代码,提高应用的整体架构层次性和可维护性。 持续关注Angular及前端领域的技术博客、论坛和GitHub项目,可以帮助开发者紧跟行业发展步伐,更好地运用过滤器这一强大工具提升应用程序的数据展示效果与用户体验。
2024-03-09 11:18:03
476
柳暗花明又一村
Shell
...态。这个状态是一个整数值,通常0表示命令成功执行,非零值则代表出现了错误或异常情况。在Shell脚本编程中,通过特殊变量$?可以获取上一条命令的退出状态,根据这个状态码,程序员可以编写条件判断语句以进行相应的错误处理。 信号(Signal) , 在操作系统层面,信号是一种进程间通信机制,用于通知进程发生了某种事件。在Shell脚本中,可以通过trap命令捕获特定的信号,如当脚本发生错误时产生的ERR信号。一旦接收到预设的信号,就会触发预先定义好的命令序列,例如进行资源清理、日志记录等操作,以实现对错误的及时响应和处理。 嵌套脚本(Nested Script) , 嵌套脚本是指在一个Shell脚本内部调用另一个Shell脚本或函数的情况。在复杂Shell脚本编写过程中,这种结构非常常见,它有助于模块化代码并提高可维护性。然而,在嵌套结构中,子脚本或函数运行时发生的错误需要正确地向父脚本传播,并在父脚本层面上得到妥善处理,否则可能会导致整个程序逻辑混乱或者资源未被正确释放的问题。为了确保这一点,可以通过检查子脚本或函数执行后的退出状态(即使用$?变量),并在必要时返回非零状态码,从而实现错误信息的有效传递和处理。
2024-03-02 10:38:18
84
半夏微凉
VUE
...PA)中。 针对这一问题,开发者社区提出了多项改进措施。其中,JWT(JSON Web Tokens)的使用得到了越来越多的关注。JWT不仅能够简化身份验证流程,还能提高系统的安全性。此外,一些开源库如express-jwt和jsonwebtoken,为开发者提供了便捷的工具来处理JWT相关的操作,从而减少因不当实现而导致的安全风险。 另外,随着微服务架构的普及,跨域资源共享(CORS)成为另一个需要关注的领域。确保正确配置CORS策略对于防止未授权访问至关重要。例如,最近Netflix公开分享了其在构建大规模微服务架构时如何处理CORS的最佳实践,其中包括详细的配置指南和常见陷阱的避免方法。 最后,持续集成/持续部署(CI/CD)流水线中的自动化安全检查也变得越来越重要。通过将安全扫描工具集成到CI/CD流程中,可以及早发现并修复潜在的安全漏洞。例如,GitHub Actions和GitLab CI等平台提供了丰富的插件和模板,帮助开发者轻松实现这一目标。 总之,通过采用最新的安全技术和最佳实践,我们可以显著提升Vue项目以及其他Web应用的安全性,从而为用户提供更加可靠的服务。
2025-01-23 15:55:50
29
灵动之光
转载文章
...a2,..,an] 表示a1,a2,..,an的最小公倍数,(a1,a2,..,an)表示a1,a2,..,an的最大公约数,其中a1,a2,..,an为非负整数。对于两个数a,b,有[a,b]=ab/(a,b),因此两个数最小公倍数可以用其最大公约数计算。但对于多个数,并没有[a1,a2,..,an]=M/(a1,a2,..,an)成立,M为a1,a2,..,an的乘积。例如:[2,3,4]并不等于24/(2,3,4)。即两个数的最大公约数和最小公倍数之间的关系不能简单扩展为n个数的情况。 本文对多个数最小公倍数和多个数最大公约数之间的关系进行了探讨。将两个数最大公约数和最小公倍数之间的关系扩展到n个数的情况。在此基础上,利用求n个数最大公约数的向量变换算法计算多个数的最小公倍数。 1. 多个数最小公倍数和多个数最大公约数之间的关系 令p为a1,a2,..,an中一个或多个数的素因子,a1,a2,..,an关于p的次数分别为r1,r2,..,rn,在r1,r2,..,rn中最大值为rc1=rc2=..=rcm=rmax,最小值为rd1=rd2=..=rdt=rmin,即r1,r2,..,rn中有m个数所含p的次数为最大值,有t个数所含p的次数为最小值。例如:4,12,16中关于素因子2的次数分别为2,2,4,有1个数所含2的次数为最大值,有2个数所含2的次数为最小值;关于素因子3的次数分别为0,1,0,有1个数所含3的次数为最大值,有2个数所含3的次数为最小值。 对最大公约数有,只包含a1,a2,..,an中含有的素因子,且每个素因子次数为a1,a2,..,an中该素因子的最低次数,最低次数为0表示不包含[1]。 对最小公倍数有,只包含a1,a2,..,an中含有的素因子,且每个素因子次数为a1,a2,..,an中该素因子的最高次数[1]。 定理1:[a1,a2,..,an]=M/(M/a1,M/a2,..,M/an),其中M为a1,a2,..,an的乘积,a1,a2,..,an为正整数。 例如:对于4,6,8,10,有[4,6,8,10]=120,而M=46810=1920,M/(M/a1,M/a2,..,M/an) =1920/(6810,4810,4610,468)=1920/16=120。 证明: M/a1,M/a2,..,M/an中p的次数都大于等于r1+r2+..+rn-rmax,且有p的次数等于r1+r2+..+rn-rmax的。这是因为 (1) M/ai中p的次数为r1+r2+..+rn-ri,因而M/a1,M/a2,..,M/an中p的次数最小为r1+r2+..+rn-rmax。 (2) 对于a1,a2,..,an中p的次数最大的项aj(1项或多项),M/aj中p的次数为r1+r2+..+rn-rmax。 或者对于a1,a2,..,an中p的次数最大的项aj,M/aj中p的次数小于等于M/ak,其中ak为a1,a2,..,an中除aj外其他的n-1个项之一,而M/aj中p的次数为r1+r2+..+rn-rmax。 因此,(M/a1,M/a2,..,M/an)中p的次数为r1+r2+..+rn-rmax,从而M/(M/a1,M/a2,..,M/an)中p的次数为rmax。 上述的p并没有做任何限制。由于a1,a2,..,an中包含的所有素因子在M/(M/a1,M/a2,..,M/an)中都为a1,a2,..,an中的最高次数,故有[a1,a2,..,an]=M/(M/a1,M/a2,..,M/an)成立。 得证。 定理1对于2个数的情况为[a,b]=ab/(ab/a,ab/b)=ab/(b,a)=ab/(a,b),即[a,b]=ab/(a,b)。因此,定理1为2个数最小公倍数公式[a,b]=ab/(a,b)的扩展。利用定理1能够把求多个数的最小公倍数转化为求多个数的最大公约数。 2.多个数最大公约数的算法实现 根据定理1,求多个数最小公倍数可以转化为求多个数的最大公约数。求多个数的最大公约数(a1,a2,..,an)的传统方法是多次求两个数的最大公约数,即 (1) 用辗转相除法[2]计算a1和a2的最大公约数(a1,a2) (2) 用辗转相除法计算(a1,a2)和a3的最大公约数,求得(a1,a2,a3) (3) 用辗转相除法计算(a1,a2,a3)和a4的最大公约数,求得(a1,a2,a3,a4) (4) 依此重复,直到求得(a1,a2,..,an) 上述方法需要n-1次辗转相除运算。 本文将两个数的辗转相除法扩展为n个数的辗转相除法,即用一次n个数的辗转相除法计算n个数的最大公约数,基本方法是采用反复用最小数模其它数的方法进行计算,依据是下面的定理2。 定理2:多个非负整数a1,a2,..,an,若aj>ai,i不等于j,则在a1,a2,..,an中用aj-ai替换aj,其最大公约数不变,即 (a1,a2,..,aj-1,aj,aj+1,..an)=(a1,a2,..,aj-1,aj-ai,aj+1,..an)。 例如:(34,24,56,68)=(34,24,56-34,68)=(34,24,22,68)。 证明: 根据最大公约数的交换律和结合率,有 (a1,a2,..,aj-1,aj,aj+1,..an)= ((ai,aj),(a1,a2,..,ai-1,ai+1,..aj-1,aj+1,..an))(i>j情况),或者 (a1,a2,..,aj-1,aj,aj+1,..an)= ((ai,aj),(a1,a2,..,aj-1,aj+1,..ai-1,ai+1,..an))(i<j情况)。 而对(a1,a2,..,aj-1,aj-ai,aj+1,..an),有 (a1,a2,..,aj-1,aj-ai,aj+1,..an)= ((ai, aj-ai),( a1,a2,..,ai-1,ai+1,.. aj-1,aj+1,..an))(i>j情况),或者 (a1,a2,..,aj-1,aj-ai,aj+1,..an)= ((ai, aj-ai),( a1,a2,..,aj-1,aj+1,.. ai-1,ai+1,..an))(i<j情况)。 因此只需证明(ai,aj)=( ai, aj-ai)即可。 由于(aj-ai)= aj-ai,因此ai,aj的任意公因子必然也是(aj-ai)的因子,即也是ai,( aj-ai)的公因子。由于aj = (aj-ai)+ai,因此ai,( aj-ai)的任意公因子必然也是aj的因子,即也是ai,aj的公因子。所以,ai,aj的最大公约数和ai,(aj-ai) 的最大公约数必须相等,即(ai,aj)=(ai,aj-ai)成立。 得证。 定理2类似于矩阵的初等变换,即 令一个向量的最大公约数为该向量各个分量的最大公约数。对于向量<a1,a2,..,an>进行变换:在一个分量中减去另一个分量,新向量和原向量的最大公约数相等。 求多个数的最大公约数采用反复用最小数模其它数的方法,即对其他数用最小数多次去减,直到剩下比最小数更小的余数。令n个正整数为a1,a2,..,an,求多个数最大共约数的算法描述为: (1) 找到a1,a2,..,an中的最小非零项aj,若有多个最小非零项则任取一个 (2) aj以外的所有其他非0项ak用ak mod aj代替;若没有除aj以外的其他非0项,则转到(4) (3) 转到(3) (4) a1,a2,..,an的最大公约数为aj 例如:对于5个数34, 56, 78, 24, 85,有 (34, 56, 78, 24, 85)=(10,8,6,24,13)=(4,2,6,0,1)=(0,0,0,0,1)=1, 对于6个数12, 24, 30, 32, 36, 42,有 (12, 24, 30, 32, 36, 42)=(12,0,6,8,0,6)=(0,0,0,2,0,6)=(0,0,0,2,0,0)=2。 3. 多个数最小共倍数的算法实现 求多个数最小共倍数的算法为: (1) 计算m=a1a2..an (2) 把a1,a2,..,an中的所有项ai用m/ai代换 (3) 找到a1,a2,..,an中的最小非零项aj,若有多个最小非零项则任取一个 (4) aj以外的所有其他非0项ak用ak mod aj代替;若没有除aj以外的其他非0项,则转到(6) (5) 转到(3) (6) 最小公倍数为m/aj 上述算法在VC环境下用高级语言进行了编程实现,通过多组求5个随机数最小公倍数的实例,与标准方法进行了比较,验证了其正确性。标准计算方法为:求5个随机数最小公倍数通过求4次两个数的最小公倍数获得,而两个数的最小公倍数通过求两个数的最大公约数获得。 5.结论 计算多个数的最小公倍数是常见的基本运算。n个数的最小公倍数可以表示成另外n个数的最大公约数,因而可以通过求多个数的最大公约数计算。求多个数最大公约数可采用向量转换算法一次性求得。 本篇文章为转载内容。原文链接:https://blog.csdn.net/u012349696/article/details/21233457。 该文由互联网用户投稿提供,文中观点代表作者本人意见,并不代表本站的立场。 作为信息平台,本站仅提供文章转载服务,并不拥有其所有权,也不对文章内容的真实性、准确性和合法性承担责任。 如发现本文存在侵权、违法、违规或事实不符的情况,请及时联系我们,我们将第一时间进行核实并删除相应内容。
2023-10-04 16:29:43
39
转载
转载文章
...进行规范化处理,将其数值范围调整到特定区间内,如0-1之间或者均值为0、标准差为1的标准正态分布区间。在Python中,可以使用sklearn库提供的StandardScaler()函数来实现这一操作。特征缩放有助于消除特征间量纲的影响,使得不同规模的特征在机器学习算法中具有可比性,从而优化模型训练效果。 独热编码 , 独热编码是一种将离散类别型特征转换为数值型特征的方法,主要用于解决分类特征在机器学习算法中的处理问题。在本文提到的场景下,Python的sklearn库提供了OneHotEncoder()函数,用于将非数值型、类别型特征转化为多维度的二进制向量表示,每个维度对应原类别特征的一个可能取值,而具体维度上的值则代表该类别的出现与否。这样处理后的特征形式更便于输入到许多基于数值计算的机器学习模型中进行训练和预测。
2024-02-09 12:42:15
704
转载
Logstash
...t types”这个问题。咱会详细地扒一扒这个错误背后的来龙去脉,再配上些实实在在的代码例子,让大家伙儿能更好地理解这问题,手把手带你把它给解决了哈! 1. Sortfilter介绍 在Logstash的众多过滤器中,Sortfilter是一个非常实用的功能组件,它可以按照指定字段对事件进行排序。比如在处理一些时间戳乱七八糟、不连贯的日志时,我们完全可以借助Sortfilter这个小帮手,把它给咱们按照时间顺序排排队、整整队。 ruby filter { sort { order => "asc" field => "@timestamp" } } 上述配置会按照@timestamp字段(通常为日志的时间戳)的升序对事件进行排序。 2. “Cannot sort array of different types”问题解析 然而,在某些情况下,当我们尝试对包含不同类型元素的数组字段进行排序时,就会遇到“Cannot sort array of different types”的错误提示。这是因为Sortfilter在内部执行排序操作时要求所有待排序的元素必须是同一类型。例如,如果某个字段是一个数组,其中包含了数字和字符串,那么就无法直接对其进行排序: json { "my_array": [1, "two", 3, "four"] } 在这种情况下,如果你试图用Sortfilter对"my_array"进行排序,Logstash将会抛出上述错误,因为数字和字符串不具备可比性,无法明确确定其排序规则。 3. 解决方案及思考过程 面对这个问题,我们需要采取一些策略来确保数组内的元素类型一致,然后再进行排序。以下是一种可能的解决方案: 3.1 类型转换 首先,我们可以通过mutate插件的convert或gsub函数,将数组内所有的元素转换为同一种类型,如全部转换为字符串或数值。 ruby filter { mutate { convert => { "[my_array]" => "string" } 将数组元素转为字符串 } sort { order => "asc" field => "[my_array]" } } 请注意,这种方式虽能解决问题,但可能会丢失原始数据的一些特性,比如数值大小关系。若数组内混有数字和字符串,且需要保留数字间的大小关系,则需谨慎使用。 3.2 分别处理并合并 另一种方法是对数组进行拆分,分别对不同类型的数据进行排序,再合并结果。不过呢,这通常意味着需要处理更复杂的逻辑,讲到对Logstash配置文件的编写,那可能会让你觉得有些烧脑,不够一目了然,就像解一个九连环谜题一样。 4. 探讨与总结 在日常使用Logstash的过程中,理解并妥善处理数据类型是非常关键的。特别是在处理像排序这种对数据类型特别依赖的任务时,咱们得确保数据的“整齐划一”和“可比性”,就像排队买票,每个人都得按照身高或者年龄排好队,这样才能顺利进行。虽然乍一看,“Sortfilter: Cannot sort array of different types”这个问题好像挺基础,但实际上它悄悄点出了我们在应对各种类型混杂的数据时,不得不面对的一个大难题——就是在确保数据本身含义不被扭曲的前提下,如何把数据收拾得整整齐齐、妥妥当当,做好有效的数据清洗和预处理工作。 因此,在设计和实施Logstash管道时,不仅要关注功能实现,更要注重对原始数据特性的深入理解和恰当处理。这样子做,咱们才能让Logstash这家伙更贴心地帮我们处理数据分析和可视化的事儿,进而从海量数据中淘出真正的金子来。
2023-03-09 18:30:41
303
秋水共长天一色
Etcd
...cd使用了Go语言的标准日志库logrus,提供了多个级别的日志输出,包括Debug、Info、Warning、Error以及Fatal五个等级。不同的日志级别对应不同的信息详细程度: - Debug:记录详细的调试信息,用于开发阶段的问题排查。 - Info:提供运行时的基本信息,如节点启动、客户端连接等。 - Warning:记录潜在错误或非预期行为,但不影响程序正常运行。 - Error:记录已发生错误,可能影响部分功能。 - Fatal:记录严重错误,导致进程终止。 2. 设置Etcd日志级别 Etcd的日志级别可以通过启动参数--log-level来设定。下面是一段启动Etcd并将其日志级别设置为info的示例代码: bash ./etcd --name my-etcd-node \ --data-dir /var/lib/etcd \ --listen-peer-urls http://localhost:2380 \ --listen-client-urls http://localhost:2379 \ --initial-cluster-token etcd-cluster-1 \ --initial-cluster=my-etcd-node=http://localhost:2380 \ --advertise-client-urls http://localhost:2379 \ --log-level=info 上述命令行中--log-level=info表示我们只关心Info及以上级别的日志信息。 3. 输出方式与格式化 Etcd默认将日志输出到标准错误(stderr),你也可以通过--log-output参数指定输出文件,例如: bash ./etcd --log-output=/var/log/etcd.log ... 此外,Etcd还支持JSON格式的日志输出,只需添加启动参数--log-format=json即可: bash ./etcd --log-format=json ... 4. 实践应用与思考 在日常运维过程中,我们可能会遇到各种场景需要调整Etcd的日志级别。比如,当我们的集群闹脾气、出现状况时,我们可以临时把日志的“放大镜”调到Debug级别,这样就能捞到更多更细枝末节的内部运行情况,像侦探一样迅速找到问题的幕后黑手。而在平时一切正常运转的日子里,为了让日志系统保持高效、易读,我们一般会把它调到Info或者Warning这个档位,就像给系统的日常表现打个合适的标签。 同时,合理地选择日志输出方式也很重要。直接输出至终端有利于实时监控,但不利于长期保存和分析。所以,在实际的生产环境里,我们通常会选择把日志稳稳地存到磁盘上,这样一来,以后想回过头来找找线索、分析问题什么的,就方便多了。 总的来说,熟练掌握Etcd日志级别的调整和输出方式,不仅能让我们更好地理解Etcd的工作状态,更能提升我们对分布式系统管理和运维的实战能力。这就像一位超级厉害的侦探大哥,他像拿着放大镜一样细致地研究Etcd日志,像读解神秘密码那样解读其中的含义。通过这种抽丝剥茧的方式,他成功揭开了集群背后那些不为人知的小秘密,确保我们的系统能够稳稳当当地运行起来。
2023-01-29 13:46:01
832
人生如戏
Apache Atlas
...格,允许软件应用通过标准HTTP方法(如GET、POST、PUT、DELETE等)与服务器进行交互,以获取或修改资源状态。在本文中,Apache Atlas的REST API是开发者与Atlas服务之间进行通信的关键桥梁,通过发送JSON格式的数据来创建、读取、更新和删除Atlas中的实体对象。 元数据 , 元数据是对数据的数据,即关于数据的信息。在大数据环境中,元数据可以描述数据集的结构、来源、所有权、生命周期、安全性要求等多种属性。例如,在Apache Atlas中,一个Hive表的元数据可能包括表名、列名、表的创建时间、所属数据库、表的所有者以及表的权限信息等,这些信息有助于用户理解和管理实际的数据内容。 实体模型 , 在Apache Atlas中,实体模型是用来描述和管理不同类型的业务对象或IT组件的抽象框架。每个实体类型都有特定的一组属性和关系,比如Hive表实体类型就包含了名称、描述、所属数据库等属性。实体模型允许用户根据实际业务需求定义和扩展不同的实体类型,并通过实体之间的关联关系构建出丰富的元数据图谱。 访问控制列表(ACL) , 访问控制列表是一种安全机制,用于指定哪些用户或角色有权访问特定的系统资源或执行特定的操作。在Apache Atlas中,ACL用于管理用户的权限,确保只有具备足够权限的用户才能成功地执行诸如创建实体之类的操作。通过调整和配置ACL,管理员可以精细地控制各个用户或角色在Atlas平台上的操作权限,从而保障系统的安全性和数据的完整性。
2023-06-25 23:23:07
561
彩虹之上
转载文章
...装饰物亮度匹配的实际问题中,快速傅里叶变换(FFT)展现出了其强大的优化能力。通过巧妙地将问题转化为求解序列卷积的最大值,我们可以借助FFT技术将原本可能需要O(n^2)时间复杂度的运算降低至O(nlogn),从而高效找到最优解。实际上,FFT的应用远不止于此,它在信号处理、图像处理、数据压缩等领域都有着广泛而深入的应用。 近日,在科学计算领域,《自然》杂志报道了一项利用FFT算法优化能源传输网络的研究成果。科研团队成功运用FFT分析了电网中各个节点间的电力波动情况,通过对大量实时数据进行快速卷积计算,精准预测并优化了电能分配策略,极大地提高了能源传输效率和稳定性,这再次验证了FFT在实际工程问题中的强大作用。 此外,深度学习领域的研究者也在探索如何结合FFT与卷积神经网络(CNN),以提升模型训练速度和推理效率。一项发表于《IEEE Transactions on Neural Networks and Learning Systems》的论文中,研究人员创新性地提出了一种基于FFT的卷积操作方法,可以显著减少CNN中的计算量,尤其在处理大规模图像识别任务时效果尤为明显。 总的来说,从日常生活中的情侣手环亮度调整问题到关乎国计民生的能源传输优化,再到前沿的人工智能技术突破,快速傅里叶变换始终以其独特的数学魅力和高效的计算性能发挥着关键作用。随着科学技术的发展,我们有理由相信FFT将在更多领域带来革命性的解决方案。
2023-01-20 17:51:37
524
转载
转载文章
...元组的概念,以高效地表示和处理多维数据。在处理大规模数据分析任务时,用户可以通过创建不同类型的元组来表达复杂的键值对或更丰富的数据结构,从而更好地适应多样化的大数据场景。 在未来,随着JDK的发展和社区对数据结构需求的深入挖掘,元组类库可能会进一步丰富和完善,提供更为灵活且高性能的API,使得开发者能够更加自如地在各类项目中运用元组这一强大的工具,解决更多类型安全和数据组合的问题。而随着Java模块化系统(JPMS)的成熟,对于元组库的依赖管理也将更加便捷,有助于推动其在更多实际项目中的落地应用。
2023-09-17 17:43:51
257
转载
转载文章
...) 概率与统计的中心问题,都是random variable, PMF与PDF PMF:probability mass function,概率质量函数,是离散型随机变量在各特定取值上的概率。与概率密度函数(PDF:probability density function)的不同之处在于:概率质量函数是对离散型随机变量定义的,本身代表该值的概率;概率密度函数是针对连续型随机变量定义的,本身不是概率(连续型随机变量单点测度为0),只有在对连续随机变量的pdf在某一给定的区间内进行积分才是概率。 notation 假设X 是一个定义在可数样本空间S 上的离散型随机变量S⊆R ,则其概率质量函数PMF为: fX(x)={Pr(X=x),0,x∈Sx∈R∖S 注意这在所有实数上,包括那些X 不可能等于的实数值上,都定义了pmf,只不过在这些X 不可能取的实数值上,fX(x) 取值为0(x∈R∖S,Pr(X=x)=0 )。 离散型随机变量概率质量函数(pmf)的不连续性决定了其累积分布函数(cdf)也不连续。 共轭先验(conjugate prior) 所谓共轭(conjugate),描述刻画的是两者之间的关系,单独的事物不构成共轭,举个通俗的例子,兄弟这一概念,只能是两者才能构成兄弟。所以,我们讲这两个人是兄弟关系,A是B的兄弟,这两个分布成共轭分布关系,A是B的共轭分布。 p(θ|X)=p(θ)p(X|θ)p(x) p(X|θ) :似然(likelihood) p(θ) :先验(prior) p(X) :归一化常数(normalizing constant) 我们定义:如果先验分布(p(θ) )和似然函数(p(X|θ) )可以使得先验分布(p(θ) )和后验分布(p(θ|X) )有相同的形式(如,Beta(a+k, b+n-k)=Beta(a, b)binom(n, k)),那么就称先验分布与似然函数是共轭的(成Beta分布与二项分布是共轭的)。 几个常见的先验分布与其共轭分布 先验分布 共轭分布 伯努利分布 beta distribution Multinomial Dirichlet Distribution Gaussian, Given variance, mean unknown Gaussian Distribution Gaussian, Given mean, variance unknown Gamma Distribution Gaussian, both mean and variance unknown Gaussian-Gamma Distribution 最大似然估计(MLE) 首先来看,大名鼎鼎的贝叶斯公式: p(θ|X)=p(θ)p(X|θ)p(X) 可将θ 看成欲估计的分布的参数,X 表示样本,p(X|θ) 则表示似然。 现给定样本集\mathcal{D}=\{x_1,x_2,\ldots,x_N\}D={x1,x2,…,xN} ,似然函数为: p(\mathcal{D}|\theta)=\prod_{n=1}^Np(x_n|\theta) p(D|θ)=∏n=1Np(xn|θ) 为便于计算,再将其转换为对数似然函数形式: \ln p(\mathcal{D}|\theta)=\sum_{n=1}^N\ln p(x_n|\theta) lnp(D|θ)=∑n=1Nlnp(xn|θ) 我们不妨以伯努利分布为例,利用最大似然估计的方式计算其分布的参数(pp ),伯努利分布其概率密度函数(pdf)为: f_X(x)=p^x(1-p)^{1-x}=\left \{ \begin{array}{ll} p,&\mathrm{x=1},\\ q\equiv1-p ,&\mathrm{x=0},\\ 0,&\mathrm{otherwise} \end{array} \right. fX(x)=px(1−p)1−x=⎧⎩⎨⎪⎪p,q≡1−p,0,x=1,x=0,otherwise 整个样本集的对数似然函数为: \ln p(\mathcal{D}|\theta)=\sum_{n=1}^N\ln p(x_n|\theta)=\sum_{n=1}^N\ln (\theta^{x_n}(1-\theta)^{1-x_n})=\sum_{n=1}^Nx_n\ln\theta+(1-x_n)\ln(1-\theta) lnp(D|θ)=∑n=1Nlnp(xn|θ)=∑n=1Nln(θxn(1−θ)1−xn)=∑n=1Nxnlnθ+(1−xn)ln(1−θ) 等式两边对\thetaθ 求导: \frac{\partial \ln(\mathcal{D}|\theta)}{\partial \theta}=\frac{\sum_{n=1}^Nx_n}{\theta}-\frac{N}{1-\theta}+\frac{\sum_{n=1}^Nx_n}{1-\theta} ∂ln(D|θ)∂θ=∑Nn=1xnθ−N1−θ+∑Nn=1xn1−θ 令其为0,得: θml=∑Nn=1xnN Beta分布 f(μ|a,b)=Γ(a+b)Γ(a)Γ(b)μa−1(1−μ)b−1=1B(a,b)μa−1(1−μ)b−1 Beta 分布的峰值在a−1b+a−2 处取得。其中Γ(x)≡∫∞0ux−1e−udu 有如下性质: Γ(x+1)=xΓ(x)Γ(1)=1andΓ(n+1)=n! 我们来看当先验分布为 Beta 分布时的后验分布: p(θ)=1B(a,b)θa−1(1−θ)b−1p(X|θ)=(nk)θk(1−θ)n−kp(θ|X)=1B(a+k,b+n−k)θa+k−1(1−θ)b+n−k−1 对应于python中的math.gamma()及matlab中的gamma()函数(matlab中beta(a, b)=gamma(a)gamma(b)/gamma(a+b))。 条件概率(conditional probability) P(X|Y) 读作: P of X given Y ,下划线读作given X :所关心事件 Y :条件(观察到的,已发生的事件),conditional 条件概率的计算 仍然从样本空间(sample space)的角度出发。此时我们需要定义新的样本空间(给定条件之下的样本空间)。所以,所谓条件(conditional),本质是对样本空间的进一步收缩,或者叫求其子空间。 比如一个人答题,有A,B,C,D 四个选项,在答题者对题目一无所知的情况下,他答对的概率自然就是 14 ,而是如果具备一定的知识,排除了 A,C 两个错误选项,此时他答对的概率简单计算就增加到了 12 。 本质是样本空间从S={A,B,C,D} ,变为了S′={B,D} 。 新样本空间下P(A|排除A/C)=0,P(C|排除A/C)=0 ,归纳出来,也即某实验结果(outcome,oi )与某条件Y 不相交,则: P(oi|Y)=0 最后我们得到条件概率的计算公式: P(oi|Y)=P(oi)P(o1)+P(o2)+⋯+P(on)=P(oi)P(Y)Y={o1,o2,…,on} 考虑某事件X={o1,o2,q1,q2} ,已知条件Y={o1,o2,o3} 发生了,则: P(X|Y)=P(o1|Y)+P(o2|Y)+0+0=P(o1)P(Y)+P(o2)P(Y)=P(X∩Y)P(Y) 条件概率与贝叶斯公式 条件概率: P(X|Y)=P(X∩Y)P(Y) 贝叶斯公式: P(X|Y)=P(X)P(Y|X)P(Y) 其实是可从条件概率推导贝叶斯公式的: P(A|B)=P(B|A)=P(A|B)P(B)===P(B|A)=P(A∩B)P(B)P(A∩B)P(A)P(A∩B)P(B)P(B)P(A∩B)P(A)P(B|A)P(A|B)P(B)P(A) 证明:P(B,p|D)=P(B|p,D)P(p|D) P(B,p|D)====P(B,p,D)P(D)P(B|p,D)P(p,D)P(D)P(B|p,D)P(p,D)P(D)P(B|p,D)P(p|D) References [1] 概率质量函数 本篇文章为转载内容。原文链接:https://blog.csdn.net/lanchunhui/article/details/49799405。 该文由互联网用户投稿提供,文中观点代表作者本人意见,并不代表本站的立场。 作为信息平台,本站仅提供文章转载服务,并不拥有其所有权,也不对文章内容的真实性、准确性和合法性承担责任。 如发现本文存在侵权、违法、违规或事实不符的情况,请及时联系我们,我们将第一时间进行核实并删除相应内容。
2024-02-26 12:45:04
517
转载
Beego
用户权限管理问题:如何在Beego中构建安全的系统 一、前言 嘿,各位开发者朋友们,今天我们要聊的是一个在Web开发中极为重要的话题——用户权限管理。特别是在用Beego框架开发的时候,怎么才能高效地搞定权限控制,保证系统安全和用户隐私不被侵犯呢?这事儿可不只是技术活儿,简直就像是在搞艺术一样!在这篇文章里,咱们一起来捣鼓一下怎么在Beego项目中整出一个既高效又安全的用户权限管理系统。 二、为什么要重视用户权限管理? 想象一下,如果你的应用程序没有良好的权限管理体系,那么可能会发生以下情况: - 普通用户能够访问到管理员级别的功能。 - 系统数据可能被恶意篡改。 - 用户的敏感信息(如密码、银行卡号等)可能泄露。 这些后果都是灾难性的,不仅损害了用户对你的信任,也可能导致法律上的麻烦。所以啊,每个开发者都得认真搞个牢靠的权限控制系统,不然麻烦可就大了。 三、Beego中的权限管理基础 Beego本身并没有内置的权限管理模块,但我们可以利用其灵活的特性来构建自己的权限控制系统。以下是几种常见的实现方式: 1. 基于角色的访问控制(RBAC) - 这是一种常用的权限管理模型,它通过将权限分配给角色,再将角色分配给用户的方式简化了权限管理。 - 示例代码: go type Role struct { ID int64 Name string } type User struct { ID int64 Username string Roles []Role // 用户可以拥有多个角色 } func (u User) HasPermission(permission string) bool { for _, role := range u.Roles { if role.Name == permission { return true } } return false } 2. JWT(JSON Web Token)认证 - JWT允许你在不依赖于服务器端会话的情况下验证用户身份,非常适合微服务架构。 - 示例代码: go package main import ( "github.com/astaxie/beego" "github.com/dgrijalva/jwt-go" "net/http" "time" ) var jwtSecret = []byte("your_secret_key") type Claims struct { Username string json:"username" jwt.StandardClaims } func loginHandler(c beego.Context) { username := c.Input().Get("username") password := c.Input().Get("password") // 这里应该有验证用户名和密码的逻辑 token := jwt.NewWithClaims(jwt.SigningMethodHS256, Claims{ Username: username, StandardClaims: jwt.StandardClaims{ ExpiresAt: time.Now().Add(time.Hour 72).Unix(), }, }) tokenString, err := token.SignedString(jwtSecret) if err != nil { c.Ctx.ResponseWriter.WriteHeader(http.StatusInternalServerError) return } c.Data[http.StatusOK] = []byte(tokenString) } func authMiddleware() beego.ControllerFunc { return func(c beego.Controller) { tokenString := c.Ctx.Request.Header.Get("Authorization") token, err := jwt.ParseWithClaims(tokenString, &Claims{}, func(token jwt.Token) (interface{}, error) { return jwtSecret, nil }) if claims, ok := token.Claims.(Claims); ok && token.Valid { // 将用户信息存储在session或者全局变量中 c.SetSession("user", claims.Username) c.Next() } else { c.Ctx.ResponseWriter.WriteHeader(http.StatusUnauthorized) } } } 3. 中间件与拦截器 - 利用Beego的中间件机制,我们可以为特定路由添加权限检查逻辑,从而避免重复编写相同的权限校验代码。 - 示例代码: go func AuthRequiredMiddleware() beego.ControllerFunc { return func(c beego.Controller) { if !c.GetSession("user").(string) { c.Redirect("/login", 302) return } c.Next() } } func init() { beego.InsertFilter("/admin/", beego.BeforeRouter, AuthRequiredMiddleware) } 四、实际应用案例分析 让我们来看一个具体的例子,假设我们正在开发一款在线教育平台,需要对不同类型的用户(学生、教师、管理员)提供不同的访问权限。例如,只有管理员才能删除课程,而学生只能查看课程内容。 1. 定义用户类型 - 我们可以通过枚举类型来表示不同的用户角色。 - 示例代码: go type UserRole int const ( Student UserRole = iota Teacher Admin ) 2. 实现权限验证逻辑 - 在每个需要权限验证的操作之前,我们都需要先判断当前登录用户是否具有相应的权限。 - 示例代码: go func deleteCourse(c beego.Controller) { if userRole := c.GetSession("role"); userRole != Admin { c.Ctx.ResponseWriter.WriteHeader(http.StatusForbidden) return } // 执行删除操作... } 五、总结与展望 通过上述讨论,我们已经了解了如何在Beego框架下实现基本的用户权限管理系统。当然,实际应用中还需要考虑更多细节,比如异常处理、日志记录等。另外,随着业务越做越大,你可能得考虑引入一些更复杂的权限管理系统了,比如可以根据不同情况灵活调整的权限分配,或者可以精细到每个小细节的权限控制。这样能让你的系统管理起来更灵活,也更安全。 最后,我想说的是,无论采用哪种方法,最重要的是始终保持对安全性的高度警惕,并不断学习最新的安全知识和技术。希望这篇文章能对你有所帮助! --- 希望这样的风格和内容符合您的期待,如果有任何具体需求或想要进一步探讨的部分,请随时告诉我!
2024-10-31 16:13:08
166
初心未变
转载文章
...PC++基于行业标准和开放规范,旨在鼓励生态系统的协作和创新。 多架构编程面临的挑战 在以数据为中心的环境中,专用工作负载的数量不断增长。专用负载通常因为没有通用的编程语言或API而需要使用不同的语言和库进行编程,这就需要维护各自独立的代码库。 由于跨平台的工具支持不一致,因此开发人员必须学习和使用一整套不同的工具。单独投入精力给每种硬件平台开发软件。 oneAPI则可以利用一种统一的编程模型以及支持并行性的库,支持包括CPU、GPU、FPGA等硬件等同于原生高级语言的开发性能,并且可以与现有的HPC编程模型交互。 SYCL SYCL支持C++数据并行编程,SYCL和OpenCL一样都是由Khronos Group管理的,SYCL是建立在OpenCL之上的跨平台抽象层,支持用C++用单源语言方式编写用于异构处理器的与设备无关的代码。 DPC++ DPC++(Data Parallel C++)是一种单源语言,可以将主机代码和异构加速器内核写在同一个文件当中,在主机中调用DPC++程序,计算由加速器执行。DPC++代码简洁且效率高,并且是开源的。现有的CUDA应用、Fortran应用、OpenCL应用都可以用不同方式很方便地迁移到DPC++当中。 下图显示了原来使用不同架构的HPC开发人员的一些推荐的转换方法。 编译和运行DPC++程序 编译和运行DPC++程序主要包括三步: 初始化环境变量 编译DPC++源代码 运行程序 例如本地运行,在本地系统上安装英特尔基础工具套件,使用以下命令编译和运行DPC++程序。 source /opt/intel/inteloneapi/setvars.shdpcpp simple.cpp -o simple./simple 编程实例 实现矢量加法 以下实例描述了使用DPC++实现矢量加法的过程和源代码。 queue类 queue类用来提交给SYCL执行的命令组,是将作业提交到运算设备的一种机制,多个queue可以映射到同一个设备。 Parallel kernel Parallel kernel允许代码并行执行,对于一个不具有相关性的循环数据操作,可以用Parallel kernel并行实现 在C++代码中的循环实现 for(int i=0; i < 1024; i++){a[i] = b[i] + c[i];}); 在Parallel kernel中的并行实现 h.parallel_for(range<1>(1024), [=](id<1> i){A[i] = B[i] + C[i];}); 通用的并行编程模板 h.parallel_for(range<1>(1024), [=](id<1> i){// CODE THAT RUNS ON DEVICE }); range用来生成一个迭代序列,1为步长,在循环体中,i表示索引。 Host Accessor Host Accessor是使用主机缓冲区访问目标的访问器,它使访问的数据可以在主机上使用。通过构建Host Accessor可以将数据同步回主机,除此之外还可以通过销毁缓冲区将数据同步回主机。 buf是存储数据的缓冲区。 host_accessor b(buf,read_only); 除此之外还可以将buf设置为局部变量,当系统超出buf生存期,buf被销毁,数据也将转移到主机中。 矢量相加源代码 根据上面的知识,这里展示了利用DPC++实现矢量相加的代码。 //第一行在jupyter中指明了该cpp文件的保存位置%%writefile lab/vector_add.cppinclude <CL/sycl.hpp>using namespace sycl;int main() {const int N = 256;// 初始化两个队列并打印std::vector<int> vector1(N, 10);std::cout<<"\nInput Vector1: "; for (int i = 0; i < N; i++) std::cout << vector1[i] << " ";std::vector<int> vector2(N, 20);std::cout<<"\nInput Vector2: "; for (int i = 0; i < N; i++) std::cout << vector2[i] << " ";// 创建缓存区buffer vector1_buffer(vector1);buffer vector2_buffer(vector2);// 提交矢量相加任务queue q;q.submit([&](handler &h) {// 为缓存区创建访问器accessor vector1_accessor (vector1_buffer,h);accessor vector2_accessor (vector2_buffer,h);h.parallel_for(range<1>(N), [=](id<1> index) {vector1_accessor[index] += vector2_accessor[index];});});// 创建主机访问器将设备中数据拷贝到主机当中host_accessor h_a(vector1_buffer,read_only);std::cout<<"\nOutput Values: ";for (int i = 0; i < N; i++) std::cout<< vector1[i] << " ";std::cout<<"\n";return 0;} 运行结果 统一共享内存 (Unified Shared Memory USM) 统一共享内存是一种基于指针的方法,是将CPU内存和GPU内存进行统一的虚拟化方法,对于C++来说,指针操作内存是很常规的方式,USM也可以最大限度的减少C++移植到DPC++的代价。 下图显示了非USM(左)和USM(右)的程序员开发视角。 类型 函数调用 说明 在主机上可访问 在设备上可访问 设备 malloc_device 在设备上分配(显式) 否 是 主机 malloc_host 在主机上分配(隐式) 是 是 共享 malloc_shared 分配可以在主机和设备之间迁移(隐式) 是 是 USM语法 初始化: int data = malloc_shared<int>(N, q); int data = static_cast<int >(malloc_shared(N sizeof(int), q)); 释放 free(data,q); 使用共享内存之后,程序将自动在主机和运算设备之间隐式移动数据。 数据依赖 使用USM时,要注意数据之间的依赖关系以及事件之间的依赖关系,如果两个线程同时修改同一个内存区,将产生不可预测的结果。 我们可以使用不同的选项管理数据依赖关系: 内核任务中的 wait() 使用 depends_on 方法 使用 in_queue 队列属性 wait() q.submit([&](handler &h) {h.parallel_for(range<1>(N), [=](id<1> i) { data[i] += 2; });}).wait(); // <--- wait() will make sure that task is complete before continuingq.submit([&](handler &h) {h.parallel_for(range<1>(N), [=](id<1> i) { data[i] += 3; });}); depends_on auto e = q.submit([&](handler &h) { // <--- e is event for kernel taskh.parallel_for(range<1>(N), [=](id<1> i) { data[i] += 2; });});q.submit([&](handler &h) {h.depends_on(e); // <--- waits until event e is completeh.parallel_for(range<1>(N), [=](id<1> i) { data[i] += 3; });}); in_order queue property queue q(property_list{property::queue::in_order()}); // <--- this will make sure all the task with q are executed sequentially 练习1:事件依赖 以下代码使用 USM,并有三个提交到设备的内核。每个内核修改相同的数据阵列。三个队列之间没有数据依赖关系 为每个队列提交添加 wait() 在第二个和第三个内核任务中实施 depends_on() 方法 使用 in_order 队列属性,而非常规队列: queue q{property::queue::in_order()}; %%writefile lab/usm_data.cppinclude <CL/sycl.hpp>using namespace sycl;static const int N = 256;int main() {queue q{property::queue::in_order()};//用队列限制执行顺序std::cout << "Device : " << q.get_device().get_info<info::device::name>() << "\n";int data = static_cast<int >(malloc_shared(N sizeof(int), q));for (int i = 0; i < N; i++) data[i] = 10;q.parallel_for(range<1>(N), [=](id<1> i) { data[i] += 2; });q.parallel_for(range<1>(N), [=](id<1> i) { data[i] += 3; });q.parallel_for(range<1>(N), [=](id<1> i) { data[i] += 5; });q.wait();//wait阻塞进程for (int i = 0; i < N; i++) std::cout << data[i] << " ";std::cout << "\n";free(data, q);return 0;} 执行结果 练习2:事件依赖 以下代码使用 USM,并有三个提交到设备的内核。前两个内核修改了两个不同的内存对象,第三个内核对前两个内核具有依赖性。三个队列之间没有数据依赖关系 %%writefile lab/usm_data2.cppinclude <CL/sycl.hpp>using namespace sycl;static const int N = 1024;int main() {queue q;std::cout << "Device : " << q.get_device().get_info<info::device::name>() << "\n";//设备选择int data1 = malloc_shared<int>(N, q);int data2 = malloc_shared<int>(N, q);for (int i = 0; i < N; i++) {data1[i] = 10;data2[i] = 10;}auto e1 = q.parallel_for(range<1>(N), [=](id<1> i) { data1[i] += 2; });auto e2 = q.parallel_for(range<1>(N), [=](id<1> i) { data2[i] += 3; });//e1,e2指向两个事件内核q.parallel_for(range<1>(N),{e1,e2}, [=](id<1> i) { data1[i] += data2[i]; }).wait();//depend on e1,e2for (int i = 0; i < N; i++) std::cout << data1[i] << " ";std::cout << "\n";free(data1, q);free(data2, q);return 0;} 运行结果 UMS实验 在主机中初始化两个vector,初始数据为25和49,在设备中初始化两个vector,将主机中的数据拷贝到设备当中,在设备当中并行计算原始数据的根号值,然后将data1_device和data2_device的数值相加,最后将数据拷贝回主机当中,检验最后相加的和是否是12,程序结束前将内存释放。 %%writefile lab/usm_lab.cppinclude <CL/sycl.hpp>include <cmath>using namespace sycl;static const int N = 1024;int main() {queue q;std::cout << "Device : " << q.get_device().get_info<info::device::name>() << "\n";//intialize 2 arrays on hostint data1 = static_cast<int >(malloc(N sizeof(int)));int data2 = static_cast<int >(malloc(N sizeof(int)));for (int i = 0; i < N; i++) {data1[i] = 25;data2[i] = 49;}// STEP 1 : Create USM device allocation for data1 and data2int data1_device = static_cast<int >(malloc_device(N sizeof(int),q));int data2_device = static_cast<int >(malloc_device(N sizeof(int),q));// STEP 2 : Copy data1 and data2 to USM device allocationq.memcpy(data1_device, data1, sizeof(int) N).wait();q.memcpy(data2_device, data2, sizeof(int) N).wait();// STEP 3 : Write kernel code to update data1 on device with sqrt of valueauto e1 = q.parallel_for(range<1>(N), [=](id<1> i) { data1_device[i] = std::sqrt(25); });auto e2 = q.parallel_for(range<1>(N), [=](id<1> i) { data2_device[i] = std::sqrt(49); });// STEP 5 : Write kernel code to add data2 on device to data1q.parallel_for(range<1>(N),{e1,e2}, [=](id<1> i) { data1_device[i] += data2_device[i]; }).wait();// STEP 6 : Copy data1 on device to hostq.memcpy(data1, data1_device, sizeof(int) N).wait();q.memcpy(data2, data2_device, sizeof(int) N).wait();// verify resultsint fail = 0;for (int i = 0; i < N; i++) if(data1[i] != 12) {fail = 1; break;}if(fail == 1) std::cout << " FAIL"; else std::cout << " PASS";std::cout << "\n";// STEP 7 : Free USM device allocationsfree(data1_device, q);free(data1);free(data2_device, q);free(data2);// STEP 8 : Add event based kernel dependency for the Steps 2 - 6return 0;} 运行结果 本篇文章为转载内容。原文链接:https://blog.csdn.net/MCKZX/article/details/127630566。 该文由互联网用户投稿提供,文中观点代表作者本人意见,并不代表本站的立场。 作为信息平台,本站仅提供文章转载服务,并不拥有其所有权,也不对文章内容的真实性、准确性和合法性承担责任。 如发现本文存在侵权、违法、违规或事实不符的情况,请及时联系我们,我们将第一时间进行核实并删除相应内容。
2023-07-22 10:28:50
321
转载
转载文章
...上传早已不是什么新鲜问题,在.net 2.0时代,HTML5也还没有问世,要实现这样的功能,要么是改web.config,要么是用flash,要么是用一些第三方控件,然而这些解决问题的方法要么很麻烦,比如改配置,要么不稳定,比如文件上G以后,上传要么死掉,要么卡住,通过设置web.config并不能很好的解决这些问题。 这是一个Html5统治浏览器的时代,在这个新的时代,这种问题已被简化并解决,我们可以利用Html5分片上传的技术,那么Plupload则是一个对此技术进行封装的前端脚本库,这个库的好处是可以自动检测浏览器是否支持html5技术,不支持再检测是否支持flash技术,甚至是sliverlight技术,如果支持,就使用检测到的技术。 那么这个库到哪里下载,怎么搭建呢,比较懒的童鞋还是用Install-Package Plupload搞定吧,一个命令搞定所有事 Plupload支持的功能这里就不细说了,什么批量上传,这里我没有用到,主要是感觉它支持的事件非常丰富,文件选取后的事件,文件上传中的事件(可获得文件的上传进度),文件上传成功的事件,文件上传失败的事件,等等 我的例子主要是上传一个单个文件,并显示上传的进度条(使用jQuery的一个进度条插件) 下面的例子主要是为文件上传交给 UploadCoursePackage.ashx 来处理 /ProgressBar/ var progressBar = $("loading").progressbar({ width: '500px', color: 'B3240E', border: '1px solid 000000' }); /Plupload/ //实例化一个plupload上传对象 var uploader = new plupload.Uploader({ browse_button: 'browse', //触发文件选择对话框的按钮,为那个元素id runtimes: 'html5,flash,silverlight,html4',//兼容的上传方式 url: "Handlers/UploadCoursePackage.ashx", //后端交互处理地址 max_retries: 3, //允许重试次数 chunk_size: '10mb', //分块大小 rename: true, //重命名 dragdrop: false, //允许拖拽文件进行上传 unique_names: true, //文件名称唯一性 filters: { //过滤器 max_file_size: '999999999mb', //文件最大尺寸 mime_types: [ //允许上传的文件类型 { title: "Zip", extensions: "zip" }, { title: "PE", extensions: "pe" } ] }, //自定义参数 (键值对形式) 此处可以定义参数 multipart_params: { type: "misoft" }, // FLASH的配置 flash_swf_url: "../Scripts/plupload/Moxie.swf", // Silverligh的配置 silverlight_xap_url: "../Scripts/plupload/Moxie.xap", multi_selection: false //true:ctrl多文件上传, false 单文件上传 }); //在实例对象上调用init()方法进行初始化 uploader.init(); uploader.bind('FilesAdded', function (uploader, files) { $("<%=fileSource.ClientID %>").val(files[0].name); $.ajax( { type: 'post', url: 'HardDiskSpace.aspx/GetHardDiskFreeSpace', data: {}, dataType: 'json', contentType: 'application/json;charset=utf-8', success: function (result) { //选择文件以后检测服务器剩余磁盘空间是否够用 if (files.length > 0) { if (parseInt(files[0].size) > parseInt(result.d)) { $('error-msg').text("文件容量大于剩余磁盘空间,请联系管理员!"); } else { $('error-msg').text(""); } } }, error: function (xhr, err, obj) { $('error-msg').text("检测服务器剩余磁盘空间失败"); } }); }); uploader.bind('UploadProgress', function (uploader, file) { var percent = file.percent; progressBar.progress(percent); }); uploader.bind('FileUploaded', function (up, file, callBack) { var data = $.parseJSON(callBack.response); if (data.statusCode === "1") { $("<%=hfPackagePath.ClientID %>").val(data.filePath); var id = $("<%=hfCourseID.ClientID %>").val(); __doPostBack("save", id); } else { hideLoading(); $('error-msg').text(data.message); } }); uploader.bind('Error', function (up, err) { alert("文件上传失败,错误信息: " + err.message); }); /Plupload/ 后台 UploadCoursePackage.ashx 的代码也重要,主要是文件分片跟不分片的处理方式不一样 using System; using System.Collections.Generic; using System.Linq; using System.Web; using System.IO; namespace WebUI.Handlers { /// <summary> /// UploadCoursePackage 的摘要说明 /// </summary> public class UploadCoursePackage : IHttpHandler { public void ProcessRequest(HttpContext context) { context.Response.ContentType = "text/plain"; int statuscode = 1; string message = string.Empty; string filepath = string.Empty; if (context.Request.Files.Count > 0) { try { string resourceDirectoryName = System.Configuration.ConfigurationManager.AppSettings["resourceDirectory"]; string path = context.Server.MapPath("~/" + resourceDirectoryName); if (!Directory.Exists(path)) Directory.CreateDirectory(path); int chunk = context.Request.Params["chunk"] != null ? int.Parse(context.Request.Params["chunk"]) : 0; //获取当前的块ID,如果不是分块上传的。chunk则为0 string fileName = context.Request.Params["name"]; //这里写的比较潦草。判断文件名是否为空。 string type = context.Request.Params["type"]; //在前面JS中不是定义了自定义参数multipart_params的值么。其中有个值是type:"misoft",此处就可以获取到这个值了。获取到的type="misoft"; string ext = Path.GetExtension(fileName); //fileName = string.Format("{0}{1}", Guid.NewGuid().ToString(), ext); filepath = resourceDirectoryName + "/" + fileName; fileName = Path.Combine(path, fileName); //对文件流进行存储 需要注意的是 files目录必须存在(此处可以做个判断) 根据上面的chunk来判断是块上传还是普通上传 上传方式不一样 ,导致的保存方式也会不一样 FileStream fs = new FileStream(fileName, chunk == 0 ? FileMode.OpenOrCreate : FileMode.Append); //write our input stream to a buffer Byte[] buffer = null; if (context.Request.ContentType == "application/octet-stream" && context.Request.ContentLength > 0) { buffer = new Byte[context.Request.InputStream.Length]; context.Request.InputStream.Read(buffer, 0, buffer.Length); } else if (context.Request.ContentType.Contains("multipart/form-data") && context.Request.Files.Count > 0 && context.Request.Files[0].ContentLength > 0) { buffer = new Byte[context.Request.Files[0].InputStream.Length]; context.Request.Files[0].InputStream.Read(buffer, 0, buffer.Length); } //write the buffer to a file. if (buffer != null) fs.Write(buffer, 0, buffer.Length); fs.Close(); statuscode = 1; message = "上传成功"; } catch (Exception ex) { statuscode = -1001; message = "保存时发生错误,请确保文件有效且格式正确"; Util.LogHelper logger = new Util.LogHelper(); string path = context.Server.MapPath("~/Logs"); logger.WriteLog(ex.Message, path); } } else { statuscode = -404; message = "上传失败,未接收到资源文件"; } string msg = "{\"statusCode\":\"" + statuscode + "\",\"message\":\"" + message + "\",\"filePath\":\"" + filepath + "\"}"; context.Response.Write(msg); } public bool IsReusable { get { return false; } } } } 再附送一个检测服务器端硬盘剩余空间的功能吧 using System; using System.Collections.Generic; using System.IO; using System.Linq; using System.Web; using System.Web.Script.Services; using System.Web.Services; using System.Web.UI; using System.Web.UI.WebControls; namespace WebUI { public partial class CheckHardDiskFreeSpace : System.Web.UI.Page { protected void Page_Load(object sender, EventArgs e) { } /// <summary> /// 获取磁盘剩余容量 /// </summary> /// <returns></returns> [WebMethod] public static string GetHardDiskFreeSpace() { const string strHardDiskName = @"F:\"; var freeSpace = string.Empty; var drives = DriveInfo.GetDrives(); var myDrive = (from drive in drives where drive.Name == strHardDiskName select drive).FirstOrDefault(); if (myDrive != null) { freeSpace = myDrive.TotalFreeSpace+""; } return freeSpace; } } } 效果展示: 详细配置信息可以参考这篇文章:http://blog.ncmem.com/wordpress/2019/08/12/plupload%e4%b8%8a%e4%bc%a0%e6%95%b4%e4%b8%aa%e6%96%87%e4%bb%b6%e5%a4%b9-2/ 本篇文章为转载内容。原文链接:https://blog.csdn.net/weixin_45525177/article/details/100654639。 该文由互联网用户投稿提供,文中观点代表作者本人意见,并不代表本站的立场。 作为信息平台,本站仅提供文章转载服务,并不拥有其所有权,也不对文章内容的真实性、准确性和合法性承担责任。 如发现本文存在侵权、违法、违规或事实不符的情况,请及时联系我们,我们将第一时间进行核实并删除相应内容。
2023-12-19 09:43:46
127
转载
转载文章
...u_class) 将数值型的变量指定为分类型anova_lm(lm, typ = 2) 可以看到不同教育水平之间的月均消费支出之间的差异是显著的,继续用多重检验来看哪些处理之间是显著的。 from statsmodels.stats.multicomp import MultiComparison 引入MultiComparison进行tukey多重检验mc = MultiComparison(raw_1['avg_exp'],raw_1['edu_class'])tukey_result = mc.tukeyhsd(alpha = 0.5)print(tukey_result) 结果是每个处理之间因变量差异的显著性,最后一列reject都为True说明各组之间均存在显著差异。 三、模型建立与诊断 3.1 一元线性回归及模型解读 以Income为自变量,以avg_exp为因变量建立一元线形回归并对模型结果进行解释 lm_1 = ols('avg_exp ~ Income', data = raw_1).fit()print(lm_1.summary()) 首先从第一部分可以看到R^2为0.454,整个模型的F检验p值小于0.05,说明模型通过显著性检验。 其次模型结果的第二块也表明自变量和截距也通过显著性检验。 最后一部分主要是对残差进行检验,左侧Omnibus、Prob(Omnibus)主要是对偏度Skew和峰度Kurtosis进行检验,正态分布的偏度为0,峰度为3,模型的Prob(Omnibus)值为0.156大于0.05,说明不能拒绝残差符合正态分布。 右侧Durbin-Watson主要是对残差的自相关性进行检(改检验可表示为,为残差之间的相关系数),Durbin-Watson的取值范围是0-4,越接近2说明残差不存在自相关性,越接近0说明存在正相关,越接近4说明存在负相关性。 右侧Jarque-Bera (JB)、Prob(JB)是对残差正态性检验,可以用来判断残差是否符合正态分布,本案例中Prob(JB)值为0.173 > 0.05,基不能拒绝残差服从正态分布。 右侧Cond. No.是多重共线性检验,该值越大,共线性越严重。 整体上看模型虽然拟合效果没那么好,但是显著性通过了检验。接下来看一下模型具体的系数,Income的系数为97.7说明模型收入越高信用卡消费越高,是符合业务预期的。 3.2 残差可视化分析 接下来对残差进一步进行可视化分析,主要看残差是否满足以下几个假定,并尝试通过对自变量、因变量进行调整来优化模型。首先来回顾一下残差需要满足的几个假定: a.残差的要服从均值为0,方差为的正态分布; b.残差之间要相互独立 c.残差和自变量没有相关性 (1)通过残差图进行模型优化 模型avg_exp ~ Income的自变量与残差分布图、残差qq图、模型拟合情况图即自变量与因变量及其预测值的图像 lm_1 = ols('avg_exp ~ Income', data = raw_1).fit() 建模raw_1['resid_1'] = lm_1.resid 模型残差raw_1['resid_1_rank'] = raw_1['resid_1'].rank(ascending = False, pct = True) 计算残差的百分位数raw_1['pred_1'] = lm_1.predict() 添加预测值plt.figure(figsize = (20, 6)) 自变量与残差分布图ax1 = plt.subplot(131)ax1.scatter('Income', 'resid', data = raw_1)ax1.set_title('Income & resid') 残差的qq图ax2 = plt.subplot(132)stats.probplot(raw_1['resid_1_rank'], dist = 'norm', plot = ax2) 模型拟合情况图,自变量与因变量以及模型预测值ax3 = plt.subplot(133)ax3.scatter('Income', 'avg_exp', data = raw_1)ax3.plot('Income', 'pred_1', data = raw_1, color = 'red')ax3.legend()ax3.text(12, 1920, 'pred func R^2: %.2f'% lm_1.rsquared)ax3.set_title('Income & avg_exp') 从第一个自变量和残差散点图可以看出,残差基本符合对称分布,但随着自变量增大,残差也在变大,存在方差不齐的情况。第二个图残差的qq图可以看出,残差近似正态分布。第三个图可以看模型的拟合效果并不是很好,R^2只有0.45。对avg_exp取对数,能够改善预测值越大残差越大的情况,但由于只对因变量取对数导致模型不好解释,对自变量Income同时取对数,代码和以上类似,只是改变因变量和自变量形式而已,以下是残差图,可以看到残差的异方差现象被有效的抑制,并且R^2也得到了提高。 (2)通过残差图发现强影响点 仔细观察以上图像结果,左下侧有两个较为异常的数据,对模型的拟和效果有较大的影响, 对于这种影响较大的可将其进行删除并重新建模: 计算学生化残差raw_1['resid_t'] = (raw_1['resid_2'] - raw_1['resid_2'].mean())/raw_1['resid_2'].std() raw_1[abs(raw_1['resid_t']) > 2] 将残差大于2的筛选出来 将强影响点删除后,得到的结果如下,模型结果更稳定。 3.3 多元线性回归 上一篇文章有说到多重共线性会对模型产生致命的影响,用方差膨胀因子来处理的话会非常繁琐。通过正则化处理如Lasso回归,能够产生某些严格等于0的系数,从而达到变量筛选的目的。接下来以Lasso为例,首先用LassoCV来找到最优的alpha。由于statsmodels中的ols的fit_regularized方法没有很好的实现,所以用sklearn中linear_model模块来进行建模 from sklearn.preprocessing import StandardScaler sklearn进行线性回归前必须要进行标准化from sklearn.linear_model import LassoCV Lasso的交叉验证方法con_xcols = ['Age', 'Income', 'dist_home_val', 'dist_avg_income']scaler = StandardScaler()X = scaler.fit_transform(raw_1[con_xcols])y = raw_1['avg_exp_ln']lasso_alphas = np.logspace(-3, 0, 100, base = 10)lcv = LassoCV(alphas = lasso_alphas, cv = 10)lcv.fit(X, y)print('best alpha %.4f' % lcv.alpha_)print('the r-square %.4f' % lcv.score(X, y)) 接下来画出不同alpha下的岭迹图,来看alpha值对系数的影响 from sklearn.linear_model import Lassocoefs = []lasso = Lasso()for i in lasso_alphas:lasso.set_params(alpha = i)lasso.fit(X, y)coefs.append(lasso.coef_)ax = plt.gca()ax.plot(lasso_alphas, coefs)ax.set_xscale('log')ax.set_xlabel('$\\alpha$')ax.set_ylabel('coefs value') 从图中可以看到随着alpha的增大,系数不断在减小,有些系数会优先收缩为0,再继续增大时所欲系数都会为0,通过该特性从而达到变量筛选的目的。将LassoCV得到的系数打印出来,可以看到用户月均信用卡支出和当地小区均价、当地人均收入成正比,当地人均收入水平的影响更大。 以上就是线形回归在应用时的注意事项。 本篇文章为转载内容。原文链接:https://blog.csdn.net/baidu_26137595/article/details/123766191。 该文由互联网用户投稿提供,文中观点代表作者本人意见,并不代表本站的立场。 作为信息平台,本站仅提供文章转载服务,并不拥有其所有权,也不对文章内容的真实性、准确性和合法性承担责任。 如发现本文存在侵权、违法、违规或事实不符的情况,请及时联系我们,我们将第一时间进行核实并删除相应内容。
2023-11-23 15:52:56
106
转载
转载文章
...式。它们可以实现模拟标准类型和重载操作符等。比如__init__()和__del__()就分别充当了构造器和析够器的功能。 这些特殊这些方法都是以双下划线(__)开始及结尾的。下表进行了总结: 基本定制型 C.__init__(self[, arg1, ...]) 构造器(带一些可选的参数) C.__new__(self[, arg1, ...]) 构造器(带一些可选的参数);通常用在设置不变数据类型的子类。 C.__del__(self) 解构器 C.__str__(self) 可打印的字符输出;内建str()及print 语句 C.__repr__(self) 运行时的字符串输出;内建repr() 和操作符 C.__unicode__(self) Unicode 字符串输出;内建unicode() C.__call__(self, args) 表示可调用的实例 C.__nonzero__(self) 为object 定义False 值;内建bool() (从2.2 版开始) C.__len__(self) “长度”(可用于类);内建len() 对象(值)比较 C.__cmp__(self, obj) 对象比较;内建cmp() C.__lt__(self, obj) C.__le__(self, obj) 小于/小于或等于;对应 C.__gt__(self, obj) C.__ge__(self, obj) 大于/大于或等于;对应>及>=操作符 C.__eq__(self, obj) C.__ne__(self, obj) 等于/不等于;对应==,!=及<>操作符 属性 C.__getattr__(self, attr) 获取属性;内建getattr();仅当属性没有找到时调用 C.__setattr__(self, attr, val) 设置属性 C.__delattr__(self, attr) 删除属性 C.__getattribute__(self, attr) 获取属性;内建getattr();总是被调用 C.__get__(self, attr) (描述符)获取属性 C.__set__(self, attr, val) (描述符)设置属性 C.__delete__(self, attr) (描述符)删除属性 数值类型:二进制操作符 C.__add__(self, obj) 加;+操作符 C.__sub__(self, obj) 减;-操作符 C.__mul__(self, obj) 乘;操作符 C.__div__(self, obj) 除;/操作符 C.__truediv__(self, obj) True 除;/操作符 C.__floordiv__(self, obj) Floor 除;//操作符 C.__mod__(self, obj) 取模/取余;%操作符 C.__divmod__(self, obj) 除和取模;内建divmod() C.__pow__(self, obj[, mod]) 乘幂;内建pow();操作符 C.__lshift__(self, obj) 左移位;< 数值类型:二进制操作符 C.__rshift__(self, obj) 右移;>>操作符 C.__and__(self, obj) 按位与;&操作符 C.__or__(self, obj) 按位或;|操作符 C.__xor__(self, obj) 按位与或;^操作符 数值类型:一元操作符 C.__neg__(self) 一元负 C.__pos__(self) 一元正 C.__abs__(self) 绝对值;内建abs() C.__invert__(self) 按位求反;~操作符 数值类型:数值转换 C.__complex__(self, com) 转为complex(复数);内建complex() C.__int__(self) 转为int;内建int() C.__long__(self) 转 .long;内建long() C.__float__(self) 转为float;内建float() 数值类型:基本表示法(String) C.__oct__(self) 八进制表示;内建oct() C.__hex__(self) 十六进制表示;内建hex() 数值类型:数值压缩 C.__coerce__(self, num) 压缩成同样的数值类型;内建coerce() C.__index__(self) 在有必要时,压缩可选的数值类型为整型(比如:用于切片索引等等) 序列类型 C.__len__(self) 序列中项的数目 C.__getitem__(self, ind) 得到单个序列元素 C.__setitem__(self, ind,val) 设置单个序列元素 C.__delitem__(self, ind) 删除单个序列元素 C.__getslice__(self, ind1,ind2) 得到序列片断 C.__setslice__(self, i1, i2,val) 设置序列片断 C.__delslice__(self, ind1,ind2) 删除序列片断 C.__contains__(self, val) 测试序列成员;内建in 关键字 C.__add__(self,obj) 串连;+操作符 C.__mul__(self,obj) 重复;操作符 C.__iter__(self) 创建迭代类;内建iter() 映射类型 C.__len__(self) mapping 中的项的数目 C.__hash__(self) 散列(hash)函数值 C.__getitem__(self,key) 得到给定键(key)的值 C.__setitem__(self,key,val) 设置给定键(key)的值 C.__delitem__(self,key) 删除给定键(key)的值 C.__missing__(self,key) 给定键如果不存在字典中,则提供一个默认值 一:简单定制 classRoundFloatManual(object):def __init__(self, val):assert isinstance(val, float), "Value must be a float!"self.value= round(val, 2)>>> rfm =RoundFloatManual(42) Traceback (mostrecent call last): File"", line 1, in? File"roundFloat2.py", line 5, in __init__assertisinstance(val, float), \ AssertionError: Value must be a float!>>> rfm =RoundFloatManual(4.2)>>>rfm >>> printrfm 它因输入非法而异常,但如果输入正确时,就没有任何输出了。在解释器中,我们得到一些信息,却不是我们想要的。print(使用str())和真正的字符串对象表示(使用repr())都没能显示更多有关我们对象的信息。这就需要实现__str__()和__repr__()二者之一,或者两者都实现。加入下面的方法: def __str__(self):return str(self.value) 现在我们得到下面的: >>> rfm = RoundFloatManual(5.590464)>>>rfm >>> printrfm5.59 >>> rfm = RoundFloatManual(5.5964)>>> printrfm5.6 但是在解释器中转储(dump)对象时,仍然显示的是默认对象符号,要修复它,只需要覆盖__repr__()。可以让__repr__()和__str__()具有相同的代码,但最好的方案是:__repr__ = __str__ 在带参数5.5964的第二个例子中,我们看到它舍入值刚好为5.6,但我们还是想显示带两位小数的数。可以这样修改: def __str__(self):return '%.2f' % self.value 这里就同时具备str()和repr()的输出了: >>> rfm =RoundFloatManual(5.5964)>>>rfm5.60 >>>printrfm5.60 所有代码如下: classRoundFloatManual(object):def __init__(self,val):assert isinstance(val, float), "Valuemust be a float!"self.value= round(val, 2)def __str__(self):return '%.2f' %self.value__repr__ = __str__ 二:数值定制 定义一个Time60,其中,将整数的小时和分钟作为输入传给构造器: classTime60(object):def __init__(self, hr, min): self.hr=hr self.min= min 1:显示 需要在显示实例的时候,得到一个有意义的输出,那么就要覆盖__str__()(如果有必要的话,__repr__()也要覆盖): def __str__(self):return '%d:%d' % (self.hr, self.min) 比如: >>> mon =Time60(10, 30)>>> tue =Time60(11, 15)>>> >>> printmon, tue10:30 11:15 2:加法 Python中的重载操作符很简单。像加号(+),只需要重载__add__()方法,如果合适,还可以用__radd__()及__iadd__()。注意,实现__add__()的时候,必须认识到它返回另一个Time60对象,而不修改原mon或tue: def __add__(self, other):return self.__class__(self.hr + other.hr, self.min + other.min) 在类中,一般不直接调用类名,而是使用self 的__class__属性,即实例化self 的那个类,并调用它。调用self.__class__()与调用Time60()是一回事。但self.__class__()的方式更好。 >>> mon = Time60(10, 30)>>> tue = Time60(11, 15)>>> mon +tue >>> print mon +tue21:45 如果没有定义相对应的特殊方法,但是却使用了该方法对应的运算,则会引起一个TypeError异常: >>> mon -tue Traceback (mostrecent call last): File"", line 1, in? TypeError:unsupported operand type(s)for -: 'Time60' and 'Time60' 3:原位加法 __iadd__(),是用来支持像mon += tue 这样的操作符,并把正确的结果赋给mon。重载一个__i__()方法的唯一秘密是它必须返回self: def __iadd__(self, other): self.hr+=other.hr self.min+=other.minreturn self 下面是结果输出: >>> mon = Time60(10,30)>>> tue = Time60(11,15)>>>mon10:30 >>>id(mon)401872 >>> mon +=tue>>>id(mon)401872 >>>mon21:45 下面是Time60的类的完全定义: classTime60(object):'Time60 - track hours and minutes' def __init__(self,hr, min):'Time60 constructor - takes hours andminutes'self.hr=hr self.min=mindef __str__(self):'Time60 - string representation' return '%d:%d' %(self.hr, self.min)__repr__ = __str__ def __add__(self, other):'Time60 - overloading the additionoperator' return self.__class__(self.hr + other.hr,self.min +other.min)def __iadd__(self,other):'Time60 - overloading in-place addition'self.hr+=other.hr self.min+=other.minreturn self 4:升华 在这个类中,还有很多需要优化和改良的地方。首先看下面的例子: >>> wed =Time60(12, 5)>>>wed12:5 正确的显示应该是:“12:05” >>> thu =Time60(10, 30)>>> fri =Time60(8, 45)>>> thu +fri18:75 正确的显示应该是:19:15 可以做出如下修改: def __str__(self):return '%02d:%02d'%(self.hr, self.min)__repr__ = __str__ def __add__(self, othertime): tmin= self.min +othertime.min thr= self.hr +othertime.hrreturn self.__class__(thr + tmin/60, tmin%60)def __iadd__(self, othertime): self.min+=othertime.min self.hr+=othertime.hr self.hr+= self.min/60self.min%= 60 return self 三:迭代器 迭代器对象本身需要支持以下两种方法,它们组合在一起形成迭代器协议: iterator.__iter__() 返回迭代器对象本身。 iterator.next() 从容器中返回下一个元素。 实现了__iter__()和next()方法的类就是一个迭代器。自定义迭代器的例子如下: RandSeq(Random Sequence),传入一个初始序列,__init__()方法执行前述的赋值操作。__iter__()仅返回self,这就是如何将一个对象声明为迭代器的方式,最后,调用next()来得到迭代器中连续的值。这个迭代器唯一的亮点是它没有终点。代码如下: classRandSeq(object):def __init__(self, seq): self.data=seqdef __iter__(self):returnselfdefnext(self):return choice(self.data) 运行它,将会看到下面的输出: >>> from randseq importRandSeq>>> for eachItem in RandSeq(('rock', 'paper', 'scissors')): ...printeachItem ... scissors scissors rock paper paper scissors ...... 四:多类型定制 现在创建另一个新类,NumStr,由一个数字-字符对组成,记为n和s,数值类型使用整型(integer)。用[n::s]来表示它,这两个数据元素构成一个整体。NumStr有下面的特征: 初始化: 类应当对数字和字符串进行初始化;如果其中一个(或两)没有初始化,则使用0和空字符串,也就是, n=0 且s=''作为默认。 加法: 定义加法操作符,功能是把数字加起来,把字符连在一起;比如,NumStr1=[n1::s1]且NumStr2=[n2::s2]。则NumStr1+NumStr2 表示[n1+n2::s1+s2],其中,+代表数字相加及字符相连接。 乘法: 类似的, 定义乘法操作符的功能为, 数字相乘,字符累积相连, 也就是,NumStr1NumStr2=[n1n::s1n]。 False 值:当数字的数值为 0 且字符串为空时,也就是当NumStr=[0::'']时,这个实体即有一个false值。 比较: 比较一对NumStr对象,比如,[n1::s1] vs. [n2::s2],有九种不同的组合。对数字和字符串,按照标准的数值和字典顺序的进行比较。 如果obj1< obj2,则cmp(obj1, obj2)的返回值是一个小于0 的整数, 当obj1 > obj2 时,比较的返回值大于0, 当两个对象有相同的值时, 比较的返回值等于0。 我们的类的解决方案是把这些值相加,然后返回结果。为了能够正确的比较对象,我们需要让__cmp__()在 (n1>n2) 且 (s1>s2)时,返回 1,在(n1s2),或相反),返回0. 反之亦然。代码如下: classNumStr(object):def __init__(self, num=0, string=''): self.__num =num self.__string =stringdef __str__(self):return '[%d :: %r]' % (self.__num, self.__string)__repr__ = __str__ def __add__(self, other):ifisinstance(other, NumStr):return self.__class__(self.__num + other.__num, self.__string + other.__string)else:raise TypeError, 'Illegal argument type for built-in operation' def __mul__(self, num):ifisinstance(num, int):return self.__class__(self.__num num, self.__string num)else:raise TypeError, 'Illegal argument type for built-inoperation' def __nonzero__(self):return self.__num or len(self.__string)def __norm_cval(self, cmpres):returncmp(cmpres, 0)def __cmp__(self, other):return self.__norm_cval(cmp(self.__num, other.__num))+\ self.__norm_cval(cmp(self.__string,other.__string)) 执行一些例子: >>> a =NumStr(3, 'foo')>>> b =NumStr(3, 'goo')>>> c =NumStr(2, 'foo')>>> d =NumStr()>>> e =NumStr(string='boo')>>> f =NumStr(1)>>>a [3 :: 'foo']>>>b [3 :: 'goo']>>>c [2 :: 'foo']>>>d [0 ::'']>>>e [0 ::'boo']>>>f [1 :: '']>>> a True>>> b False>>> a ==a True>>> b 2[6 :: 'googoo']>>> a 3[9 :: 'foofoofoo']>>> b +e [3 :: 'gooboo']>>> e +b [3 :: 'boogoo']>>> if d: 'not false'...>>> if e: 'not false'...'not false' >>>cmp(a, b)-1 >>>cmp(a, c)1 >>>cmp(a, a) 0 如果在__str__中使用“%s”,将导致字符串没有引号: return '[%d :: %s]' % (self.__num, self.__string)>>> printa [3 :: foo] 第二个元素是一个字符串,如果用户看到由引号标记的字符串时,会更加直观。要做到这点,使用“repr()”表示法对代码进行转换,把“%s”替换成“%r”。这相当于调用repr()或者使用单反引号来给出字符串的可求值版本--可求值版本的确要有引号: >>> printa [3 :: 'foo'] __norm_cval()不是一个特殊方法。它是一个帮助我们重载__cmp__()的助手函数:唯一的目的就是把cmp()返回的正值转为1,负值转为-1。cmp()基于比较的结果,通常返回任意的正数或负数(或0),但为了我们的目的,需要严格规定返回值为-1,0 和1。 对整数调用cmp()及与 0 比较,结果即是我们所需要的,相当于如下代码片断: def __norm_cval(self, cmpres):if cmpres<0:return -1 elif cmpres>0:return 1 else:return 0 两个相似对象的实际比较是比较数字,比较字符串,然后返回这两个比较结果的和。 本篇文章为转载内容。原文链接:https://blog.csdn.net/weixin_30849865/article/details/112989450。 该文由互联网用户投稿提供,文中观点代表作者本人意见,并不代表本站的立场。 作为信息平台,本站仅提供文章转载服务,并不拥有其所有权,也不对文章内容的真实性、准确性和合法性承担责任。 如发现本文存在侵权、违法、违规或事实不符的情况,请及时联系我们,我们将第一时间进行核实并删除相应内容。
2023-04-19 14:30:42
132
转载
站内搜索
用于搜索本网站内部文章,支持栏目切换。
知识学习
实践的时候请根据实际情况谨慎操作。
随机学习一条linux命令:
groups user
- 显示指定用户的所属组。
推荐内容
推荐本栏目内的其它文章,看看还有哪些文章让你感兴趣。
2023-04-28
2023-08-09
2023-06-18
2023-04-14
2023-02-18
2023-04-17
2024-01-11
2023-10-03
2023-09-09
2023-06-13
2023-08-07
2023-03-11
历史内容
快速导航到对应月份的历史文章列表。
随便看看
拉到页底了吧,随便看看还有哪些文章你可能感兴趣。
时光飞逝
"流光容易把人抛,红了樱桃,绿了芭蕉。"