clc,clear;
%% Parameters related to households
gamm=2;%跨期替代弹性的相反数（风险厌恶系数）
bet=0.985;%主观贴现因子
varpi=4.060;%劳动负效用系数
varphi=1.455;%劳动曲率参数
%% Parameters related to open economies
iot=1.5;%本国商品与进口商品之间的替代弹性
nu=0.20;%开放水平
B_Y=0.15;%外债占GDP的稳态比率
psi=0.05;%国别利率溢价参数
rho_Rstar=0.99;%国外利率冲击持续性参数
sigma_Rstar=0.06;%国外利率冲击标准差
%% Parameters related to banks
chi=0.5;%金融摩擦程度
%theta=0.475;%可挪用资产的份额
%kxi=8.16/1000;%家庭转移到银行的份额
sigm=0.96;%银行生存率
%% Other parameters
alph=0.33;%有效资本份额
eta=1.5;%资本价格对投资的弹性
delta=0.025;%资本折旧率
G_Y=0.3;%政府支出占GDP的稳态比率
Sk=0.0005;%风险溢价
Lev=6.5;%银行杠杆
taub=0.01;
tauk=0.01;
rho_Z=0.90;%技术冲击持续性参数
sigma_Z=0.05;%技术冲击标准差
%% Steady State
Z_ss=1;%(1)
Lambda_ss=bet;%(2)
R_ss=1/bet;%(3)
Q_ss=1;%(4)
q_ss=1;%外生假定
gq_ss=1;%(5)
e_ss=1;%(6)
CA_ss=0;%(7)
Rk_ss=R_ss+Sk;%外生假定（8）
Y_K=Q_ss*(Rk_ss-1+delta)/(gq_ss^(-1)*alph); %gq=Ph_ss/P_ss
I_K=delta;
I_Y=I_K*(Y_K)^(-1);
Rb_ss=R_ss-chi*(Rk_ss-R_ss);%外生假定（40）
TB_Y=(Rb_ss-1)*B_Y;%外生假定
K_L=(Y_K)^(-1/(1-alph));
L_ss=((1-alph)/chi*gq_ss^(-1)*(K_L)^alph)^(1/(varphi-1));%gq_ss=P_ss/Ph_ss（10）
K_ss=(K_L)*L_ss;%（11）
S_ss=K_ss;%（12）
Y_ss=Z_ss*K_ss^(alph)*L_ss^(1-alph);%（13）
I_ss=(I_Y)*Y_ss;%（14）
C_Y=1/gq_ss-I_Y-G_Y-(0.5*eta*(I_ss/I_ss-1)^(2)*I_ss)/Y_ss-TB_Y;
B_ss=(B_Y)*Y_ss;%（15）
G_ss=(G_Y)*Y_ss;%（16）
Th_ss=G_ss;%（17）
C_ss=(C_Y)*Y_ss;%（18）
Ch_ss=(1-nu)*(1/gq_ss)^(-iot)*C_ss;%（19）
Cf_ss=nu*e_ss^(-iot)*C_ss;%（20）
TB_ss=(TB_Y)*Y_ss;%（21）
EX_ss=q_ss^(-1)*(Y_ss-(1-nu)*gq_ss^(iot)*(C_ss+I_ss+G_ss+0.5*eta*(I_ss/I_ss-1)^(2)*I_ss));%（22）
varrho_ss=(C_ss-chi/varphi*L_ss^(varphi))^(-gamm);%（23）
N_ss=K_ss*(Lev)^(-1);%Lev=K_N=3 银行杠杆（24）
Rbstar_ss=Rb_ss;%（25）
Rstar_ss=Rb_ss;%（26）
phia_ss=Q_ss*S_ss/N_ss + chi*e_ss*B_ss/N_ss;%（27）
D_ss=Q_ss*S_ss-e_ss*B_ss-N_ss;%（28）
Ne_ss=sigm*(Rk_ss*Q_ss*S_ss-Rb_ss*e_ss*B_ss-R_ss*D_ss);%（29）
kxi=1/(Rk_ss*Q_ss*S_ss)*(N_ss+sigm*Rb_ss*e_ss*B_ss+sigm*R_ss*D_ss-sigm*Rk_ss*Q_ss*S_ss);
Nn_ss=kxi*Rk_ss*Q_ss*S_ss;%（30）
phib_ss=phia_ss/chi;%（31）
Omegaa_ss=(1-sigm)/(1-sigm*(phia_ss*Lambda_ss*(Rk_ss-R_ss)+Lambda_ss*R_ss));%（32）
Nu_ss=Lambda_ss*Omegaa_ss*R_ss;%（33）
mu_ss=Lambda_ss*Omegaa_ss*(Rk_ss-R_ss);%（34）
mub_ss=Lambda_ss*Omegaa_ss*(R_ss-Rb_ss);%（35）
theta=mu_ss+Nu_ss/phia_ss;
%w_ss=(1-alph)*gq_ss^(-1)*(Y_K)*(K_L);%（36）
w_ss=chi*L_ss^(varphi-1);
Taub_ss=taub*(log((e_ss*B_ss)/(e_ss*B_ss)));%（37）
Tauk_ss=tauk*(log((Q_ss*S_ss)/(Q_ss*S_ss)));%（38）
Zeta_ss=Tauk_ss*Rk_ss*Q_ss*S_ss+Taub_ss*Rb_ss*e_ss*B_ss;%（39）
Nus_ss=mu_ss*Q_ss;
Nub_ss=(Nu_ss-mub_ss)*e_ss;
save _param;
%% dyanre program
dynare Copy_of_soe_bank

% model_diagnostics(M_, options_, oo_)

