• 大小: 6.81MB
    文件类型: .zip
    金币: 1
    下载: 0 次
    发布日期: 2023-08-12
  • 语言: Matlab
  • 标签: 数学建模  matlab  

资源简介

copula 理论及应用 Frim算法求最小生成树问题 Kruskal算法求最小生成树问题 参数估计与检验 层次分析法 非线性方程组求解 二次指数平滑法 加权移动平滑法 解线性方程组的迭代法 利用MATLAB生成Word和Excel文档 聚类分析程序 矩阵特征值计算 模拟退火算法 神经网络算法 图像处理程序 主成分分析程序最短路径算法等、

资源截图

代码片段和文件信息

%--------------------------------------------------------------------------
%                         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个文件信息

评论

共有 条评论