專利名稱：頻率成分分析裝置的制作方法
技術領域：
本發(fā)明涉及將針對內(nèi)燃機的運轉(zhuǎn)參數(shù)的檢測值的頻率成分分析與內(nèi)燃機的旋轉(zhuǎn)同步來進行的頻率成分分析裝置。
背景技術：
在專利文獻1中公開了這樣的信號處理裝置使用離散傅立葉變換(Discrete Fourier Transform)算法(以下稱為“DFT算法”)和高速傅立葉變換(Fast Fourier Transform)算法(以下稱為“FFT算法”)，進行安裝在內(nèi)燃機上的爆震傳感器的輸出信號的頻率成分分析。FFT算法為了提高運算速度而采用了特殊的算法，因而存在不能獲得需要的頻率成分(具體地說是爆震傳感器輸出信號的中心頻率成分)的強度的情況。因此，在上述信號處理裝置中，針對爆震傳感器輸出信號的中心頻率成分(強度最大的成分)，使用 DFT算法來計算強度，針對中心頻率成分以外的頻率成分，使用FFT算法來計算強度。專利文獻1 日本特開平11-303673號公報爆震傳感器輸出信號的中心頻率成分分析有必要基本上與內(nèi)燃機的旋轉(zhuǎn)同步進行，因而在內(nèi)燃機旋轉(zhuǎn)速度高的狀態(tài)下，需要進一步縮短頻率成分分析的運算時間。然而，由于在專利文獻1公開的裝置中沒有考慮到這一點，因而可確保用于使用通過頻率成分分析獲得的頻率成分強度的爆震判別處理的時間縮短，準確的判定很有可能變得困難。

發(fā)明內(nèi)容
本發(fā)明是著眼于上述這一點而作成的，本發(fā)明的目的是提供這樣的頻率成分分析裝置在與內(nèi)燃機的旋轉(zhuǎn)同步進行檢測參數(shù)的頻率成分分析時，可更適當?shù)貓?zhí)行頻率成分強度運算，特別是可提高在內(nèi)燃機旋轉(zhuǎn)速度高的狀態(tài)下的運算速度。為了達到上述目的，本發(fā)明提供一種與所述內(nèi)燃機的旋轉(zhuǎn)同步來進行針對內(nèi)燃機的運轉(zhuǎn)參數(shù)的檢測值(VKNK)的頻率成分分析的頻率成分分析裝置。該頻率成分分析裝置具有采樣單元，要素強度計算單元，以及頻率成分強度計算單元。所述采樣單元以預定時間(TSMP)間隔對所述運轉(zhuǎn)參數(shù)進行采樣，將采樣值轉(zhuǎn)換為數(shù)字值。所述要素強度計算單元針對預定數(shù)(ND)的采樣值計算與所述檢測值(VKNK)內(nèi)包含的多個頻率成分對應的第1要素的強度(DMFTS)和相對于所述第1要素相位偏移了 90度的第2要素的強度(DMFTC)。頻率成分強度計算單元使用所述第1要素強度(DMFTQ和所述第2要素強度(DMFTC)，與所述內(nèi)燃機的旋轉(zhuǎn)同步地計算所述多個頻率的頻率成分強度。而且，當所述內(nèi)燃機的旋轉(zhuǎn)速度 (NE)在設定閾值(NETH)以上時，所述頻率成分強度計算單元將所述第1要素強度的累計值的一部分和所述第2要素強度的累計值的一部分分別置換為上次計算值(TRMS1，TRMC1) 來進行所述頻率成分強度(STFT)的計算。根據(jù)該結(jié)構(gòu)，以預定時間間隔對內(nèi)燃機的運轉(zhuǎn)參數(shù)進行采樣，將采樣值轉(zhuǎn)換為數(shù)字值。根據(jù)獲得的數(shù)字值，針對預定數(shù)的采樣值計算與多個頻率成分對應的第1要素的強度和相對于第1要素相位偏移了 90度的第2要素的強度，使用第1要素強度和第2要素強度，與內(nèi)燃機的旋轉(zhuǎn)同步地計算多個頻率的頻率成分強度。當內(nèi)燃機的旋轉(zhuǎn)速度在設定閾值以上時，將第1要素強度的累計值的一部分和第2要素強度的累計值的一部分分別置換為上次計算值來進行頻率成分強度的計算。因此，在計算頻率成分強度時，無需再次進行第 1和第2要素強度的累計運算的一部分，可提高運算速度。其結(jié)果，能充分確保使用計算出的頻率成分強度的后處理的時間。期望的是，所述設定閾值(NETH)是根據(jù)獲得所述預定數(shù)(ND)的采樣值所需要的時間(TS)來設定的。根據(jù)該結(jié)構(gòu)，內(nèi)燃機旋轉(zhuǎn)速度的設定閾值是根據(jù)獲得預定數(shù)的采樣值所需要的時間來設定的。具體地說，當內(nèi)燃機旋轉(zhuǎn)速度增高且頻率成分強度的運算周期比獲得預定數(shù)的采樣值所需要的時間短時，在頻率成分強度的計算中應用的第1和第2要素的強度的上次計算值和本次計算值的一部分重復，因而能針對該重復的強度使用上次計算值。期望的是，所述設定閾值(NE)是根據(jù)所述頻率成分強度(STFT)的計算所需要的每單位時間的處理步驟數(shù)(賂)來設定的。根據(jù)該結(jié)構(gòu)，內(nèi)燃機旋轉(zhuǎn)速度的設定閾值是根據(jù)頻率成分強度的計算所需要的每單位時間的處理步驟數(shù)來設定的。如上所述，當頻率成分強度的運算周期比獲得預定數(shù)的采樣值所需要的時間短時，在頻率成分強度的計算中應用的第1和第2要素的強度的上次計算值和本次計算值的一部分重復，因而能針對該重復的強度使用上次計算值。然而，在實際的運算裝置中，存在這樣的情況通過進行將要素強度的上次計算值應用于頻率成分強度的運算的處理，處理步驟數(shù)增加，反而招致運算速度的降低。因此，通過將即使進行應用上次計算值的處理、每單位時間的步驟數(shù)也不增加的內(nèi)燃機旋轉(zhuǎn)速度作為設定閾值，可進行與運算裝置的性能(特性)對應的適當?shù)膹姸冗\算。期望的是，所述頻率成分強度計算單元使用在以與所述內(nèi)燃機的旋轉(zhuǎn)同步的觸發(fā)信號(CRK中斷)的產(chǎn)生定時為中心的預定期間(TS)內(nèi)獲得的采樣值來進行所述頻率成分強度(STFT)的計算。根據(jù)該結(jié)構(gòu)，由于使用在以與內(nèi)燃機的旋轉(zhuǎn)同步的觸發(fā)信號的產(chǎn)生定時為中心的預定期間內(nèi)所采樣的數(shù)字值來進行頻率成分強度的計算，因而可進行以成為對象的內(nèi)燃機的旋轉(zhuǎn)角度為中心的頻率成分分析。其結(jié)果，與采用使用在以產(chǎn)生觸發(fā)信號的時刻為起點的預定期間內(nèi)或者在以產(chǎn)生觸發(fā)信號的時刻為終點的預定期間內(nèi)所采樣的數(shù)字值的方法相比，可獲得難以受到內(nèi)燃機旋轉(zhuǎn)變動的影響的效果。


圖1是示出本發(fā)明的一實施方式涉及的內(nèi)燃機及其控制裝置的結(jié)構(gòu)的圖。圖2是用于說明爆震傳感器輸出的采樣和頻率成分分析的圖。圖3是用于說明計算頻率成分強度(STFT)的順序和定時關系的時序圖。圖4是用于說明頻率成分強度(STFT)的計算期間(TS)和CRK中斷的產(chǎn)生周期 (CRME)的關系的時序圖。圖5是示出發(fā)動機轉(zhuǎn)速(NE)和每單位時間的運算步驟數(shù)的關系的圖。圖6是用于說明計算頻率成分強度(STFT)的順序的流程圖。
圖7是將通過頻率成分分析獲得的強度數(shù)據(jù)表示為頻譜時序映射圖和二值化頻譜時序映射圖的圖。圖8是示出進氣門的落位噪聲的二值化頻譜時序映射圖的圖。圖9是為了說明噪聲去除處理而示出二值化頻譜時序映射圖的圖。圖10是為了說明使用主模式映射圖的爆震判定而示出二值化頻譜時序映射圖的圖。圖11是示出加權映射圖的一例的圖。圖12是爆震判定處理的流程圖。圖13是示出在圖12的處理中參照的映射圖的圖。圖14是在圖12所示的處理中執(zhí)行的二值化數(shù)據(jù)映射圖計算處理的流程圖。圖15是在圖14所示的處理中執(zhí)行的二值化處理的流程圖。圖16是示出在圖15的處理中參照的映射圖的圖。圖17是在圖12所示的處理中執(zhí)行的噪聲去除處理的流程圖。圖18是在圖12所示的處理中執(zhí)行的適合率計算處理的流程圖。圖19是用于說明在圖18的處理中參照的映射圖的圖。圖20是在圖12所示的處理中執(zhí)行的噪聲學習處理的流程圖。圖21是在圖20所示的處理中執(zhí)行的噪聲映射圖更新處理的流程圖。圖22是示出由爆震傳感器檢測出的振動的頻率成分強度的圖。
具體實施例方式以下，參照

本發(fā)明的實施方式。圖1是本發(fā)明的一實施方式涉及的內(nèi)燃機(以下稱為“發(fā)動機”)及其控制裝置的整體結(jié)構(gòu)圖，例如在4氣缸的發(fā)動機1的進氣管2的中途配設有節(jié)氣門3。節(jié)氣門3與檢測節(jié)氣門開度TH的節(jié)氣門開度傳感器4連接，傳感器4的檢測信號被提供給電子控制單元 (以下稱為“ECU”)5。針對每個氣缸，在位于發(fā)動機1與節(jié)氣門3之間、且處于進氣管2中的未圖示的進氣門的稍稍上游側(cè)位置設置有燃料噴射閥6，各噴射閥與未圖示的燃料泵連接并與ECU5電連接，根據(jù)來自該ECU5的信號來控制燃料噴射閥6的開啟時間。在發(fā)動機1的各氣缸內(nèi)設有火花塞7，火花塞7與E⑶5連接。E⑶5將點火信號提供給火花塞7。在節(jié)氣門3的下游側(cè)設有檢測進氣壓PBA的進氣壓傳感器8以及檢測進氣溫度TA 的進氣溫度傳感器9。在發(fā)動機1的主體上安裝有檢測發(fā)動機冷卻水溫度TW的冷卻水溫度傳感器10和非共振型的爆震傳感器11。傳感器8 11的檢測信號被提供給ECU5。作為爆震傳感器11，例如使用可檢測從5kHz到25kHz的頻帶的振動的傳感器。在進氣管2的節(jié)氣門3的上游側(cè)設有檢測吸入空氣流量GA的吸入空氣流量傳感器13，其檢測信號被提供給ECU5。E⑶5與檢測發(fā)動機1的曲軸(未圖示)的旋轉(zhuǎn)角度的曲軸角度位置傳感器12連接，與曲軸的旋轉(zhuǎn)角度對應的信號被提供給ECU5。曲軸角度位置傳感器12由以下傳感器構(gòu)成在發(fā)動機1的特定氣缸的預定曲軸角度位置輸出脈沖(以下稱為“CYL脈沖”)的氣缸判別傳感器；針對各氣缸的吸入行程開始時的上止點(TDC)，在預定曲軸角度之前的曲軸角度位置(在4氣缸發(fā)動機中是每180度曲軸角度)輸出TDC脈沖的TDC傳感器；以及以比TDC脈沖短的一定曲軸角度周期(例如6度的周期)產(chǎn)生1個脈沖(以下稱為“CRK脈沖”)的CRK傳感器，CYL脈沖、TDC脈沖以及CRK脈沖被提供給E⑶5。這些脈沖用于燃料噴射正時、點火正時等的各種正時控制、發(fā)動機轉(zhuǎn)速(發(fā)動機旋轉(zhuǎn)速度)NE的檢測。發(fā)動機1具有氣門工作特性改變裝置20，氣門工作特性改變裝置20具有第1氣門工作特性改變機構(gòu)和第2氣門工作特性改變機構(gòu)，第1氣門工作特性改變機構(gòu)連續(xù)變更進氣門(未圖示)的氣門升程量和開啟角(氣門開啟期間)，第2氣門工作特性改變機構(gòu)連續(xù)變更驅(qū)動進氣門的凸輪的以曲軸旋轉(zhuǎn)角度為基準的工作相位。ECU5將升程量控制信號和工作相位控制信號提供給氣門工作特性改變裝置20，進行進氣門的工作控制。第1和第2氣門工作特性改變機構(gòu)的結(jié)構(gòu)分別在例如日本特開2008-2M18號公報和日本特開 2000-227013號公報中示出。ECU5由以下電路等構(gòu)成具有對來自各種傳感器的輸入信號波形進行整形，將電壓電平修正為預定電平，將模擬信號值轉(zhuǎn)換為數(shù)字信號值等的功能的輸入電路；中央運算處理單元(以下稱為“CPU”)；存儲由該CPU執(zhí)行的各種運算程序和運算結(jié)果等的存儲電路 (存儲器)；以及將驅(qū)動信號提供給燃料噴射閥6和火花塞7的輸出電路。在本實施方式中，進行爆震傳感器11的輸出信號的頻率成分分析，根據(jù)該分析結(jié)果獲得的頻率成分強度進行爆震判定。因此首先，說明頻率成分分析的概要。圖2(a)示出爆震傳感器11的輸出信號波形，圖2(b)是將該圖2(a)的期間TS的波形進行放大示出的圖。在本實施方式中，將采樣周期TSMP設定為20微秒，以連續(xù)檢測出的50個數(shù)據(jù)作為對象，進行基于離散傅立葉變換(DFT)的頻率成分分析。該頻率成分分析的結(jié)果如圖2(c)所示。圖2(c)的縱軸是頻率成分強度STFT，在本實施方式中，與從5kHz到25kHz的頻帶中的每IkHz的頻率(5，6，7，…，24，25kHz)對應的頻率成分強度STFT是針對每6度曲軸角度計算出的(發(fā)動機1的曲軸每旋轉(zhuǎn)6度)。圖3是用于說明上述頻率成分分析的時序圖，圖3 (a)示出與依次存儲有通過每20 微秒對爆震傳感器11的輸出信號進行采樣獲得的50個檢測數(shù)據(jù)VKNK (數(shù)字值)的存儲器的地址對應的地址編號，圖3 (b)示出使用檢測數(shù)據(jù)VKNK所計算的正弦波成分強度DMFTS 和余弦波成分強度DMFTC對應于檢測數(shù)據(jù)VKNK被存儲在存儲器內(nèi)的狀態(tài)。地址編號與圖 3(a)相同，對應于1個檢測數(shù)據(jù)VKNK的采樣值，每1毫秒計算與從頻率5kHz到25kHz的 21個頻率對應的正弦波成分強度DMFTS和余弦波成分強度DMFTC，將它們存儲在存儲器內(nèi)。圖3(c)示出從曲軸角度位置傳感器12輸出的CRK脈沖，在本實施方式中，CRK脈沖的下降定時被用作運算執(zhí)行定時的基準(以下稱為“CRK中斷”)。在本實施方式中，存儲有發(fā)生CRK中斷的定時，使用與在以該定時為中心的采樣期間TS中獲得的50個檢測數(shù)據(jù)對應的正弦波成分強度DMFTS和余弦波成分強度DMFTC來計算各頻率(5 25kHz)的頻率成分強度STFT (j，i)(圖3(d))。這里，“j”是表示頻率的指數(shù)參數(shù)(以下稱為“頻率指數(shù)”)，j = 0、l、2、…、20對應于頻率5、6、7、…、25kHz?！癷”是表示發(fā)生CRK中斷的曲軸角度CA(將活塞位于上止點時的角度位置設定為0度)的指數(shù)參數(shù)(以下稱為“曲軸角度指數(shù)”)，i = 0、1、2、…、14對應于曲軸角度6、12、18、…、90度。在圖3所示的例子中，指數(shù)參數(shù)m(= 1 50)表示在頻率成分強度STFT (j，i)的計算中應用的數(shù)據(jù)。在最開始的CRK中斷定時，使用以地址編號“5”的檢測數(shù)據(jù)(對應于m = 25的檢測數(shù)據(jù))為中心的50個數(shù)據(jù)來計算頻率成分強度STFT(j，i)，在下一個CRK中斷定時，使用以地址編號“45”的檢測數(shù)據(jù)(對應于m = 25的檢測數(shù)據(jù))為中心的50個數(shù)據(jù)來計算頻率成分強度STFT(j，i+1)。另外，頻率成分強度STFT(j，i)是表示相對強度的
無量綱量。通過將在頻率成分強度STFT(j，i)的計算中應用的數(shù)據(jù)作為以與m = 25對應的檢測數(shù)據(jù)為中心的50個數(shù)據(jù)，換句話說，在以CRK中斷發(fā)生定時為中心的采樣期間TS中獲得的50個檢測數(shù)據(jù)，可進行以作為對象的曲軸角度為中心的頻率成分分析。其結(jié)果，與采用使用在以發(fā)生CRK中斷的時刻為起點的采樣期間內(nèi)或者在以發(fā)生CRK中斷的時刻為終點的采樣期間內(nèi)所采樣的檢測數(shù)據(jù)的方法相比，可獲得難以受到內(nèi)燃機旋轉(zhuǎn)變動的影響的效^ ο圖4是示出CRK中斷發(fā)生周期CRME和在頻率成分強度STFT的計算中應用的檢測數(shù)據(jù)的采樣期間(表示為“STFT計算”)TS的關系的圖。圖4(a)示出發(fā)動機轉(zhuǎn)速NE是 IOOOrpm的例子，采樣期間TS和中斷產(chǎn)生周期CRME —致。因此，當發(fā)動機轉(zhuǎn)速NE超過 IOOOrpm時，采樣期間TS—部分重復。圖4(b)示出發(fā)動機轉(zhuǎn)速NE稍高于2000rpm的例子， 采樣期間TS的重復部分變長。在本實施方式中，在這樣采樣期間TS重復的狀態(tài)下，通過使用上次計算值作為與重復部分相關的正弦波成分強度DMFTS的累計值和余弦波成分強度DMFTC的累計值，縮短頻率成分強度STFT的計算所需要的運算時間。不過，在使用上次計算值的運算中，當可使用的重復數(shù)據(jù)數(shù)的比率小時，依賴于運算裝置的特性，每單位時間的處理步驟數(shù)增加。圖5所示的實線Lll表示在本實施方式中使用的CPU中的在使用了上次計算值的情況下的每單位時間的處理步驟數(shù)NS和發(fā)動機轉(zhuǎn)速NE的關系，在達到轉(zhuǎn)速NESAT (以下稱為“步驟數(shù)飽和轉(zhuǎn)速NESAT”)之前，處理步驟數(shù)NS 增加，在NE ^ NESAT的范圍內(nèi)為一定值NSSAT。因此在本實施方式中，將閾值NETH設定為步驟數(shù)飽和轉(zhuǎn)速NESAT，當發(fā)動機轉(zhuǎn)速NE在閾值NETH以上時，進行使用上次計算值的運算。圖5所示的虛線L12表示在使用FFT(高速傅立葉變換)算法的情況下的發(fā)動機轉(zhuǎn)速NE和處理步驟數(shù)NS的關系。從圖5可以看出，在使用FFT算法的情況下，相對于發(fā)動機轉(zhuǎn)速NE的增加，處理步驟數(shù)NS線性增加，與此相對，通過進行使用DFT算法、且在發(fā)動機轉(zhuǎn)速NE高的區(qū)域(NE > NESAT)使用上次計算值的運算，能抑制處理步驟數(shù)NS的增加。圖6是示出頻率成分強度STFT計算順序的流程圖。在步驟Sl中，取得檢測數(shù)據(jù) VKNK，在步驟S2中，根據(jù)下述式(1)和(2)計算正弦波成分強度DMFTS (j，k)和余弦波成分強度DMFTC(j，k)。這里，指數(shù)參數(shù)k是圖3(b)所示的地址編號，取從“1”到ND(在本實施方式中是“50”)的值。并且At相當于(采樣周期XND)，在本實施方式中是1毫秒。DMFTS (j, k) = VKNK(k) X sin {2 π X (j+5) X1000X AtXk/ND} (1)DMFTC(j，k) = VKNK(k) X cos {2 π X (j+5) X 1000XAtXk/ND} (2)每個采樣周期(20微秒)執(zhí)行步驟Sl和S2的處理。在步驟S3中，判定是否發(fā)生了 CRK中斷，當未發(fā)生時回到步驟Si，在發(fā)生了 CRK中斷的定時，判別發(fā)動機轉(zhuǎn)速NE是否在閾值NETH以上。當NE < NETH時，根據(jù)不使用上次使用數(shù)據(jù)的下述第1式(A)計算頻率成分強度STFT，另一方面，當NE > NETH時，根據(jù)使用上次使用數(shù)據(jù)的下述第2式(B)計算頻率成分強度STFT (j，i)。指數(shù)參數(shù)m是被設定成使發(fā)生了 CRK中斷時的地址編號k為“25”的修正地址編號(參照圖3(d))。指數(shù)參數(shù)mp是在上次計算時應用的修正地址編號，mx是上次邊界指數(shù)，my是本次邊界指數(shù)，分別由下述式 (3)和⑷給出。
權利要求
1.一種頻率成分分析裝置，該頻率成分分析裝置與內(nèi)燃機的旋轉(zhuǎn)同步地進行針對所述內(nèi)燃機的運轉(zhuǎn)參數(shù)的檢測值的頻率成分分析，所述頻率成分分析裝置具有采樣單元，其以預定時間間隔對所述運轉(zhuǎn)參數(shù)進行采樣，將采樣值轉(zhuǎn)換為數(shù)字值；要素強度計算單元，其針對預定數(shù)目的采樣值計算與所述檢測值中包含的多個頻率成分對應的第1要素的強度以及相對于所述第1要素相位偏移了 90度的第2要素的強度；以及頻率成分強度計算單元，其使用所述第1要素強度和所述第2要素強度，與所述內(nèi)燃機的旋轉(zhuǎn)同步地計算所述多個頻率的頻率成分強度，當所述內(nèi)燃機的旋轉(zhuǎn)速度大于等于設定閾值時，所述頻率成分強度計算單元將所述第 1要素強度的累計值的一部分和所述第2要素強度的累計值的一部分分別置換為上次的計算值來進行所述頻率成分強度的計算。
2.根據(jù)權利要求1所述的頻率成分分析裝置，其中，所述設定閾值是根據(jù)獲得所述預定數(shù)目的采樣值所需要的時間而設定的。
3.根據(jù)權利要求1所述的頻率成分分析裝置，其中，所述設定閾值是根據(jù)所述頻率成分強度的計算所需要的每單位時間的處理步驟數(shù)而設定的。
4.根據(jù)權利要求1至3中任一項所述的頻率成分分析裝置，其中，所述頻率成分強度計算單元使用在以與所述內(nèi)燃機的旋轉(zhuǎn)同步的觸發(fā)信號的產(chǎn)生定時為中心的預定期間內(nèi)獲得的采樣值來進行所述頻率成分強度的計算。
5.一種頻率成分分析方法，該頻率成分分析方法與內(nèi)燃機的旋轉(zhuǎn)同步地進行針對所述內(nèi)燃機的運轉(zhuǎn)參數(shù)的檢測值的頻率成分分析，所述頻率成分分析方法具有以下步驟a)以預定時間間隔對所述運轉(zhuǎn)參數(shù)進行采樣，將采樣值轉(zhuǎn)換為數(shù)字值；b)針對預定數(shù)目的采樣值計算與所述檢測值中包含的多個頻率成分對應的第1要素的強度以及相對于所述第1要素相位偏移了 90度的第2要素的強度；以及c)使用所述第1要素強度和所述第2要素強度，與所述內(nèi)燃機的旋轉(zhuǎn)同步地計算所述多個頻率的頻率成分強度，當所述內(nèi)燃機的旋轉(zhuǎn)速度大于等于設定閾值時，將所述第1要素強度的累計值的一部分和所述第2要素強度的累計值的一部分分別置換為上次的計算值來進行所述頻率成分強度的計算。
6.根據(jù)權利要求5所述的頻率成分分析方法，其中，所述設定閾值是根據(jù)獲得所述預定數(shù)目的采樣值所需要的時間而設定的。
7.根據(jù)權利要求5所述的頻率成分分析方法，其中，所述設定閾值是根據(jù)所述頻率成分強度的計算所需要的每單位時間的處理步驟數(shù)而設定的。
8.根據(jù)權利要求5至7中任一項所述的頻率成分分析方法，其中，所述頻率成分強度計算單元使用在以與所述內(nèi)燃機的旋轉(zhuǎn)同步的觸發(fā)信號的產(chǎn)生定時為中心的預定期間內(nèi)獲得的采樣值來進行所述頻率成分強度的計算。
全文摘要
本發(fā)明提供一種與內(nèi)燃機的旋轉(zhuǎn)同步來進行針對內(nèi)燃機的運轉(zhuǎn)參數(shù)的檢測值的頻率成分分析的頻率成分分析裝置。以預定時間間隔對運轉(zhuǎn)參數(shù)進行采樣，將采樣值轉(zhuǎn)換為數(shù)字值。針對預定數(shù)的采樣值計算與檢測值內(nèi)包含的多個頻率成分對應的第1要素的強度和相對于第1要素相位偏移了90度的第2要素的強度。使用第1要素強度和第2要素強度，與內(nèi)燃機的旋轉(zhuǎn)同步地計算多個頻率的頻率成分強度。當內(nèi)燃機的旋轉(zhuǎn)速度在設定閾值以上時，將第1要素強度的累計值的一部分和第2要素強度的累計值的一部分分別置換為上次計算值來進行頻率成分強度的計算。
文檔編號G01M15/00GK102165172SQ20098013779
公開日2011年8月24日 申請日期2009年9月9日 優(yōu)先權日2008年9月26日
發(fā)明者義煎將之, 倉內(nèi)淳史, 北村夏子, 小松弘崇, 廣田和彥, 明本禧洙, 村上真, 高木治郎 申請人:本田技研工業(yè)株式會社