基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型

曹勇 ，李孟剛 ，李剛，常友玲

（東北大學(xué) 秦皇島分校1a.經(jīng)濟學(xué)院；1b.管理學(xué)院，河北 秦皇島 066004；2.北京交通大學(xué) 中國產(chǎn)業(yè)安全研究中心，北京 100044）

商業(yè)銀行信貸資金的優(yōu)化配置是金融資產(chǎn)優(yōu)化配置的一個特定領(lǐng)域，從金融資產(chǎn)優(yōu)化配置決策所采用的方法看，非線性數(shù)學(xué)規(guī)劃方法被經(jīng)常采用。不同金融資產(chǎn)優(yōu)化配置模型的差別，主要在所建立的目標(biāo)函數(shù)或約束條件上具有特色。Alexander等[1]通過最小化投資組合收益率的方差建立了股票投資組合優(yōu)化模型。Quaranta等[2]通過最小化投資組合的CVaR 值建立了股票投資組合的優(yōu)化配置模型。劉艷萍等[3]通過最大化貸款組合的效用函數(shù)建立了銀行貸款組合的優(yōu)化配置模型。Brandtner[4]對比研究了最小化CVaR 與最小化組合收益率方差對投資組合優(yōu)化配置決策的影響。

由于貸款組合的效用函數(shù)是一個關(guān)于貸款組合期望收益率的增函數(shù)和關(guān)于收益率方差的減函數(shù)，又由于無論是投資組合收益率的方差，還是投資組合的VaR 或CVaR 值都是反映投資組合收益率的波動幅度或分布特征，也就是投資組合的風(fēng)險，所以計量貸款組合的收益率和風(fēng)險是構(gòu)建銀行貸款組合優(yōu)化配置模型的基礎(chǔ)。Frey等[5]利用Bernoulli混合模型研究了信貸組合VaR 和CVaR 值的計算。Dietsch等[6]采用次序Probit模型和Gamma分布研究了法國中小企業(yè)貸款組合的風(fēng)險計量。Lucas等[7]采用狀態(tài)轉(zhuǎn)換模型探討了經(jīng)濟擴張和衰退情況下的信用轉(zhuǎn)移矩陣及其對信貸組合風(fēng)險的影響。Rosen等[8]研究了信貸組合風(fēng)險產(chǎn)生的原因及其對風(fēng)險的邊際貢獻。Tsaig等[9]探討了信用遷移對銀行貸款組合和債券組合價值方差的影響。

上述關(guān)于貸款組合風(fēng)險計量方法的研究為銀行貸款組合優(yōu)化配置模型的構(gòu)建奠定了基礎(chǔ)。遲國泰等[10]在貸款組合收益率既定及VaR 約束條件下，通過最小化貸款組合收益率的方差建立銀行貸款組合優(yōu)化決策模型。郭戰(zhàn)琴等[11]在VaR 約束下，通過最小化貸款組合收益率的方差和最大化貸款組合期望收益率建立貸款組合優(yōu)化的多目標(biāo)規(guī)劃模型。劉艷萍等[3]在VaR 約束下，通過最大化貸款組合的效用函數(shù)建立銀行貸款組合的優(yōu)化配置模型。洪忠誠等[12]利用信用風(fēng)險遷移矩陣，通過最小化貸款組合的CVaR 值建立銀行貸款組合的優(yōu)化模型。

現(xiàn)有貸款組合優(yōu)化配置模型的研究取得了長足進展，但仍然存在需要進一步深入研究的問題?，F(xiàn)有商業(yè)銀行信貸組合優(yōu)化配置模型多是探討銀行信貸資金在公司間的優(yōu)化配置，這類研究對于商業(yè)銀行同時對多家公司進行貸款的決策問題具有重要參考意義，但現(xiàn)有研究很少探討商業(yè)銀行信貸資金如何在行業(yè)之間進行優(yōu)化配置。

本文在現(xiàn)有研究基礎(chǔ)上，進一步探討銀行信貸資金在行業(yè)間的優(yōu)化配置問題，從而使商業(yè)銀行信貸資金在整體上達到最優(yōu)。本文的創(chuàng)新與特色在于：以行業(yè)內(nèi)代表性上市公司市值加權(quán)的股票價格作為行業(yè)平均股價，以市值加權(quán)的每股負(fù)債作為行業(yè)平均每股負(fù)債。根據(jù)行業(yè)平均股價與行業(yè)平均每股負(fù)債，應(yīng)用KMV 模型計算行業(yè)違約概率，解決了行業(yè)違約概率不能直接計算的難題。在應(yīng)用KMV模型計算行業(yè)違約概率的基礎(chǔ)上，利用正態(tài)Copula函數(shù)計算m個行業(yè)2m種違約狀態(tài)的聯(lián)合違約概率，解決了行業(yè)違約狀態(tài)的聯(lián)合概率因不可直接觀測而難以計量的問題。以對m個行業(yè)的貸款權(quán)重為優(yōu)化變量，以銀行整體信貸資金的差異系數(shù)最小化為目標(biāo)函數(shù)，建立基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型，求解出每個行業(yè)的最優(yōu)貸款權(quán)重，能夠使得銀行整體信貸資金獲取單位收益的風(fēng)險最小。本文構(gòu)建的基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型，對于商業(yè)銀行信貸資金在行業(yè)間的優(yōu)化配置與風(fēng)險管理具有重要參考意義。

1 基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型構(gòu)建原理

1.1 研究目標(biāo)

商業(yè)銀行在貸款發(fā)放過程中，既需要在微觀上考慮對哪些公司貸款及貸款比例為多少，能夠使得貸款風(fēng)險一定時收益最大，或收益一定時風(fēng)險最小；也需要在宏觀上控制信貸資金在各行業(yè)間的配置比例，從而分散貸款的集中度，在保證信貸資金收益率的同時降低風(fēng)險，使得商業(yè)銀行信貸資金在整體上達到最優(yōu)。

本文構(gòu)建了基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型，為銀行信貸資金在行業(yè)間的優(yōu)化配置提供決策參考。

1.2 開展研究的難點

（1）行業(yè)違約風(fēng)險難以計量?，F(xiàn)有違約風(fēng)險計量模型，如Logistic、KMV、Credit Risk+、Credit Metrics等模型，均是針對公司或上市公司，其違約風(fēng)險計量的依據(jù)是公司財務(wù)指標(biāo)或上市公司的股票價格與每股負(fù)債等。而對于行業(yè)而言，既不存在確定的財務(wù)指標(biāo)值，也無對應(yīng)的股票價格與每股負(fù)債。因此，行業(yè)違約風(fēng)險難以計量，而違約風(fēng)險的計量又是貸款定價與信貸資金優(yōu)化配置的基礎(chǔ)參數(shù)。

（2）行業(yè)間違約狀態(tài)的聯(lián)合概率不可觀測。違約狀態(tài)指在m筆貸款到期時，每筆貸款是否如期償付本息的狀況，由于每筆貸款均存在到期違約的可能性，即每筆貸款到期時均存在違約或不違約兩種狀態(tài)，故理論上在m筆貸款到期時共有2m種違約狀態(tài)。將每一種違約狀態(tài)可能發(fā)生的概率稱為違約狀態(tài)的聯(lián)合概率。由于行業(yè)間違約狀態(tài)的聯(lián)合概率實際不可觀測，故違約狀態(tài)聯(lián)合概率的測算存在困難。

（3）銀行信貸資金收益率波動的分布特征難以確定。在上述不同的違約狀態(tài)下，銀行信貸資金的收益率是變化的，如果能夠確定信貸資金收益率的分布特征，則可據(jù)以計算銀行信貸資金收益率的方差、VaR 或CVaR 值等，進而可構(gòu)建銀行貸款資金優(yōu)化配置模型。但是由于行業(yè)間違約狀態(tài)的聯(lián)合概率不可直接觀測，故銀行信貸資金收益率的期望及其方差難以計算，亦即銀行信貸資金收益率的分布特征難以確定。

1.3 解決難點的思路

（1）行業(yè)違約風(fēng)險計量的思路。將行業(yè)視為行業(yè)內(nèi)代表性上市公司的加總，以行業(yè)內(nèi)代表性上市公司市值加權(quán)的股票價格作為行業(yè)平均股價，以市值加權(quán)的每股負(fù)債作為行業(yè)平均每股負(fù)債。根據(jù)行業(yè)平均股價與行業(yè)平均每股負(fù)債數(shù)據(jù)，應(yīng)用KMV模型計算行業(yè)違約概率，解決行業(yè)違約風(fēng)險難以計量的問題。

（2）行業(yè)間違約狀態(tài)聯(lián)合概率測算的思路。將商業(yè)銀行對m個行業(yè)的貸款在到期時區(qū)分為2m種違約狀態(tài)。例如銀行對A 和B兩個行業(yè)貸款，則在貸款到期時，有2m=22=4種違約狀態(tài)，分別為：A、B行業(yè)均不違約、A 行業(yè)違約而B行業(yè)不違約、A 行業(yè)不違約但B行業(yè)違約、A、B行業(yè)均違約。在上述應(yīng)用KMV模型計算行業(yè)違約概率的基礎(chǔ)上，再應(yīng)用正態(tài)Copula函數(shù)計算各違約狀態(tài)的聯(lián)合概率，從而解決行業(yè)間違約狀態(tài)聯(lián)合概率不可直接觀測的問題。

（3）銀行信貸資金收益率分布特征確定的思路。將行業(yè)違約概率與貸款違約損失率的乘積作為對該行業(yè)貸款的信用風(fēng)險溢價。在銀行貸款基準(zhǔn)利率的基礎(chǔ)上再加上行業(yè)貸款信用風(fēng)險溢價，作為對該行業(yè)貸款的貸款利率。在貸款到期時，如行業(yè)不違約，則對該行業(yè)貸款的收益率為貸款利率；如違約，則貸款的收益率為負(fù)且其值等于違約損失率。根據(jù)2m種違約狀態(tài)下整體信貸資金的收益率及各違約狀態(tài)的聯(lián)合概率，確定銀行整體信貸資金收益率波動的分布特征，從而解決銀行整體信貸資金收益率分布特征難以確定的問題。

在上述解決難點思路的基礎(chǔ)上，構(gòu)建基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型，方法為：根據(jù)2m種違約狀態(tài)下整體信貸資金的收益率及各違約狀態(tài)的聯(lián)合概率，計算銀行整體信貸資金收益率的期望及其方差，進一步計算銀行整體信貸資金的差異系數(shù)，即信貸資金獲取單位收益所承擔(dān)的風(fēng)險。將銀行全部信貸資金視為單位1，以對m個行業(yè)的貸款權(quán)重為優(yōu)化變量，以銀行整體信貸資金的差異系數(shù)最小化為目標(biāo)函數(shù)，建立基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型，為商業(yè)銀行信貸資金在行業(yè)間的優(yōu)化配置提供決策參考。

基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型構(gòu)建原理如圖1所示。

圖1 基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型構(gòu)建原理圖

2 基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型

2.1 行業(yè)平均股價與平均每股負(fù)債的計算

2.1.1 行業(yè)平均股價的計算假設(shè)某行業(yè)有n家代表性上市公司。記Pi為第i個代表性上市公司的股票價格，Qi為第i個代表性上市公司的股票數(shù)量，w i為第i個上市公司的市值在n家代表性上市公司總市值中所占的權(quán)重（i=1，2，…，n），則有

根據(jù)行業(yè)內(nèi)n家代表性上市公司的股票價格P i及以市值計算的權(quán)重w i，可計算得到行業(yè)平均股票價格。記P為該行業(yè)的平均股票價格，則有

式（2）的經(jīng)濟學(xué)含義：行業(yè)的平均股票價格等于行業(yè)內(nèi)代表性上市公司的股票價格與其根據(jù)市值計算的權(quán)重的乘積之和。

2.1.2 行業(yè)平均每股負(fù)債的計算 記STDi為行業(yè)內(nèi)第i家代表性上市公司的每股短期負(fù)債，STD為該行業(yè)的平均每股短期負(fù)債，則有

式（3）的經(jīng)濟學(xué)含義：行業(yè)平均每股短期負(fù)債等于行業(yè)內(nèi)代表性上市公司的每股短期負(fù)債與其根據(jù)市值計算的權(quán)重的乘積之和。

類似地，記LTDi為行業(yè)內(nèi)第i家上市公司的每股長期負(fù)債，LTD 為該行業(yè)的平均每股長期負(fù)債，則有

式（4）的經(jīng)濟學(xué)含義：行業(yè)平均每股長期負(fù)債等于行業(yè)內(nèi)代表性上市公司的每股長期負(fù)債與其根據(jù)市值計算的權(quán)重的乘積之和。

2.2 行業(yè)平均每股資產(chǎn)價值及資產(chǎn)收益波動率的計算

2.2.1 行業(yè)平均每股收益波動率的計算 行業(yè)平均每股收益率是指根據(jù)行業(yè)平均股票價格計算的每股收益率。記Pt為某行業(yè)在t時期的平均股票價格，rt為該行業(yè)在t時期的平均每股收益率，則有[13]

式（5）的經(jīng)濟學(xué)含義：t時期的行業(yè)平均每股收益率r t等于t時期的行業(yè)平均股價P t取自然對數(shù)后的1階差分。

記T為考察的行業(yè)平均每股收益率的期間數(shù)為該期間行業(yè)平均每股收益率的均值，σS為某行業(yè)平均每股收益波動率，則有[14]

式（6）的經(jīng)濟學(xué)含義：行業(yè)平均每股收益波動率σS為行業(yè)平均每股收益率r t的標(biāo)準(zhǔn)差。

計算行業(yè)平均每股收益波動率的作用是通過行業(yè)平均股票價格和行業(yè)平均每股收益波動率計算行業(yè)平均每股資產(chǎn)價值及每股資產(chǎn)收益波動率。

2.2.2 計算行業(yè)平均每股資產(chǎn)價值及資產(chǎn)收益波動率的聯(lián)立方程 由于上市公司財務(wù)報表中資產(chǎn)價值核算的是公司資產(chǎn)的歷史成本，并不能準(zhǔn)確地反映公司資產(chǎn)的市場價值。為解決該問題，KMV 模型通過上市公司的股票價格及股票收益波動率來測算公司的資產(chǎn)價值及資產(chǎn)收益波動率。其測算原理是將上市公司股票價值視為以公司資產(chǎn)價值為標(biāo)的資產(chǎn)，以公司負(fù)債為執(zhí)行價格的看漲期權(quán)價值，從而構(gòu)建計算上市公司每股資產(chǎn)價值及資產(chǎn)收益波動率的聯(lián)立方程[13]。

本文將行業(yè)視為行業(yè)內(nèi)代表性上市公司的加總，因此，將行業(yè)平均股票價格替代KMV 模型中的上市公司股票價格，將行業(yè)平均每股負(fù)債替代KMV 模型中的公司每股負(fù)債，應(yīng)用KMV 模型的原理和方法，測算行業(yè)違約概率，從而解決行業(yè)違約概率不可觀測的難題。

記S為行業(yè)平均股票價格，σS為行業(yè)平均每股收益波動率，V為行業(yè)平均每股資產(chǎn)價值，σV為行業(yè)平均每股資產(chǎn)收益波動率，D為行業(yè)平均每股負(fù)債，rf為無風(fēng)險利率，N（·）為標(biāo)準(zhǔn)正態(tài)分布的累積概率分布函數(shù)，則有[13]

式中：

式（7）為Black-Scholes期權(quán)定價公式，其經(jīng)濟學(xué)含義為：行業(yè)平均股票價格S等于行業(yè)平均每股資產(chǎn)價值與以d1為參數(shù)的標(biāo)準(zhǔn)正態(tài)累積概率的乘積VN（d1）減去行業(yè)平均每股負(fù)債D與以d2為參數(shù)的標(biāo)準(zhǔn)正態(tài)累積概率乘積DN（d2）以無風(fēng)險利率rf折現(xiàn)后的差。

為求解式（7）中V及σV這兩個未知變量，KMV模型引入股票收益波動率σS與資產(chǎn)收益波動率σV，兩者之間存在如下函數(shù)關(guān)系[13]：

式（8）中各符號的含義同前。式（7）、（8）聯(lián)立構(gòu)成非線性方程組，解此非線性方程組可計算得到行業(yè)平均每股資產(chǎn)價值及每股資產(chǎn)收益波動率。

2.2.3 求解聯(lián)立方程的算法 利用Newton迭代法可求解式（7）、（8）聯(lián)立構(gòu)成的非線性方程組，為此構(gòu)造函數(shù)分別為[14]：

求解行業(yè)平均每股資產(chǎn)價值V與每股資產(chǎn)收益波動率σV的迭代公式為[14]

式中：V0為行業(yè)平均每股資產(chǎn)價值的給定初值；σV0為每股資產(chǎn)收益波動率的給定初值。應(yīng)用式（11）多次迭代可計算行業(yè)平均每股資產(chǎn)價值V及每股資產(chǎn)收益波動率σV。迭代停止的條件是相鄰兩次迭代結(jié)果的誤差足夠小，如小于10-5。

利用Newton迭代算法求解式（7）、（8）聯(lián)立構(gòu)成的非線性方程組，計算得到行業(yè)平均每股資產(chǎn)價值V及每股資產(chǎn)收益波動率σV，為進一步采用KMV 模型計算行業(yè)違約概率奠定了基礎(chǔ)。

2.3 基于KMV 模型與正態(tài)Copula函數(shù)的行業(yè)違約狀態(tài)聯(lián)合概率的計算

2.3.1 基于KMV 模型的行業(yè)違約概率測算

KMV 模型根據(jù)上市公司資產(chǎn)價值、公司負(fù)債和公司資產(chǎn)收益波動率計算公司違約距離，再通過違約距離計算公司違約概率[13]。

本文將行業(yè)平均每股資產(chǎn)價值替代KMV 模型中的上市公司每股資產(chǎn)價值，將行業(yè)平均每股負(fù)債替代公司每股負(fù)債，將行業(yè)資產(chǎn)收益波動率替代上市公司資產(chǎn)收益波動率，再利用KMV 模型的方法計算行業(yè)違約概率。

記STD 為行業(yè)平均每股短期負(fù)債，LTD 為行業(yè)平均每股長期負(fù)債，D為行業(yè)違約臨界點，KMV模型將違約臨界點D設(shè)定為短期負(fù)債STD 加γ倍的長期負(fù)債LTD[13]，即

式（12）的經(jīng)濟學(xué)含義：行業(yè)違約臨界點D等于行業(yè)平均每股短期負(fù)債STD 加上行業(yè)平均每股長期負(fù)債LTD 乘以折扣系數(shù)γ值后的和。

根據(jù)穆迪KMV 公司的統(tǒng)計數(shù)據(jù)，當(dāng)公司資產(chǎn)價值低于其短期負(fù)債加0.5倍的長期負(fù)債時，公司違約最為頻繁，故取γ=0.5[13]。由于本文將行業(yè)視為行業(yè)內(nèi)代表性上市公司的加總，故對行業(yè)違約臨界點的測算也取γ=0.5。

記V為行業(yè)平均每股資產(chǎn)價值，σV為行業(yè)平均每股資產(chǎn)收益波動率，DD為行業(yè)違約距離。根據(jù)KMV 模型，行業(yè)違約距離[13]

式（13）的經(jīng)濟學(xué)含義：行業(yè)違約距離DD等于行業(yè)平均每股資產(chǎn)價值V減去行業(yè)違約臨界點D的差，再除以行業(yè)平均每股資產(chǎn)價值V與行業(yè)平均每股資產(chǎn)收益波動率σV之積的商。

在假定行業(yè)平均每股資產(chǎn)價值V服從正態(tài)分布的條件下，行業(yè)違約概率為[13]

式中，N（·）為標(biāo)準(zhǔn)正態(tài)分布的累計概率分布函數(shù)。

式（14）的經(jīng)濟學(xué)含義：行業(yè)違約概率PD等于以負(fù)的行業(yè)違約距離DD為參數(shù)的標(biāo)準(zhǔn)正態(tài)分布的累計概率。查標(biāo)準(zhǔn)正態(tài)累積概率分布表，可得到給定參數(shù)-DD的概率值。由式（14）可知，DD與PD呈反向變化關(guān)系，DD越大，PD越小。

根據(jù)行業(yè)平均每股資產(chǎn)價值、行業(yè)平均每股負(fù)債和行業(yè)每股資產(chǎn)收益波動率計算行業(yè)違約距離，再通過行業(yè)違約距離計算行業(yè)違約概率，解決了行業(yè)違約概率不可觀測的問題。

2.3.2 行業(yè)違約距離之間相關(guān)系數(shù)矩陣的計算 在所考察的某段時期內(nèi)取T個時間點，在每個時間點上，應(yīng)用式（13）計算m個行業(yè)的違約距離。

記DD kt為第k個行業(yè)在第t個時點上的違約距離，DD jt為第j個行業(yè)在第t個時點上的違約距離，其中，k=1，2，…，m，j=1，2，…，m，t=1，2，…，T。記ρkj為第k、j個行業(yè)違約距離的相關(guān)系數(shù)，則有

應(yīng)用式（15）可計算出任意兩個行業(yè)違約距離之間的相關(guān)系數(shù)，從而得到m個行業(yè)違約距離之間的相關(guān)系數(shù)矩陣，并記R=[ρkj]m×m為m個行業(yè)違約距離的相關(guān)系數(shù)矩陣。

2.3.3 基于正態(tài)Copula函數(shù)的行業(yè)違約狀態(tài)聯(lián)合概率的計算 由于對每一個行業(yè)的貸款在到期時均有違約和不違約兩種狀態(tài)，故對m個行業(yè)貸款的違約狀態(tài)共有2m種。

對于單個行業(yè)的違約概率可用式（12）～（14）中的KMV 模型計算。

由于KMV 模型采用正態(tài)分布完成從違約距離到違約概率的計算（見式（14）），故本文采用正態(tài)Copula函數(shù)計算m個行業(yè)2m種違約狀態(tài)的聯(lián)合概率。

用s=1，2，…，2m分別標(biāo)記m個行業(yè)的2m種違約狀態(tài)；記JPDs為第s種違約狀態(tài)的聯(lián)合概率，則違約狀態(tài)聯(lián)合概率為[15]

式中：R=[ρkj]m×m為m個行業(yè)違約距離的相關(guān)系數(shù)矩陣；ΦR（·）為以相關(guān)系數(shù)矩陣R為參數(shù)的m維正態(tài)分布的聯(lián)合概率分布函數(shù)；Φ-1（·）為1維標(biāo)準(zhǔn)正態(tài)分布函數(shù)的逆函數(shù)。

由式（14）可知，-DD=Φ-1（PD），故式（16）可變形為

調(diào)用Matlab 軟件中的多元正態(tài)分布函數(shù)mvncdf（）可方便地計算出第s種違約狀態(tài)的聯(lián)合違約概率JPDs。

計算m個行業(yè)違約狀態(tài)聯(lián)合概率的作用是為進一步計算對m個行業(yè)貸款的貸款組合收益率均值和方差奠定基礎(chǔ)。

2.4 基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型

2.4.1 基于信用風(fēng)險溢價的行業(yè)貸款利率的計算 記rk為銀行對第k個行業(yè)的貸款利率，rb為商業(yè)銀行貸款的基準(zhǔn)利率，PD k為第k個行業(yè)根據(jù)KMV 模型計算的行業(yè)違約概率，LGD 為銀行貸款的違約損失率，則有

式中：右端第1項為銀行貸款的基準(zhǔn)利率rb；第2項為第k個行業(yè)的違約概率PD k乘以銀行貸款的違約損失率LGD，稱為銀行對第k個行業(yè)貸款的信用風(fēng)險溢價CS k。

式（18）的經(jīng)濟學(xué)含義：對第k個行業(yè)貸款的利率等于銀行貸款的基準(zhǔn)利率加上對第k個行業(yè)貸款的信用風(fēng)險溢價。根據(jù)式（18）確定行業(yè)貸款利率的好處是能夠使行業(yè)違約風(fēng)險導(dǎo)致的期望損失在貸款利率中得到補償。

2.4.2 貸款組合收益率均值、方差與差異系數(shù)的計算 設(shè)銀行同時對m個行業(yè)貸款，記rk為銀行對第k個行業(yè)貸款的利率，由式（18）計算得到。將銀行整體信貸資金設(shè)為1，并設(shè)w k為銀行對第k個行業(yè)貸款的權(quán)重，顯然，有約束成立。

為簡單起見，假設(shè)對m個行業(yè)貸款的到期日相同，且期限均為1年。用示性函數(shù)I（s，k）表示在第s種違約狀態(tài)下第k個行業(yè)是否違約，如果在第s種違約狀態(tài)下第k個行業(yè)違約，則I（s，k）=1；否則，I（s，k）=0。其中，s=1，2，…，2m，k=1，2，…，m。

記rps為銀行整體信貸資金在第s種違約狀態(tài)下的收益率，則有

如果在第s種違約狀態(tài)下，第k個行業(yè)不違約，則I（s，k）=0，且1-I（s，k）=11，即式（19）右端大括號中起作用的項變?yōu)閞k，即銀行對第k個行業(yè)的貸款到期時獲得r k的收益率。

如果在第s種違約狀態(tài)下，第k個行業(yè)違約，則I（s，k）=1，且1-I（s，k）=0，即式（19）右端大括號中起作用的項變?yōu)?LGD，即銀行對第k個行業(yè)的貸款到期時損失LGD。

綜上所述，式（19）的經(jīng)濟學(xué)含義為：在第s種違約狀態(tài)下，銀行整體信貸資金的收益率r ps等于對m個行業(yè)貸款的收益率rk或損失率-LGD的加權(quán)和。

記E（rp）為銀行整體信貸資金在2m種違約狀態(tài)下的期望收益率，JPDs為第s種違約狀態(tài)的聯(lián)合概率（見式（17）），則有

式（20）的經(jīng)濟學(xué)含義為：銀行整體信貸資金的期望收益率E（rp）等于2m種違約狀態(tài)的聯(lián)合概率JPDs與該違約狀態(tài)下銀行整體信貸資金收益率r ps的乘積和。

記var（rp）為銀行整體信貸資金收益率的方差，則由統(tǒng)計學(xué)中方差的計算公式易得

從而，銀行整體信貸資金收益率的標(biāo)準(zhǔn)差為

記θ為銀行整體信貸資金收益率的差異系數(shù)，由差異系數(shù)的定義[3]，有

式（23）的經(jīng)濟學(xué)含義：銀行整體信貸資金的差異系數(shù)θ等于銀行整體信貸資金收益率的標(biāo)準(zhǔn)差std（rp）除以銀行整體信貸資金的期望收益率E（rp）。銀行整體信貸資金的差異系數(shù)θ反映了銀行通過對m個行業(yè)發(fā)放貸款，獲得單位收益率所承擔(dān)的風(fēng)險大小。差異系數(shù)越大，說明銀行獲得單位收益率所承擔(dān)的風(fēng)險也越大。

2.4.3 基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型 當(dāng)商業(yè)銀行整體信貸資金中對各行業(yè)的貸款權(quán)重w k變化時，銀行整體信貸資金的期望收益率E（rp）及其標(biāo)準(zhǔn)差std（rp）、銀行整體信貸資金的差異系數(shù)θ也相應(yīng)變化。

以銀行整體信貸資金中對m個行業(yè)的貸款權(quán)重w k為決策變量；以最小化銀行整體信貸資金的差異系數(shù)θ為目標(biāo)函數(shù)；以銀行整體信貸資金的期望收益率E（rp）大于或等于目標(biāo)收益率ro、貸款權(quán)重大于或等于0、且貸款權(quán)重之和等于1為約束條件，構(gòu)建銀行信貸資金行業(yè)間優(yōu)化配置模型：

目標(biāo)函數(shù)式（24）通過最小化商業(yè)銀行整體信貸資金的差異系數(shù)θ，能夠使銀行整體信貸資金獲取單位收益率所承擔(dān)的風(fēng)險達到最小。約束條件式（25）使銀行整體信貸資金的期望收益率大于或等于目標(biāo)收益率；式（26）使銀行信貸資金全部貸放出去；式（27）使銀行對任一行業(yè)的貸款為非負(fù)值。

式（24）～（27）可使用Matlab軟件的優(yōu)化工具箱或Excel軟件的規(guī)劃求解等多種軟件求解。通過求解式（24）～（27），可得到對每個行業(yè)貸款的最優(yōu)權(quán)重，為商業(yè)銀行信貸資金在行業(yè)間的優(yōu)化配置提供決策參考。

3 基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型算例

3.1 樣本選取

選取機械設(shè)備業(yè)、建筑業(yè)、金屬非金屬業(yè)、運輸倉儲業(yè)和批發(fā)零售業(yè)5個行業(yè)作為樣本行業(yè)，每個樣本行業(yè)均選取4家代表性樣本公司，共20家樣本公司，如表1所示。

3.2 行業(yè)平均股價與平均每股負(fù)債的計算

3.2.1 行業(yè)平均股價的計算 20家樣本上市公司2011年的股票周收盤價格如表2所示。表2 中樣本上市公司股票周收盤價格來自Wind資訊金融終端。鑒于篇幅，表2僅列出了部分樣本上市公司股票的部分周收盤價格。

表1 樣本行業(yè)與樣本公司名稱

表2 上市公司股票周收盤價格

根據(jù)上市公司市值計算的樣本公司在行業(yè)中的權(quán)重如表3所示。以表3 許繼電器平均股價14.586 1 元/股為例，說明樣本公司平均股價的計算過程。將表2許繼電器51個周的周收盤價格取平均值，可計算得到許繼電器平均股價為14.586 1元/股。表3第3列其他數(shù)據(jù)的計算類推，不贅述。

表3 根據(jù)上市公司市值計算的樣本公司權(quán)重

以表3 許繼電器在機械設(shè)備業(yè)中的權(quán)重等于0.063 0為例，說明樣本公司在行業(yè)中所占權(quán)重的計算過程。將表3機械設(shè)備業(yè)4家樣本公司的平均股價及其股本總數(shù)代入式（1），可計算得到許繼電器在機械設(shè)備業(yè)中的權(quán)重等于0.063 0。表3第5列其他數(shù)據(jù)的計算類推。

根據(jù)行業(yè)內(nèi)代表性上市公司的股價及其權(quán)重，可計算得到行業(yè)平均股價、及行業(yè)平均每股收益率，計算結(jié)果如表4所示。

表4 行業(yè)平均股價及行業(yè)平均每股收益率

以表4中機械設(shè)備業(yè)在2011年1月7日的平均周收盤股價10.365 3元/股為例，說明行業(yè)平均股價的計算過程。將表2中許繼電器、云內(nèi)動力、特變電工和中國北車的周收盤價格，及表3中上述4家上市公司的公司權(quán)重代入式（2），可計算得到機械設(shè)備業(yè)在2011 年1 月7 日的平均周收盤股價為10.365 3元/股。表4第3～7列的行業(yè)平均股價的計算過程類推。

由于每股收益率的計算需要用到本期及上一期的股票價格，故第1周2011年1月7日對應(yīng)的每股收益率沒有計算，即表4第1行第8～12列的數(shù)據(jù)為空。

以表4中機械設(shè)備業(yè)在2011年1月14日的平均每股收益率-0.033 4為例，說明行業(yè)平均每股收益率的計算過程。將表4 中機械設(shè)備業(yè)在2011年1月14日和2011年1月7日的周收盤價格代入式（5），可計算得到機械設(shè)備業(yè)的平均每股收益率為-0.033 4。表4第8～12列的行業(yè)平均每股收益率的計算過程類推。

3.2.2 行業(yè)平均每股負(fù)債的計算20家樣本上市公司的每股短期負(fù)債與每股長期負(fù)債數(shù)據(jù)如表5所示。表5中第3、4列的數(shù)據(jù)來自Wind資訊金融終端。以表5中許繼電器違約臨界點D=9.696 0元/股為例，說明上市公司違約臨界點的計算過程。將表5中許繼電器每股短期負(fù)債、每股長期負(fù)債及參數(shù)γ=0.5[13]代入式（12），計算得到許繼電器違約臨界點D=9.696 0元/股。表5第5列其他數(shù)據(jù)的計算類推。

表5 樣本上市公司的每股負(fù)債

3.3 行業(yè)平均每股資產(chǎn)價值及資產(chǎn)收益波動率的計算

3.3.1 行業(yè)平均每股收益波動率的計算根據(jù)行業(yè)平均股價S、行業(yè)平均每股收益波動率σS求解式（7）、（8），可計算得到行業(yè)每股資產(chǎn)價值V及行業(yè)資產(chǎn)收益波動率σV，計算結(jié)果如表6所示。

以表6 中機械設(shè)備業(yè)加權(quán)股價的均值為8.484 5元/股為例，說明行業(yè)加權(quán)股價均值的計算過程。對表4 中機械設(shè)備業(yè)51 個周的加權(quán)股價取平均值，得到機械設(shè)備業(yè)加權(quán)股價的均值為8.484 5元/股。

表6 樣本行業(yè)違約距離與違約概率

以表6 中機械設(shè)備業(yè)股票收益波動率σS=0.272 1為例，說明行業(yè)股票收益波動率的計算過程。將表4中機械設(shè)備業(yè)股票收益率及其均值、時期數(shù)T=50代入式（6），可計算得到機械設(shè)備業(yè)股票收益波動率σS=0.272 1。

以表6 中機械設(shè)備業(yè)違約臨界點D=7.350 5元/股為例，說明行業(yè)違約臨界點的計算過程。將表5中機械設(shè)備業(yè)代表性上市公司許繼電器、云內(nèi)動力、特變電工和中國北車的違約臨界點分別與表3中上述4家上市公司的權(quán)重相乘，然后相加，計算得到機械設(shè)備業(yè)違約臨界點D=7.350 5元/股。

3.3.2 行業(yè)平均每股資產(chǎn)價值及資產(chǎn)收益波動率的計算通過求解由式（7）、（8）構(gòu)成的聯(lián)立方程組，計算得到行業(yè)每股資產(chǎn)價值和資產(chǎn)收益波動率。

由于表6第6、7列的數(shù)據(jù)通過求解聯(lián)立方程組同時得到，故一并說明其計算過程。以列機械設(shè)備業(yè)每股資產(chǎn)價值V=15.632 0元/股、資產(chǎn)收益波動率σV=0.147 7為例，說明行業(yè)每股資產(chǎn)價值和資產(chǎn)收益波動率的計算過程。

將表6中機械設(shè)備業(yè)的加權(quán)股價均值、股票收益波動率、違約臨界點及無風(fēng)險利率rf=0.028 0[16]分別代入式（7）、（8），利用式（9）～（11）所示的Newton迭代法，計算得到機械設(shè)備業(yè)每股資產(chǎn)價值V=15.632 0元/股、資產(chǎn)收益波動率σV=0.147 7。上述求解過程可利用Matlab軟件求解非線性方程組的fsolve函數(shù)實現(xiàn)。

3.4 基于KMV 模型與正態(tài)Copula函數(shù)的行業(yè)違約狀態(tài)聯(lián)合概率的計算

3.4.1 基于KMV 模型的行業(yè)違約概率測算以表6中機械設(shè)備業(yè)違約距離DD=3.587 4為例，說明行業(yè)違約距離的計算過程。將表6中機械設(shè)備業(yè)的每股資產(chǎn)價值、資產(chǎn)收益波動率及違約臨界點代入式（13），計算得到機械設(shè)備業(yè)違約距離DD=3.587 4。

以表6中機械設(shè)備業(yè)違約概率PD=0.000 167為例，說明行業(yè)違約概率的計算過程。將表6中機械設(shè)備業(yè)違約距離代入式（14），計算得到PD=0.000 167。

3.4.2 行業(yè)違約距離之間相關(guān)系數(shù)矩陣的計算應(yīng)用正態(tài)Copula函數(shù)計算m個行業(yè)2m種違約狀態(tài)的聯(lián)合概率時，需用到m個行業(yè)違約距離之間的相關(guān)系數(shù)矩陣，為此，本文計算了5 個樣本行業(yè)2011年4個季度的違約距離。計算結(jié)果如表7所示。

表7各列數(shù)據(jù)的計算過程與表6相應(yīng)列數(shù)據(jù)的計算過程類似，只是計算的時間區(qū)間由年度變?yōu)榧径?。因此，?各列數(shù)據(jù)的計算過程不贅述。

根據(jù)表7中5 個樣本行業(yè)4 個季度的違約距離，利用式（15）可計算出違約距離的相關(guān)系數(shù)矩陣。但直接計算出的行業(yè)違約距離相關(guān)系數(shù)矩陣近似為奇異矩陣，不能直接應(yīng)用于正態(tài)Copula函數(shù)計算m個行業(yè)2m種違約狀態(tài)的聯(lián)合概率。本文根據(jù)行業(yè)4 個季度的違約距離，首先計算出行業(yè)違約距離的均值與方差，然后根據(jù)行業(yè)違約距離的均值和方差并利用正態(tài)分布，每個行業(yè)隨機生成20個違約距離，再應(yīng)用式（15）計算違約距離的相關(guān)系數(shù)矩陣，可解決直接計算相關(guān)系數(shù)矩陣近似為奇異矩陣的問題。

鑒于篇幅，每個樣本行業(yè)隨機生成的20個違約距離略。5個樣本行業(yè)違約距離之間的相關(guān)系數(shù)矩陣如表8所示。表8 所示的相關(guān)系數(shù)矩陣可利用Excel軟件的correl（）函數(shù)方便得到。

根據(jù)樣本行業(yè)的違約距離及違約距離之間的相關(guān)系數(shù)矩陣，利用正態(tài)Copula函數(shù)可計算得到各違約狀態(tài)的聯(lián)合概率，結(jié)果如表9所示。

表7 樣本行業(yè)季度違約距離

表8 行業(yè)違約距離的相關(guān)系數(shù)矩陣

表9 樣本行業(yè)的違約狀態(tài)及其聯(lián)合概率

3.4.3 樣本行業(yè)的違約狀態(tài)本文算例中共有5個樣本行業(yè)，由于在貸款到期時，每個行業(yè)均存在違約或非違約兩種可能的狀態(tài)，故樣本行業(yè)的違約狀態(tài)共有25=32種，其中，“0”表示非違約，“1”表示違約。

表9第1行第2～6列表示在第1種違約狀態(tài)下，5個樣本行業(yè)的違約狀況，由于表9 第1 行第2～6列的數(shù)據(jù)均為“0”，表示在第1種違約狀態(tài)下，5個樣本行業(yè)均不違約。表9第2行第2～6列為第2種違約狀態(tài)下，5個樣本行業(yè)的違約狀況，由于表9第2行第2列的數(shù)字為“1”，表示在第2種違約狀態(tài)下，機械設(shè)備在貸款到期時違約；而表9第2行第3～6列的數(shù)字為“0”，表示其他4個樣本行業(yè)非違約，即在貸款到期時正常還本付息。

表9其他各行第2～6列表示在對應(yīng)違約狀態(tài)下，5個樣本行業(yè)是否違約的“0”“1”標(biāo)識的含義類推。特別地，在第32種違約狀態(tài)下，表9第32行第2～6列數(shù)字均為“1”，表示在第32種違約狀態(tài)下，5個樣本行業(yè)均違約。限于篇幅，表9中僅列出了部分違約狀態(tài)下各樣本行業(yè)是否違約的標(biāo)識。

3.4.4 各違約狀態(tài)下樣本行業(yè)的違約距離當(dāng)某行業(yè)貸款到期時，根據(jù)KMV 模型，其違約概率PD=N（-DD），即式（14），其中，N（·）為標(biāo)準(zhǔn)正態(tài)分布的累積概率分布函數(shù)。因此，該行業(yè)違約概率對應(yīng)的違約距離為-DD；相應(yīng)地，該行業(yè)不違約的概率為1-PD，由標(biāo)準(zhǔn)正態(tài)分布累積概率分布函數(shù)的性質(zhì)，有1-PD=1-N（-DD）=N（DD）。因此，該行業(yè)非違約概率對應(yīng)的違約距離為DD。

表9第1行第7～11列表示在第1種違約狀態(tài)下，5個樣本行業(yè)違約概率或非違約概率對應(yīng)的違約距離。由表9第1行第2～6列數(shù)據(jù)知，在第1種違約狀態(tài)下，5個樣本行業(yè)均不違約，故表9第1行第7～11列數(shù)據(jù)分別為5個樣本行業(yè)非違約概率對應(yīng)的違約距離DD i，而5個樣本行業(yè)非違約概率對應(yīng)的違約距離DD i來自于表6第1～5行第8列數(shù)據(jù)。

表9第2行第7～11列表示在第2種違約狀態(tài)下，5個樣本行業(yè)違約概率或非違約概率對應(yīng)的違約距離。由表9第2行第2列數(shù)據(jù)知，在第2種違約狀態(tài)下，機械設(shè)備業(yè)違約，故表9第2行第7列數(shù)據(jù)為機械設(shè)備業(yè)違約概率對應(yīng)的違約距離-DD i，而機械設(shè)備業(yè)違約概率對應(yīng)的違約距離來自于表6第1行第8列數(shù)據(jù)取負(fù)號。

由表9第2行第3～6列數(shù)據(jù)可知，在第2種違約狀態(tài)下，除機械設(shè)備業(yè)以外的其余4個樣本行業(yè)非違約，故表9 第2 行第8～11 列數(shù)據(jù)分別為該4個樣本行業(yè)非違約概率對應(yīng)的違約距離DD i，而該4個樣本行業(yè)非違約概率對應(yīng)的違約距離來自于表6第2～5行第8列。

表9其他各行第7～11列數(shù)據(jù)的含義、來源及其計算過程類推。特別地，在第32種違約狀態(tài)下，5個樣本行業(yè)均違約。因此，表9第32行第7～11列數(shù)據(jù)分別為5個樣本行業(yè)違約概率對應(yīng)的違約距離-DD i，而5個樣本行業(yè)違約概率對應(yīng)的違約距離來自于表6第1～5行第8列數(shù)據(jù)取負(fù)號。

限于篇幅，表9中僅列出了部分違約狀態(tài)下各樣本行業(yè)是否違約對應(yīng)的違約距離。

3.4.5 基于正態(tài)Copula函數(shù)的行業(yè)違約狀態(tài)聯(lián)合概率的計算以表9中第1種違約狀態(tài)的聯(lián)合概率JPDs=0.871 709 為例，說明違約狀態(tài)聯(lián)合概率JPDs的計算過程。

將表9中5個樣本行業(yè)的違約距離DD i、表8中5個樣本行業(yè)違約距離的相關(guān)系數(shù)矩陣代入式（17），可計算得到第1 種違約狀態(tài)的聯(lián)合概率JPDs=0.871 709。

3.5 基于違約狀態(tài)聯(lián)合概率的銀行信貸資金行業(yè)間優(yōu)化配置模型應(yīng)用算例

3.5.1 基于信用風(fēng)險溢價的行業(yè)貸款利率的計算表10所示為根據(jù)銀行貸款基準(zhǔn)利率及行業(yè)信用風(fēng)險溢價計算出的行業(yè)貸款利率。

表10 樣本行業(yè)的貸款利率

表10 第3 列為1 年期銀行貸款的基準(zhǔn)利率rb=0.065 6[17]，數(shù)據(jù)來自中國人民銀行網(wǎng)站。第4列為樣本行業(yè)的違約概率，數(shù)據(jù)來自表6第9列。第5列為銀行貸款的違約損失率LGD=0.598 0[18]，數(shù)據(jù)來自穆迪公司的技術(shù)報告。第6列為行業(yè)貸款的信用風(fēng)險溢價CS i，等于表10第4列的行業(yè)違約概率PD i乘以第5列的銀行貸款違約損失率LGD。第7列為行業(yè)貸款利率ri，等于表10第3列的銀行貸款基準(zhǔn)利率rb加上第6列的信用利差或信用風(fēng)險溢價CSi，如式（18）所示。

3.5.2 銀行整體信貸資金收益率均值、方差與差異系數(shù)的計算表11所示為在各違約狀態(tài)下，對5個樣本行業(yè)貸款的收益率或損失率及整體貸款組合的收益率或損失率。

表11第1列為32種違約狀態(tài)的序號，該序號與表9第1列違約狀態(tài)的序號相一致。表11第2～6列為各違約狀態(tài)下，對5個樣本行業(yè)貸款的收益率或損失率。在貸款到期時如行業(yè)不違約，則為收益率，其值等于貸款利率ri；如行業(yè)違約，則為損失率，其值等于貸款違約損失率LGD。

表11第1 行第2～6列表示在第1 種違約狀態(tài)下，5個樣本行業(yè)貸款的收益率。由表9 第1 行第2～6列數(shù)據(jù)可知，在第1種違約狀態(tài)下，5個樣本行業(yè)均不違約。因此，在第1 種違約狀態(tài)下，對5個樣本行業(yè)貸款的收益率即為行業(yè)貸款利率ri，即表11第1行第2～6列的數(shù)據(jù)來自表10第7列。

表11 各違約狀態(tài)下銀行整體信貸資金的收益率

表11第2行第2～6列表示在第2種違約狀態(tài)下，5個樣本行業(yè)貸款的收益率或損失率。由表9第2行第2～6列數(shù)據(jù)可知，在第2種違約狀態(tài)下，機械設(shè)備業(yè)違約，而其余4個樣本行業(yè)不違約。因此，在第2種違約狀態(tài)下，對機械設(shè)備業(yè)的貸款收益率為-0.598 0，亦即損失率為0.598 0，如表11第2行第2列所示，該數(shù)據(jù)來自表10第1行第5列。對其余4個樣本行業(yè)貸款的收益率即為該行業(yè)的貸款利率，如表11第2行第3～6列所示，數(shù)據(jù)分別來自表10第2～5行第7列。

表11其他各行第2～6列數(shù)據(jù)的含義、來源及其計算過程類推。特別地，在第32種違約狀態(tài)下，5個樣本行業(yè)均違約。因此，對5 個樣本行業(yè)貸款的收益率均為-0.598 0，亦即損失率均為0.598 0，如表11第32 行第2～6 列數(shù)據(jù)所示，數(shù)據(jù)來自表10第5列。

限于篇幅，表11中僅列出了部分違約狀態(tài)下各樣本行業(yè)貸款的收益率或損失率。

將銀行整體信貸資金看作單位1，表11 第7列為在5個樣本行業(yè)初始貸款權(quán)重均為0.200 0的條件下，各違約狀態(tài)下銀行整體信貸資金的收益率。5個樣本行業(yè)初始貸款權(quán)重如表12第3列所示。以表11第1行第7列第1 種違約狀態(tài)下，銀行整體信貸資金收益率=0.065 7 為例，說明各違約狀態(tài)下銀行整體信貸資金收益率rps的計算過程。

將表11第1行第2～6列5個樣本行業(yè)貸款收益率ri分別乘以對5個樣本行業(yè)貸款的權(quán)重w i，然后相加，可計算得到第1種違約狀態(tài)下銀行整體信貸資金的收益率為

將上述計算結(jié)果列于表11第1行第7列。表11第7列其他數(shù)據(jù)的計算過程類推。表11第8列為各違約狀態(tài)的聯(lián)合概率JPDs，數(shù)據(jù)來自表9 第12列。

表12所示為將銀行整體信貸資金看作單位1的情況下，5個樣本行業(yè)的初始貸款權(quán)重與最優(yōu)貸款權(quán)重。

表12 樣本行業(yè)的初始貸款權(quán)重與最優(yōu)貸款權(quán)重

其中，表12 第3 列的初始貸款權(quán)重為主觀設(shè)定，但初始權(quán)重的設(shè)置并不影響最優(yōu)貸款權(quán)重的求解結(jié)果。第4 列的最優(yōu)貸款權(quán)重，通過求解式（24）～（27）得到。

根據(jù)各違約狀態(tài)下銀行整體信貸資金的收益率及各違約狀態(tài)的聯(lián)合概率，可計算得到銀行整體信貸資金收益率的均值、方差及其差異系數(shù)，計算結(jié)果如表13所示。

表13 銀行整體信貸資金收益率均值、方差與差異系數(shù)

表13第1行第3列表示在5個樣本行業(yè)貸款權(quán)重均為0.200 0的情況下，銀行整體信貸資金收益率的均值。其計算過程為：將表11中各違約狀態(tài)下的銀行整體信貸資金收益率r ps及其聯(lián)合概率JPDs代入式（20），可計算得到銀行整體信貸資金收益率的均值E（rp）=0.045 8。

表13第2行第3列為在5個樣本行業(yè)貸款權(quán)重均為0.200 0的情況下，銀行整體信貸資金收益率的標(biāo)準(zhǔn)差。其計算過程為：將表11中各違約狀態(tài)下的銀行整體信貸資金收益率rps及其聯(lián)合概率JPDs，以及表13中銀行整體信貸資金收益率的均值E（rp）代入式（21），可計算得到銀行整體信貸資金收益率的方差，對其開平方得到銀行整體信貸資金收益率的標(biāo)準(zhǔn)差std（rp）=0.053 8。

表13中第3行第3列為等權(quán)重情況下銀行整體信貸資金的差異系數(shù)θ。其計算過程為：將表13中第2行第3列等權(quán)重情況下銀行整體信貸資金收益率的標(biāo)準(zhǔn)差std（rp）=0.053 8、第1行第3列等權(quán)重情況下銀行整體信貸資金收益率的均值E（rp）=0.045 8代入式（23），計算得到等權(quán)重情況下銀行整體信貸資金的差異系數(shù)θ=1.173 5，列于表13第3行第3列。

3.5.3 基于違約狀態(tài)聯(lián)合概率的銀行信貸資金行業(yè)間配置最優(yōu)貸款權(quán)重的計算表12給出了銀行信貸資金在行業(yè)間優(yōu)化配置的行業(yè)最優(yōu)貸款權(quán)重。為簡單起見，設(shè)初始權(quán)重均為0.200 0，初始權(quán)重的設(shè)置并不影響最優(yōu)貸款權(quán)重的求解結(jié)果。

根據(jù)初始貸款權(quán)重計算的各違約狀態(tài)下銀行整體信貸資金的收益率r ps列于表11第7列；根據(jù)初始貸款權(quán)重計算的銀行整體信貸資金的期望收益率E（rp）、銀行整體信貸資金收益率的標(biāo)準(zhǔn)差std（r p）及銀行整體信貸資金的差異系數(shù)θ分別列于表13第1～3行第3列。

當(dāng)貸款權(quán)重發(fā)生變化時，各違約狀態(tài)下銀行整體信貸資金收益率、銀行整體信貸資金收益率的均值和標(biāo)準(zhǔn)差及銀行整體信貸資金的差異系數(shù)也相應(yīng)發(fā)生變化。

在各樣本行業(yè)貸款權(quán)重之和為1，且貸款權(quán)重均大于或等于0的約束條件下，總存在一組貸款權(quán)重，使得銀行整體信貸資金的差異系數(shù)達到最小，即銀行整體信貸資金獲得單位收益率所承擔(dān)的風(fēng)險最小，該組貸款權(quán)重即為最優(yōu)貸款權(quán)重。

上述最優(yōu)貸款權(quán)重求解問題實質(zhì)是一個非線性數(shù)學(xué)規(guī)劃問題（見式（24）～（27））。該非線性數(shù)學(xué)規(guī)劃問題以銀行整體信貸資金中樣本行業(yè)的貸款權(quán)重為決策變量；以銀行整體信貸資金的期望收益率大于或等于目標(biāo)收益率、貸款權(quán)重大于或等于0、且貸款權(quán)重之和等于1為約束條件；以最小化商業(yè)銀行整體信貸資金的差異系數(shù)θ為目標(biāo)函數(shù)，來確定對各樣本行業(yè)的最優(yōu)貸款權(quán)重，從而建立銀行整體信貸資金行業(yè)間最優(yōu)配置決策方案。

上述數(shù)學(xué)規(guī)劃問題可利用Excel軟件中的規(guī)劃求解功能方便地得到對各樣本行業(yè)的最優(yōu)貸款權(quán)重。

3.5.4 最優(yōu)權(quán)重與等權(quán)重下銀行整體信貸資金收益率與風(fēng)險的比較根據(jù)各樣本行業(yè)的最優(yōu)貸款權(quán)重，可重新計算各違約狀態(tài)下銀行整體信貸資金的收益率，從而進一步計算最優(yōu)貸款權(quán)重下銀行整體信貸資金收益率的均值和標(biāo)準(zhǔn)差，及最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金的差異系數(shù)。

最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金收益率的均值E（r p）、標(biāo)準(zhǔn)差std（rp）及最優(yōu)貸款組合的差異系數(shù)θ分別列于表13第1～3行第4列。其計算過程與初始等權(quán)重情況下銀行整體信貸資金收益率均值、標(biāo)準(zhǔn)差及差異系數(shù)的計算過程類似，即與表13第1～3行第3列數(shù)據(jù)的計算過程類似。

比較表13第1行第3、4列數(shù)據(jù)可知，初始等權(quán)重情況下銀行整體信貸資金的平均收益率為0.045 8，而最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金的平均收益率為0.054 0，最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金的平均收益率高于初始等權(quán)重下銀行整體信貸資金的平均收益率。

比較表13第2行第3、4列數(shù)據(jù)可知，初始等權(quán)重情況下銀行整體信貸資金收益率的標(biāo)準(zhǔn)差為0.053 8，而最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金收益率的標(biāo)準(zhǔn)差為0.046 6，最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金收益率的標(biāo)準(zhǔn)差小于初始等權(quán)重下銀行整體信貸資金收益率的標(biāo)準(zhǔn)差，亦即最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金的風(fēng)險低于初始等權(quán)重下銀行整體信貸資金的風(fēng)險。

比較表13第3行第3、4列數(shù)據(jù)可知，初始等權(quán)重情況下銀行整體信貸資金的差異系數(shù)為1.173 5，而最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金的差異系數(shù)為0.863 3，最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金的差異系數(shù)小于初始等權(quán)重下銀行整體信貸資金的差異系數(shù)，亦即最優(yōu)行業(yè)貸款權(quán)重下銀行整體信貸資金獲得單位收益率所承擔(dān)的風(fēng)險小于初始等權(quán)重下銀行整體信貸資金獲得單位收益率所承擔(dān)的風(fēng)險。

4 結(jié)語

本文以行業(yè)內(nèi)代表性上市公司市值加權(quán)的股票價格作為行業(yè)平均股價，以市值加權(quán)的每股負(fù)債作為行業(yè)平均每股負(fù)債，根據(jù)行業(yè)平均股價與行業(yè)平均每股負(fù)債，應(yīng)用KMV 模型計算行業(yè)違約概率，進而將商業(yè)銀行對m個行業(yè)的貸款在到期時區(qū)分為2m種違約狀態(tài)，應(yīng)用正態(tài)Copula函數(shù)計算各違約狀態(tài)的聯(lián)合概率。根據(jù)2m種違約狀態(tài)下信貸資金的收益率及各違約狀態(tài)的聯(lián)合概率，計算銀行整體信貸資金的差異系數(shù)，即信貸資金獲取單位收益所承擔(dān)的風(fēng)險。以對m個行業(yè)的貸款權(quán)重為優(yōu)化變量，以銀行整體信貸資金的差異系數(shù)最小化為目標(biāo)函數(shù)，建立了基于違約狀態(tài)聯(lián)合概率的商業(yè)銀行信貸資金行業(yè)間優(yōu)化配置模型，為商業(yè)銀行信貸資金在行業(yè)間的優(yōu)化配置提供決策參考。