基于最小二乘法對(duì)艾煙濃度的光學(xué)濃度和PM10質(zhì)量濃度的數(shù)據(jù)擬合

韓　麗　趙　遷　劉　平　楊　佳　黃　暢　哈　略　和　蕊　王　昊　惠　鑫　潘秀頡　趙百孝

(1 北京中醫(yī)藥大學(xué)養(yǎng)生學(xué)研究所，北京，100029； 2 北京中醫(yī)藥大學(xué)針灸推拿學(xué)院，北京，100029； 3 北京新能源汽車有限公司，北京，102606； 4 華北電網(wǎng)有限公司北京電力醫(yī)院，北京，100073； 5 中國(guó)人民解放軍軍事醫(yī)學(xué)科學(xué)院，北京，100039)

?

基于最小二乘法對(duì)艾煙濃度的光學(xué)濃度和PM10質(zhì)量濃度的數(shù)據(jù)擬合

韓麗1,2趙遷3劉平2,4楊佳2黃暢2哈略2和蕊2王昊2惠鑫2潘秀頡5趙百孝2

(1 北京中醫(yī)藥大學(xué)養(yǎng)生學(xué)研究所，北京，100029； 2 北京中醫(yī)藥大學(xué)針灸推拿學(xué)院，北京，100029； 3 北京新能源汽車有限公司，北京，102606； 4 華北電網(wǎng)有限公司北京電力醫(yī)院，北京，100073； 5 中國(guó)人民解放軍軍事醫(yī)學(xué)科學(xué)院，北京，100039)

目的：基于最小二乘法對(duì)艾煙濃度檢測(cè)中的光學(xué)濃度和PM10質(zhì)量濃度進(jìn)行數(shù)據(jù)擬合，實(shí)現(xiàn)數(shù)據(jù)間的折算。方法：以CSED-A氣體染毒儀和HOPE-MED 8050系動(dòng)式染毒柜進(jìn)行光學(xué)濃度檢測(cè)，以遮光率(%)表示；以P5L2C便攜式微電腦粉塵儀檢測(cè)艾煙中可吸入顆粒物PM10濃度，以質(zhì)量濃度(mg/m3)表示。采用P5L2C光散射式數(shù)字微電腦粉塵測(cè)試儀分別監(jiān)測(cè)不同染毒柜染毒腔內(nèi)不同遮光率水平的艾煙的質(zhì)量濃度，監(jiān)測(cè)時(shí)采樣口原則上與鼠的高度一致，進(jìn)行多個(gè)遮光率的多次質(zhì)量-體積濃度的監(jiān)測(cè)，每個(gè)濃度重復(fù)3次。采用最小二乘法，運(yùn)用MATLAB分別對(duì)CSED-A氣體染毒儀和HOPE-MED 8050系列動(dòng)式染毒柜進(jìn)行曲線擬合，以實(shí)現(xiàn)數(shù)據(jù)間的折算。結(jié)果：CSED-A氣體染毒儀的擬合公式f(t)=414.4403*t^2+263.8095*t-3.0725 r2=85.8315；HOPE-MED 8050系列動(dòng)式染毒柜的擬合公式f(t)=182.3730*t^2+141.2080*t+4.4667 r2=3.7984。結(jié)論:實(shí)現(xiàn)通過擬合公式對(duì)艾煙光學(xué)濃度和PM10質(zhì)量濃度進(jìn)行折算，對(duì)今后艾煙實(shí)驗(yàn)的進(jìn)行起到指導(dǎo)作用。

艾煙；艾灸；遮光率；可吸入顆粒物；針灸

灸法是中醫(yī)學(xué)的特色療法，有著悠久的歷史，曾在數(shù)個(gè)朝代作為中醫(yī)學(xué)的主流療法。艾灸在防病治病[1-2]方面有著顯著的療效，適應(yīng)病種也日益增多。中國(guó)科學(xué)院院士承淡安先生認(rèn)為針灸的獨(dú)特作用不僅在于熱，更在于其特具的芳香氣味，艾煙被認(rèn)為[3-5]是艾灸發(fā)揮療效的一個(gè)途徑。近代也有研究發(fā)現(xiàn)艾煙能夠消毒空氣[6]、抑菌殺菌[7]，也有報(bào)道艾灸過敏[8-9]現(xiàn)象，過敏者會(huì)出現(xiàn)扁桃體腫大、咽癢、咽干、目痛、咽鼓管癢和昏睡等癥狀。

隨著艾灸安全性的研究的深入，艾煙濃度的檢測(cè)，也由簡(jiǎn)單的以艾條數(shù)目[10]描述演變?yōu)橘|(zhì)量濃度和光學(xué)濃度[11-12]表示方法，本文針對(duì)兩款可用于艾煙暴露的全自動(dòng)染毒儀器的濃度表示方法進(jìn)行了艾煙質(zhì)量濃度和光學(xué)濃度的折算，為今后艾煙實(shí)驗(yàn)的開展提供方便。

1　儀器與材料

1.1CSED-A氣體染毒儀CSED-A氣體染毒儀[13-14]，見圖1。為軍事醫(yī)學(xué)科學(xué)院研制全自動(dòng)動(dòng)態(tài)染毒裝置，北京慧榮和科技有限公司生產(chǎn)，目前已經(jīng)廣泛應(yīng)用于動(dòng)物毒理實(shí)驗(yàn)。煙霧濃度用激光的遮光率(%)來表示。

圖1　CSED-A

圖2　HOPE-MED 8050

1.2HOPE-MED 8050系列動(dòng)式染毒柜HOPE-MED 8050系列動(dòng)式染毒柜[11-12,15-18]，見圖2。為完全按照國(guó)際和國(guó)內(nèi)權(quán)威機(jī)構(gòu)標(biāo)準(zhǔn)設(shè)計(jì)生產(chǎn)，天津合普公司生產(chǎn)。煙霧濃度用激光的遮光率(%)來表示。

1.3P5L2C便攜式微電腦粉塵儀P5L2C便攜式微電腦粉塵儀[5,19-20]，北京賓達(dá)綠創(chuàng)科技有限公司生產(chǎn)，適用于環(huán)境中可吸入顆粒物PM10濃度及作業(yè)場(chǎng)所呼吸性粉塵、總粉塵濃度的快速測(cè)定煙霧濃度用質(zhì)量濃度(mg/m3)來表示。

1.4其他實(shí)驗(yàn)材料市售南陽(yáng)漢醫(yī)有限責(zé)任公司生產(chǎn)的三年陳艾灸條，規(guī)格18 mm×200 mm×10支，動(dòng)式染毒柜(CSED-A)，燃燒爐，酒精燈。

2　實(shí)驗(yàn)方法

2.1操作方法設(shè)置相關(guān)參數(shù)及時(shí)間。啟動(dòng)艾煙發(fā)生裝置，煙氣導(dǎo)入染毒腔，倉(cāng)內(nèi)艾煙煙氣濃度穩(wěn)定，各項(xiàng)參數(shù)達(dá)到預(yù)設(shè)值光學(xué)濃度透光度，煙氣動(dòng)態(tài)平穩(wěn)后，采用P5L2C光散射式數(shù)字微電腦粉塵測(cè)試儀(北京賓達(dá)綠創(chuàng)科技有限公司)監(jiān)測(cè)該染毒柜染毒腔內(nèi)的艾煙的質(zhì)量濃度。粉塵儀的位置應(yīng)均勻分布，采樣口原則上與大鼠的高度一致，進(jìn)行多個(gè)遮光率的多次質(zhì)量-體積濃度的監(jiān)測(cè)，每個(gè)濃度重復(fù)3次。記錄不同系統(tǒng)的光學(xué)濃度遮光率與質(zhì)量濃度的對(duì)應(yīng)數(shù)據(jù)，進(jìn)行數(shù)據(jù)擬合。

2.2數(shù)據(jù)曲線擬合構(gòu)造一個(gè)比較簡(jiǎn)單的函數(shù)y=φ(x)，使函數(shù)在觀測(cè)點(diǎn)的值等于已知的數(shù)值或?qū)?shù)值，尋找這樣的函數(shù)φ(x)。

2.2.1曲線擬合-最小二乘法已知離散點(diǎn)上的數(shù)據(jù)集，求得一解析函數(shù)y=f(x)，使f(x)在原離散點(diǎn)xi上盡可能接近給定yi的值，這一過程叫曲線擬合。

最常用的曲線擬合是最小二乘法曲線擬合[21-22]，擬合結(jié)果可使誤差的平方和最小，即找出使下式最小的f(x)：

最小二乘法曲線擬合的幾何意義見下圖，即諸點(diǎn)到曲線的距離平方和最小。

圖3　最小二乘法的幾何意義

通常，在解決實(shí)際問題時(shí)先將已知數(shù)據(jù)的散點(diǎn)圖畫出，然后設(shè)計(jì)擬合的曲線類型，最后根據(jù)某種準(zhǔn)則選定最佳的曲線。多項(xiàng)式擬合就是選擇適當(dāng)?shù)亩囗?xiàng)式對(duì)數(shù)據(jù)集進(jìn)行擬合。

2.2.2調(diào)用MATLAB命令實(shí)現(xiàn)擬合應(yīng)用MATLAB進(jìn)行多項(xiàng)式擬合的命令為：P=polyfit(x,y,m)

其中，x、y為已知數(shù)據(jù)點(diǎn)向量，分別表示橫、縱坐標(biāo)，m為擬合多項(xiàng)式的次數(shù)，結(jié)果返回m次擬合多項(xiàng)式系數(shù)，從高次到低次存放在向量p中。

y0=polyval(p,x0)，可求得多項(xiàng)式在x0處的值y0。

3　結(jié)果

3.1CSED-A氣體染毒儀用2次多項(xiàng)式曲線擬合數(shù)據(jù)點(diǎn)，編寫Matlab程序如下:

x=[0.04,0.06,0.08,0.1,0.2,0.4,0.5];

y=[6,14,21,28,71,162,236];

p=polyfit(x,y,2);

t=[0.04,0.06,0.08,0.1,0.2,0.4,0.5];

s=polyval(p,t);

plot(x,y,t,s)

得到

p=414.4403263.8095-3.0725

即擬合函數(shù)為f(t)=414.4403*t^2+263.8095*t-3.0725

招聘辦公室設(shè)在二樓，主考官把應(yīng)試者集中在一間會(huì)議室里，告訴大家，應(yīng)聘裝修工的等下去做現(xiàn)場(chǎng)考試，應(yīng)聘管理人員的到經(jīng)理辦公室分別面試。應(yīng)聘裝修工的被帶走了，就剩下三名管理人員應(yīng)試者坐在那，一個(gè)一個(gè)被叫到經(jīng)理室面試，竹韻是最后一個(gè)被叫進(jìn)去的。

得到試驗(yàn)數(shù)據(jù)及擬合曲線如圖4所示。

誤差的平方和為85.8315。

圖4　CSED-A氣體染毒儀數(shù)據(jù)及擬合曲線

圖5　HOPE-MED 8050系列動(dòng)式染毒柜數(shù)據(jù)擬合曲線

3.2HOPE-MED 8050系列動(dòng)式染毒柜用2次多項(xiàng)式曲線擬合數(shù)據(jù)點(diǎn)，編寫Matlab程序如下:

x=[0.004,0.008,0.02,0.05,0.11,0.15,0.2];

y=[3.882,5.82,8.408,12.08,22.57,28.75,40.44];

P=polyfit(x,y,2);

t=[0.004,0.008,0.02,0.05,0.11,0.15,0.2];

s=polyval(p,t);

plot(x,y,t,s)

得到p=182.3730141.20804.4667

即擬合函數(shù)為f(t)=182.3730*t^2+141.2080*t+4.4667

得到試驗(yàn)數(shù)據(jù)及擬合曲線如圖5所示。

誤差的平方和為3.7984。

4　討論與結(jié)論

實(shí)驗(yàn)中所使用的染毒系統(tǒng)均為全自動(dòng)動(dòng)態(tài)染毒裝置，監(jiān)測(cè)染毒柜內(nèi)的煙氣、氧氣等各項(xiàng)指標(biāo)的濃度，實(shí)時(shí)準(zhǔn)確，為理想中煙氣染毒染毒裝置，并已經(jīng)廣泛應(yīng)用于艾煙的毒理學(xué)實(shí)驗(yàn)[5，13-14，19-20]。實(shí)驗(yàn)通過對(duì)不同的染毒系統(tǒng)進(jìn)行不同遮光率條件下的質(zhì)量濃度的監(jiān)測(cè)取得光學(xué)濃度和質(zhì)量濃度的對(duì)應(yīng)數(shù)據(jù)，經(jīng)METLAB基于最小二乘法原理進(jìn)行多項(xiàng)式擬合遮光率與質(zhì)量濃度間的轉(zhuǎn)換公式，實(shí)現(xiàn)了2個(gè)艾煙染毒儀器中不同濃度的不同單位間的轉(zhuǎn)換，使艾煙濃度的描述更加多元化易于接受。

CSED-A氣體染毒儀[13-14，23-25]，為軍事醫(yī)學(xué)科學(xué)院研制全自動(dòng)動(dòng)態(tài)染毒裝置，煙霧濃度以激光的遮光率(%)來表示。本系統(tǒng)由內(nèi)腔、外腔、加熱制冷系統(tǒng)、氣體系統(tǒng)、控制系統(tǒng)、上位機(jī)控制系統(tǒng)六部分組成。艾煙在煙氣發(fā)生裝置產(chǎn)生之后導(dǎo)入該染毒儀，并在整個(gè)的實(shí)驗(yàn)過程中保證該濃度的穩(wěn)定性，同時(shí)能檢測(cè)和控制染毒腔內(nèi)的溫度、濕度、壓差、氧氣濃度以及受試煙霧的濃度以及其濃度累積值，并實(shí)時(shí)、動(dòng)態(tài)顯示其數(shù)值與實(shí)時(shí)曲線，以保障實(shí)驗(yàn)條件的單一性和可重復(fù)性。目前本儀器已廣泛應(yīng)用于動(dòng)物毒理實(shí)驗(yàn)。

HOPE-MED 8050系列動(dòng)式染毒柜[5，19-20，26-30]，為天津合普公司完全按照國(guó)際和國(guó)內(nèi)權(quán)威機(jī)構(gòu)標(biāo)準(zhǔn)設(shè)計(jì)生產(chǎn)。系統(tǒng)由動(dòng)物染毒柜0.3 m3的染毒倉(cāng)，風(fēng)速調(diào)節(jié)裝置，氧氣濃度、溫度、濕度、微壓力檢測(cè)裝置，鼠糞收集清洗系統(tǒng)和壓縮氣體凈化干燥器，動(dòng)物行為視頻監(jiān)視系統(tǒng)，以及I，II，III級(jí)廢氣處理系統(tǒng)組成。可用計(jì)算機(jī)控制或主機(jī)控制自動(dòng)完成小動(dòng)物單一或多種濃度情況下急性、亞急性、亞慢性和慢性的吸入染毒試驗(yàn)。動(dòng)靜洗三合一式染毒功能，具有控速度控時(shí)間控制燃燒量功能。煙霧濃度用激光的遮光率(%)來表示。

P5L2C便攜式微電腦粉塵儀[31-36]適用于環(huán)境中可吸入顆粒物PM10濃度及作業(yè)場(chǎng)所呼吸性粉塵、總粉塵濃度的快速測(cè)定。微電腦粉塵儀是在引進(jìn)日本柴田P-5L2型便攜式數(shù)字粉塵儀基礎(chǔ)上加入單片機(jī)控制的新一代快速測(cè)塵儀器，山東省曲阜華強(qiáng)礦業(yè)設(shè)備有限公司適用于公共場(chǎng)所可吸入顆粒物PM10濃度的快速測(cè)定、工礦企業(yè)生產(chǎn)現(xiàn)場(chǎng)勞動(dòng)衛(wèi)生方面呼吸性粉塵、總粉塵的測(cè)定。近年來有千臺(tái)以上P-5型儀器在全國(guó)各地使用，是同類測(cè)試儀器中最暢銷的產(chǎn)品。符合衛(wèi)生部標(biāo)準(zhǔn)《公共場(chǎng)所空氣中可吸入顆粒物(PM10)測(cè)定方法——光散射法》(99年唯一通過標(biāo)委會(huì)審定)；符合勞動(dòng)行業(yè)標(biāo)準(zhǔn)《空氣中粉塵濃度的光散射測(cè)定法》LD98-1996；符合鐵路行業(yè)標(biāo)準(zhǔn)《鐵路作業(yè)現(xiàn)場(chǎng)粉塵測(cè)定相對(duì)質(zhì)量濃度與質(zhì)量濃度轉(zhuǎn)換方法》TB/T2323-92。

CSED-A氣體染毒儀和HOPE-MED 8050系列動(dòng)式染毒柜分別被檢測(cè)，其艾煙的質(zhì)量濃度與光學(xué)濃度分別為f(t)=414.4403*t^2+263.8095*t-3.0725和f(t)=182.3730*t^2+141.2080*t+4.4667。分析公式不同的原因在于CSED-A氣體染毒儀和HOPE-MED 8050系列動(dòng)式染毒柜中檢測(cè)激光透光度的探頭的距離不同，距離越遠(yuǎn)，可檢測(cè)的艾煙濃度的范圍越小；探頭的距離越近，可檢測(cè)到的艾煙濃度的范圍越大。同一光學(xué)濃度透光度，探頭的距離越近，艾煙濃度越高，探頭的距離越遠(yuǎn)，艾煙濃度越小。光學(xué)濃度的使用應(yīng)當(dāng)與激光探頭的距離同時(shí)標(biāo)準(zhǔn)，方能準(zhǔn)確描述艾煙的濃度。

基于以上，實(shí)驗(yàn)基于最小二乘法原理對(duì)常用的艾煙暴露系統(tǒng)的艾煙濃度進(jìn)行了光學(xué)濃度和質(zhì)量濃度的數(shù)據(jù)分析擬合，為今后艾煙的毒理學(xué)實(shí)驗(yàn)提供了依據(jù)。

[1]韓麗，劉銅華，趙百孝，等.艾灸用于養(yǎng)生延緩衰老的研究進(jìn)展[J].世界中醫(yī)藥，2014，9(12)：1693-1696，1700.

[2]韓麗，趙百孝，劉銅華，等.艾灸防病保健與“長(zhǎng)壽基因”sirt1[J].中華中醫(yī)藥雜志，2015，30(6)：2023-2025.

[3]Wang Lei，Han Li，Liu Ping，etc.Moxa smoke：Is it a beneficial or hazardous factor for infertility？ A preclinical study on sperm parameters and sex hormones in male rats[J].European Journal of Integrative Medicine，2016，8(1)：28-35.

[4]Xu Huanfang，Zhao Baixiao，Cui Yingxue，etc.Effects of moxa smoke on monoamine neurotransmitters in samp8 mice[J].Evidence-Based Complementary and Alternative Medicine，2013.

[5]Cui Y.X.，Zhao B.X.，Huang Y.H.，etc.Effects of moxa(folium artemisiae argyi)smoke exposure on heart rate and heart rate variability in healthy young adults：A randomized，controlled human study[J].Evidence-Based Complementary And Alternative Medicine，2013.

[6]張國(guó)山，蘭蕾，常小榮，等.艾煙空氣消毒的研究進(jìn)展[J].世界中西醫(yī)結(jié)合雜志，2011，6(11)：1006-1009.

[7]呂熒，葉春枚，高建芳，等.熏灸治療外科感染性疾病575例(附艾煙抑菌試驗(yàn))[J].安徽中醫(yī)學(xué)院學(xué)報(bào)，1988，7(4)：36-37.

[8]Park J.E.，Lee S.S.，Lee M.S.，etc.Adverse events of moxibustion：A systematic review[J].Complement Ther Med，2010，18(5)：215-223.

[9]馮禎鈺，安靜.艾灸過敏一例[J].中國(guó)針灸，2005，25(12)：899.

[10]蘭蕾，常小榮，譚靜，等.艾煙的急性毒理試驗(yàn)[J].光明中醫(yī)，2011，26(10)：1992-1995.

[11]黃劍，韓麗，林岷瑜，等.不同濃度艾煙干預(yù)12周對(duì)大鼠肺與骨骼肌的影響[J].環(huán)球中醫(yī)藥，2014，7(9)：661-665.

[12]楊佳，趙百孝，韓麗，等.艾煙可控式動(dòng)態(tài)吸入染毒對(duì)大鼠血清中l(wèi)dl-r、icam-1及心臟組織形態(tài)的影響[J].中國(guó)針灸，2014,34(6)：573-577.

[13]吳彥，楊鴻，劉紅，等.中藥組分清除煙氣有害物質(zhì)及降低細(xì)胞毒性的研究[J].中國(guó)中藥雜志，2011，36(22)：3184-3188.

[14]劉永，胡旺順，齊紹武，等.亞慢性和長(zhǎng)期主流煙氣暴露小鼠血清乳酸脫氫酶和ɑ-羥丁酸脫氫酶活力的變化研究[J].中國(guó)煙草學(xué)報(bào)，2008，14(4)：55-57.

[15]易均鳳.自噬溶酶體途徑在熱習(xí)服保護(hù)熱射病腦損傷中的作用[D].重慶：第三軍醫(yī)大學(xué),2015.

[16]周次利，趙繼夢(mèng)，吳璐一，等.艾灸生成物對(duì)大鼠肺形態(tài)與臟器系數(shù)的影響[J].環(huán)球中醫(yī)藥，2014，7(6)：406-410.

[17]劉平，潘秀頡，韓麗，等.長(zhǎng)期艾煙干預(yù)對(duì)wistar大鼠外周血t淋巴細(xì)胞亞群及 cd4+ cd25+trek的影響[J].中國(guó)針灸，2013，33(2)：145-148.

[18]韓麗，趙百孝，劉銅華，等.不同濃度艾煙干預(yù)對(duì)大鼠肺臟及血清sod，mda表達(dá)的影響[J].中國(guó)針灸，2015，35(7)：687-690.

[19]黃茶熙，趙百孝，劉平，等.京津地區(qū)艾灸場(chǎng)所夏季可吸入顆粒物(pm10)的質(zhì)量濃度及微觀形貌分析[J].中華中醫(yī)藥雜志，2012，27(12)：3104-3108.

[20]Xu H.F.，Zhao B.X.，Cui Y.X.，etc.Effects of moxa smoke on monoamine neurotransmitters in samp8 mice[J].Evidence-Based Complementary And Alternative Medicine，2013.

[21]龔循強(qiáng)，李通，陳西江.總體最小二乘法在曲線擬合中的應(yīng)用[J].地礦測(cè)繪，2012，28(3)：4-6.

[22]陳良波，鄭亞青.基于最小二乘法的曲線擬合研究[J].無(wú)錫職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2012，11(5)：52-55.

[23]韓麗.艾煙的毒理學(xué)實(shí)驗(yàn)研究[D].北京：北京中醫(yī)藥大學(xué),2013.

[24]吳彥，楊鴻，楊陟華，等.中藥組分清除煙氣有害物質(zhì)及降低細(xì)胞毒性的研究[A].廣西煙草學(xué)會(huì)2011年學(xué)術(shù)年會(huì)論文集[C].廣西，2011:7.

[25]易軍艷.基于界面行為的多孔瀝青混合料凍融損傷特性研究[D].哈爾濱：哈爾濱工業(yè)大學(xué),2012.

[26]譚鐵紅.大氣混合污染物對(duì)大鼠肺組織及血清肺表面活性蛋白d表達(dá)的影響[D].長(zhǎng)春：吉林大學(xué),2013.

[27]張玉超.甲醛暴露對(duì)小鼠造血系統(tǒng)轉(zhuǎn)錄因子及骨髓組織毒性效應(yīng)的研究[D].武漢：華中師范大學(xué),2014.

[28]易均鳳，何根林，楊舉，等.熱習(xí)服對(duì)熱射病小鼠中樞神經(jīng)細(xì)胞自噬及凋亡的影響[J].第三軍醫(yī)大學(xué)學(xué)報(bào)，2015，37(11)：1059-1063.

[29]徐雄利，潘滬湘，任小孟，等.用于微生物氣溶膠動(dòng)態(tài)吸入實(shí)驗(yàn)暴露濃度監(jiān)測(cè)的采檢方法[J].海軍醫(yī)學(xué)雜志，2014，35(4)：324-325.[30]殷霄，賈肖輝，程浩，等.1,2-二氯乙烷對(duì)sd大鼠亞急性吸入毒性研究[J].中國(guó)職業(yè)醫(yī)學(xué)，2015，42(4)：378-384,391.

[31]朱鈺.基于現(xiàn)實(shí)生態(tài)位的道路植物景觀調(diào)查與分析[D].長(zhǎng)沙：中南林業(yè)科技大學(xué),2015.

[32]梁紅，陳曉雙，達(dá)良俊.上海佘山國(guó)家森林公園空氣負(fù)離子動(dòng)態(tài)及其主要影響因子[J].城市環(huán)境與城市生態(tài)，2014，7(1)：7-11.

[33]劉志云，賈艷艷，張玉潔.焊接煙塵濃度分布規(guī)律實(shí)驗(yàn)研究[J].科學(xué)技術(shù)與工程，2014，14(17)：178-181.

[34]黃玉海，李軍，崔瑩雪，等.艾煙對(duì)健康成年人血壓、呼吸頻率、心率、心電、血氧飽和度的影響[J].世界中醫(yī)藥，2014，9(6)：784-787.

[35]童海兵，邵丹，張珊，等.不同飼養(yǎng)密度對(duì)肉雞舍內(nèi)有害氣體、粉塵和微生物的影響[J].中國(guó)家禽，2014，36(20)：30-33.

[36]艾麗皎，馮義龍，譚興晏.重慶地區(qū)屋頂綠化現(xiàn)狀調(diào)查與分析[J].中國(guó)園林，2015，31(11)：27-30.

(2016-06-28收稿責(zé)任編輯：洪志強(qiáng))

Data-fitting on the Smoke Density, Optical Density and PM 10 Density of Moxa based on Least Square Method

Han Li1,2, Zhao Qian3, Liu Ping2,4, Yang Jia2, Huang Chang2, Ha Lue2, He Rui2,Wang Hao2, Hui Xin2, Pan Xiujie5, Zhao Baixiao2

(1InstituteofHealthPreservationofTraditionalChineseMedicine,BeijingUniversityofChineseMedicine,Beijing100029,China; 2SchoolofAcupuncture-MoxibustionandTuina,BeijingUniversityofChineseMedicine,Beijing100029,China; 3BeijingElectricVehicleCO.,LTD,Beijing102606,China; 4BeijingElectricPowerHospital,Beijing100073,China; 5ChineseMilitaryAcademyofMedicalSciences,Beijing100039,China)

Objective:Based on the least square method, to data fit the optical density and PM10 mass concentration of moxa smoke, and to realize the conversion between the optical density and PM10 mass concentration. Methods: The optical density were tested by CSED-A automatic exposure device and HOPE-MED 8050 automatic exposure device, and expressed as a shading rate (%); and the concentration of particulates of PM10 in moxa smoke was detected by P5L2C portable microcomputer instrument, and expressed as a mass concentration (mg/m3). When different shading rates were set up separately in chambers, the mass concentrations of PM10 were detected 3 mins for three times. Therefore, every shading rate could be expressed in PM10. The optical density in both device and mass concentration of MP10 were curve fitted using the least square method. Results: The fitting formula in CSED-A was f (t)=t ^ 2 + 263.8095 * 414.4403 * t-3.0725 r2=85.8315 and the fitting formula in HOPE-MED 8050 was f (t)=t + t ^ 2 + 182.3730 * 182.3730 * 4.4667 r2=3.7984. Conclusion: The optical density and PM10 mass concentration of moxa smoke could be converted to each other, which would play a guiding role in further research.

Moxa Smoke; Moxibustion; Shading Rate; PM; Acupunture

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(“973”計(jì)劃)項(xiàng)目(編號(hào)：2009CB522906)；國(guó)家自然科學(xué)基金項(xiàng)目(編號(hào)：81373730，81403449，81574068)；國(guó)家國(guó)際科技合作專項(xiàng)項(xiàng)目(編號(hào)：2011DFA31370)

韓麗(1981.04—)，女，博士，助理研究員，研究方向：保健灸的機(jī)制及安全性，E-mail;doctorhanli@163.com

趙百孝(1963.03—)，博士，教授，博士研究生導(dǎo)師，北京中醫(yī)藥大學(xué)針灸推拿學(xué)院院長(zhǎng)，研究方向：艾灸作用與原理，Tel：(010)64286737，E-mail:baixiao100@vip.sina.com

R245.81

A

10.3969/j.issn.1673-7202.2016.08.008