stepan
/
quest-iii-estimation
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Use strcmp instead of strmatch.

master
Stéphane Adjemian (Argos) 2023-10-12 17:47:26 +02:00
parent 9c06a79bb6
commit c4161c4605
Signed by: stepan
GPG Key ID: A6D44CB9C64CE77B
1 changed files with 10 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
Q3_steadystate.m
View File

@ -25,9 +25,9 @@ GYW0 = GY0;
% DELTAKGE = (DELTAGE+GPOP0)/(1+GPOP0);
DELTAKE = (DELTAE+GPOP0);
DELTAKGE = (DELTAGE+GPOP0);
M_.params(strmatch('E_EX_RW',M_.param_names,'exact')) = E_EX_RW;
M_.params(strmatch('GYW0',M_.param_names,'exact')) = GYW0;
M_.params(strmatch('GTFP0',M_.param_names,'exact')) = GTFP0;
M_.params(strcmp('E_EX_RW', M_.param_names)) = E_EX_RW;
M_.params(strcmp('GYW0', M_.param_names)) = GYW0;
M_.params(strcmp('GTFP0', M_.param_names)) = GTFP0;
% end;
%SE=0.5*(SESUM+SEDIFF);
SWE=SE;
@ -73,7 +73,7 @@ E_BWRY = (E_EX_RW - E_EX_R)/RPREME;
TBYN = (-E_INOM+GPOP0+GP0+GY0)*E_BWRY;
E_TBYN=TBYN;
E_LYWY = LYWY0;
M_.params(strmatch('LYWY0',M_.param_names,'exact')) = LYWY0;
M_.params(strcmp('LYWY0', M_.param_names)) = LYWY0;
YWY = exp(LYWY0);
% E_LIK = log(DELTAKE+GY0+GPCPI0)-(GY0+GPCPI0);
@ -89,7 +89,7 @@ else
YKN = (E_EX_R+RPREMK+DELTAE)/(1-TP)/(1-TAUE)/(1-ALPHAE);
E_LYKPPI=log(YKN);
ISN=IK/YKN;
M_.params(strmatch('ISN',M_.param_names,'exact')) = ISN;
M_.params(strcmp('ISN', M_.param_names)) = ISN;
assignin('base','ISN',ISN);
end
@ -169,7 +169,7 @@ E_ZEPS_IG = 0;
E_LGSN = log(GSN);
E_LISN = log(ISN);
CSN = 1-(IGSN+GSN+ISN+TBYN);
M_.params(strmatch('CSN',M_.param_names,'exact')) = CSN;
M_.params(strcmp('CSN', M_.param_names)) = CSN;
E_LCSN = log(CSN);
E_LBGYN=log(BGTAR);
YWR = 1/((1-TAUE)*ALPHAE/L0*(1+E_LOL)) ;
@ -180,7 +180,7 @@ E_TW = 0.2;
% TW0=E_TW/(E_WPTU)^TW1;
M_.params(strmatch('TW0',M_.param_names,'exact')) = TW0;
M_.params(strcmp('TW0', M_.param_names)) = TW0;
E_TRW = TRSN;
E_WS = exp(E_LL-E_LYWR);
E_TAXYN = (E_INOM-GP0-GY0-GPOP0)*BGTAR+IGSN+GSN+TRSN*E_WS ...
@ -272,13 +272,13 @@ end
A1E = (1-TAUE)*(1-ALPHAE)*exp(E_LYKPPI)/UCAP0;
assignin('base','A1E',A1E);
M_.params(strmatch('A1E',M_.param_names,'exact')) = A1E;
M_.params(strcmp('A1E', M_.param_names)) = A1E;
M_.params(strmatch('OMEGE',M_.param_names,'exact')) = OMEGE;
M_.params(strcmp('OMEGE', M_.param_names)) = OMEGE;
assignin('base','OMEGE',OMEGE);
if DDYN==0
RPREMK = (1-TAUE)*(1-ALPHAE)*(1-TP)*exp(E_LYKPPI) - DELTAE-E_EX_R;
M_.params(strmatch('RPREMK',M_.param_names,'exact')) = RPREMK;
M_.params(strcmp('RPREMK', M_.param_names)) = RPREMK;
assignin('base','RPREMK',RPREMK);
end