电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

Copula理论及MATLAB应用实例

Copula理论及MATLAB应用实例_第1页
1/6
Copula理论及MATLAB应用实例_第2页
2/6
Copula理论及MATLAB应用实例_第3页
3/6
%-------------------------------------------------------------------------- % 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 进行正态性检验 [h,p] = jbtest(X) % Jarque-Bera 检验 [h,p] = kstest(X,[X,normcdf(X,mean(X),std(X))]) % Kolmogorov-Smirnov 检验 [h, p] = lillietest(X) % Lilliefors 检验 % 分别调用jbtest、kstest 和 lillietest 函数对 Y 进行正态性检验 [h,p] = jbtest(Y) % Jarqu e-Bera 检验 [h,p] = kstest(Y,[Y,normcdf(Y,mean(Y),std(Y))]) % Kolmogorov -Smirnov 检验 [h, p] = lillietest(Y) % Lilliefors 检验 %****************************求经验分布函数值******************************* % 调用ecdf 函数求 X 和 Y 的经验分布函数 [fx , Xsort] = ecdf(X); [fy , Ysort] = ecdf(Y); % 调用spline 函数,利用样条插值法求原始样本点处的经验分布函数值 U1...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

Copula理论及MATLAB应用实例

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部