diary on ; % 使用当前 diary 日志文件名启用记录。
disp(['*********开始时间:',datestr(now,'yyyy-mm-dd HH:MM:SS'),'**********************************']);
*********开始时间:2025-07-08 21:07:11**********************************

%% ------全面对标--------- %%
clear
clc
load Data_Frechet ; 

% 反事实的初始参数设置
N          = 32;
vbarrier   = 0.9;  % gamma_hat 外商壁垒降低幅度, 1(不调整) ,0.9(降低10%）
barrier    = -1 ;  % 3(不调整) or 0(中国对其他国家) or 1(其他国家对中国) or -1(双边) or 2(多边) 
preference = 0.9 ; % 0   or 0.9 :beta_hat is 1(不调整) or 0.9(降低偏向性)
SOEratio   = 0 ;   % 
alpha_hat  = 1.1 ; % 1   or 1.1 :alpha_hat is 1(不调整) or 1.1(提高商业考虑) 
par        = fun_parament(N,barrier,vbarrier,preference,alpha_hat,SOEratio);

% 反事实求解(请按说明文件要求运行)
x0         = ones(N,1);
[w_hat,M_hat,K_hat] = fun_Iteration_Frechet(x0,data,par,N);
m.w_hat    = w_hat;
m.K_hat    = K_hat;
m.M_hat    = M_hat;
[Result,X_iln,spsi_lambda] = get_result_Frechet(m,data,par,N);
Result_benchmarking        = Result;
CPTPP_benchmarking         = fun_CPTPP(Result_benchmarking);

writematrix(Result_benchmarking,'Benchmarking.xls','Sheet',1 ,'Range','A1:H32');
writematrix(CPTPP_benchmarking ,'Benchmarking.xls','Sheet',1 ,'Range','I1:P32');
writematrix(spsi_lambda,'Result_Benchmarking_decompose.xls','Sheet',1);
%%
disp(['*********结束时间:',datestr(now,'yyyy-mm-dd HH:MM:SS'),'**********************************']);
*********结束时间:2025-07-08 21:08:52**********************************
diary off; %结束运行
