fprintf('\n======================Robust_non-manufacturing======================');
fprintf('\n======================Robust_non-manufacturing======================');
fprintf('\n======================Robust_non-manufacturing======================');
fprintf('\n==================================================\n\n\n\n\n\n\n\n\n');


%% Load Parameters
clear;  clc;  close

load("alpha7\data\Baseline.mat");
Case1_alpha7 = Case1;
Case2_alpha7 = Case2;
Case3_alpha7 = Case3;
Case4_alpha7 = Case4;

load("alpha8\data\Baseline.mat");
Case1_alpha8 = Case1;
Case2_alpha8 = Case2;
Case3_alpha8 = Case3;
Case4_alpha8 = Case4;

load("alpha9\data\Baseline.mat");
Case1_alpha9 = Case1;
Case2_alpha9 = Case2;
Case3_alpha9 = Case3;
Case4_alpha9 = Case4;

wel = [Case1_alpha7.wel_c; Case1_alpha7.wel_u; ...
       Case2_alpha7.wel_c; Case2_alpha7.wel_u; ...
       Case3_alpha7.wel_c; Case3_alpha7.wel_u; ...
       Case4_alpha7.wel_c; Case4_alpha7.wel_u; ...
       Case1_alpha8.wel_c; Case1_alpha8.wel_u; ...
       Case2_alpha8.wel_c; Case2_alpha8.wel_u; ...
       Case3_alpha8.wel_c; Case3_alpha8.wel_u; ...
       Case4_alpha8.wel_c; Case4_alpha8.wel_u; ...
       Case1_alpha9.wel_c; Case1_alpha9.wel_u; ...
       Case2_alpha9.wel_c; Case2_alpha9.wel_u; ...
       Case3_alpha9.wel_c; Case3_alpha9.wel_u; ...
       Case4_alpha9.wel_c; Case4_alpha9.wel_u];
Pi = [Case1_alpha7.Pi_c; Case1_alpha7.Pi_u; ...
      Case2_alpha7.Pi_c; Case2_alpha7.Pi_u; ...
      Case3_alpha7.Pi_c; Case3_alpha7.Pi_u; ...
      Case4_alpha7.Pi_c; Case4_alpha7.Pi_u; ...
      Case1_alpha8.Pi_c; Case1_alpha8.Pi_u; ...
      Case2_alpha8.Pi_c; Case2_alpha8.Pi_u; ...
      Case3_alpha8.Pi_c; Case3_alpha8.Pi_u; ...
      Case4_alpha8.Pi_c; Case4_alpha8.Pi_u; ...
      Case1_alpha9.Pi_c; Case1_alpha9.Pi_u; ...
      Case2_alpha9.Pi_c; Case2_alpha9.Pi_u; ...
      Case3_alpha9.Pi_c; Case3_alpha9.Pi_u; ...
      Case4_alpha9.Pi_c; Case4_alpha9.Pi_u];
M = [Case1_alpha7.M_c; Case1_alpha7.M_u; ...
     Case2_alpha7.M_c; Case2_alpha7.M_u; ...
     Case3_alpha7.M_c; Case3_alpha7.M_u; ...
     Case4_alpha7.M_c; Case4_alpha7.M_u; ...
     Case1_alpha8.M_c; Case1_alpha8.M_u; ...
     Case2_alpha8.M_c; Case2_alpha8.M_u; ...
     Case3_alpha8.M_c; Case3_alpha8.M_u; ...
     Case4_alpha8.M_c; Case4_alpha8.M_u; ...
     Case1_alpha9.M_c; Case1_alpha9.M_u; ...
     Case2_alpha9.M_c; Case2_alpha9.M_u; ...
     Case3_alpha9.M_c; Case3_alpha9.M_u; ...
     Case4_alpha9.M_c; Case4_alpha9.M_u];
w = [Case1_alpha7.w_c; Case1_alpha7.w_u; ...
     Case2_alpha7.w_c; Case2_alpha7.w_u; ...
     Case3_alpha7.w_c; Case3_alpha7.w_u; ...
     Case4_alpha7.w_c; Case4_alpha7.w_u; ...
     Case1_alpha8.w_c; Case1_alpha8.w_u; ...
     Case2_alpha8.w_c; Case2_alpha8.w_u; ...
     Case3_alpha8.w_c; Case3_alpha8.w_u; ...
     Case4_alpha8.w_c; Case4_alpha8.w_u; ...
     Case1_alpha9.w_c; Case1_alpha9.w_u; ...
     Case2_alpha9.w_c; Case2_alpha9.w_u; ...
     Case3_alpha9.w_c; Case3_alpha9.w_u; ...
     Case4_alpha9.w_c; Case4_alpha9.w_u];
P = [Case1_alpha7.P_c; Case1_alpha7.P_u; ...
     Case2_alpha7.P_c; Case2_alpha7.P_u; ...
     Case3_alpha7.P_c; Case3_alpha7.P_u; ...
     Case4_alpha7.P_c; Case4_alpha7.P_u; ...
     Case1_alpha8.P_c; Case1_alpha8.P_u; ...
     Case2_alpha8.P_c; Case2_alpha8.P_u; ...
     Case3_alpha8.P_c; Case3_alpha8.P_u; ...
     Case4_alpha8.P_c; Case4_alpha8.P_u; ...
     Case1_alpha9.P_c; Case1_alpha9.P_u; ...
     Case2_alpha9.P_c; Case2_alpha9.P_u; ...
     Case3_alpha9.P_c; Case3_alpha9.P_u; ...
     Case4_alpha9.P_c; Case4_alpha9.P_u];

result = reshape([round(wel,2),round(Pi,2),round(M,2),round(w,2),round(P,2)],[24,5]);
Table1 = array2table(result, 'VariableNames', {'福利', '商品利润', '数据利润', '工资', '价格指数'});
writetable(Table1, '表XI1.xlsx', 'WriteVariableNames', true);

fprintf('\n=================================================================\n');
fprintf('         表XI1 不同数据中间商市场结构以及不同数据政策下的福利\n');
fprintf('=================================================================\n\n');
disp(Table1);

%% Replicate Section 3 of the Paper
clear;  clc;  close
load("alpha7\data\Counterfactual_TariffUnilateral_UC.mat");
Case1_alpha7_UC = Case1;
Case2_alpha7_UC = Case2;
Case3_alpha7_UC = Case3;
Case4_alpha7_UC = Case4;

load("alpha7\data\Counterfactual_TariffUnilateral_CU.mat");
Case1_alpha7_CU = Case1;
Case2_alpha7_CU = Case2;
Case3_alpha7_CU = Case3;
Case4_alpha7_CU = Case4;

load("alpha8\data\Counterfactual_TariffUnilateral_UC.mat");
Case1_alpha8_UC = Case1;
Case2_alpha8_UC = Case2;
Case3_alpha8_UC = Case3;
Case4_alpha8_UC = Case4;

load("alpha8\data\Counterfactual_TariffUnilateral_CU.mat");
Case1_alpha8_CU = Case1;
Case2_alpha8_CU = Case2;
Case3_alpha8_CU = Case3;
Case4_alpha8_CU = Case4;

load("alpha9\data\Counterfactual_TariffUnilateral_UC.mat");
Case1_alpha9_UC = Case1;
Case2_alpha9_UC = Case2;
Case3_alpha9_UC = Case3;
Case4_alpha9_UC = Case4;

load("alpha9\data\Counterfactual_TariffUnilateral_CU.mat");
Case1_alpha9_CU = Case1;
Case2_alpha9_CU = Case2;
Case3_alpha9_CU = Case3;
Case4_alpha9_CU = Case4;

opt_alpha7 = [Case1_alpha7_CU.t_cu_opt; Case1_alpha7_UC.t_uc_opt; ...
       Case2_alpha7_CU.t_cu_opt; Case2_alpha7_UC.t_uc_opt; ...
       Case3_alpha7_CU.t_cu_opt; Case3_alpha7_UC.t_uc_opt; ...
       Case4_alpha7_CU.t_cu_opt; Case4_alpha7_UC.t_uc_opt];
wel_alpha7 = [Case1_alpha7_CU.wel_c_opt; Case1_alpha7_UC.wel_u_opt; ...
       Case2_alpha7_CU.wel_c_opt; Case2_alpha7_UC.wel_u_opt; ...
       Case3_alpha7_CU.wel_c_opt; Case3_alpha7_UC.wel_u_opt; ...
       Case4_alpha7_CU.wel_c_opt; Case4_alpha7_UC.wel_u_opt];

opt_alpha8 = [Case1_alpha8_CU.t_cu_opt; Case1_alpha8_UC.t_uc_opt; ...
       Case2_alpha8_CU.t_cu_opt; Case2_alpha8_UC.t_uc_opt; ...
       Case3_alpha8_CU.t_cu_opt; Case3_alpha8_UC.t_uc_opt; ...
       Case4_alpha8_CU.t_cu_opt; Case4_alpha8_UC.t_uc_opt];
wel_alpha8 = [Case1_alpha8_CU.wel_c_opt; Case1_alpha8_UC.wel_u_opt; ...
       Case2_alpha8_CU.wel_c_opt; Case2_alpha8_UC.wel_u_opt; ...
       Case3_alpha8_CU.wel_c_opt; Case3_alpha8_UC.wel_u_opt; ...
       Case4_alpha8_CU.wel_c_opt; Case4_alpha8_UC.wel_u_opt];

opt_alpha9 = [Case1_alpha9_CU.t_cu_opt; Case1_alpha9_UC.t_uc_opt; ...
       Case2_alpha9_CU.t_cu_opt; Case2_alpha9_UC.t_uc_opt; ...
       Case3_alpha9_CU.t_cu_opt; Case3_alpha9_UC.t_uc_opt; ...
       Case4_alpha9_CU.t_cu_opt; Case4_alpha9_UC.t_uc_opt];
wel_alpha9 = [Case1_alpha9_CU.wel_c_opt; Case1_alpha9_UC.wel_u_opt; ...
       Case2_alpha9_CU.wel_c_opt; Case2_alpha9_UC.wel_u_opt; ...
       Case3_alpha9_CU.wel_c_opt; Case3_alpha9_UC.wel_u_opt; ...
       Case4_alpha9_CU.wel_c_opt; Case4_alpha9_UC.wel_u_opt];

result = reshape([round(opt_alpha7,2),round(wel_alpha7,2),round(opt_alpha8,2),round(wel_alpha8,2),round(opt_alpha9,2),round(wel_alpha9,2)],[8,6]);
Table2 = array2table(result, 'VariableNames', {'关税_alpha7', '福利_alpha7', '关税_alpha8', '福利_alpha8', '关税_alpha9', '福利_alpha9'});
writetable(Table2, '表XI2.xlsx', 'WriteVariableNames', true);

fprintf('\n=================================================================\n');
fprintf('         表XI2 稳健性检验：不同原始数据权重下的最优关税与福利\n');
fprintf('=================================================================\n\n');
disp(Table2);

%% Table 3：纳什关税和福利
clear;  clc;  close
load("alpha7\data\Counterfactual_TariffNash.mat");
Case1_alpha7_Nash = Case1;
Case2_alpha7_Nash = Case2;
Case3_alpha7_Nash = Case3;
Case4_alpha7_Nash = Case4;

load("alpha8\data\Counterfactual_TariffNash.mat");
Case1_alpha8_Nash = Case1;
Case2_alpha8_Nash = Case2;
Case3_alpha8_Nash = Case3;
Case4_alpha8_Nash = Case4;

load("alpha9\data\Counterfactual_TariffNash.mat");
Case1_alpha9_Nash = Case1;
Case2_alpha9_Nash = Case2;
Case3_alpha9_Nash = Case3;
Case4_alpha9_Nash = Case4;


opt_alpha7 = [Case1_alpha7_Nash.t_cu_opt; Case1_alpha7_Nash.t_uc_opt; ...
       Case2_alpha7_Nash.t_cu_opt; Case2_alpha7_Nash.t_uc_opt; ...
       Case3_alpha7_Nash.t_cu_opt; Case3_alpha7_Nash.t_uc_opt; ...
       Case4_alpha7_Nash.t_cu_opt; Case4_alpha7_Nash.t_uc_opt];
wel_alpha7 = [Case1_alpha7_Nash.wel_c; Case1_alpha7_Nash.wel_u; ...
       Case2_alpha7_Nash.wel_c; Case2_alpha7_Nash.wel_u; ...
       Case3_alpha7_Nash.wel_c; Case3_alpha7_Nash.wel_u; ...
       Case4_alpha7_Nash.wel_c; Case4_alpha7_Nash.wel_u];

opt_alpha8 = [Case1_alpha8_Nash.t_cu_opt; Case1_alpha8_Nash.t_uc_opt; ...
       Case2_alpha8_Nash.t_cu_opt; Case2_alpha8_Nash.t_uc_opt; ...
       Case3_alpha8_Nash.t_cu_opt; Case3_alpha8_Nash.t_uc_opt; ...
       Case4_alpha8_Nash.t_cu_opt; Case4_alpha8_Nash.t_uc_opt];
wel_alpha8 = [Case1_alpha8_Nash.wel_c; Case1_alpha8_Nash.wel_u; ...
       Case2_alpha8_Nash.wel_c; Case2_alpha8_Nash.wel_u; ...
       Case3_alpha8_Nash.wel_c; Case3_alpha8_Nash.wel_u; ...
       Case4_alpha8_Nash.wel_c; Case4_alpha8_Nash.wel_u];

opt_alpha9 = [Case1_alpha8_Nash.t_cu_opt; Case1_alpha8_Nash.t_uc_opt; ...
       Case2_alpha8_Nash.t_cu_opt; Case2_alpha8_Nash.t_uc_opt; ...
       Case3_alpha8_Nash.t_cu_opt; Case3_alpha8_Nash.t_uc_opt; ...
       Case4_alpha8_Nash.t_cu_opt; Case4_alpha8_Nash.t_uc_opt];
wel_alpha9 = [Case1_alpha8_Nash.wel_c; Case1_alpha8_Nash.wel_u; ...
       Case2_alpha8_Nash.wel_c; Case2_alpha8_Nash.wel_u; ...
       Case3_alpha8_Nash.wel_c; Case3_alpha8_Nash.wel_u; ...
       Case4_alpha8_Nash.wel_c; Case4_alpha8_Nash.wel_u];


result = reshape([round(opt_alpha7,2),round(wel_alpha7,2),round(opt_alpha8,2),round(wel_alpha8,2),round(opt_alpha9,2),round(wel_alpha9,2)],[8,6]);
Table3 = array2table(result, 'VariableNames', {'关税_alpha7', '福利_alpha7', '关税_alpha8', '福利_alpha8', '关税_alpha9', '福利_alpha9'});
writetable(Table3, '表XI3.xlsx', 'WriteVariableNames', true);

fprintf('\n=================================================================\n');
fprintf('         表XI3 不同数据中间商市场结构以及不同数据政策下的纳什关税与福利\n');
fprintf('=================================================================\n\n');
disp(Table3);
