%--------------------------------------------------------------------------
%?????????????????????????Copula理論及應用實例
%--------------------------------------------------------------------------

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


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


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

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


%******************************正態性檢驗***********************************
%?分別調用jbtest、kstest和lillietest函數對X進行正態性檢驗
[hp]?=?jbtest（X）??%?Jarque-Bera檢驗
[hp]?=?kstest（X[Xnormcdf（Xmean（X）std（X））]）??%?Kolmogorov-Smirnov檢驗
[h?p]?=?lillietest（X）??%?Lilliefors檢驗

%?分別調用jbtest、kstest和lillietest函數對Y進行正態性檢驗
[hp]?=?jbtest（Y）??%?Jarque-Bera檢驗
[hp]?=?kstest（Y[Ynormcdf（Ymean（Y）std（Y））]）??%?Kolmogorov-Smirnov檢驗
[h?p]?=?lillietest（Y）??%?Lilliefors檢驗


%****************************求經驗分布函數值*******************************
%?調用ecdf函數求X和Y的經驗分布函數
[fx?Xsort]?=?ecdf（X）;
[fy?Ysort]?=?ecdf（Y）;
%?調用spline函數，利用樣條插值法求原始樣本點處的經驗分布函數值
U1?=?spline（Xsort（2:end）fx（2:end）X）;
V1?=?spline（Ysort（2:end）fy（2:end）Y）;

%?調用ecdf函數求X和Y的經驗分布函數
[fx?Xsort]?=?ecdf（X）;
[fy?Ysort]?=?ecdf（Y）;
%?提取fx和fy的第2個至最后一個元素，即排序后樣本點處的經驗分布函數值
fx?=?fx（2:end）;
fy?=?fy（2:end）;

%?通過排序和反排序恢復原始樣本點處的經驗分布函數值U1和V1
[Xsortid]?=?sort（X）;
[idsortid]?=?sort（id）;
U1?=?fx（id）;
[Ysortid]?=?sort（Y）;
[idsortid]?=?sort（id）;
V1?=?fy（id）;


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


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

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


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

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