%函数设定
function F=kmv01(EtoD,r,T,EquityTheta,x)
d1=(log(x(1)*EtoD)+(r+0.5*x(2)^2)*T)/(x(2)*sqrt(T));
d2=d1-x(2)*sqrt(T);
F=[x(1)*normcdf(d1)-exp(-r*T)*normcdf(d2)/EtoD-1;normcdf(d1)*x(1)*x(2)-EquityTheta];
% 注 EquityTheta 企业股权价值波动率年化
--------------------------------------------------------------------------------------------------------------
function [Va,AssetTheta]=kmv02(E,D,r,T,EquityTheta)
EtoD=E/D;
x0=[1,1];
VaThetaX=fsolve(@(x) kmv01(EtoD,r,T,EquityTheta,x),x0);
AssetTheta=VaThetaX(2);
Va=VaThetaX(1)*E;
% 注:Va企业资产价值;AssetTheta企业资产价值的波动率
--------------------------------------------------------------------------------------------------------------
%条件设定
clc
clear all
% 同济堂(600090)
r=0.015;% 央行一年期定期存款利率
T=1;
SD=21.41e8; %短期债务
LD=25.87e8;%长期债务
E=69e8;%12个月的平均股权价值
priceTheta=0.2913;% 年化股权价值波动率=8.41%*sqrt(12)
EquityTheta=priceTheta;
DP=SD+0.5*LD;
D=DP;%违约点
%解方程组的
[Va,AssetTheta]=kmv02(E,D,r,T,EquityTheta);
DD=(Va-DP)/(Va*AssetTheta);%违约距离计算
----------------------------------------------------------------------------------------------------------------- %计算结果
disp(['市场价值',num2str(Va)])
disp(['资产波动率',num2str(AssetTheta)])
disp(['违约距离',num2str(DD)])
- 1
- 2
前往页