CN103678921A - 一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法 - Google Patents
一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法 Download PDFInfo
- Publication number
- CN103678921A CN103678921A CN201310697949.7A CN201310697949A CN103678921A CN 103678921 A CN103678921 A CN 103678921A CN 201310697949 A CN201310697949 A CN 201310697949A CN 103678921 A CN103678921 A CN 103678921A
- Authority
- CN
- China
- Prior art keywords
- texture
- silicon steel
- regression
- component
- orientation silicon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Soft Magnetic Materials (AREA)
Abstract
本发明涉及一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,包括下列步骤:测定无取向硅钢中有益、有害织构组分含量;所有数据的标准化处理;不同织构含量数据的降维处理;计算特征值,确定主成分及其表达式;作回归分析并对回归方程作显著性检验;利用标准差标准化法的逆运算将回归方程转换成不同织构组分含量与无取向硅钢磁性能间的多元线性关系。本发明能够有效地对多变量问题进行分析,使原有变量所代表的信息更集中、更典型的体现出来,消除了变量间相关性带来的影响,从定量的角度揭示不同织构组分含量影响无取向硅钢磁性能的规律,为实际生产提供指导,以实现更高磁感、低铁损电工钢产品生产工艺的优化和改进。
Description
技术领域
本发明涉及无取向硅钢性能的控制技术领域,尤其涉及一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法。
背景技术
近年来,随着人们对节能、减耗、环保的日益重视,各类电机、发电机、压缩机等产品正朝高效率化、高精度化、小型化方向发展,故对无取向硅钢的磁性能的要求越来越高,极力追求更高磁感、低铁损的优异磁性能。有研究表明,影响无取向硅钢磁感的主要因素是织构,分析不同织构组分含量对无取向硅钢磁性能的影响,并找出织构对无取向硅钢磁性能的影响规律,能够为实际生产提供指导,以实现更高磁感、低铁损电工钢产品生产工艺的优化和改进。
目前,国内外关于织构对无取向硅钢磁性能的影响已有大量的研究,但只是定性的从机理、织构含量统计等方面初步考查了不同织构组分与磁性能的关系,未能从定量的角度描述无取向硅钢磁性能与不同织构组分间的关系,且国内外关于建立数学模型来分析不同织构组分对无取向硅钢磁性能的影响的研究并不多见。
发明内容
本发明旨在提供一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,运用主成分回归分析法,结合无取向硅钢的磁性能及试验测定的不同织构组分含量的数据进行分析,建立不同织构组分含量与无取向硅钢磁性能间的多元线性关系,并对其进行显著性检验,从而分析织构对无取向硅钢磁性能影响的规律,为实际生产更高磁感、低铁损电工钢产品提供方向性的指导。
为此,本发明采用的技术方案是:
一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其包括如下步骤:
(a)测定无取向硅钢中的不同织构含量,得到有益、有害织构组分含量;
(b)对无取向硅钢的磁性能(铁损P1.5/50、磁感B50)及测定的不同的织构含量数据采用标准差标准化法进行标准化处理;
(c)标准化后的不同织构含量数据的降维处理;
(d)计算特征方程中所有特征值,根据特征值确定主成分的数量,再计算主成分的成分矩阵及单位正交特征向量,得到主成分表达式;
(e)分别以标准化后的铁损、磁感值作为因变量,以提取的主成分为自变量作回归分析;
(f)对得到的回归方程作显著性检验;
(g)利用标准差标准化法的逆变换运算将线性回归方程转换成关于原n个相关变量的多元线性方程;
(h)由上述方法能够从定量的角度分析不同织构组分含量对无取向硅钢磁性能的影响。
进一步地,所述步骤(a)中试验试样的观察面分为轧面和纵截面,采用EBSD技术,在100~200倍视场下,扫描步长选为2.5~5μm,获取试验数据。
进一步地,所述步骤(a)中利用Channel5取向分析软件测定不同织构组分含量,将试验数据导入Project Manager软件分析,选择正交坐标系,观察面为轧面样品的旋转角度参数为0°,0°,0°,观察面为纵截面样品的旋转角度为参数0°,90°,0°,将数据用Tango程序打开,选定需标注出的织构类型,得到该样品的欧拉角衬度图及取向分布图等,并得定量统计计算出各织构的含量,其中有益织构组分主要统计{100}面织构、{110}<001>Goss织构,有害织构组分主要统计{111}<110>织构、{111}<112>织构。
进一步地,所述步骤(b)中数据的标准差标准化法如下:
进一步地,所述步骤(c)中的降维处理具体操作是:将n个相关变量组合成n个独立变量,变换前后保持变量的方差和不变,选择前m(m<n)个独立变量的作用代替原n个相关变量的作用。
进一步地,所述步骤(d)中提取的主成分所满足的条件是:成分特征值大于1。
进一步地,所述步骤(d)中计算单位正交特征向量公式如下:
其中,i表示主成分数量,p表示自变量个数,li表示主成分向量,λi表示第i个主成分的特征值。得到主成分表达式Zi=ti T·[x1 *,x2 *…,xp *]T。
进一步地,所述步骤(f)中多元回归方程线性关系的显著性检验利用统计量F来实现,涉及的参数的计算方法统计于方差分析表中,如表1所示,表中数据除N外均指标准化后的数据。
表1方差分析表
其中,S回表示因变量对各自变量作线性回归的回归平方和,S剩表示残差平方和,S总表示总离差平方和,为由非标准化系数构成回归方程计算出的回归值,即yi的回归拟合值,为所有因变量yi的平均值,N为统计的样品数。
令Fsig(p,N-p-1)=F,即当sig小于0.05时,差异显著,所检验的回归方程呈线性是可靠的,反之,当sig大于0.05时,差异不显著,则所检验的回归方程的线性关系是不可靠的。
本发明的技术优点及效果如下:
1.主成分回归分析法能够有效地对多变量问题进行分析,使原有变量所代表的信息更集中、更典型的体现出来,各变量均能以相同的权重参加运算,消除了变量间相关性带来的影响,简化了回归方程的结构;对大量的原始数据进行主成分回归分析,能够从错综复杂的影响因素中获取主要的信息来定量分析变量间的多元线性关系。
2.利用本发明的方法分析不同织构组分含量对无取向硅钢磁性能的影响时,能够可靠的得到它们之间的多元线性关系,通过比较多元线性回归方程中各自变量的系数,可以从定量的角度揭示不同织构组分含量影响无取向硅钢磁性能的规律。
附图说明
图1为本发明基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法流程示意图。
图2为本发明实施例的不同织构组分含量对无取向硅钢磁感影响程度的柱状图。
具体实施方式
下面结合实施例和附图对本发明作进一步说明。
实施例采用某钢厂提供的经连铸、热轧(2.6mm厚)、冷轧(0.5mm厚)、连退及表面涂层后的无取向硅钢成品检测试样,选取10组磁性能各不相同的试样进行研究,各组试样的磁性能见表2。
表2实施例试样的磁性能
利用ZEISS ULTRA55场发射扫描电镜的EBSD***及Channel5取向分析软件测定试样不同织构含量,试样的观察面分为轧面和纵截面,在100~200倍,优选100倍视场下,扫描步长选为2.5~5μm,优选2.5μm,统计计算出各织构的含量,其中有益织构组分主要统计{100}面织构、{110}<001>Goss织构,有害织构组分主要统计{111}<110>织构、{111}<112>织构,不同织构组分含量的统计数据如表3所示。
表3不同织构组分含量统计数据
本发明借助IBM SPSS Statistics19.0软件对不同织构组分含量统计数据进行主成分回归分析。
首先采用标准差标准化法对无取向硅钢的磁性能(铁损P1.5/50、磁感B50)及统计的不同织构组分含量进行标准化处理,数据分析中,不同织构组分含量对应自变量编号如表4,标准化后的相关数据如表5所示,量纲为1。
表4不同织构组分含量对应自变量编号
表5标准化后的铁损、磁感及不同织构含量统计数据
对标准化后的不同织构含量数据进行降维处理后,计算得到特征根和方差贡献率如表6所示。
表6各成分的特征根和方差贡献率
第一成分的特征值大于1,且方差累积贡献率已达到76.752%,基本上包含了以上自变量x1、x2…x4的所有信息,为了以尽可能少的指标反映尽量多的信息,故选取第一个成分作为主成分,其成分矩阵如表7所示。
表7主成分的成分矩阵
计算出主成分的单位正交特征向量为:t1=[0.55,-0.54,0.46,-0.45]T
相应的主成分表达式为:
Z1=0.55x1 *-0.54x2 *+0.46x3 *-0.45x4 * (3)
将标准化后的铁损YP *、磁感YB *与主成分Z1作多元线性回归分析,表8为分析过程中得到的方差分析表,计算出的回归方程为:
YP *=-0.466Z1 (4)
YB *=0.786Z1 (5)
表8多元线性回归分析的方差分析表
表8中,在分析不同织构组分含量与铁损的关系时,sig值为0.174,大于0.05,所以方程4中Z1为不显著变量,由于Z1为唯一自变量,此时不宜作自变量的显著性检验,只能说明不同织构组分含量不是线性的影响无取向硅钢的铁损。
分析不同织构组分含量与磁感的关系时,由表8可知sig值为0.007,小于0.05,说明方程5的线性关系可靠,利用公式1进行逆运算,可得到磁感YB关于不同织构组分含量的多元线性回归方程,如式6所示。
YB=1.74146+0.00275x1-0.00102x2+0.00224x3-0.00049x4 (6)
由式6及附图2可知,轧面、纵截面上有益织构组分含量的增加能改善无取向硅钢的磁感,而有害织构组分则起到负面作用。在优化磁感时,增加有益织构组分含量比降低有害织构组分含量效果要好。
Claims (8)
1.一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其特征在于,包括如下步骤:
1.1测定无取向硅钢中的不同织构含量,得到有益、有害织构组分含量;
1.2对无取向硅钢的磁性能(铁损P1.5/50和磁感B50)及测定的不同的织构含量数据采用标准差标准化法进行标准化处理;
1.3标准化后的不同织构含量数据的降维处理;
1.4计算特征方程中所有特征值,根据特征值确定主成分的数量,再计算主成分的成分矩阵及单位正交特征向量,得到主成分表达式;
1.5分别以标准化后的铁损、磁感值作为因变量,以提取的主成分为自变量作回归分析;
1.6对得到的回归方程作显著性检验;
1.7利用标准差标准化法的逆变换运算将线性回归方程转换成关于原n个相关变量的多元线性方程;
1.8由上述方法能够从定量的角度分析不同织构组分含量对无取向硅钢磁性能的影响。
2.如权利要求1所述的基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其特征在于,所述步骤1.1中试验试样的观察面分为轧面和纵截面,采用EBSD技术,在100~200倍视场下,扫描步长选为2.5~5μm,获取试验数据。
3.如权利要求1所述的基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其特征在于,所述步骤1.1中利用Channel5取向分析软件测定不同织构组分含量,将试验数据导入ProjectManager软件分析,选择正交坐标系,观察面为轧面样品的旋转角度参数为0°,0°,0°,观察面为纵截面样品的旋转角度为参数0°,90°,0°,将数据用Tango程序打开,选定需标注出的织构类型,得到该样品的欧拉角衬度图及取向分布图等,并得定量统计计算出各织构的含量,其中有益织构组分主要统计{100}面织构、{110}<001>Goss织构,有害织构组分主要统计{111}<110>织构、{111}<112>织构。
5.如权利要求4所述的基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其特征在于,所述步骤1.3中的降维处理具体操作是:将n个相关变量组合成n个独立变量,变换前后保持变量的方差和不变,选择前m个独立变量的作用代替原n个相关变量的作用,其中m<n。
6.如权利要求5所述的基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其特征在于,所述步骤1.4中提取的主成分所满足的条件是:成分特征值大于1。
7.如权利要求6所述的基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其特征在于,所述步骤1.4中计算单位正交特征向量公式如下:
其中,i表示主成分数量,p表示自变量个数,li表示主成分向量,λi表示第i个主成分的特征值,得到主成分表达式Z i =t i T·[X1 *,X2 *...,X p *]T。
8.如权利要求7所述的基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法,其特征在于,所述步骤1.6中多元回归方程线性关系的显著性检验利用统计量F来实现,涉及的参数的计算方法统计于方差分析表中,如表1所示,表中数据除N外均指标准化后的数据;
表1方差分析表
其中,S回表示因变量对各自变量作线性回归的回归平方和,S剩表示残差平方和,S总表示总离差平方和,为由非标准化系数构成回归方程计算出的回归值,即yi的回归拟合值,为所有因变量yi的平均值,N为统计的样品数;
令Fsig(p,N-p-1)=F,即当sig小于0.05时,差异显著,所检验的回归方程呈线性是可靠的,反之,当sig大于0.05时,差异不显著,则所检验的回归方程的线性关系是不可靠的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310697949.7A CN103678921B (zh) | 2013-12-18 | 2013-12-18 | 一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310697949.7A CN103678921B (zh) | 2013-12-18 | 2013-12-18 | 一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103678921A true CN103678921A (zh) | 2014-03-26 |
CN103678921B CN103678921B (zh) | 2016-06-15 |
Family
ID=50316451
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310697949.7A Active CN103678921B (zh) | 2013-12-18 | 2013-12-18 | 一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103678921B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825288A (zh) * | 2015-12-07 | 2016-08-03 | 北京师范大学 | 一种消除复杂***中回归数据共线性问题的优化分析方法 |
CN105956274A (zh) * | 2016-05-04 | 2016-09-21 | 武汉钢铁股份有限公司 | 通过织构指数评价无取向硅钢磁性能的方法 |
CN109754164A (zh) * | 2018-12-18 | 2019-05-14 | 内蒙古鄂尔多斯资源股份有限公司东昊厂 | 一种羊绒纺织生产计算投入量的动态制成率模型 |
CN112461735A (zh) * | 2020-11-02 | 2021-03-09 | 航天科工防御技术研究试验中心 | 一种加速腐蚀的试验方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090292473A1 (en) * | 2008-05-23 | 2009-11-26 | Baker Hughes Incorporated | Real-Time NMR Distribution While Drilling |
CN101707919A (zh) * | 2007-06-15 | 2010-05-12 | 英国石油化学品有限公司 | 在线分析蒸气相工艺流的方法 |
-
2013
- 2013-12-18 CN CN201310697949.7A patent/CN103678921B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101707919A (zh) * | 2007-06-15 | 2010-05-12 | 英国石油化学品有限公司 | 在线分析蒸气相工艺流的方法 |
US20090292473A1 (en) * | 2008-05-23 | 2009-11-26 | Baker Hughes Incorporated | Real-Time NMR Distribution While Drilling |
Non-Patent Citations (4)
Title |
---|
WANG J.ET AL.: "Effect of heating rate on micostructure evolution and magnetic properties of cold rolled non-oriented electrical steel", 《JOURNAL OF ICON AND STTEL RESEARCH》 * |
党耀国等人编著: "《应用多元统计分析》", 31 May 2012 * |
桂福生等人: "成分与工艺对冷轧无取向低碳和低硅电工钢织构与磁性能的影响", 《钢铁研究总院学报》 * |
许令峰等人: "晶体取向对冷轧无取向硅钢磁时效的影响", 《功能材料》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825288A (zh) * | 2015-12-07 | 2016-08-03 | 北京师范大学 | 一种消除复杂***中回归数据共线性问题的优化分析方法 |
CN105825288B (zh) * | 2015-12-07 | 2019-12-10 | 北京师范大学 | 一种消除复杂***中回归数据共线性问题的优化分析方法 |
CN105956274A (zh) * | 2016-05-04 | 2016-09-21 | 武汉钢铁股份有限公司 | 通过织构指数评价无取向硅钢磁性能的方法 |
CN109754164A (zh) * | 2018-12-18 | 2019-05-14 | 内蒙古鄂尔多斯资源股份有限公司东昊厂 | 一种羊绒纺织生产计算投入量的动态制成率模型 |
CN112461735A (zh) * | 2020-11-02 | 2021-03-09 | 航天科工防御技术研究试验中心 | 一种加速腐蚀的试验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103678921B (zh) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104820873A (zh) | 一种基于金属定量构效关系的淡水急性基准预测方法 | |
CN103678921A (zh) | 一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法 | |
CN103823974A (zh) | 无取向硅钢磁性能影响因素的主成分回归分析法 | |
CN104316660A (zh) | 快速检测豆瓣酱风味、等级及产品溯源的方法 | |
Yang et al. | Estimating the regional eco-efficiency in China based on bootstrapping by-production technologies | |
CN103632013A (zh) | 基于主成分回归分析的晶粒尺寸影响无取向硅钢磁性能的分析方法 | |
CN105021564A (zh) | 基于近红外光谱分析技术测定烟叶中麦角甾醇含量的方法 | |
CN106844208A (zh) | 用于反应堆物理计算程序适用性验证的方法 | |
Shiu et al. | Technical change and total factor productivity growth for Chinese provinces: A panel data analysis | |
CN105808910B (zh) | 锆铌压力管在candu反应堆中的辐照变形计算方法 | |
CN105243608A (zh) | 基于输变电工程造价设计要素的线路工程造价估算方法 | |
CN102994736B (zh) | 罩式炉退火模型修正方法 | |
Xie et al. | Measuring total‐factor energy environment efficiency, energy‐saving and carbon emission‐reduction potential in China's food industry: Based on a meta‐frontier slacks‐based measure model | |
CN103823975A (zh) | 织构组分对无取向硅钢磁感影响的主成分回归分析法 | |
CN106649202B (zh) | 基于多样性变量加权plsr模型的工业过程软测量方法 | |
CN107871215A (zh) | 一种标准对经济发展影响评估模型及验证方法 | |
Li et al. | Big data oriented macro-quality index based on customer satisfaction index and PLS-SEM for manufacturing industry | |
Ding et al. | Controllable points setting method for plan view pattern control in plate rolling process | |
CN105021532B (zh) | 一种快速检测x70金相组织及组织含量的方法 | |
CN103678922B (zh) | 基于主成分回归分析的夹杂物影响无取向硅钢磁性能的分析方法 | |
CN109033125A (zh) | 时序数据支配集信息提取方法 | |
Yafei et al. | Temperature and carbon content dependence of electrical resistivity of carbon steel | |
CN113791131A (zh) | 一种构建锆石二次离子质谱年龄三维分布图像的方法 | |
Parteka | Employment and export specialization patterns versus gdp per capita performance-unifying approach | |
CN105700475A (zh) | 一种实现大范围环境温度的机床稳健性热误差补偿的数据处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |