• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于類龍伯格預測器的財務(wù)困境預警研究

    2023-03-02 01:03:15任嘉嵩伍胡星宇任煒杰
    會計之友 2023年5期

    任嘉嵩 伍胡星宇 任煒杰

    【摘 要】 預測財務(wù)困境一直是財務(wù)管理中的重要環(huán)節(jié),而傳統(tǒng)基于卡爾曼濾波的動態(tài)預警算法在應對噪聲分布不確定問題時難以保證預測結(jié)果的準確率,一定程度上降低了該方法的實際價值。針對其缺陷,文章創(chuàng)新性地提出了基于龍伯格觀測器理論和線性矩陣不等式技術(shù)的預測算法,在噪聲分布情況不能有效辨識或辨識誤差較大時預測效果更優(yōu),同時簡化了預測步驟,降低了計算復雜度。文章以滬深兩市電力企業(yè)為例,結(jié)合龍伯格觀測器理論設(shè)計預警模型,用線性矩陣不等式計算增益,應用數(shù)據(jù)分析和參數(shù)估計對各家電力企業(yè)財務(wù)狀態(tài)進行跟蹤和預測。同時與傳統(tǒng)基于卡爾曼濾波的算法橫向?qū)Ρ龋辛Φ赝怀鲱慅埌馗耦A測器的預警效果和實用價值。

    【關(guān)鍵詞】 財務(wù)困境; 類龍伯格預測器; 線性矩陣不等式; 狀態(tài)空間方程

    【中圖分類號】 F224.9? 【文獻標識碼】 A? 【文章編號】 1004-5937(2023)05-0038-08

    一、引言

    財務(wù)困境理論是公司管理中一個重要的研究領(lǐng)域[1]。當出現(xiàn)內(nèi)部管理不善、行業(yè)局勢動蕩、產(chǎn)品競爭力不足等情況,可能導致公司財務(wù)狀態(tài)由健康轉(zhuǎn)入困境甚至導致破產(chǎn)。不過公司陷入財務(wù)困境是一個動態(tài)的、連續(xù)的、可預測的變化過程[2]。20世紀60年代以來,西方學者對財務(wù)困境預警展開了廣泛的研究。從線性概率[3]、多元判別分析[4]、邏輯回歸[5]等靜態(tài)預警模型,到人工神經(jīng)網(wǎng)絡(luò)[6]、卡爾曼濾波[7]等動態(tài)預警模型,相關(guān)研究成果層見疊出。靜態(tài)預測多基于發(fā)生財務(wù)困境前的單一時間樣本數(shù)據(jù),其優(yōu)點是不需要考慮資金時間價值,計算簡便;缺點則是忽視了不同時間點之間的有機聯(lián)系。因為企業(yè)從健康轉(zhuǎn)入財務(wù)困境是一個累積變異特性的過程,財務(wù)狀態(tài)所表現(xiàn)的突變性往往是這種累積變異特性達到臨界值造成的。因此,在對財務(wù)困境設(shè)計動態(tài)預警模型時應具備兩個條件:第一,所選財務(wù)指標數(shù)據(jù)應該追溯數(shù)期,可反映財務(wù)陷入困境前較長時間的財務(wù)變動;第二,預警系統(tǒng)的設(shè)計應充分考慮到財務(wù)狀態(tài)的累積變異特性。

    在動態(tài)預警模型中,卡爾曼濾波作為當前較為主流的預測方法得到了廣大學者的青睞[7-15]。孫曉琳等[8]將廣泛應用于慣性導航、定位系統(tǒng)等領(lǐng)域的卡爾曼濾波引入財務(wù)系統(tǒng),研究結(jié)果表明對?觹ST公司的財務(wù)預測擁有較高的辨識能力和跟蹤精度。莊倩[9]在前人的基礎(chǔ)上對財務(wù)預警模型做出改進,設(shè)計了超前n步預測??柭鼮V波通過輸入輸出觀測數(shù)據(jù),以實時遞推的方式不斷修正系統(tǒng)狀態(tài)變量估計,且無需存儲大量數(shù)據(jù),通過觀測新數(shù)據(jù)即可得到新濾波值,但在設(shè)計上存在一定的缺陷:卡爾曼濾波的增益矩陣是在線更新的,存在因計算負擔過大而導致系統(tǒng)不穩(wěn)定以致預測失敗的情況[10]。同時,卡爾曼濾波需知道系統(tǒng)噪聲和測量噪聲的協(xié)方差矩陣,即噪聲的統(tǒng)計學特征。當噪聲的分布情況不能有效辨識或辨識誤差較大時,預測結(jié)果的準確性將受到影響。

    針對上述缺陷,本文以龍伯格(Luenberger)觀測器[16]代替卡爾曼濾波作為模型設(shè)計對象。龍伯格觀測器是一種有效解決動態(tài)系統(tǒng)控制率問題的方法,設(shè)計龍伯格觀測器的初衷是為了獲取模型所需的系統(tǒng)狀態(tài)變量。這些變量本身很難直接得到,但在實際控制系統(tǒng)中又有明顯的利用價值,基于狀態(tài)空間方程建立的龍伯格觀測器可將系統(tǒng)的輸入值與輸出值轉(zhuǎn)化為系統(tǒng)狀態(tài)變量以滿足模型需求。龍伯格觀測器結(jié)構(gòu)簡易,性能表現(xiàn)良好,控制效果優(yōu)異,在工業(yè)控制工程、電液控制系統(tǒng)和前饋控制策略中,選擇合適的反饋增益實現(xiàn)估計值和系統(tǒng)被估計值的零誤差,從而得到理想的系統(tǒng)狀態(tài)量[17]。與卡爾曼濾波方法相比,龍伯格觀測器離線計算增益,既能降低計算復雜度,又能規(guī)避在線更新增益的不穩(wěn)定情況。

    借鑒龍伯格觀測器在控制領(lǐng)域中采取合適的增益反饋實現(xiàn)系統(tǒng)零誤差,本文采用線性矩陣不等式(Linear Matrix Inequality,LMI)技術(shù)計算增益,同時在計算過程中引入抑制噪聲的設(shè)計條件,盡可能減少噪聲對預警結(jié)果精確度的影響。有學者研究發(fā)現(xiàn),線性矩陣不等式可用于解決非線性時滯的穩(wěn)定性等工程控制問題[18]。后期,學者提出多種線性矩陣不等式數(shù)值解法,包括替代凸投影算法、內(nèi)點法等。內(nèi)點法又分為中心點法、投影法、原始對偶法,其共同思路都是把線性矩陣不等式問題看成凸優(yōu)化問題處理。1995年矩陣實驗室(Matlab)集成了線性矩陣不等式工具箱(LMI Toolbox),這使得高維的線性矩陣不等式計算變成可能,同時也推動了運用線性矩陣不等式解決系統(tǒng)與控制問題的熱潮。

    基于這兩個理論在數(shù)學與控制領(lǐng)域的優(yōu)秀表現(xiàn),本文創(chuàng)新性地將龍伯格觀測器與線性矩陣不等式技術(shù)相結(jié)合設(shè)計了類龍伯格預測器,探究其在財務(wù)領(lǐng)域的適用性。相比傳統(tǒng)的卡爾曼濾波預測設(shè)計,本文的創(chuàng)新點如下:第一,類龍伯格預測器采用時不變的增益矩陣L,不需要實時在線更新計算,一方面降低計算復雜度,另一方面能有效保證預測系統(tǒng)的穩(wěn)定性,規(guī)避因在線計算而造成系統(tǒng)不穩(wěn)定的情況;第二,類龍伯格預測器無需知道系統(tǒng)噪聲和測量噪聲的協(xié)方差矩陣,因此當噪聲辨識誤差較大時,本文的方法依然具有較高的精確性;第三,預測系統(tǒng)的增益矩陣可利用Matlab中的線性矩陣不等式工具箱(LMI Toolbox)進行高效求解。

    本文設(shè)計類龍伯格預測器應用于財務(wù)困境預警體系,基于龍伯格觀測器理論構(gòu)建類龍伯格預測系統(tǒng),并運用線性矩陣不等式技術(shù)計算增益;將我國電力企業(yè)財務(wù)指標與模型相結(jié)合,同時對數(shù)據(jù)進行預處理和主成分分析;對比卡爾曼濾波方法,驗證類龍伯格預測器對財務(wù)狀態(tài)的跟蹤和預測效果。

    對全文所用到的符號做出如下解釋。Rn表示n維的歐幾里德空間,x∈Rn和A∈Rn×n分別代表一個n維的向量x和一個n×n維的矩陣A。矩陣P?酆0(P?芻0)表示P是一個正定(負定)的矩陣。AT為矩陣A的矩陣轉(zhuǎn)置。0和I分別為合適維度的全零矩陣和單位矩陣。在一個對稱矩陣中,?觹表示對應位置的元素可以由矩陣的對稱性得到,例如:

    ■?圳■

    二、預警體系設(shè)計

    類龍伯格預測器是一種時域方法,它基于系統(tǒng)狀態(tài)空間方程的構(gòu)造,利用系統(tǒng)中的可測量和觀測量誤差作為系統(tǒng)反饋,通過選擇合適的反饋增益使得反饋誤差迅速逼近零,以此獲得待觀測量。類龍伯格預測器的關(guān)鍵技術(shù)之一是狀態(tài)空間方程。在動態(tài)數(shù)學中,狀態(tài)空間方程描述動態(tài)系統(tǒng)可分為三個步驟:引出狀態(tài)變量、建立狀態(tài)方程、建立用于觀測各狀態(tài)的測量方程。

    本文用狀態(tài)方程和測量方程構(gòu)造財務(wù)困境動態(tài)預測方程組,系統(tǒng)中用于表征目標各時間域的最小內(nèi)部變量組xk∈Rn以列向量呈現(xiàn),式中x1,x2,…,xn分別對應x的n維度狀態(tài)變量:

    xk=x1,kx2,k■xn,k

    狀態(tài)空間即狀態(tài)變量的集合,狀態(tài)變量的維數(shù)等同于狀態(tài)空間的維數(shù)。其中狀態(tài)方程用于表征目標各時期財務(wù)指標間的聯(lián)系,測量方程用于表征目標財務(wù)指標與財務(wù)狀態(tài)的關(guān)聯(lián)。引入時間序列的狀態(tài)空間方程可列式為:

    狀態(tài)方程:

    xk=Axk-1+wk-1? ? ? ? ? (1)

    測量方程:

    Zk=Hxk+vk? ? ? ? ? ?(2)

    在式1和式2中,xk∈Rn是狀態(tài)向量;zk∈Rp是測量向量;wk∈Rn是系統(tǒng)噪聲向量;vk∈Rp是測量噪聲向量;A∈Rn×n是狀態(tài)轉(zhuǎn)移矩陣;H∈Rp×n是測量矩陣。

    (一)預警系統(tǒng)設(shè)計

    1.類龍伯格預測器設(shè)計

    不同于傳統(tǒng)的基于卡爾曼濾波技術(shù)的動態(tài)系統(tǒng)預警方法[8-15],本文設(shè)計全新的財務(wù)困境預測系統(tǒng),即類龍伯格預測系統(tǒng)。該系統(tǒng)由如下兩個子模型構(gòu)成,即式3、式4。

    財務(wù)過程狀態(tài)預測模型:

    ■k+1=A■k+L(zk-H■k)? ? ? ? (3)

    財務(wù)困境判別模型:

    ■k+1=H■k+1? ? ? ? ? ? ?(4)

    其中,■k+1為財務(wù)狀態(tài)的預測值,■k+1為根據(jù)狀態(tài)預測值得到的困境判別變量,L∈Rn×p為類龍伯格預測器增益矩陣。

    為了保證財務(wù)過程狀態(tài)向量的估計值■k不斷趨近于真實的狀態(tài)值(■(xk-■k)=0),進而能夠?qū)崿F(xiàn)判別財務(wù)困境的目的,式3需要滿足穩(wěn)定且不斷收斂的設(shè)計條件。定義ek=xk-■k為預測系統(tǒng)預測誤差,那么期望的設(shè)計條件也等價于系統(tǒng)的動態(tài)誤差系統(tǒng)穩(wěn)定收斂(■ek=0),由定義可知:

    ek+1=xk+1-■k+1? ? ? ? ? (5)

    在式5中帶入式1和式3可得:

    ek+1=Axk+wk-(A■k+L(zk-H■k))

    整理可得到預警系統(tǒng)的動態(tài)誤差表達式:

    ek+1=(A-LH)ek+wk-Lvk? ? ? ? (6)

    2.增益矩陣計算

    為了實時有效地預警財務(wù)困境,需要尋找合適的預警系統(tǒng)增益矩陣,使得預測的誤差系統(tǒng)穩(wěn)定且收斂。本節(jié)將尋找增益矩陣的問題轉(zhuǎn)化為線性矩陣不等式的求解問題,在實現(xiàn)預測要求的同時,有效克服了傳統(tǒng)基于卡爾曼濾波方法需要在每個時刻k計算預測增益矩陣的缺陷。接下來將提出基于線性矩陣不等式技術(shù)[18]計算財務(wù)過程狀態(tài)預測模型式3中增益矩陣L的定理并給出相應證明過程。

    定理1:若存在對稱正定矩陣P∈Rn×n、矩陣W∈Rn×p以及標量γ>0,使得線性矩陣不等式7和式8成立。

    P?酆0? ? ? ? ? (7)

    ■?芻0? (8)

    那么,所設(shè)計的財務(wù)困境預警系統(tǒng)模型式3、式4是穩(wěn)定收斂的,并且系統(tǒng)的動態(tài)誤差滿足如下歐幾里德范數(shù)條件。

    ■<γ■? ? ? ? ? ?(9)

    預警系統(tǒng)模型式3、式4的增益可由式10計算得到:

    L=P-1W? ? ? ? ? ? ?(10)

    為了簡化證明過程,式6可簡寫為:

    ek+1=(A-LH)ek+Fδk? ? (11)

    其中,F(xiàn)=[In-L],δk=wkvk。

    分析系統(tǒng)穩(wěn)定性時,本文首先選取能量函數(shù)Vk=e■■Pek,其中,P為合適維度的正定對稱權(quán)重矩陣。若誤差系統(tǒng)式11的能量不斷衰減,則可以認為該系統(tǒng)穩(wěn)定且收斂的,因為本文所考慮的財務(wù)困境預警系統(tǒng)模型式3、式4為離散數(shù)學模型,那么誤差系統(tǒng)能量不斷衰減用數(shù)學形式可表示為:

    ΔV=Vk+1-Vk<0? ? (12)

    即動態(tài)系統(tǒng)能量的差分函數(shù)始終小于零。考慮到Vk=e■■Pek和Vk+1=e■■Pek+1,同時將式11帶入式12可得:

    ΔV=(e■■(A-LH)T+δ■■FT)P×((A-LH))ek+Fδk)-e■■

    Pek? ?(13)

    若式13對于k≥0始終滿足小于零的條件,則動態(tài)誤差系統(tǒng)式11是穩(wěn)定的。但是考慮到狀態(tài)空間系統(tǒng)模型式1、式2中存在未知噪聲(系統(tǒng)噪聲wk和測量噪聲vk)的干擾,在設(shè)計類龍伯格預測系統(tǒng)時需要考慮抑制噪聲的設(shè)計條件。若要滿足式9所提出的范數(shù)條件,則需要能量函數(shù)的差分(ΔV)滿足:

    ΔV<-e■■ek+γ2δ■■δk

    將式13結(jié)果帶入可得到:

    (e■■(A-LH)T+δ■■FT)P((A-LH)ek-1+Fδk-1)-e■■Pek-1+e■■ek-1-γ2δ■■δk-1<0

    將上式寫為矩陣的形式,即:

    ek-1vk-1■■ek-1■k-1<0? (14)

    其中,

    ?漬1=(A-LH)TP(A-LH)-P+I

    ?漬2=(A-LH)TPF

    ?漬3=FTP(A-LH)

    ?漬4=FTPF-γ2I

    觀察可知,式14為二次型形式,故式14成立等價于:

    ■?芻0? ? (15)

    根據(jù)矩陣Schur補引理[19],由式15繼續(xù)推導可以得到:

    ■?芻0

    由定義可知,F(xiàn)[In -L],對上式進一步推導可得

    ■?芻0

    不難發(fā)現(xiàn),當前所得到的不等式為非線性的矩陣不等式(即存在待求得權(quán)重矩陣P和增益矩陣L以乘積PL的形式出現(xiàn)),求解上存在一定難度,故令W=PL∈Rn×p,可得到線性矩陣不等式8。

    (二)卡爾曼濾波算法

    作為類龍伯格預測器的對照組,基于卡爾曼濾波的預警方法同樣依賴于狀態(tài)空間模型式1和式2,并且在應用卡爾曼優(yōu)化迭代算法時需要提供系統(tǒng)噪聲(wk)和測量噪聲(vk)的協(xié)方差矩陣,故需假定wk和vk都是均值為零、相互獨立并且符合正態(tài)分布的高斯白噪聲序列,則其協(xié)方差矩陣分別定義為Q和R,用公式可表示為:

    E[wk]=0,E[wkw■■]=Q

    E[vk]=0,E[vkv■■]=R

    E[wkv■■]=E[vkw■■]=0

    參照張東等[14]設(shè)計的財務(wù)危機預警機制,本文基于卡爾曼濾波的算法可概括為如下五個步驟。

    第一步,假定k-1時最優(yōu)估計■■已知,進行一步預測:

    ■■=A■■■

    第二步,計算一步預測誤差方差矩陣:

    P■=A■P■A■■+Qk-1

    第三步,計算卡爾曼增益矩陣:

    Kk=P■H■■[HkP■H■■+Rk]-1

    其中,K為基于卡爾曼濾波的財務(wù)困境預警所需的增益,本文類龍伯格預測器所用增益為L。

    第四步,基于卡爾曼增益,結(jié)合最小均方誤差原則修正第一步的預測結(jié)果■■,得到最優(yōu)估計值:

    ■■=■■+Kk[Zk-Hk■■]

    第五步,計算最優(yōu)估計值■■的誤差方差矩陣:

    P■=[I-KkHk]P■

    卡爾曼濾波在隨機線性離散系統(tǒng)中的濾波計算是不斷重復預測與修正的遞歸過程,類龍伯格預測器的預測過程與之類似,但兩者仍存在以下區(qū)別:(1)本文所用方法以時不變的增益矩陣計算結(jié)果,而基于卡爾曼濾波方法的增益矩陣需要不斷更新[14];(2)本文方法在預測時不涉及協(xié)方差矩陣Q和R的計算,而基于卡爾曼濾波方法的預測效果與噪聲分布情況(即Q、R)能否有效辨識高度相關(guān)。這兩點都會對最終的測量結(jié)果造成影響,因此,可將卡爾曼濾波預測結(jié)果設(shè)置為對照組對比兩者的預測精準度。

    三、實證分析

    利用我國電力企業(yè)的財務(wù)指標對上文所設(shè)計的類龍伯格預測器展開印證。首先進行樣本數(shù)據(jù)的篩選;其次對數(shù)據(jù)做指標一致性和無量綱化處理,通過主成分分析得到公司綜合得分;最后界定模型參數(shù)和閾值,為財務(wù)困境預警分析做準備[12]。

    (一)樣本篩選

    電力行業(yè)的碳排放強度一直遠高于其他行業(yè),習近平總書記多次強調(diào)要加速綠色環(huán)保發(fā)展,革新高碳產(chǎn)業(yè)技術(shù),提倡資源高效利用。未來市場局勢變幻莫測,技術(shù)革新的同時往往伴隨著市場經(jīng)濟下滑等風險。因此,基于財務(wù)指標對電力企業(yè)展開財務(wù)困境預警將具有較高的現(xiàn)實意義。電力企業(yè)本身具有技術(shù)成熟度高、知識體系全面、發(fā)展趨勢穩(wěn)定的行業(yè)特征,利用公開數(shù)據(jù)可以較為準確地衡量不同時期電力企業(yè)的發(fā)展變化。本文選取滬深兩市電力企業(yè)為研究對象,依據(jù)現(xiàn)有財務(wù)指標進行財務(wù)困境預警研究。

    本文參考莊倩等[11]從償債能力、發(fā)展能力、經(jīng)營能力、現(xiàn)金流量、盈利能力5個財務(wù)類別中篩選資產(chǎn)負債率、可持續(xù)增長率、營業(yè)凈利率等25項指標,從電力企業(yè)中挑選?觹ST公司和健康公司作為研究對象,考慮季節(jié)波動性大、半年度指標參考價值不高等因素,以一年為一個周期,將公司被?觹ST當期記為k,向前追溯10期,分別對應k-1,k-2,…,k-10,用于各年財務(wù)狀態(tài)的記錄和對比。

    本文數(shù)據(jù)主要取自國泰安數(shù)據(jù)庫,對于缺失數(shù)據(jù)通過查閱公司年報等資料進行了核實和補充。最終選取32家?觹ST公司和32家健康公司,將每家公司為期11期的25項財務(wù)指標(指標體系見表1)作為樣本數(shù)據(jù)。

    (二)數(shù)據(jù)分析

    1.指標預處理

    財務(wù)困境預警結(jié)果由各指標所表達的信息決定,但不同指標數(shù)值大小與其表達信息的效果并不一致,即存在正負指標混雜的情況,需對其進行指標一致化處理。同時,指標的單位不同,則數(shù)據(jù)之間不具備可比性,可通過無量綱化處理,在保留數(shù)據(jù)差異性的同時將數(shù)值確定在[0,1]之間。

    (1)正指標預處理公式

    X'j=■

    (2)負指標預處理公式

    X'j=■

    式中,Xmax和Xmin為該指標對應公司各年份財務(wù)狀態(tài)的最高值與最低值,Xj和X'j為該指標預處理前后的數(shù)值。

    2.主成分分析

    為反映電力企業(yè)的真實情況,指標的選擇需全面且有代表性。但若處理所有的數(shù)據(jù)不僅會導致信息重疊,指標間的多重線性關(guān)系,還會引起較大誤差,增加預測結(jié)果的不確定性和復雜程度。因此,本文選擇主成分分析法,不僅能減少區(qū)別度不高的冗余指標個數(shù),還能盡可能地保留原有數(shù)據(jù)信息。主成分辨別時可依據(jù)兩項原則:(1)特征值>1;(2)累計貢獻率>97%,可較好地實現(xiàn)數(shù)據(jù)降維目的。最終,本文使用Matlab R2018b提取每家公司相應的主成分。受篇幅限制,僅將?觹ST東方主成分分析結(jié)果(表2)呈現(xiàn)于文中。

    在確認?觹ST東方主成分個數(shù)的同時,還可通過主成分分析得出基于這6個主成分11期的Z得分。依據(jù)主成分貢獻度數(shù)值和Z得分對應權(quán)重計算最終的綜合得分,得到公司不同時期的財務(wù)狀態(tài)。

    (三)參數(shù)估計

    若要進行有效預警,需要對式3和式4中的參數(shù)A和H進行參數(shù)辨識。另外,雖然本文提出的預測方法無需知道噪聲的統(tǒng)計學特征,但為了對比張東等[14]所用方法的預測效果,系統(tǒng)噪聲協(xié)方差矩陣(Q)、測量噪聲協(xié)方差矩陣(R)也需要通過辨識得到,同時本文沿用其極大似然估計方法辨識系統(tǒng)參數(shù)。以?觹ST東方公司的綜合得分為例,用Matlab R2018b計算可得出A=-0.153,H=1,對卡爾曼濾波預測模型中的協(xié)方差矩陣可辨識為Q=0.638,R=0。

    (四)閾值確定

    財務(wù)困境預警是依據(jù)財務(wù)判別閾值對不同時刻的財務(wù)狀態(tài)做出相應預警信號的,因此如何設(shè)置合理的閾值對預測準確度而言至關(guān)重要。本文依照中國證券監(jiān)督管理委員會發(fā)布的“針對財務(wù)狀況異常”的界定標準,借鑒莊倩等[13]的統(tǒng)計分析法設(shè)立本文財務(wù)困境判斷閾值,在95%置信水平計算財務(wù)困境的置信上下限。

    置信上限=XNST-za/2■

    置信下限=XST+za/2■

    式中,XNST和XST分別指樣本中所有?觹ST公司和健康公司財務(wù)狀態(tài)的平均值,za/2為置信系數(shù),■和■分別指樣本中所有?觹ST公司和健康公司財務(wù)狀態(tài)的標準差。

    利用Matlab R2018b計算得到財務(wù)困境判斷閾值:置信上限為-1.079,置信下限為-2.232。當預測值高于-1.079時,可判定該時刻財務(wù)狀態(tài)健康;當預測值介于-1.079和-2.232之間時,可判定該時刻公司面臨輕度財務(wù)困境;當預測值低于-2.232時,可判定該時刻公司陷入重度財務(wù)困境,公司必須及時管控風險,否則將面臨破產(chǎn)。

    四、財務(wù)困境預警效果分析

    (一)預測結(jié)果動態(tài)分析

    本文對64家樣本公司進行動態(tài)預測,參照張東等[14]文獻中的模型參數(shù)設(shè)定卡爾曼濾波對照組,并與本文所用方法對照呈現(xiàn)于圖中,受篇幅限制僅呈現(xiàn)4家財務(wù)狀態(tài)預測結(jié)果(圖1—圖4),并在本節(jié)最后展示全部公司的財務(wù)預警正確率。

    本文的方法是針對基于卡爾曼濾波方法在噪聲分布問題上存在的缺陷而設(shè)計的,為突出本文設(shè)計方法的預測優(yōu)勢,對?觹ST華電進行特殊處理(對噪聲協(xié)方差矩陣設(shè)定50%的不確定性)。同時以曲線圖展現(xiàn)預警結(jié)果,可更直觀地展現(xiàn)數(shù)據(jù)跟蹤和預警結(jié)果。

    1.?觹ST公司預測結(jié)果

    圖1為?觹ST東方公司,該公司在2008年(k-2)以前財務(wù)狀況一直呈良好波動的發(fā)展趨勢,但在2009年(k-1)財務(wù)變化趨勢陡降,數(shù)值由-0.786跌落至-3.689,陷入重度財務(wù)困境;2010年(k)狀態(tài)略微好轉(zhuǎn),但仍低于置信下限值。相比之下,基于卡爾曼濾波的預測值雖整體趨勢與真實值保持一致,但當財務(wù)狀態(tài)出現(xiàn)較大變動時,預測效果變差。該公司2010年被滬深證券標記?觹ST,而通過模型可以提早一年給出預警,也證明了類龍伯格預測器的良好表現(xiàn)。

    2.健康公司預測結(jié)果

    圖2和圖3是對皖能電力和寶新能源的預測,雖都為健康公司,但這兩家公司的財務(wù)狀態(tài)并不相似。皖能電力自2002年(k-10)以來,雖在2008年(k-4)存在財務(wù)狀態(tài)下滑的情況,但整體趨勢都是穩(wěn)步上升;而寶新能源雖然還未被標記成?觹ST,但在2017年(k-2)和2018年(k-1)數(shù)值都處于重度財務(wù)困境的邊緣,亟須得到企業(yè)的關(guān)注。寶新能源應當攻克現(xiàn)有難題,努力調(diào)整公司運營和財務(wù)狀態(tài),規(guī)避財務(wù)困境。

    從圖2和圖3中不難發(fā)現(xiàn),兩個模型對企業(yè)的財務(wù)狀況都進行了較好的跟蹤和預測。對皖能電力而言,兩種預測方法都能精準地體現(xiàn)財務(wù)變化趨勢,辨識出企業(yè)不斷上升的財務(wù)狀態(tài)。對寶新能源公司而言,雖并未標記?觹ST,但預測結(jié)果顯示其財務(wù)狀態(tài)并不樂觀,這說明本文所用方法對一些尚處在“健康”狀態(tài)的公司同樣具有一定的警示作用。

    3.特殊處理后預測結(jié)果

    圖4和圖5為?觹ST華電的特殊處理前后對比圖,此處的特殊處理為參數(shù)A和H保持不變,系統(tǒng)噪聲(wk)和測量噪聲(vk)對應的協(xié)方差矩陣Q和R各設(shè)定50%的不確定性,即模擬噪聲誤差對預警效果帶來的影響。如圖所示,兩種預測方法原先都有很好的預警效果,特殊處理后,類龍伯格預測器的預測曲線幾乎沒有變化,而卡爾曼濾波方法無法辨識出重度財務(wù)困境,失去了預警效果。

    噪聲誤差是一種客觀存在的未知干擾,數(shù)學建模無法保證每次都能得到最優(yōu)結(jié)果,當系統(tǒng)中的噪聲存在較大誤差時,輕則干擾預測結(jié)果精確度,重則嚴重偏離真實的模型動態(tài)。鑒于類龍伯格預測器無需系統(tǒng)噪聲和測量噪聲的協(xié)方差矩陣,因此當噪聲的分布情況不能有效辨識或辨識誤差較大時,該方法的預測結(jié)果依然具有較高精確性。

    (二)預警正確率分析

    本文選取樣本32家?觹ST公司和32家健康公司應用分類識別錯誤率對本次預警系統(tǒng)正確率進行測試。參考孫曉琳[15]采用的錯誤分類識別規(guī)定如下:當一家?觹ST公司被錯誤識別為健康公司時,可認為屬于型I類錯誤;當一家健康公司被錯誤識別為?觹ST公司時,可認為屬于型II類錯誤。兩種預警模型的識別正確率見表3。

    對比預警模型的識別正確率,可發(fā)現(xiàn)兩種模型都較好地實現(xiàn)了對財務(wù)狀態(tài)實際得分的追蹤和預警。而本文方法的優(yōu)勢在于:當噪聲分布情況不能有效辨識或辨識誤差較大時,依舊擁有較好的預測水平。建模、測量手段,甚至是財務(wù)數(shù)據(jù)真實性、市場波動等因素都會使得系統(tǒng)存在較大誤差,故本文提出的預警方法在兼顧了預測性能和設(shè)計復雜度的同時,還具有相對較強的可適用性。

    五、結(jié)論

    針對財務(wù)動態(tài)預警中存在的噪聲誤差、在線計算量大等問題,本文創(chuàng)新性地設(shè)計了類龍伯格預測器模型,對電力企業(yè)財務(wù)狀態(tài)展開預測,并設(shè)置卡爾曼濾波方法作為對照組。依據(jù)分析結(jié)果可知,類龍伯格預測器在財務(wù)領(lǐng)域具有較高的實踐價值,同時在處理噪聲和增益方面比卡爾曼濾波更具優(yōu)勢,預測成本小,預警正確率也更高。

    總的來看,實驗結(jié)果大致符合建模預期??紤]到本文首次基于龍伯格觀測器與線性矩陣不等式技術(shù)設(shè)計類龍伯格預測器進行財務(wù)困境預警,缺少與現(xiàn)有管理系統(tǒng)的有機結(jié)合,在未來可作為學習研究方向繼續(xù)挖掘類龍伯格預測器在財務(wù)領(lǐng)域的潛在價值。

    【參考文獻】

    [1] GILSON S C.Management turnover and financial distress[J].Journal of Financial Economics,1989,25(2):241-262.

    [2] JOHN T A.Accounting measures of corporate liquidity,leverage,and costs of financial distress[J].Financial Management,1993,22(3):91-100.

    [3] PURI M,ROCHOLL J,STEFFEN S.Global retail lending in the aftermath of the US financial crisis:distinguishing between supply and demand effects[J].Journal of Financial Economics,2011,100(3):556-578.

    [4] LIN T H.A cross model study of corporate financial distress prediction in Taiwan:multiple discriminant analysis,logit,probit and neural networks models[J].Neurocomputing,2009,72(16-18):3507-3516.

    [5] CHEN M Y.Predicting corporate financial distress based on integration of decision tree classification and logistic regression[J].Expert Systems with Applications,2011,38(9):11261-11272.

    [6] LI Y,MA W.Applications of artificial neural networks in financial economics:a survey[C].2010 International Symposium on Computational Intelligence and Design,2010.

    [7] OZBEK L,OZLALE U.Employing the extended Kalman filter in measuring the output gap[J].Journal of Economic Dynamics and Control,2005,29(9):1611-1622.

    [8] 孫曉琳,田也壯,王文彬.基于Kalman濾波的企業(yè)財務(wù)危機動態(tài)預警模型[J].系統(tǒng)管理學報,2010(4):408-414.

    [9] 莊倩.基于卡爾曼濾波的財務(wù)困境預測動態(tài)性研究[M].南京:東南大學出版社,2016.

    [10] 潘雅瓊,劉艷.科技型企業(yè)財務(wù)危機預警研究——基于生存分析法[J].會計之友,2021(8):64-70.

    [11] 莊倩,陳良華.基于卡爾曼濾波的企業(yè)財務(wù)困境動態(tài)預警模型[J].統(tǒng)計與決策,2015(24):190-192.

    [12] 吳樹暢,張雪,于靜.經(jīng)營風險與財務(wù)杠桿關(guān)系研究——基于異質(zhì)性負債視角[J].會計之友,2021(19):59-64.

    [13] 莊倩,陳良華.財務(wù)困境理論的預測動態(tài)性改進與應用[J].審計與經(jīng)濟研究,2014,29(5):69-76.

    [14] 張東,姚喬茜,王樂,等.基于卡爾曼濾波的財務(wù)困境預測方法與應用研究——以電力企業(yè)為例[J].東南大學學報(哲學社會科學版),2017(5):132-140.

    [15] 孫曉琳.基于狀態(tài)空間的財務(wù)危機動態(tài)預警模型在中國的實證研究[J].中國軟科學,2013(4):140-147.

    [16] LUENBERGER D.An introduction to observers[J].IEEE Transactions on Automatic Control,1971,16(6):596-602.

    [17] 任煒杰,郭勝輝,朱樹先.無人駕駛轉(zhuǎn)向預瞄控制系統(tǒng)狀態(tài)區(qū)間估計[J].控制工程,2021,28(9):1773-1779.

    [18] 楊曉光,張慶靈,靖新.內(nèi)點法在解線性矩陣不等式問題中的應用[C].中國數(shù)學力學物理學高新技術(shù)交叉研究學會第十二屆學術(shù)年會論文集,2008.

    [19] 俞立.魯棒控制:線性矩陣不等式處理方法[M].北京:清華大學出版社,2002.

    亚洲精品粉嫩美女一区| 国产一区二区亚洲精品在线观看| 精品乱码久久久久久99久播| 国内揄拍国产精品人妻在线| 久久亚洲精品不卡| 亚洲18禁久久av| 亚洲,欧美精品.| 波多野结衣高清作品| a级毛片免费高清观看在线播放| 成年版毛片免费区| 国产精品99久久久久久久久| 特级一级黄色大片| 少妇的逼水好多| 国产探花在线观看一区二区| 国产激情偷乱视频一区二区| 国产精品野战在线观看| a级一级毛片免费在线观看| 人人妻人人看人人澡| 一个人看视频在线观看www免费| 熟妇人妻久久中文字幕3abv| 亚洲av五月六月丁香网| 麻豆一二三区av精品| 色尼玛亚洲综合影院| 中文字幕免费在线视频6| 欧美国产日韩亚洲一区| 高潮久久久久久久久久久不卡| 国产伦精品一区二区三区四那| 欧美高清成人免费视频www| netflix在线观看网站| 日韩国内少妇激情av| 久久久久免费精品人妻一区二区| 中文字幕av在线有码专区| 国产成+人综合+亚洲专区| 国产乱人伦免费视频| 国内精品美女久久久久久| 12—13女人毛片做爰片一| 人人妻人人澡欧美一区二区| 亚洲片人在线观看| АⅤ资源中文在线天堂| 一级毛片久久久久久久久女| 精品久久久久久久末码| 日本 欧美在线| 最近最新免费中文字幕在线| 亚洲最大成人手机在线| 51国产日韩欧美| 老司机午夜福利在线观看视频| 久久伊人香网站| 脱女人内裤的视频| 欧美乱色亚洲激情| 又紧又爽又黄一区二区| 久久热精品热| 欧美bdsm另类| 日本精品一区二区三区蜜桃| 国产精品爽爽va在线观看网站| 偷拍熟女少妇极品色| 亚洲黑人精品在线| 久久午夜亚洲精品久久| 桃色一区二区三区在线观看| 如何舔出高潮| 淫妇啪啪啪对白视频| 91狼人影院| 亚洲精品一卡2卡三卡4卡5卡| 麻豆久久精品国产亚洲av| 午夜免费成人在线视频| 搡女人真爽免费视频火全软件 | 国产精品嫩草影院av在线观看 | 天天躁日日操中文字幕| 日韩 亚洲 欧美在线| 亚洲第一欧美日韩一区二区三区| 精品一区二区三区av网在线观看| 搞女人的毛片| 亚洲电影在线观看av| 国产av在哪里看| 色吧在线观看| 亚洲性夜色夜夜综合| 最近最新免费中文字幕在线| 午夜激情欧美在线| 国产av一区在线观看免费| 日韩欧美在线乱码| 国产精品亚洲一级av第二区| av天堂在线播放| 香蕉av资源在线| 俄罗斯特黄特色一大片| 禁无遮挡网站| 两人在一起打扑克的视频| 老女人水多毛片| 91狼人影院| 国产午夜福利久久久久久| 美女免费视频网站| 亚洲国产欧洲综合997久久,| 日本a在线网址| 国产欧美日韩一区二区精品| 国产精品伦人一区二区| 亚洲在线观看片| 伊人久久精品亚洲午夜| 色哟哟哟哟哟哟| av国产免费在线观看| www.999成人在线观看| 99国产精品一区二区蜜桃av| 麻豆成人午夜福利视频| 亚洲一区高清亚洲精品| 久久久精品大字幕| 亚洲欧美日韩东京热| 三级国产精品欧美在线观看| 色噜噜av男人的天堂激情| 欧美午夜高清在线| 在线观看一区二区三区| 久久这里只有精品中国| 久久精品影院6| 国产伦人伦偷精品视频| 午夜福利在线观看吧| 免费在线观看亚洲国产| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久香蕉精品热| bbb黄色大片| 精品久久久久久成人av| 真人做人爱边吃奶动态| 国产老妇女一区| 午夜影院日韩av| 国产欧美日韩精品亚洲av| 欧美激情在线99| 非洲黑人性xxxx精品又粗又长| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 91字幕亚洲| 国产精品一及| 一个人免费在线观看电影| 欧美成狂野欧美在线观看| 成熟少妇高潮喷水视频| 欧美最黄视频在线播放免费| 伊人久久精品亚洲午夜| 成人av一区二区三区在线看| 欧美bdsm另类| 午夜福利免费观看在线| 日本一二三区视频观看| 熟女人妻精品中文字幕| 丰满人妻熟妇乱又伦精品不卡| 99精品在免费线老司机午夜| 又黄又爽又免费观看的视频| 国产一区二区三区视频了| 中文字幕人成人乱码亚洲影| 变态另类丝袜制服| 免费人成视频x8x8入口观看| 国产人妻一区二区三区在| 久久精品国产亚洲av天美| 日韩欧美在线乱码| 变态另类成人亚洲欧美熟女| 婷婷精品国产亚洲av在线| 国产精品免费一区二区三区在线| 亚州av有码| 欧美xxxx黑人xx丫x性爽| 99热这里只有精品一区| 亚洲最大成人中文| 最好的美女福利视频网| 淫秽高清视频在线观看| 男人的好看免费观看在线视频| 精品人妻偷拍中文字幕| 亚洲片人在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久大精品| 日韩有码中文字幕| 夜夜夜夜夜久久久久| 有码 亚洲区| 国产成人影院久久av| 一卡2卡三卡四卡精品乱码亚洲| 亚洲午夜理论影院| 国产爱豆传媒在线观看| 国产精品av视频在线免费观看| 最近视频中文字幕2019在线8| 亚洲欧美日韩高清专用| 国内久久婷婷六月综合欲色啪| 亚洲成av人片在线播放无| 国产视频一区二区在线看| 麻豆成人av在线观看| 免费人成在线观看视频色| 国产精品98久久久久久宅男小说| 色噜噜av男人的天堂激情| 国产黄色小视频在线观看| 在线看三级毛片| 国产精品不卡视频一区二区 | 最近最新免费中文字幕在线| 国产久久久一区二区三区| 国产精品亚洲美女久久久| 亚洲电影在线观看av| 国产主播在线观看一区二区| 免费观看人在逋| 亚洲人成伊人成综合网2020| 国产精品野战在线观看| 三级国产精品欧美在线观看| 黄片小视频在线播放| 一区二区三区免费毛片| 亚洲avbb在线观看| 又黄又爽又免费观看的视频| 少妇人妻一区二区三区视频| 18禁在线播放成人免费| 中文字幕精品亚洲无线码一区| 国产精品野战在线观看| 久久久久久九九精品二区国产| 99在线人妻在线中文字幕| 在现免费观看毛片| 十八禁国产超污无遮挡网站| 熟女电影av网| 黄片小视频在线播放| 99国产精品一区二区蜜桃av| 不卡一级毛片| 国产伦精品一区二区三区视频9| 听说在线观看完整版免费高清| 欧美色欧美亚洲另类二区| 99国产精品一区二区三区| 一级a爱片免费观看的视频| 在线观看美女被高潮喷水网站 | 亚洲欧美清纯卡通| 欧美色欧美亚洲另类二区| 中文字幕av在线有码专区| 12—13女人毛片做爰片一| 欧美+亚洲+日韩+国产| АⅤ资源中文在线天堂| 国产精品自产拍在线观看55亚洲| 午夜免费激情av| 直男gayav资源| 精品福利观看| 亚洲美女视频黄频| 欧美xxxx黑人xx丫x性爽| av黄色大香蕉| 国产精品亚洲一级av第二区| 我的女老师完整版在线观看| 亚洲精品久久国产高清桃花| 十八禁国产超污无遮挡网站| 亚洲精品在线美女| 长腿黑丝高跟| 尤物成人国产欧美一区二区三区| 欧美又色又爽又黄视频| 日韩亚洲欧美综合| 午夜福利欧美成人| 国产成人aa在线观看| 老熟妇乱子伦视频在线观看| 波多野结衣高清作品| 99久国产av精品| 99久久久亚洲精品蜜臀av| 亚洲欧美日韩卡通动漫| 91狼人影院| 色吧在线观看| 国产成人欧美在线观看| av福利片在线观看| 久久久成人免费电影| 亚洲五月婷婷丁香| 少妇高潮的动态图| 乱人视频在线观看| 99久久精品一区二区三区| 搡女人真爽免费视频火全软件 | 精品午夜福利在线看| 久久草成人影院| 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| 成人精品一区二区免费| 欧美黑人欧美精品刺激| av黄色大香蕉| 日本 av在线| 老司机午夜十八禁免费视频| 我的老师免费观看完整版| 91久久精品国产一区二区成人| 赤兔流量卡办理| 我的女老师完整版在线观看| 日韩欧美在线二视频| 国产精品影院久久| 国产真实乱freesex| 亚洲色图av天堂| or卡值多少钱| 欧美潮喷喷水| 亚洲激情在线av| 一级av片app| 精品人妻偷拍中文字幕| 日韩亚洲欧美综合| 欧美三级亚洲精品| 欧美激情在线99| 在线国产一区二区在线| 亚洲内射少妇av| 国内久久婷婷六月综合欲色啪| 国产野战对白在线观看| 午夜激情欧美在线| 久久香蕉精品热| 国产成人福利小说| 一级作爱视频免费观看| 美女被艹到高潮喷水动态| 国产精品综合久久久久久久免费| 国产黄a三级三级三级人| 日韩 亚洲 欧美在线| 麻豆一二三区av精品| 国产午夜精品久久久久久一区二区三区 | 桃红色精品国产亚洲av| 九色国产91popny在线| 99久久久亚洲精品蜜臀av| 国产精品女同一区二区软件 | 国产成人欧美在线观看| 好看av亚洲va欧美ⅴa在| 自拍偷自拍亚洲精品老妇| 我要搜黄色片| 特大巨黑吊av在线直播| 黄色日韩在线| 国产精品久久久久久精品电影| 性色avwww在线观看| 九色国产91popny在线| 成人毛片a级毛片在线播放| 久久精品综合一区二区三区| 欧美黄色淫秽网站| 国产毛片a区久久久久| www.色视频.com| 麻豆国产av国片精品| 人人妻人人看人人澡| 3wmmmm亚洲av在线观看| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 深夜a级毛片| 欧美xxxx黑人xx丫x性爽| 日韩大尺度精品在线看网址| av天堂在线播放| 欧美在线一区亚洲| 午夜久久久久精精品| ponron亚洲| 免费观看人在逋| 精品久久久久久久久久久久久| av在线蜜桃| 舔av片在线| 内地一区二区视频在线| 亚洲 欧美 日韩 在线 免费| 久久久久久久午夜电影| 91麻豆av在线| 国产精品一区二区性色av| 听说在线观看完整版免费高清| 中文字幕免费在线视频6| 亚洲一区二区三区色噜噜| a级毛片免费高清观看在线播放| 亚洲第一欧美日韩一区二区三区| 变态另类成人亚洲欧美熟女| 99在线人妻在线中文字幕| 欧美丝袜亚洲另类 | 国产乱人伦免费视频| 久久人人爽人人爽人人片va | 在线十欧美十亚洲十日本专区| 国产私拍福利视频在线观看| 美女 人体艺术 gogo| 老司机深夜福利视频在线观看| 国产伦精品一区二区三区视频9| 国产欧美日韩一区二区精品| 国产成+人综合+亚洲专区| 欧美午夜高清在线| 精品无人区乱码1区二区| 免费无遮挡裸体视频| 哪里可以看免费的av片| 亚洲第一欧美日韩一区二区三区| 一区二区三区激情视频| or卡值多少钱| 欧美极品一区二区三区四区| 亚洲 欧美 日韩 在线 免费| 最近视频中文字幕2019在线8| 日本在线视频免费播放| 国内少妇人妻偷人精品xxx网站| 天美传媒精品一区二区| 国产在线精品亚洲第一网站| 人人妻,人人澡人人爽秒播| 中亚洲国语对白在线视频| 国产探花在线观看一区二区| 91av网一区二区| 男人和女人高潮做爰伦理| 亚洲一区二区三区不卡视频| 国产欧美日韩一区二区三| 中文字幕精品亚洲无线码一区| 国产白丝娇喘喷水9色精品| 久久人人爽人人爽人人片va | 婷婷亚洲欧美| 长腿黑丝高跟| 精品日产1卡2卡| 999久久久精品免费观看国产| 在线观看美女被高潮喷水网站 | 深夜a级毛片| 老鸭窝网址在线观看| 成人一区二区视频在线观看| 久久热精品热| 欧美成人性av电影在线观看| 免费在线观看日本一区| 3wmmmm亚洲av在线观看| 国产淫片久久久久久久久 | 热99在线观看视频| bbb黄色大片| 久久国产精品影院| 最近视频中文字幕2019在线8| 精品人妻一区二区三区麻豆 | 亚洲国产色片| 人妻制服诱惑在线中文字幕| 国产欧美日韩精品亚洲av| 亚洲欧美精品综合久久99| 我的老师免费观看完整版| 女人十人毛片免费观看3o分钟| av专区在线播放| 草草在线视频免费看| 欧美xxxx性猛交bbbb| 亚洲精品一卡2卡三卡4卡5卡| 噜噜噜噜噜久久久久久91| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 毛片女人毛片| 黄色配什么色好看| 久久精品综合一区二区三区| 亚洲国产精品久久男人天堂| 好男人电影高清在线观看| 97热精品久久久久久| 国产亚洲精品av在线| 深夜精品福利| 精品久久久久久成人av| 国产三级中文精品| 午夜老司机福利剧场| 内射极品少妇av片p| 少妇的逼好多水| 久久国产精品人妻蜜桃| 色噜噜av男人的天堂激情| 99热6这里只有精品| 小说图片视频综合网站| 在线观看美女被高潮喷水网站 | 黄色日韩在线| 最好的美女福利视频网| 亚洲五月婷婷丁香| 人妻制服诱惑在线中文字幕| 琪琪午夜伦伦电影理论片6080| 又爽又黄a免费视频| 美女大奶头视频| 国产人妻一区二区三区在| 精品不卡国产一区二区三区| 亚洲七黄色美女视频| 中文字幕高清在线视频| eeuss影院久久| x7x7x7水蜜桃| 不卡一级毛片| 日本 欧美在线| 久久久久久久久久成人| 色吧在线观看| 特大巨黑吊av在线直播| 嫩草影视91久久| 美女cb高潮喷水在线观看| 老熟妇乱子伦视频在线观看| 我要看日韩黄色一级片| 看免费av毛片| 亚洲激情在线av| 成熟少妇高潮喷水视频| 亚洲精品在线观看二区| 好男人在线观看高清免费视频| 一个人看的www免费观看视频| 国内毛片毛片毛片毛片毛片| 成人欧美大片| 久久久久久久久久黄片| 婷婷丁香在线五月| 人人妻,人人澡人人爽秒播| 18美女黄网站色大片免费观看| 亚洲无线观看免费| 琪琪午夜伦伦电影理论片6080| 久久人妻av系列| 国产精品av视频在线免费观看| 首页视频小说图片口味搜索| 国产免费一级a男人的天堂| 国产成人影院久久av| 成年版毛片免费区| 色综合站精品国产| 又黄又爽又刺激的免费视频.| 欧美精品国产亚洲| 亚洲无线观看免费| 久久天躁狠狠躁夜夜2o2o| 亚洲精品在线美女| 国产日本99.免费观看| www日本黄色视频网| 亚洲av电影不卡..在线观看| 少妇人妻一区二区三区视频| 亚洲成人久久爱视频| 国内精品美女久久久久久| 午夜福利在线观看吧| 成年版毛片免费区| 日日干狠狠操夜夜爽| 一本综合久久免费| 乱人视频在线观看| 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 看免费av毛片| 波多野结衣高清作品| 国产在线精品亚洲第一网站| 免费av不卡在线播放| 91麻豆精品激情在线观看国产| 日本黄色片子视频| 亚洲 国产 在线| 欧美色视频一区免费| 欧美三级亚洲精品| 亚洲五月婷婷丁香| 99视频精品全部免费 在线| 久久天躁狠狠躁夜夜2o2o| 国产白丝娇喘喷水9色精品| 久久久久久久精品吃奶| 黄色配什么色好看| 日日干狠狠操夜夜爽| 搞女人的毛片| 午夜亚洲福利在线播放| 亚洲无线观看免费| 亚洲性夜色夜夜综合| 韩国av一区二区三区四区| 九九热线精品视视频播放| 人人妻人人看人人澡| 国产精品爽爽va在线观看网站| 天堂av国产一区二区熟女人妻| 国产私拍福利视频在线观看| 很黄的视频免费| 亚洲人成网站在线播放欧美日韩| 性色avwww在线观看| 亚洲成人免费电影在线观看| 最近最新免费中文字幕在线| 国产亚洲精品av在线| 九九久久精品国产亚洲av麻豆| 中出人妻视频一区二区| 午夜福利高清视频| 久久久久久久精品吃奶| 久久久久久久久久成人| 免费观看人在逋| 国产免费男女视频| а√天堂www在线а√下载| 人妻制服诱惑在线中文字幕| av视频在线观看入口| 婷婷亚洲欧美| 亚洲精品影视一区二区三区av| 伦理电影大哥的女人| 色播亚洲综合网| 欧美日韩瑟瑟在线播放| 欧美成人性av电影在线观看| 欧美日韩瑟瑟在线播放| 色播亚洲综合网| 变态另类成人亚洲欧美熟女| 国模一区二区三区四区视频| 九九在线视频观看精品| 国产精品久久视频播放| .国产精品久久| 国产精品久久久久久人妻精品电影| 在线十欧美十亚洲十日本专区| 精品久久久久久久久亚洲 | 午夜福利在线观看吧| 波野结衣二区三区在线| 欧美zozozo另类| 国产精品1区2区在线观看.| 成人无遮挡网站| 一级黄色大片毛片| 一本久久中文字幕| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲精品久久久com| 国产精品一区二区性色av| 成人av在线播放网站| 国产乱人视频| 日本撒尿小便嘘嘘汇集6| 色在线成人网| 国产精品伦人一区二区| 欧美bdsm另类| 欧美国产日韩亚洲一区| 久久精品国产清高在天天线| 日本一本二区三区精品| 国产欧美日韩精品亚洲av| a级一级毛片免费在线观看| 我要看日韩黄色一级片| 内射极品少妇av片p| av视频在线观看入口| 亚洲av五月六月丁香网| 久久久久国内视频| 日本黄色视频三级网站网址| 欧美色视频一区免费| 国产乱人视频| 国产亚洲欧美98| 一级a爱片免费观看的视频| 国产成人啪精品午夜网站| 亚洲综合色惰| 精品一区二区三区av网在线观看| 少妇的逼好多水| 亚洲午夜理论影院| 搡老妇女老女人老熟妇| 伦理电影大哥的女人| 国产成人影院久久av| 精品一区二区三区人妻视频| 三级毛片av免费| 久久国产乱子伦精品免费另类| 一卡2卡三卡四卡精品乱码亚洲| 黄片小视频在线播放| 赤兔流量卡办理| 一夜夜www| 亚洲经典国产精华液单 | 日本成人三级电影网站| 男人和女人高潮做爰伦理| 人人妻,人人澡人人爽秒播| 国产v大片淫在线免费观看| 99久久精品热视频| 啦啦啦韩国在线观看视频| 国产精品女同一区二区软件 | 国产亚洲精品久久久久久毛片| 国产真实乱freesex| 自拍偷自拍亚洲精品老妇| 亚洲av电影不卡..在线观看| 又粗又爽又猛毛片免费看| 国产精品亚洲av一区麻豆| 欧美bdsm另类| 乱码一卡2卡4卡精品| 亚洲内射少妇av| 国产熟女xx| 丝袜美腿在线中文| 日本黄色视频三级网站网址| 精品人妻熟女av久视频| 神马国产精品三级电影在线观看| 国产精品伦人一区二区| 久久99热这里只有精品18| 淫妇啪啪啪对白视频| 淫秽高清视频在线观看| 黄色一级大片看看| 亚洲五月婷婷丁香| 午夜福利高清视频| 99热这里只有是精品在线观看 | 国产单亲对白刺激| 精品国产三级普通话版| 又黄又爽又刺激的免费视频.|