采用窗函數的棉纖維線密度分布表征

蘇玉恒， 孔繁榮

(河南工程學院 紡織新產品開發河南省工程實驗室， 河南 鄭州 450007)

采用窗函數的棉纖維線密度分布表征

蘇玉恒， 孔繁榮

(河南工程學院 紡織新產品開發河南省工程實驗室， 河南 鄭州 450007)

為完整表征棉纖維的線密度分布，在引入表達棉纖維縱向幾何形態的窗函數的基礎上，推導出棉纖維線密度分布的密度函數表達式。對選用的5種籽棉試樣進行去籽和堿處理后，采用顯微鏡攝影和圖像處理的方法獲得棉纖維直徑沿長度方向的變化數據，運用最小二乘法擬合出了不同棉纖維試樣的窗函數。在規定以棉纖維中點為特定標記位的情況下，利用測試獲得的標記位直徑數據擬合了纖維間直徑的分布密度函數，從而構建了完整表征棉纖維線密度的分布密度函數。結果表明，棉纖維的窗函數為不對稱的二次多項式，而纖維間的直徑變化則服從正態分布，經對所獲得的分布函數期望的計算和數值模擬，在合理修正的基礎上，與實測數據有較高的一致性。

棉纖維； 線密度分布； 窗函數； 模擬

棉纖維線密度的分布與成紗的強力、不勻等指標有密切的關聯性[1-2]。由于棉纖維橫截面為腰圓形，且縱向中段粗，根端較細，梢端最細[3]，因此，棉纖維線密度的表達不僅是纖維間粗細變化，還與纖維的形態和長度相關。邱燕煒[4]在研究理想紗條的統計特征時，為了弱化理想紗條的條件，提出用窗函數的概念來反映纖維沿長度方向的線密度分布。Muneem Shahriar等[5]提出用計算機視覺觀察棉纖維縱向分布變化，由于棉纖維橫截面不是標準圓形，用纖維投影寬度作為直徑折算纖維橫截面時，需乘以一個纖維面積折算系數。林倩等[6]在假設每根纖維線密度相同的情況下，引入長度分布的非參數核密度估計，計算了不同長度下棉纖維之間的線密度分布，較為全面地給出了棉纖維線密度不勻及直徑不勻的計算方法。王曉東等[3]通過將梳理好的棉纖維須叢沿長度方向切斷稱重的方法，測量了棉纖維線密度沿其長度方向的分布，給出了當根梢顛倒合并后線密度分布的拋物線方程式，可以看作是纖維形態窗函數的一種表達。

上述研究結果為線密度的表征奠定了良好的基礎，但這些研究并沒有給出完整的表達方法和驗證。特別是在纖維形態和纖維之間線密度的表征方面并沒有給出合適的方法。本文在對棉纖維進行適當處理的基礎上，消除轉曲和非標準圓截面對測量的影響，采用顯微圖像觀察法獲得纖維縱向直徑的變化規律，在采用特定歸一化處理的基礎上，擬合出纖維的縱向直徑變化規律，即纖維形態窗函數。通過對纖維歸一化標準位置直徑的測試，擬合其分布的密度函數，用于表達纖維間線密度的分布特征。這樣將纖維間標記位特征分布和窗函數結合，便可完整表征棉纖維的線密度。

1 棉纖維線密度表征模型

纖維的線密度通常用長度為1 km的纖維的質量來表示，在采用中斷稱重法進行測量時所表示的是纖維的平均線密度。對于每根纖維來說，則是各不相同的，是一個隨機變量；對于一批纖維而言，其中的纖維線密度形成一種分布。對于從一個纖維集合中抽出的任一根纖維，設其線密度為Nt，則有：

(1)

式中：Wt為纖維的質量，g；L為纖維的長度，mm。

式(1)可進一步轉變為

(2)

式中：Vt為纖維的體積，mm3；φ為纖維的直徑，mm；ρ為纖維的密度，g/mm3；m為棉纖維的充實度系數(標準值為0.577[7])。

對于一根纖維而言，顯然φ是沿長度方向變化的一個變量，而對一個纖維集合而言，不同纖維間的φ是隨機變化的，是一個隨機變量，因此φ可以表示為

φ(x)=f(x)·ξ

(3)

式中：f(x)稱為窗函數，表示1根纖維延長度方向直徑變化的函數，反映了纖維的縱向形態，對1批纖維可假設其形態是一樣的；ξ表示纖維上某一特定位置的直徑，ξ反映了纖維之間直徑的變化，是一個隨機變量。那么對纖維集合體中的1根纖維，其線密度可進一步表示為

(4)

式中，ξ和L均為隨機變量。不同長度的纖維將會有不同的窗函數，為排除長度的影響，可將窗函數進行歸一化處理，這樣x的值域變為[0，1]，式(4)變為

(5)

式中：

(6)

可以看出，Nt與ξ2的分布相關，而與L的分布無關，即線密分布與長度分布相互獨立。設Nt的密度函數為n(r)，ξ具有密度函數q(t)，令Nt=kξ2，則棉纖維線密度的分布函數為

N(Nt)=P(Nt≤r)=P(kξ2≤r)=

(7)

所以有：

(8)

2 棉纖維線密度分布密度函數的擬合

2.1 窗函數的擬合

窗函數不僅有截斷作用，而且能夠起到平滑的作用[8]，可以更有效地保留原信號的頻率信息，在很多領域得到了應用。在對棉纖維的研究中，通常把棉纖維當作均勻體，邱燕煒[4]在研究理想紗條的統計特征時，提出了纖維線密度的窗函數的概念，用來表示棉纖維沿長度方向的形態變化。

選取5個來自河南的細絨棉籽棉試樣，由籽棉上取下纖維，分清纖維頭端和尾端，將取下的纖維浸泡在氫氧化鈉溶液中10 min，以消除纖維卷曲,取出洗掉堿液，用鑷子夾取纖維平直地放在載玻片上，在顯微鏡下進行觀察。轉動橫向手輪移動載玻片，每隔1 mm拍照1次，直至拍攝完整根纖維。圖1示出1～5號試樣的圖像。

圖1 試樣的纖維圖像(×400)Fig.1 Fiber images for samples(×400). (a) Sample 1; (b) Sample 2; (c) Sample 3; (d) Sample 4; (e) Sample 5

圖2 試樣的窗函數擬合曲線Fig.2 Fitted curves of windows functions for samples. (a) Sample 1; (b) Sample 2; (c) Sample 3; (d) Sample 4; (e) Sample 5

以棉纖維根端點為原點，向梢端點為正方向，建立坐標。將所拍攝的整根纖維照片經繪圖軟件處理，獲得纖維寬度的像素數，表示纖維在該處的直徑；以纖維中間位置(總長度的0.5處)作為基準，稱為標記位；其他各處直徑與標記位處的直徑相除，得到各處直徑相對于標記位的比率作為縱坐標。以纖維總長度為1，以觀測位置距纖維根端距離與總長度的比值表示觀測點的相對位置，作為橫坐標，以排除長度變化對窗函數的影響，繪制纖維縱向形態的散點圖。采用最小二乘法對纖維縱向直徑變化數據進行擬合，5個測試試樣的實測數據散點及擬合方程曲線如圖2所示，擬合方程見式(9)～(13)。

f1(x)=-1.332 0x2+1.055 2x+0.834 2

(9)

f2(x)=-0.490 4x2+0.154 3x+0.951 5

(10)

f3(x)=-0.266 7x2+0.156 0x+0.889 7

(11)

f4(x)=-0.605 0x2+0.452 2x+0.841 2

(12)

f5(x)=-1.093 5x2+0.850 1x+0.774 9

(13)

以上5個擬合方程采用F檢驗，結果如表1所示。表中數據表明5個函數均具有很好的顯著性。

表1 試樣1～5的窗函數檢驗結果Tab.1 Test results of fitted functions for samples 1 to 5

注：α為顯著性水平；v1、v2為自由度；F為統計量；Fα為臨界值。

從式(9)～(13)可見，棉纖維縱向形態變化，基本呈現為不對稱的二次多項式，且各系數符號相同。

其中試樣1和試樣2的各系數較為接近，而試樣3、試樣4及試樣5的系數較為接近，這表明棉纖維形態具有一定的趨同性，因此，采用窗函數來描述棉纖維的縱向形態變化將具有良好的通用性。

2.2 棉纖維間線密度分布的估計

纖維間的線密度分布，即ξ的分布，采用標記位置的直徑分布來表示。在顯微鏡下對5個試樣標記位置的纖維直徑進行測試，并對像素點所表示的實際尺寸進行標定，獲得標記位纖維直徑的實測值。圖3示出5種試樣標記位直徑分布的頻數直方圖及正態密度函數擬合。從擬合結果來看，纖維間標記位直徑分布基本服從正態分布，各試樣分布的密度函數如式(14)～(18)所示。

圖3 試樣標記位的直徑分布Fig.3 Diameters distribution of maker positions for samples.(a) Sample 1; (b) Sample 2; (c) Sample 3; (d) Sample 4; (e) Sample 5

(14)

(15)

(16)

(17)

(18)

擬合方程采用擬合優度檢驗的結果如表2所示。

2.3 擬合結果的對比

將式(14)～(18)代入式(8)，可得到各試樣的線密度分布函數；將式(9)～(12)代入式(6)，計算各試樣的k值。試樣1線密度分布的完整密度函數表達式如式(19)所示。

r∈(0,∞)

(19)

其中：

表2 試樣1～5正態密度函數的擬合優度檢驗Tab.2 Test results of goodness of fit for samples 1 to 5

取m=0.577，ρ=1.58×10-3g/mm3，根據線密度分布密度函數可計算出線密度分布的期望值，也可采用試樣的密度函數隨機生成一定量的線密度值計算平均線密度，得到模擬線密度值。表3示出5個試樣的線密度期望值、模擬值和實測值的對比。

表3 試樣1～5的線密度期望值、模擬值和實測值Tab.3 Expect, simulation and measured values of samples 1 to 5

表3結果中5個試樣的k值表現出與其對應的纖維截面直徑期望值的偏離，表明纖維形態變化確實存在著較多不同。根據線密度的密度函數計算的線密度期望值及數字模擬值與實測線密度也存在較大的差異，理論線密度較實測值偏大，其原因可能來自2個方面：一方面是棉纖維并非圓形截面，本文研究中盡管采用堿溶液對纖維進行了處理，但與圓形截面的假設仍存在差異，面積折算系數為0.8～0.9[9]；另一方面是纖維在直徑測試時經過了堿處理，產生了一定膨脹，一般在水中截面積增加21%～34%，而在燒堿溶液中截面積增加45%，長度增加6%[10]。表3中的修正值是按截面增加40%的情況下修正的理論密度函數的期望值，結果與實測值有較好的一致性。

3 結 語

本文在引入纖維縱向尺度截取窗函數的基礎上，提出了一個綜合表征棉纖維線密度的方法，推導出了棉纖維線密分布的密度函數，該函數與纖維的長度無關，與纖維標定位置直徑分布的平方相關。在此基礎上，選用5種籽棉進行去籽處理，得到5種棉纖維試樣。采用顯微鏡法測試棉纖維直徑沿長度方向變化，運用最小二乘法進行回歸分析，擬合出纖維的窗函數，5種試樣的窗函數均為不對稱的二次多項式，這與棉纖維外形中間粗、兩端細，且不對稱的形態是相符的。

在明確特定標記位的情況下，采用顯微鏡法測試了纖維間標記直徑的變化，并估計了5種棉纖維的纖維間幾何尺度分布密度函數，擬合優度檢驗的結果表明其較好地服從正態分布。在考慮合理的修正系數的情況下，根據對5個試樣密度函數的估計和測試結果進行對比，密度函數的期望值、模擬值均與實測線密度具有較好的一致性，但修正系數的選擇方法與標準仍需進一步研究。

引用窗函數獲得的棉纖維線密度分布密度函數，排除了長度分布的干擾，建立了綜合反映纖維線密度分布特征的密度函數，可精確表征棉纖維幾何特征，為纖維集合體的結構與性能研究及對紗線品質的精確預測和推斷提供有效的支持。FZXB

[1] 張朝輝，李梅.原棉性能與紗線強力的關系[J].紡織學報，2005，26(1)：52-53. ZHANG Chaohui, LI Mei. Relationship between raw cotton property and yarn strength[J]. Journal of Textile Research, 2005,26(1)：52-53.

[2] 嚴廣松，蘇玉恒，朱進忠.紗線條干不勻的隨機模擬[J].紡織學報，2012，33(1)：33-37. YAN Guangsong, SU Yuheng, ZHU Jinzhong. Stochastic simulation of yarn unevenness[J].Journal of Textile Research,2012, 33(1)：33-37.

[3] 王曉東，姚穆.棉纖維沿長度方向線密度之分布[J].西北紡織工學院學報，1990，13(1)：1-5. WANG Xiaodong, YAO Mu.The variance of linear density of cotton fiber along its longitudinal direc-tion[J]. Journal of Northwest Institute of Textile Science and Technology, 1990,13(1):1-5.

[4] 邱燕煒.理想紗條的統計特性[J].中國紡織大學學報，1999，25(4)：68-70. QIU Yanwei. Statistical analyses for the ideal sliver[J].Journal of China Textile University,1999,25(4):68-70.

[5] MUNEEM Shahriar, IAN Scott Fleming, Hamed Sari Sarraf, et al. A machine vision system to estimate cotton fiber maturity from longitudinal view using a transfer learning approach[J].Machine Vision and Applications, 2013,24(3)：1661-1683．

[6] 林倩，郁崇文.棉纖維線密度分布表征的探討[J].東華大學學報，2010，36(2)：124-128. LIN Qian, YU Chongwen. Discussion on characherizing the linear density distribution of cotton fiber[J]. Journal of Donghua University(Natural Science),2010,36(2)：124-128.

[7] 張冶，崔玉梅，穆征.棉纖維成熟度的測試與分析[J].棉紡織技術，2000，28(10)：33-35. ZHANG Ye, CUI Yumei, MU Zheng.Text and analysis on maturity of cotton fiber[J].Cotton Textile Technology, 2000,28(10):33-35.

[8] 鄭星亮，程潔，魏任之.數字信號處理中的窗效應及窗函數的應用原則[J].北京聯合大學學報，1997，11(2)：33-37. ZHENG Xingliang, CHENG Jie, WEI Renzhi. Window effects and the principle of using window functions in digital signal processing[J]. Journal of Beijing Union University, 1997,11(2):33-37.

[9] 趙寶元，吳雄英，丁雪梅，等.縱向測定棉纖維含量方法中棉纖維面積折算系數的實驗樣本量研究[J].中國纖檢，2010(4)：79-81. ZHAO Baoyuan, WU Xiongying, DING Xuemei, et al. Study on the minimum sample size in determination cross-sectional area multiplier of cotton fibers in longitudinal direction[J].China Fiber Inspection, 2010(4)：79-81.

[10] 袁霞.棉在燒堿溶液中溶脹的結構變化研究[J].上海紡織科技，2013，41(3)：17-19. YUAN Xia. Research on structure change of cotton fiber when swelling and dissolution in caustic soda solution[J]. Shanghai Textile Science and Technology, 2013,41(3):17-19.

Characterization of cotton linear density distribution using window function

SU Yuheng, KONG Fanrong

(HenanEngineeringLaboratoryofNewTextileDevelopment,HenanInstituteofEngineering,Zhengzhou,Henan450007,China)

In order to describe the cotton linear density completely, the function of linear density distribution was derived based on the window function expressing longitudinal geometry of cotton. After deseeding and alkali treating, diameter variations along length of five cotton samples were tested by photomicrography and image processing and the window functions of five samples were fitted by least square method using the test data. Marking the midpoint of fiber as the specific point, the distribution density function which could characterize the cotton linear density completely was fitted according to the diameter data of specific point. The results show that the window function of cotton approaches to asymmetrical quadratic polynomial, the diameter variation of fibers is normal distribution, and the linear density of cotton simulated and tested has high consistency based on calculating expectation of functions and reasonable correcting.

cotton fiber; linear density distribution; window function; simulation

10.13475/j.fzxb.20160200806

2016-02-01

2016-07-22

蘇玉恒(1971—)，男，教授，碩士。主要研究方向為棉紡紗結構及加工理論與實踐。E-mail:fzsyhsyh@126.com。

TS 102.2

A