diary on ; % 使用当前 diary 日志文件名启用记录。
disp(['*********开始时间:',datestr(now,'yyyy-mm-dd HH:MM:SS'),'**********************************']);
*********开始时间:2025-07-08 15:11:37**********************************
%% ------1.中国单边降低外商投资运营壁垒--------- %%
clear
clc
load Data ;

% 反事实的初始参数设置
N          = 32;
vbarrier   = 0.9;  % gamma_hat 外商壁垒降低幅度, 1(不调整) ,0.9(降低10%)
barrier    = 0 ;   % 3(不调整) or 0(中国对其他国家) or 1(其他国家对中国) or -1(双边) or 2(多边)
preference = 0 ;   % 0   or 0.9 :beta_hat is 1(不调整) or 0.9(降低偏向性)
SOEratio   = 0 ;    
alpha_hat  = 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(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(m,data,par,N);
Result_barrier_OpenofCHN   = Result;
CPTPP_barrier_OpenofCHN    = fun_CPTPP(Result_barrier_OpenofCHN);

writematrix(Result_barrier_OpenofCHN,'Barrier.xls','Sheet',2 ,'Range','A1:H32');
writematrix(CPTPP_barrier_OpenofCHN ,'Barrier.xls','Sheet',2 ,'Range','I1:P32');
%% ------2.中国仅消除针对外资的政策歧视--------- %%
clear
clc
load Data ; 

% 反事实的初始参数设置
N          = 32;
vbarrier   = 1;    % gamma_hat 外商壁垒降低幅度, 1(不调整) ,0.9(降低10%)
barrier    = 3;    % 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   or 1.1 :alpha_hat is 1(不调整) or 1.1(提高商业考虑) 
par        = fun_parament(N,barrier,vbarrier,preference,alpha_hat,SOEratio);
par.zeta   = 1 ;   % 仅针对外资

% 反事实求解(请按说明文件要求运行)
x0         = ones(N,1);
[w_hat,M_hat,K_hat] = fun_Iteration_prefer(x0,data,par,N);

zeta =

     1

m.w_hat    = w_hat;
m.K_hat    = K_hat;
m.M_hat    = M_hat;
[Result,X_iln,spsi_lambda]     = get_result_prefer(m,data,par,N);
Result_PolicyPrefer_asymmetric = Result;
CPTPP_PolicyPrefer_asymmetric  = fun_CPTPP(Result_PolicyPrefer_asymmetric);

writematrix(Result_PolicyPrefer_asymmetric,'PolicyPrefer.xls','Sheet',2 ,'Range','A1:H32');
[警告: 已添加指定的工作表。] 
writematrix(CPTPP_PolicyPrefer_asymmetric ,'PolicyPrefer.xls','Sheet',2 ,'Range','I1:P32');
%% ------3.全面对标的特殊情况（缺失特定方面）--------- %%
%--a.缺少投资政策对标
clear
clc
load Data ; 

% 反事实的初始参数设置
N          = 32;
vbarrier   = 1;    % gamma_hat 外商壁垒降低幅度, 1(不调整) ,0.9(降低10%)
barrier    = 3;    % 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(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(m,data,par,N);
Result_partial_nogamma     = Result;
CPTPP_partial_nogamma      = fun_CPTPP(Result_partial_nogamma);

writematrix(Result_partial_nogamma,'Partial.xls','Sheet',1 ,'Range','A1:H32');
writematrix(CPTPP_partial_nogamma ,'Partial.xls','Sheet',1 ,'Range','I1:P32');

%--b.缺少偏向性政策对标
clear
clc
load Data ; 

% 反事实的初始参数设置
N          = 32;
vbarrier   = 0.9;  % gamma_hat 外商壁垒降低幅度, 1(不调整) ,0.9(降低10%)
barrier    = -1 ;  % 3(不调整) or 0(中国对其他国家) or 1(其他国家对中国) or -1(双边) or 2(多边)
preference = 0 ;   % 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(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(m,data,par,N);
Result_partial_noprefer    = Result;
CPTPP_partial_noprefer     = fun_CPTPP(Result_partial_noprefer);

writematrix(Result_partial_noprefer,'Partial.xls','Sheet',2 ,'Range','A1:H32');
[警告: 已添加指定的工作表。] 
writematrix(CPTPP_partial_noprefer ,'Partial.xls','Sheet',2 ,'Range','I1:P32');

%--c.缺少国企商业考虑对标
clear
clc
load Data ; 

% 反事实的初始参数设置(请按说明文件要求运行)
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   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(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(m,data,par,N);
Result_partial_nodistort   = Result;
CPTPP_partial_nodistort    = fun_CPTPP(Result_partial_nodistort);

writematrix(Result_partial_nodistort,'Partial.xls','Sheet',3 ,'Range','A1:H32');
[警告: 已添加指定的工作表。] 
writematrix(CPTPP_partial_nodistort ,'Partial.xls','Sheet',3 ,'Range','I1:P32');

%% 
disp(['*********结束时间:',datestr(now,'yyyy-mm-dd HH:MM:SS'),'**********************************']);
*********结束时间:2025-07-08 15:19:03**********************************
diary off; %结束运行
