%%clear;clc;%獲得購買美團(tuán)股票、京東股票、百度股票的月均盈利率數(shù)據(jù)priceuprate=[-1.46 +4.06 +58.60 -11.00 -10.08 -15.90 +16.23 -7.00 -15.00 +5.00 +8.00 -5.00 +14.56 -8.00;13.20 0.29 -1.29 -0.21 -0.64 -20.67 -1.03, 2.90 17.45 -8.35 -7.98  1.55  48.15  2.31;-7.01 +27.40 -13.35 -6.77 -0.87 +7.46 -4.58 +6.45 -4.57 -4.57 +2.84 +4.70 +2.25 +38.27]/100;r_f = 0.002;%無風(fēng)險利率priceuprate = priceuprate';size_priceuprate = size(priceuprate);average_priceuprate = sum(priceuprate)/size_priceuprate(1);%取平均利率risk_matrix = cov(priceuprate);%取方差-協(xié)方差矩陣%%num_dots = 5000;%設(shè)置二次規(guī)劃點(diǎn)數(shù)target_priceuprates = linspace(min(average_priceuprate) ,...max(average_priceuprate),num_dots);%設(shè)定固定利率weights = zeros(size_priceuprate(2) , num_dots);risks = zeros(1 , num_dots);Aeq = [ones(1 , size_priceuprate(2)) ; average_priceuprate];lb = zeros(size_priceuprate(2) , 1);%設(shè)置下界(不允許賣空)ub = ones(size_priceuprate(2) , 1);%設(shè)置上界(不允許其他為負(fù))H = 2 * risk_matrix;p = zeros(size_priceuprate(2) , 1);%給出求解函數(shù)options = optimoptions('quadprog','Algorithm','interior-point-convex',...'Display','iter','MaxIterations',200);%給出二次規(guī)劃設(shè)置%進(jìn)行二次規(guī)劃for i = 1 : num_dotstarget_priceuprate = target_priceuprates(i);beq = [1 ; target_priceuprate];[weight , risk_min , exitflag , output , lambda]=...quadprog(H , p , [] , [] , Aeq , beq , lb , ub , [] , options);if exitflag > 0weights(: , i) = weight;risks(1 , i) =sqrt( weight'*risk_matrix*weight );elserisk(1 , i) = NaN;end%檢驗(yàn)規(guī)劃結(jié)果并保存end%%excess_priceuprate = average_priceuprate - r_f;weight_tan =risk_matrix \ excess_priceuprate';weight_tan = weight_tan / sum(weight_tan);risk_tan = sqrt(weight_tan' * risk_matrix *weight_tan);priceuprate_tan = average_priceuprate * weight_tan;k_tan = (priceuprate_tan - r_f) / (risk_tan);%求解資本市場線權(quán)重、風(fēng)險以及利率，并據(jù)此求解斜率options = optimoptions('quadprog','Algorithm','interior-point-convex'...,'Display','iter','MaxIterations',200);beq = [1 , priceuprate_tan];[weight_tan1 , risk_min , exitflag , output , lambda]=...quadprog(H , p , [] , [] , Aeq , beq , lb , ub , [] , options);%再次求解權(quán)重value1 = 0;for i = 1 : size_priceuprate(2)if weight_tan(i,1) ~= weight_tan1(i,1)weight_tan = weight_tan1;break;endend%保證非負(fù)結(jié)果%%% 可視化結(jié)果figure('Position', [100, 100, 1200, 500]);% 子圖1：有效前沿和資本市場線subplot(1,2,1);plot(risks , target_priceuprates, 'b-', 'LineWidth', 2);hold on;% 繪制資本市場線cml_risks = linspace(0, max(risks)*1.2, 100);cml_priceuprate = r_f + k_tan * cml_risks;plot(cml_risks, cml_priceuprate, 'r--', 'LineWidth', 2);% 標(biāo)記切點(diǎn)plot(risk_tan, priceuprate_tan, 'ro', 'MarkerSize', 8, 'MarkerFaceColor', 'red');% 標(biāo)記無風(fēng)險資產(chǎn)plot(0, r_f, 'go', 'MarkerSize', 8, 'MarkerFaceColor', 'green');xlabel('風(fēng)險（標(biāo)準(zhǔn)差）');ylabel('預(yù)期收益率');title('有效前沿和資本市場線');legend('有效前沿', '資本市場線', '切點(diǎn)組合', '無風(fēng)險資產(chǎn)', 'Location', 'southeast');grid on;% 子圖2：切點(diǎn)組合權(quán)重subplot(1,2,2);bar(weight_tan);xlabel('資產(chǎn)');ylabel('權(quán)重');title('切點(diǎn)組合權(quán)重分配');set(gca, 'XTickLabel', {'美團(tuán)股票', '京東股票', '百度股票'});grid on;

上面這段代碼在北太天元中會一直運(yùn)行但無結(jié)果輸出，但在Matlab中可以正常運(yùn)行并輸出結(jié)果