CN2743834Y - 基于木材结构参数确定木材物理、力学特性的分析仪器 - Google Patents

基于木材结构参数确定木材物理、力学特性的分析仪器 Download PDF

Info

Publication number
CN2743834Y
CN2743834Y CN 200420096149 CN200420096149U CN2743834Y CN 2743834 Y CN2743834 Y CN 2743834Y CN 200420096149 CN200420096149 CN 200420096149 CN 200420096149 U CN200420096149 U CN 200420096149U CN 2743834 Y CN2743834 Y CN 2743834Y
Authority
CN
China
Prior art keywords
wood
module
sigma
test
tube baby
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.)
Expired - Lifetime
Application number
CN 200420096149
Other languages
English (en)
Inventor
江泽慧
周玉成
姜笑梅
井元传
吕建雄
赵亮
费本华
秦特夫
黄洛华
任海青
赵有科
殷亚方
刘君良
余雁
Original Assignee
周玉成
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 周玉成 filed Critical 周玉成
Priority to CN 200420096149 priority Critical patent/CN2743834Y/zh
Application granted granted Critical
Publication of CN2743834Y publication Critical patent/CN2743834Y/zh
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本实用新型提出一种基于木材结构参数确定木材物理、力学特性的分析仪器。该仪器包括:输入键盘,计算机主机,显示器或/和打印输出装置,所述计算机主机内包括:1)输入单元,包括:用于接收木材结构参数的输入模块,参数辨识、校对及控制器,加法器;2)接口引擎单元,包括:开启Matlab引擎的模块,数据缓冲区,实现木材结构参数数据格式转换的模块,神经元网络模型文件的调用模块;3)神经元网络非线性逼进单元,包括:用于确定木材结构参数影响木材物理力学特性的权值函数的计算模块,用于通过自学习神经元网络层权值修正模块并产生权值的线性函数的模块;4)输出单元,包括:材性量化指标的输出模块和相关输出曲线的输出模块。

Description

基于木材结构参数确定木材物理、力学特性的分析仪器
技术领域
本实用新型涉及一种由木材结构参数确定木材物理、力学特性的分析仪器。这些分析可以为木材科学领域的木材性质研究、探讨木材性质形成机理、树木优质种质资源培育、树木转基因工程、定向培育材质改良的树木新品种等方面提供科学依据及量化指标。
背景技术
从森林的培育到木材的综合利用,都与树木及木材的内部结构息息相关。为此,国内外专家学者都在木材内部结构与物理及力学关系上问题上开展了大量的研究。克恩(Cown)在(Cown et al.1991)给出了木材的螺旋纹理与木材的力学强度的关系模型,通过该模型可知螺旋纹理对木材的干燥程度和力学强度的关系,并给出根据螺旋纹理的分布状态、方向和角度确定其木材加工方法,已部分地缓解强度降低、开裂等问题。哈里斯(Harris)研究了木材纹理与木材种类的关系,并得出大部分木材呈螺旋纹理仅有少部分木材呈直纹理状态的结论(Harris 1973,1989)。狄普(Draper)、史密斯(Smith)用线性或非线性的方法对木材结构的单一特征与密度、强度等建立了回归模型。
就目前国内外对木材结构与力学、强度的关系研究存在着如下问题:
●现有揭示木材结构与力学强度的关系不能反映其全部的特征。因为就其树木与木材本身就是一类复杂的生命体,仅用其内部的某一结构特征是不能全面而准确的反映其规律的。
●现有的回归方法所得到的精度一般都在70~80%左右,这只能说成是一个大致的趋势,而无法准确地把握真正的内在性质与力学间的关系;
●现有的成果很难推广应用到森林种质资源培育、转基因工程、树木定向培育、木材材质改良以及提供木材新品种的研究和应用上来。
发明内容
本实用新型要解决的技术问题是针对国内外对树木内部结构与木材物理、力学关系的研究开发中存在的问题和不足,提出一种基于木材结构参数确定木材物理、力学特性的分析仪器。该仪器提供了对木材性质形成机理的成因分析,描述了木材这一类复杂生命体内部结构与其物理、力学间的相互关系,并保证精度在95%以上。其目的是为森林种质资源培育、转基因工程、树木定向培育、木材性质研究、木材性质及形成机理定向培育木材材质改良、木材新品种培育等领域提供科学的依据及方法,使我国能在短期内,在上述领域内位居世界前列。
本实用新型提出一种基于木材结构参数确定木材物理、力学特性的分析仪器,它包括:输入键盘,计算机主机,显示器或/和打印输出装置,所述计算机主机内包括并依次连接:
1)输入单元,其中依次包括:用于接收木材结构参数的输入模块,参数辨识、校对及控制器,加法器;
2)接口引擎单元,其中包括:用于开启Matlab引擎的模块,以及并列联接的数据缓冲区、用于实现木材结构参数数据格式转换的模块、用于神经元网络模型文件的调用模块;
3)神经元网络非线性逼进单元,其中依次包括:用于确定木材结构参数影响木材物理力学特性的权值的权值函数的计算模块,用于通过自学习神经元网络层权值修正并产生权值的线性函数的模块;
4)输出单元,其中包括:材性量化指标的输出模块和相关输出曲线的输出模块。
本实用新型的基于木材结构参数确定木材物理、力学特性的分析仪器与现有技术相比较,具有以下优点:
1、本实用新型中的木材解剖微观参数包括:微纤丝角、胞壁率、射线比量、管胞比量、早材管胞长度、晚材管胞长度、早材管胞宽度、晚材管胞宽度、早材管胞T-壁厚度、晚材管胞T-壁厚度、早材管胞R-壁厚度、晚材管胞R-壁厚度、早材管胞T-直径、晚材管胞T-直径、早材管胞壁腔比、晚材管胞壁腔比、早材管胞腔径比、晚材管胞腔径比等18个参数,它们都对该木材的密度、径向干缩、轴向干缩、弦向干缩、顺纹抗压强度、抗弯弹性模量、抗弯强度等7个重要材性指标产生不同的影响。这种不同程度的影响,在本实用新型中体现在权值函数中,也就是说,本实用新型全面的把握了上述木材结构参数与木材物理力学特性的关系。而现有的对这种关系的把握仅是上述的一个或几个参数(如微纤丝角与木材物理力学特性的关系)对木材物理力学之间的关系。因此,本实用新型可以更完全、更真实地表现出木材结构参数与其物理力学特性之间的关系。显示直观,一目了然。
2、现有的对这种复杂非线性关系的描述是基于线性或非线性回归的方法,而这种方法本身就存在着不精确或者说先天不足。本实用新型所采用的神经元网络方法,可以任意地逼进非线性,即只要神经元层数足够多,便可以达到任意理想的结果。
3、现有的算法是通过观察输入与输出之间点的分布情况采取曲线拟合方法,其准确度只能达到70%~80%之间。而本实用新型通过木材每一个结构参数对木材物理、力学特性影响的权值及所求出的权值函数,更加真实地、切合实际地描述了木材结构参数与木材物理力学特性之间关系,其准确度可达95%以上。
附图说明
图1是基于木材结构参数确定木材物理、力学特性分析仪器的结构示意图;
图2是本实用新型分析仪器的工作流程图;
图3是输入模块的示意图;
图4是接口引擎的工作流程图;
图5是标准化处理器的示意图;
图6是标准化处理器的运算模块的示意图;
图7是神经元网络的运算模块的示意图;
图8是神经元网络性能评估的线性回归图。
具体实施方式
在说明本实用新型的实施例之前,首先说明本实用新型的原理。
由于木材本身就是一个复杂的生命体或生物材料,其内部的结构与结构之间、这些结构与其自身的物理力学之间存在着相当复杂的奇异非线性关系,一般说来不存在或很困难找到显式的函数关系,即使免强找到了函数关系,那也会是一类非常复杂的非线性偏微分方程,处理起来相当难以把握。事实上,非线性奇异系统都不能用子空间来描述(它们是属于一些低维子流形),直接讨论低维子流形是比较困难的。因此,本实用新型利用了神经元网络能够逼进任意的非线性这一特性。即,通过输入模块接收木材结构参数,经过参数辨识及校对经加法器输出给接口引擎。在接口引擎中首先开辟一个数据缓冲区,在此区内实现木材结构参数数据的格式转换、神经元网络模型文件的调用。经过初步处理的数据传递给神经网络的第一层,在这一层须要标准化处理,再经神经元网络初步的确定后计算出欧氏距离,最后通过自学习模型的修正后产生标准的权值函数,这个函数已将木材结构参数对木材物理力学的影响的权值确定。在神经元网络的第二层,将收到第一层权值函数通过自学习神经网络层权值模块修正后进行正规化处理产生权值的线性函数,最后达到木材物理力学特性的理想输出,因此,只要神经元层数足够多,就可以达到任意精度的结果。
下面结合附图和实施例,进一步说明本实用新型。
本实用新型工作流程如图2所示。分析仪器的结构如图1所示,包括:输入键盘,计算机主机,显示器或打印输出装置,或者显示器和打印输出装置,所述计算机主机还包括以下四大部分并依次联接:
I.输入单元1,其中包括:木材结构参数的输入模块5,参数辨识、校对及控制器6,加法器7;本实用新型在木材结构参数的输入模块5中提供以下三种输入模块,可以选择一种或其组合:
a.基于整体模型的宏观分析的微观参数输入模块8,在该模块中共有18个木材解剖的微观参数,包括微纤丝角、胞壁率、射线比量、管胞比量、早材管胞长度、晚材管胞长度、早材管胞宽度、晚材管胞宽度、早材管胞T-壁厚度、晚材管胞T-壁厚度、早材管胞R-壁厚度、晚材管胞R-壁厚度、早材管胞T-直径、晚材管胞T-直径、早材管胞壁腔比、晚材管胞壁腔比、早材管胞腔径比、晚材管胞腔径比,根据以上18个输入参数即可确定该木材的7个重要材性指标,其中包括基本密度、径向干缩、轴向干缩、弦向干缩、顺纹抗压强度、抗弯弹性模量、抗弯强度。
b.分年龄段木材材性指标分析的参数输入模块9,在宏观分析的基础上,反应木材不同年龄段的性质变异和边材心材的差异,给出分年龄段木材材性指标分析的分模型,在输入参数时给出其木材年龄即可得出精度更高的材性指标预测结果。本实用新型能分析和预测1-15年的木材。
c.重要解剖参数的输入模块10,该模块保证基本功能前提下简化操作,减少对大量解剖参数的依赖,提供了不完全查询,即可从18个微观解剖数据中自行选择输入参数或定义输入参数个数,得出木材的7个材性指标,特别指出的是,此项不完全分析能保证其最大精度,在不完全查询的参数组合中采用了逐步回归分析法将参数归类组合,采用主成分分析法可以降低输入参数的维数。
图3所示为输入模块结构图,无论选择上述哪种分析方式,终端与界面接口都交互有数据信息及与数据相关的控制信息,即参数定义域控制、参数步长控制、数据格式控制。其中参数定义域控制主要用于校正参数输入,在误操作时给出告警,各定义域的初始设置由实验统计数据确定,仪器工作前可专门定制其所需要的参数区间。参数步长控制用于更改参数微调按钮的步长,便于使用。数据格式控制是在输入数据时选择几位有效数字和何种科学计数方法。控制信息和数据信息在输入模块内部的通信是透明的,可以自由选择和分别输入控制信息和数据信息,二者永远保持同步。
II.接口引擎单元2,其中包括:用于开启Matlab引擎的模块11,以及并列联接的数据缓冲区13、用于实现木材结构参数数据格式转换的模块14、用于神经元网络模型文件的调用模块12。本实用新型采用专门的数学运算软件Matlab作为其后台运算和分析的工具,在该单元要建立起前端和后台Matlab的通信,即Matlab接口引擎。具体说来,就是将输入的参数信息通过Matlab与界面开发程序的接口送入到Matlab工作空间,参与计算,运算结束后又能把运算结果以合适的数据格式返回,同时,提供一系列神经元网络性能和训练结果分析的图表。其工作流程如图4所示。
III.神经元网络非线性逼进单元3,其中至少包括如下两个单元:
a.神经元网络的第一层单元15,它包括:标准化处理器17,神经元网络隐含层权矩阵模块18,欧氏距离计算模块19,自学习模型的修正模块20,点乘模块21,高斯函数发生器22;所述标准化处理器17和神经元网络隐含层权矩阵模块18分别与欧氏距离计算模块19联接,然后欧氏距离计算模块19和自学习模型的修正模块20分别与点乘模块21联接,最后点乘模块21与高斯函数发生器22联接。
b.神经元网络的第二层单元16,它又包括并依次联接:神经元网络线性层权矩阵模块23,正规化处理器24,产生权值的线性函数发生器25。
在进入真正的运算网络之前,输入样本要进入标准化处理器,进行事先、事后的标准化处理,处理方式如图5所示,对于输入样本数据0-104的值域区间进行标准化处理。通过标准化处理后将输入向量和目标输出向量量化为零均值和偏差为1的标准向量。下面是通过零均值和偏差型函数的实现过程。
[pn,meanp,stdp,tn,meant,stdt]=prestd(p,t)
参数意义:p       网络输入向量
          t       目标输出向量
          pn      量化后的输入向量
          meanp   输入向量的均值
          stdp    输入向量的偏差
          tn      量化后的目标输出向量
          meant   目标输出的均值
          stdt    目标输出的偏差
具体的标准化过程如图6所示。下面给出标准化处理器的运算环节:
mean ( P ′ ) = mean ( p 11 p 21 Λ p R 1 p 12 p 22 Λ p R 2 M M M p 1 Q p 2 Q Λ p RQ )
= Σ i = 1 i = Q p 1 i Q Σ i = 1 i = Q p 2 i Q Λ Σ i = 1 i = Q p Ri Q 1 × R
std ( P ′ ) = std ( p 11 p 21 Λ p R 1 p 12 p 22 Λ p R 2 M M M p 1 Q p 2 Q Λ p RQ )
= ( 1 Q - 1 Σ i = 1 Q ( p 1 i - 1 Q Σ i = 1 Q p 1 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2 1 × R
Pn=(P-meanp×oneQ)·/(stdp×oneQ)
= ( p 11 p 12 Λ p 1 Q p 21 p 22 Λ p 2 Q M M M p R 1 p R 2 Λ p RQ - Σ i = 1 i = Q p 1 i Q Σ i = 1 i = Q p 2 i Q M Σ i = 1 i = Q p Ri Q R × 1 × 1 1 Λ 1 1 × Q ) · / ( ( 1 Q - 1 Σ i = 1 Q ( p 1 i - 1 Q Σ i = 1 Q p 1 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 M ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2 R × 1 1 1 Λ 1 1 × Q )
= p 11 - Σ i = 1 i = Q p 1 i Q p 12 - Σ i = 1 i = Q p 1 i Q Λ p 1 Q - Σ i = 1 i = Q p 1 i Q p 21 - Σ i = 1 i = Q p 2 i Q p 22 - Σ i = 1 i = Q p 2 i Q Λ p 2 Q - Σ i = 1 i = Q p 2 i Q M M M p R 1 - Σ i = 1 i = Q p Ri Q p R 2 - Σ i = 1 i = Q p Ri Q Λ p RQ - Σ i = 1 i = Q p Ri Q · / ( 1 Q - 1 Σ i = 1 Q ( p 12 - 1 Q Σ i = 1 Q p 12 ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( p 12 - 1 Q Σ i = 1 Q p 12 ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( p 12 - 1 Q Σ i = 1 Q p 12 ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 M M M ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2
p 11 - Σ i = 1 i = Q p 1 i Q ( 1 Q - 1 Σ i = 1 Q ( p 1 i - 1 Q Σ i = 1 Q p 1 i ) 2 ) 1 2 p 12 - Σ i = 1 i = Q p 1 i Q ( 1 Q - 1 Σ i = 1 Q ( p 1 i - 1 Q Σ i = 1 Q p 1 i ) 2 ) 1 2 Λ p 1 Q - Σ i = 1 i = Q p 1 i Q ( 1 Q - 1 Σ i = 1 Q ( p 1 i - 1 Q Σ i = 1 Q p 1 i ) 2 ) 1 2 p 21 - Σ i = 1 i = Q p 2 i Q ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 p 22 - Σ i = 1 i = Q p 2 i Q ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 Λ p 2 Q - Σ i = 1 i = Q p 2 i Q ( 1 Q - 1 Σ i = 1 Q ( p 2 i - 1 Q Σ i = 1 Q p 2 i ) 2 ) 1 2 M M M p R 1 - Σ i = 1 i = Q p Ri Q ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2 p R 2 - Σ i = 1 i = Q p Ri Q ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2 Λ p RQ - Σ i = 1 i = Q p Ri Q ( 1 Q - 1 Σ i = 1 Q ( p Ri - 1 Q Σ i = 1 Q p Ri ) 2 ) 1 2
mean ( T ′ ) = mean ( t 11 t 21 Λ t S 1 t 12 t 22 Λ t S 2 M M M t 1 Q t 2 Q Λ t SQ )
= Σ i = 1 i = Q t 1 i Q Σ i = 1 i = Q t 2 i Q Λ Σ i = 1 i = Q t Si Q 1 × S
std ( T ′ ) = std ( t 11 t 21 Λ t S 1 t 21 t 22 Λ t S 2 M M M t 1 Q t 2 Q Λ t SQ )
( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( t Si - 1 Q Σ i = 1 Q t Si ) 2 ) 1 2 1 × S
Tn=(T-meant×oneQ)·/(stdt×oneQ)
= ( p 11 p 12 Λ p 1 Q p 21 p 22 Λ p 2 Q M M M p S 1 p S 2 Λ p SQ - Σ i = 1 i = Q t 1 i Q Σ i = 1 i = Q t 2 i Q M Σ i = 1 i = Q t Si Q S × 1 × 1 1 Λ 1 1 × Q ) · / ( ( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 M ( 1 Q - 1 Σ i = 1 Q ( t Si - 1 Q Σ i = 1 Q t Si ) 2 ) 1 2 S × 1 1 1 Λ 1 1 × Q )
= t 12 - Σ i = 1 i = Q t 1 i Q t 12 - Σ i = 1 i = Q t 1 i Q Λ t 1 Q - Σ i = 1 i = Q t 1 i Q t 2 i - Σ i = 1 i = Q t 2 i Q t 22 - Σ i = 1 i = Q t 2 i Q Λ t 2 Q - Σ i = 1 i = Q t 2 i Q M M M t S 1 - Σ i = 1 i = Q t Si Q t S 2 - Σ i = 1 i = Q t Si Q Λ t SQ - Σ i = 1 i = Q t Si Q · / ( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 M M M ( 1 Q - 1 Σ i = 1 Q ( t Si - 1 Q Σ i = 1 Q t Si ) 2 ) 1 2 ( 1 Q - 1 Σ i = 1 Q ( t Si - 1 Q Σ i = 1 Q t Si ) 2 ) 1 2 Λ ( 1 Q - 1 Σ i = 1 Q ( t Si - 1 Q Σ i = 1 Q t Si ) 2 ) 1 2
t 11 - Σ i = 1 i = Q t 1 i Q ( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 t 12 - Σ i = 1 i = Q t 1 i Q ( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 Λ t 1 Q - Σ i = 1 i = Q t 1 i Q ( 1 Q - 1 Σ i = 1 Q ( t 1 i - 1 Q Σ i = 1 Q t 1 i ) 2 ) 1 2 t 21 - Σ i = 1 i = Q t 2 i Q ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 t 22 - Σ i = 1 i = Q t 2 i Q ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 Λ t 2 Q - Σ i = 1 i = Q t 2 i Q ( 1 Q - 1 Σ i = 1 Q ( t 2 i - 1 Q Σ i = 1 Q t 2 i ) 2 ) 1 2 M M M t R 1 - Σ i = 1 i = Q t Ri Q ( 1 Q - 1 Σ i = 1 Q ( t Ri - 1 Q Σ i = 1 Q t Ri ) 2 ) 1 2 t R 2 - Σ i = 1 i = Q t Ri Q ( 1 Q - 1 Σ i = 1 Q ( t Si - 1 Q Σ i = 1 Q t Si ) 2 ) 1 2 Λ p RQ - Σ i = 1 i = Q p Ri Q ( 1 Q - 1 Σ i = 1 Q ( t Si - 1 Q Σ i = 1 Q t Si ) 2 ) 1 2
经过标准化处理器后,进入正式的神经元网络运算模块,如图7所示:
该神经元网络模块中有Q个节点,R维输入,S维输出。神经元网络输入向量为P,神经元网络隐含层的权值设为P’,该层每个神经元节点的带权输入,即欧几里德距离‖dist‖。‖dist‖是欧几里德距离权函数,权函数把权重加到输入矩阵上以得到带权输入矩阵。对于dist(W,P),W为S×R权矩阵,P为Q维输入列向量矩阵,dist(W,P)返回S×Q维向量距离矩阵。
在该模型中将权矩阵W定义为P′,则
dist ( P Q × R ′ , P R × Q )
= dist ( p 11 p 12 K p 1 R p 21 p 22 K p 2 R M M M p Q 1 p Q 2 K p QR , p 11 p 21 K p Q 1 p 12 p 22 K p Q 2 M M M p 1 R p 2 R K p QR )
= 0 d 12 d 13 Λ d 1 Q d 21 0 d 23 Λ d 2 Q d 31 d 32 0 Λ d 3 Q M M M M d Q 1 d Q 2 d Q 3 Λ 0 Q × Q
式中dij表示矩阵P′第i个行向量与矩阵P第j个列向量间的距离,因此对角线上的元素均为0。进而将dist(P′,P)与b1作点乘,即
dist ( P ′ , P ) · * b 1 = 0 d 12 d 13 Λ d 1 Q d 21 0 d 23 Λ d 2 Q d 31 d 32 0 Λ d 3 Q M M M M d Q 1 d Q 2 d Q 3 Λ 0 · * b 11 b 12 Λ b 1 Q b 21 b 22 Λ b 2 Q M M M b Q 1 b Q 2 Λ b QQ
= 0 b 12 * d 12 b 12 * d 13 Λ b 1 Q * d 1 Q b 21 * d 21 0 b 23 * d 23 Λ b 2 Q * d 2 Q b 31 * d 31 b 32 * d 32 0 Λ b 3 Q * d 3 Q M M M M b Q 1 * d Q 1 b Q 2 * d Q 2 b Q 3 * d Q 3 Λ 0
在网络隐含层每个神经元网络输入是其带权输入与其偏值之积,如上式所示,而每个神经元的输出是网络输入的径向基函数。
利用高斯核函数(Gaussian kernel function)作为基函数的形式,如下式所示:
u j = exp [ - ( X - C j ) T ( X - C j ) 2 δ j 2 ] , j = 1,2 , K , N h
其中,uj是第j个隐层节点的输出,X=(x1,x2,K,xn)T是输入样本,Cj是高斯函数的中心值,δj是标准化常数,Nh是隐层节点数。其隐含层节点中的作用函数(核函数)对输入信号将在局部产生响应,也就是说,当输入信号靠近核函数的中央范围时,隐层节点将产生较大的输出,由此,这种神经元网络具有局部逼近能力,所以径向基函数网络也成为局部感知场网络。由上式可知,节点的输出范围在0和1之间,如果一个神经元的权向量与其输入向量相等(转置),其带权输入将为0,当其网络输入为0,则输出为1,且输入样本愈靠近节点的中心,输出值愈大。
采用高斯基函数,具备如下优点:
1、表示形式简单,即使对于多变量输入也不增加太多的复杂性;
2、径向对称;
3、光滑性好,任意阶导数存在;
4、由于该基函数表示简单且解析性好,因而便于进行理论分析。
经过基函数发生器,神经元节点输出为a1,开始进入网络线性层。在网络线性层中首先要经过一个正规化处理器,然后再进入普通的线性神经元。在正规化处理器中,采用normprod函数来计算网络的输出向量n2。normprod是一个权函数,权函数将权重加到输入矩阵上得到带权矩阵。对于normprod(W,P),W为S×R权矩阵,P为Q维输入列向量矩阵,normprod(W,P)返回S×Q维正规化点积。
在该网络中,将网络线性层的权矩阵设为网络的目标输出TS×Q,即normprod(T,a1)
= normorid ( t 11 t 12 Λ t 1 Q t 21 t 22 Λ t 2 Q M M M t S 1 t S 2 Λ t SQ , a 11 a 12 Λ a 1 Q a 21 a 22 Λ a 2 Q M M M a Q 1 a Q 2 Λ a QQ )
= Σ j = 1 j = Q t 1 j * a j 1 Σ j = 1 j = Q a j 1 Σ j = 1 j = Q t 1 j * a j 2 Σ j = 1 j = Q a j 2 Λ Σ j = 1 j = Q t 1 j * a jQ Σ j = 1 j = Q a jQ Σ j = 1 j = Q t 2 j * a j 1 Σ j = 1 j = Q a j 1 Σ j = 1 j = Q t 2 j * a j 2 Σ j = 1 j = Q a j 2 Λ Σ j = 1 j = Q t 2 j * a jQ Σ j = 1 j = Q a jQ M M M Σ j = 1 j = Q t Sj * a j 1 Σ j = 1 j = Q a j 1 Σ j = 1 j = Q t Sj * a j 2 Σ j = 1 j = Q a j 2 Λ Σ j = 1 j = Q t Sj * a jQ Σ j = 1 j = Q a jQ
本实用新型所采用的神经元网络是一个逼近器,只要隐单元足够多,它就可以逼近任意M元连续函数且对任一未知的非线性函数,总存在一组权值使网络对该函数的逼近效果最好。网络第二层也有与网络输入和目标向量相同的神经元节点数,这里将第二层的权值矩阵设为目标向量矩阵T。
IV.输出模块4,其中包括:材性量化指标的输出模块26和相关输出曲线的输出模块27。
当网络训练结束之后,用sim函数来仿真神经元网络的输出,从而与目标输出进行比较,来检验神经元网络的性能。函数postreg利用了线形回归的方法分析了神经元网络输出和目标输出的关系,即神经元网络输出变化相对于目标输出变化的变化率,从而评估了神经元网络的训练结果。
a=sim(net,p)
[m,b,r]=postreg(a,t)
函数postreg返回了3个值,m和b分别表示最优回归直线的斜率和y轴截距,当m等于1,b等于0的时候,神经元网络输出和目标输出完全相同,此时的神经元网络具有最优的性能。r表示网络输出与目标输出的相关系数,它越接近于1,表示网络输出与目标输出越接近,神经元网络性能越好。函数postreg显示的图形中,横坐标为目标输出,纵坐标为网络输出,“o”表示数据,理想回归直线(神经元网络输出等于目标输出时的直线)由实线表示,最优回归直线由虚线表示。在神经元网络输出模块中,分别给出各个输出材性指标的用来衡量神经元网络性能线性回归图形,输出图形如图8所示,从仿真图中可以看出其精度达到了98.999%。
以下列表是采用本实用新型分析仪器对杉木和72杨的分析结果。
杉木微观解剖结构数据(输入I)
  微纤丝角   12.5   晚材管胞T-壁厚度   8.80   晚材管胞壁腔比   0.5195
  早材管胞长度   2481.39   早材管胞R-壁厚度   6.29   早材管胞腔径比   0.8456
  晚材管胞长度   2121.32   晚材管胞R-壁厚度   11.32   晚材管胞腔径比   0.6562
  早材管胞宽度   45.63   早材管胞T-直径   39.93   射线比量   7.11
  晚材管胞宽度   30.46   晚材管胞T-直径   30.18   管胞比量   92.89
  早材管胞T-壁厚度   5.48   早材管胞壁腔比   0.1869   胞壁率   41.79
72杨微观解剖结构数据(输入II)
  微纤丝角   29.22   腔径比   0.70764
  纤维长度   1191   纤维比量   58.72
  纤维宽度   24.19   射线比量   13.86
  导管长度   584.19   导管比量   27.42
  纤维壁厚   5.13   胞壁率   44.84
  纤维直径   17.55   导管直径   71.5
  壁腔比   0.4142   导管个数   51.9
杉木、72杨化学组成表(输入)
  杉木化学组成   72杨化学组成
  综纤维素   60.79   综纤维素   76.51
  木质素   33.62   木质素   19.29
  α-纤维素   41.55   α-纤维素   39.35
  半-纤维素   19.24   半-纤维素   37.16
  酸溶木质素   0.32   酸溶木质素   3.94
杉木微观解剖所确定的物理力学特性输出表
  目标输出   实际输出   相对误差(%)
  轴向干缩   0.2601   0.2508   3.57
  弦向干缩   7.0225   7.0683   0.65
  径向干缩   2.6540   2.6893   1.33
  基本密度   0.2491   0.2544   2.13
  抗弯强度   52.0   52.4   0.77
  顺纹抗压强度   30.7   30.784   0.27
  抗弯弹性模量   5673   5828   2.73
72杨微观解剖所确定的物理力学特性输出表
  目标输出   实际输出   相对误差(%)
  轴向干缩   0.5624   0.5528   1.71
  弦向干缩   7.0111   7.2179   2.95
  径向干缩   0.7813   1.0936   39.97
  基本密度   0.2830   0.2933   3.64
  抗弯强度   56.2   56.881   1.21
  顺纹抗压强度   25.2   25.448   0.98
  抗弯弹性模量   5260   5285.8   0.49
杉木化学组成所确定的物理力学特性输出表
  目标输出   实际输出   相对误差(%)
  轴向干缩   0.2601   0.24974   4.16
  弦向干缩   7.0225   6.896   1.8
  径向干缩   2.6540   2.569   3.2
  基本密度   0.2491   0.27146   8.98
  抗弯强度   52.0   51.288   1.37
  顺纹抗压强度   30.7   30.451   0.81
  抗弯弹性模量   5673   5577.5   1.68
72杨化学组成所确定的物理力学特性输出表
  目标输出   实际输出   相对误差(%)
  轴向干缩   0.5624   0.5726   1.81
  弦向干缩   7.0111   7.6769   9.5
  径向干缩   0.7813   1.0772   37.87
  基本密度   0.2830   0.3432   21.27
  抗弯强度   56.2   63.1   12.28
  顺纹抗压强度   25.2   28.651   13.69
  抗弯弹性模量   5260   5929   12.72

Claims (4)

1、一种基于木材结构参数确定木材物理、力学特性的分析仪器,它包括:输入键盘,计算机主机,显示器或/和打印输出装置,其特征是,所述计算机主机内包括并依次连接:
1)输入单元,其中依次包括:用于接收木材结构参数的输入模块,参数辨识、校对及控制器,加法器;
2)接口引擎单元,其中包括:用于开启Matlab引擎的模块,以及并列联接的数据缓冲区、用于实现木材结构参数数据格式转换的模块、用于神经元网络模型文件的调用模块;
3)神经元网络非线性逼进单元,其中依次包括:用于确定木材结构参数影响木材物理力学特性的权值的权值函数的计算模块,用于通过自学习神经元网络层权值修正并产生权值的线性函数的模块;
4)输出单元,其中包括:材性量化指标的输出模块和相关输出曲线的输出模块。
2、根据权利要求1所述的分析仪器,其特征是,所述用于接收木材结构参数的输入模块包括下列参数输入模块中的一种或其组合:
a.基于对木材整体模型进行宏观分析的木材解剖微观参数的输入模块;
b.分开年龄段进行木材材性指标分析的参数输入模块;
c.重要解剖参数的输入模块。
3、根据权利要求1所述的分析仪器,其特征是,所述3)神经元网络非线性逼进单元中至少包括如下两个单元:
a.神经元网络的第一层单元,它包括:标准化处理器,神经元网络隐含层权矩阵模块,欧氏距离计算模块,自学习模型的修正模块,点乘模块,高斯函数发生器;所述标准化处理器和神经元网络隐含层权矩阵模块分别与欧氏距离计算模块联接,然后欧氏距离计算模块和自学习模型的修正模块分别与点乘模块联接,最后点乘模块与高斯函数发生器联接;
b.神经元网络的第二层单元,它包括并依次联接:神经元网络线性层权矩阵模块,正规化处理器,产生权值的线性函数发生器。
4、根据权利要求2所述的分析仪器,其特征是,所述a.基于对木材整体模型进行宏观分析的木材解剖微观参数的输入模块中包括下列参数:微纤丝角、胞壁率、射线比量、管胞比量、早材管胞长度、晚材管胞长度、早材管胞宽度、晚材管胞宽度、早材管胞T-壁厚度、晚材管胞T-壁厚度、早材管胞R-壁厚度、晚材管胞R-壁厚度、早材管胞T-直径、晚材管胞T-直径、早材管胞壁腔比、晚材管胞壁腔比、早材管胞腔径比、晚材管胞腔径比。
CN 200420096149 2004-09-27 2004-09-27 基于木材结构参数确定木材物理、力学特性的分析仪器 Expired - Lifetime CN2743834Y (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200420096149 CN2743834Y (zh) 2004-09-27 2004-09-27 基于木材结构参数确定木材物理、力学特性的分析仪器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200420096149 CN2743834Y (zh) 2004-09-27 2004-09-27 基于木材结构参数确定木材物理、力学特性的分析仪器

Publications (1)

Publication Number Publication Date
CN2743834Y true CN2743834Y (zh) 2005-11-30

Family

ID=35633852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200420096149 Expired - Lifetime CN2743834Y (zh) 2004-09-27 2004-09-27 基于木材结构参数确定木材物理、力学特性的分析仪器

Country Status (1)

Country Link
CN (1) CN2743834Y (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624114A (zh) * 2020-06-05 2020-09-04 内蒙古农业大学 一种基于木材微观形态特征评估木材物理力学性质的方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624114A (zh) * 2020-06-05 2020-09-04 内蒙古农业大学 一种基于木材微观形态特征评估木材物理力学性质的方法

Similar Documents

Publication Publication Date Title
CN1275042C (zh) 一种对中药生产工艺进行实时监控的方法
CN106446600A (zh) 一种基于CRISPR/Cas9的sgRNA的设计方法
CN101832922B (zh) 一种有机肥产品近红外模型转移方法
Fan et al. To how many simultaneous hypothesis tests can normal, student's t or bootstrap calibration be applied?
CN111157698B (zh) 一种利用发射率数据获取黑土土壤全钾含量的反演方法
CN106755441B (zh) 一种基于多性状的基因组选择进行林木多性状聚合育种的方法
CN115240772B (zh) 一种基于图神经网络的解析单细胞通路活性的方法
CN106845767A (zh) 一种行业科技发展力量化评估方法与评估系统
CN106127242A (zh) 基于集成学习的年极端降水预测系统及其预测方法
CN106886793B (zh) 基于判别信息和流形信息的高光谱图像波段选择方法
CN113420795B (zh) 一种基于空洞卷积神经网络的矿物光谱分类方法
CN115881232A (zh) 一种基于图神经网络和特征融合的scRNA-seq细胞类型注释方法
CN113570161A (zh) 基于宽度迁移学习的搅拌釜反应物浓度预测模型构建方法
CN2743834Y (zh) 基于木材结构参数确定木材物理、力学特性的分析仪器
CN114034684A (zh) 一种基于一站式光谱模型的土壤属性快速测定方法
CN109358185A (zh) 基于极限学习机的橡胶配方性能预测模型及预测方法
CN110555566B (zh) 一种基于b样条分位数回归的光电概率密度预测方法
CN117057464A (zh) 基于非线性记忆学习的土壤有机质光谱预测方法和装置
CN112149260A (zh) 一种三维抗冲击负泊松比结构的设计方法
CN1755362B (zh) 基于木材结构参数确定木材物理、力学特性的分析方法及其仪器
CN1699635A (zh) 金属电沉积过程中基于神经元网络进行实时控制的方法
CN112464541B (zh) 一种考虑多尺度不确定性的混合复合材料铺层方法
CN115728290A (zh) 一种土壤中铬元素的检测方法、系统、设备及存储介质
CN114999579A (zh) 一种筛选化合物分子描述符并确定其取值范围的方法
CN108182549B (zh) 一种考虑环境成本的区域能源效率测算方法

Legal Events

Date Code Title Description
C14 Grant of patent or utility model
GR01 Patent grant
AV01 Patent right actively abandoned

Granted publication date: 20051130

Effective date of abandoning: 20101124