- idx = idxys.idxdre(idxys.year == 2015 & idxys.season == 3); %提取年份t、季度k的股指日收益率
- [~,n] = size(stk2015_3);
- A = zeros(2,2493);
- for i = 1:5
- stk = stk2015_3(:,i);
-
- u = ecdf(idx); %设定股指经验分布函数
- v = ecdf(stk); %设定个股经验分布函数
- T = length(u); %设定T
-
- options = optimset('Display','iter','TolCon',10^-12,'TolFun',10^-4,'TolX',10^-6);
-
- lower = [0,0]; %设定估计的上界与下界
- upper = [1,1];
- theta0 = [0.25;0.25]; %设定起始估计值
- [kappa,LL] = fmincon('sym_jc_CL',theta0,[],[],[],[],lower,upper,[],options,[u,v]); %估计SJC copula的两个参数
-
- A(:,i) = kappa; %A是一个空矩阵,多次迭代合并以保存输出结果
-
- end
执行代码后matlab报错:
错误使用horzcat
要串联的数组维度不一致
debug了一下,发现问题出现在 [kappa,LL]=... 这一语句上,求助高手解决