%--------------------------------------------------------------------------
%?????????????????????????Copula理論及應(yīng)用實(shí)例
%--------------------------------------------------------------------------

%******************************讀取數(shù)據(jù)*************************************
%?從文件hushi.xls中讀取數(shù)據(jù)
hushi?=?xlsread（‘hushi.xls‘）;
%?提取矩陣hushi的第5列數(shù)據(jù)，即滬市的日收益率數(shù)據(jù)
X?=?hushi（:5）;
%?從文件shenshi.xls中讀取數(shù)據(jù)
shenshi?=?xlsread（‘shenshi.xls‘）;
%?提取矩陣shenshi的第5列數(shù)據(jù)，即深市的日收益率數(shù)據(jù)
Y?=?shenshi（:5）;


%****************************繪制頻率直方圖*********************************
%?調(diào)用ecdf函數(shù)和ecdfhist函數(shù)繪制滬、深兩市日收益率的頻率直方圖
[fx?xc]?=?ecdf（X）;
figure;
ecdfhist（fx?xc?30）;
xlabel（‘滬市日收益率‘）;??%?為X軸加標(biāo)簽
ylabel（‘f（x）‘）;??%?為Y軸加標(biāo)簽
[fy?yc]?=?ecdf（Y）;
figure;
ecdfhist（fy?yc?30）;
xlabel（‘深市日收益率‘）;??%?為X軸加標(biāo)簽
ylabel（‘f（y）‘）;??%?為Y軸加標(biāo)簽


%****************************計(jì)算偏度和峰度*********************************
%?計(jì)算X和Y的偏度
xs?=?skewness（X）
ys?=?skewness（Y）

%?計(jì)算X和Y的峰度
kx?=?kurtosis（X）
ky?=?kurtosis（Y）


%******************************正態(tài)性檢驗(yàn)***********************************
%?分別調(diào)用jbtest、kstest和lillietest函數(shù)對(duì)X進(jìn)行正態(tài)性檢驗(yàn)
[hp]?=?jbtest（X）??%?Jarque-Bera檢驗(yàn)
[hp]?=?kstest（X[Xnormcdf（Xmean（X）std（X））]）??%?Kolmogorov-Smirnov檢驗(yàn)
[h?p]?=?lillietest（X）??%?Lilliefors檢驗(yàn)

%?分別調(diào)用jbtest、kstest和lillietest函數(shù)對(duì)Y進(jìn)行正態(tài)性檢驗(yàn)
[hp]?=?jbtest（Y）??%?Jarque-Bera檢驗(yàn)
[hp]?=?kstest（Y[Ynormcdf（Ymean（Y）std（Y））]）??%?Kolmogorov-Smirnov檢驗(yàn)
[h?p]?=?lillietest（Y）??%?Lilliefors檢驗(yàn)


%****************************求經(jīng)驗(yàn)分布函數(shù)值*******************************
%?調(diào)用ecdf函數(shù)求X和Y的經(jīng)驗(yàn)分布函數(shù)
[fx?Xsort]?=?ecdf（X）;
[fy?Ysort]?=?ecdf（Y）;
%?調(diào)用spline函數(shù)，利用樣條插值法求原始樣本點(diǎn)處的經(jīng)驗(yàn)分布函數(shù)值
U1?=?spline（Xsort（2:end）fx（2:end）X）;
V1?=?spline（Ysort（2:end）fy（2:end）Y）;

%?調(diào)用ecdf函數(shù)求X和Y的經(jīng)驗(yàn)分布函數(shù)
[fx?Xsort]?=?ecdf（X）;
[fy?Ysort]?=?ecdf（Y）;
%?提取fx和fy的第2個(gè)至最后一個(gè)元素，即排序后樣本點(diǎn)處的經(jīng)驗(yàn)分布函數(shù)值
fx?=?fx（2:end）;
fy?=?fy（2:end）;

%?通過(guò)排序和反排序恢復(fù)原始樣本點(diǎn)處的經(jīng)驗(yàn)分布函數(shù)值U1和V1
[Xsortid]?=?sort（X）;
[idsortid]?=?sort（id）;
U1?=?fx（id）;
[Ysortid]?=?sort（Y）;
[idsortid]?=?sort（id）;
V1?=?fy（id）;


%*******************************核分布估計(jì)**********************************
%?調(diào)用ksdensity函數(shù)分別計(jì)算原始樣本X和Y處的核分布估計(jì)值
U2?=?ksdensity（XX‘function‘‘cdf‘）;
V2?=?ksdensity（YY‘function‘‘cdf‘）;


%?**********************繪制經(jīng)驗(yàn)分布函數(shù)圖和核分布估計(jì)圖**********************
[Xsortid]?=?sort（X）;??%?為了作圖的需要，對(duì)X進(jìn)行排序
figure;??%?新建一個(gè)圖形窗口
plot（XsortU1（id）‘c‘‘LineWidth‘5）;?%?繪制滬市日收益率的經(jīng)驗(yàn)分布函數(shù)圖
hold?on
plot（XsortU2（id）‘k-.‘‘LineWidth‘2）;?%?繪制滬市日收益率的核分布估計(jì)圖
legend（‘經(jīng)驗(yàn)分布函數(shù)‘‘核分布估計(jì)‘?‘Location‘‘NorthWest‘）;?%?加標(biāo)注框
xlabel（‘滬市日收益率‘）;??%?為X軸加標(biāo)簽
ylabel（‘F（x）‘）;??%?為Y軸加標(biāo)簽

[Ysortid]?=?sort（Y）;??%?為了作圖的需要，對(duì)Y進(jìn)行排序
figure;??%?新建一個(gè)圖形窗口
plot（YsortV1（id）‘c‘‘LineWidth‘5）;?%?繪制深市日收益率的經(jīng)驗(yàn)分布函數(shù)圖
hold?on
plot（YsortV2（id）‘k-.‘‘LineWidth‘2）;?%?繪制深市日收益率的核分布估計(jì)圖
legend（‘經(jīng)驗(yàn)分布函數(shù)‘‘核分布估計(jì)‘?‘Location‘‘NorthWest‘）;?%?加標(biāo)注框
xlabel（‘深市日收益率‘）;??%?為X軸加標(biāo)簽
ylabel（‘F（x）‘）;??%?為Y軸加標(biāo)簽


%****************************繪制二元頻數(shù)直方圖*****************************
%?調(diào)用ksdensity函數(shù)分別計(jì)算原始樣本X和Y處的核分布估計(jì)值
U?=?ksdensi

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\Copula理論及應(yīng)用\
?????文件????????8161??2010-05-13?01:39??matlab程序附集\Copula理論及應(yīng)用\example06_01.m
?????文件??????365056??2009-11-08?12:14??matlab程序附集\Copula理論及應(yīng)用\hushi.xls
?????文件??????365056??2009-11-08?12:15??matlab程序附集\Copula理論及應(yīng)用\shenshi.xls
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\Frim算法求最小生成樹問(wèn)題\
?????文件?????????461??2013-08-22?22:35??matlab程序附集\Frim算法求最小生成樹問(wèn)題\prim.m
?????文件???????11287??2013-08-22?22:39??matlab程序附集\Frim算法求最小生成樹問(wèn)題\最小生成樹問(wèn)題.jpg
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\Kruskal算法求最小生成樹問(wèn)題\
?????文件?????????580??2013-08-22?22:41??matlab程序附集\Kruskal算法求最小生成樹問(wèn)題\Kruskal.m
?????文件???????11287??2013-08-22?22:39??matlab程序附集\Kruskal算法求最小生成樹問(wèn)題\最小生成樹問(wèn)題.jpg
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\一次指數(shù)平滑法\
?????文件?????????388??2013-08-19?20:34??matlab程序附集\一次指數(shù)平滑法\yczsphf.m
?????文件?????????437??2013-08-19?20:32??matlab程序附集\一次指數(shù)平滑法\一次指數(shù)平滑法.txt
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\三次指數(shù)平滑法\
?????文件?????????871??2013-08-19?20:45??matlab程序附集\三次指數(shù)平滑法\sczsphf.m
?????文件?????????422??2013-08-19?20:44??matlab程序附集\三次指數(shù)平滑法\三次指數(shù)平滑法.txt
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\主成分分析程序\
?????文件???????19456??2004-10-05?10:47??matlab程序附集\主成分分析程序\examp11_02.xls
?????文件????????1315??2010-01-22?08:15??matlab程序附集\主成分分析程序\example11_01.m
?????文件????????3923??2010-01-21?08:26??matlab程序附集\主成分分析程序\example11_02.m
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\二次指數(shù)平滑法\
?????文件?????????494??2013-08-19?20:39??matlab程序附集\二次指數(shù)平滑法\eczsphf.m
?????文件?????????437??2013-08-19?20:32??matlab程序附集\二次指數(shù)平滑法\一次指數(shù)平滑法.txt
?????目錄???????????0??2018-01-09?11:13??matlab程序附集\偏微分方程的數(shù)值解法\
?????文件??????????48??2007-10-23?19:09??matlab程序附集\偏微分方程的數(shù)值解法\DKIniU.m
?????文件??????????56??2007-10-23?20:38??matlab程序附集\偏微分方程的數(shù)值解法\EllIni2Uxl.m
?????文件??????????60??2007-10-23?20:38??matlab程序附集\偏微分方程的數(shù)值解法\EllIni2Uxr.m
?????文件??????????54??2007-10-23?20:38??matlab程序附集\偏微分方程的數(shù)值解法\EllIni2Uyl.m
?????文件??????????99??2007-10-23?20:48??matlab程序附集\偏微分方程的數(shù)值解法\EllIni2Uyr.m
?????文件?????????253??2007-09-11?19:55??matlab程序附集\偏微分方程的數(shù)值解法\Funval.m
............此處省略556個(gè)文件信息