Contributions to Game Theory and Management, XVI, 61—73
Analysis and Control of Macroeconomic Trends Based on
the Leontief Model
Zhao Guo and Dan Wang
St.Petersburg State University, 7/9, Universitetskaya nab., St. Petersburg, 199034, Russia E-mail: 13709267195@163.com, wd789564@163.com
Abstract In the present study, we consolidate the forty-five industrial sectors delineated in the U.S. input-output tables, as disseminated by the OECD, into three overarching sectors: advanced manufacturing, modern services, and a residual category termed 'others.' We adopt gross fixed capital formation as a proxy for the proportion of net profit allocated to investment, positing that the investment requisite for augmenting output is commensurate with the requisite capital intensity. This framework enables us to forecast the trajectory of total output and GDP, taking into account the interplay of multiple determinants. In addressing the inherent linear control dynamics of the input-output model, we apply classical control theory to regulate the advanced manufacturing sector. By deriving control equations that accommodate multifactorial influences, we substantiate the efficacy of this control mechanism through rigorous numerical analysis. Moreover, we reconceptualize the dynamic input-output system as a game-theoretic model characterized by a saddle-point equilibrium. By leveraging the saddle-point equilibrium strategy, we pioneer an innovative approach to resolving the complexities of dynamic input-output analysis. This methodological innovation not only enhances the precision of our predictions but also contributes a novel perspective to the literature on economic modeling and control theory.
Keywords: dynamic input-output model, program control, saddle point equilibrium strategy, differential game.
1. Introduction
The application of static input-output models and the theory of dynamic input-output models are the main focuses of current input-output research. Few researchers have employed dynamic input-output models in conjunction with game theory to investigate the optimal control of an unbalanced economy.
In the era of the digital economy, it is crucial to consider the time component in the input-output table and to be able to control the parameters to achieve the desired goals. In particular, in the case of unbalanced input-output, providing optimal control under the influence of multiple factors can help enterprises to optimize their benefits. Similarly, the government can alter tax policy based on the predicted result of the input-output table to guide resource allocation and encourage the development of each sector.
Leontief (1956) proposed the dynamic inverse model, laying the foundation of dynamic input-output modeling. Miller (2009) included the time factor in the input-output analysis. Smirnov (2021) considered government consumption as the fourth quadrant of the input-output table in the input-output equilibrium equation, and defined total output as the derivative of time. In this paper, based on (Miller, 2009) and (Smirnov, 2021), we remove Taxes fewer subsidies on intermediate and final https://doi.org/10.21638/11701/spbu31.2023.04
domestic products, which correspond directly to each sector, as the fourth quadrant according to the input-output table published by OECD. The time series is utilized to forecast current economic trends in the United States.
Given the linear control characteristic in the dynamic input-output model, we combine the classical program control theory proposed by Tamasyan (2008) and based on the viewpoint of the advanced manufacturing industry, i.e., the industry that utilizes the emerging technology as the fundamental means, proposed by Singhry (2016), Mourtzis et al. (2018) and Jin et al. (2017), we classify what sectors are included in the advanced manufacturing industry. The equations for the control of multifactor influence are obtained after applying complete control to the advanced manufacturing industry in the U.S. The effectiveness of the program control is verified by numerical computation.
Qu and Huang (1999) described optimal production strategies for consumption tracking in dynamic input-output systems. Kang et al. (1992) discussed the dynamic input-output optimal control model. And Mao (1992) indicated the dynamic input-output optimal control model with constraints. Based on the dynamic input-output model proposed by Leontief, we study the optimal policy design problem of the continuous dynamic input-output model. The dynamic input-output system is abstracted as a saddle-point equilibrium game model in the current paper using ideas from Qu, Kang, and Mao. The saddle-point equilibrium strategy is used to design a new method for solving the dynamic input-output problem, which provides a basis for decision-making by macroeconomic policymakers. The current paper also uses the saddle-point equilibrium strategy to design a new method for solving the dynamic input-output problem, which provides a basis for decision-making by macroeconomic policymakers.
2. Dynamic Input-output Modeling of a Balanced Economy 2.1. Relationship between input-output table
The input-output tables published by the OECD cover 45 economic sectors, and all values are represented in current dollars. For the input-output table, the rows upward represent the destination of the allocated uses of the product of a sector, the sum of which equals total output Xj; the columns upward represent the various inputs to the production of the product, and the sum of the values of these inputs is the total inputs, i.e., Xj. Each element of the matrix Q = {xj}j=1 represents, from the row-wise view, the amount of i product allocated to the production of j product, and from the column-wise direction, the amount of i product consumed in the production of j product. The row vector Ej represents the value added in each j sector. It contains labor wages Wj , net production taxes ONTj, and net profits Prhj. TF denotes Taxes less subsidies on intermediate and final imported products (TXS — IMP — FNL). TD denotes Taxes fewer subsidies on intermediate and final domestic products (TXS — INT — FNL). The column vector Y denotes the final expenditure on output in the industrial sector, which includes Gross Fixed Capital Formation(GFCF), Changes in inventories(III), Direct purchases abroad by residents (imports), Direct purchases abroad, Exports (cross border), and Imports (cross border).
Fig. 1 has the following equilibrium relations.
Industry Final expenditure
Industry Q Y X
TXS_IMP_FNL TF
TXS_INT_FNL TD
E
X
Fig. 1. Schematic diagram of OECD input-output table
Xj = £ Xij + TFj + TDj + Ej = £ Xi3 + Yi = Xi, i=i j=i
(1)
where i,j = 1,..., 45.
GDP is the sum of the value added Ej and TXS — INT — FNL, that is,
GDP = £ Ej + £ TDj = £ Yi + £ TFj + £ TDj, j=i j=i
i=i
jj j=m-n j=m-n
where n = 45, m = 54.
Similarly, we can obtain:
aij
xij
— Ej — TDj Xj
TDj
td ■ =_j
Xi
(2)
(3)
a^j denotes the quantity consumed by sector i per unit of output in sector j, i.e., the direct consumption coefficient. tdj denotes the sum of taxes fewer subsidies on intermediate and final domestic products in sector j as a share of total output Xi.
rWn =
Wj
Ej
rtj =
ONTj Ej
apj = £ aij. j=i
(4)
rwj represents sector j's labor wages as a share of its value added, rtj represents sector j's taxes on production as a share of its value added, and apj is the total share of intermediate consumption in sector j in sector i's annual output Xi.
Y TD,
Yr , Ej = (1 — apj ) Xj, td =^j=m-n j
GDP
GDP
(5)
Yri denotes the Yi share of GDP for final consumption in sector i, and td indicates the GDP share of the sum of taxes fewer subsidies on intermediate and final domestic products in sector m corresponding to final consumption to m — n.
By combining Eqs. (1)-(5), we can obtain the equation for GDP as affected by each relative variable:
GDP = £ Ej + £ TDj = (1 — apj ) X + td • GDP. j=i j=i
2.2. Dynamic Input-output Modeling
The characteristic properties of the necessary state variables, which provide the basis for simulating the anticipated evolution of the dynamic economy, are listed in the equation above. The economic sector is taken into account while creating the input-output table, and to treat it like a dynamic system, we must first create a set of differential equations that describe the course of economic development.
Definition 1 (Smirnov, 2021). At a constant level of technology, the change in total output Xt over time is defined as its derivative at time t, that is X(t) = , where,Xt = (X1 (t),... ,Xn(t)) is a vector of accrual-based outputs of the economic sector in terms of inputs at time t. The vector X(t) describes the acceleration of production in all sectors of the economy.
At the same level of technology and initial output X0, the relative increase in output requires a proportional increase in the comparative gross fixed asset
formation △iFcFF and the change in inventories △jjjI, where GFCF0 and III0 are the initial gross fixed asset formation and the initial inventory change, respectively. The amount of investment Cp required to increase output is proportional to the requisite acceleration Xj(t). As a result, when combined with the definition 1, the ratio that defines the amount of investment in each sector of the economy will take the following form:
Cpi(t) = Fei ■ Xj(t), (7)
where the Fei value represents the capital density of each sector of the economy. It is the coefficient of proportionality between the growth of output Xi(t) and the amount of investment Cpi(t) required to ensure it.
Fei indicates the output per unit time ^st = ^^-Ti^7^. Consequently, Fei leads to an acceleration in the production of goods and services. Its value to each manufacturer is determined by the ratio, i.e:
= Cpi(t) = Cpi(t) (T2 — Ti) = Cpi(t) (T2 — Ti) ei Xi (T2) — Xi (Ti) Xi (Ti) (— 1) . ()
For the input-output tables published by the OECD, capital intensity can be written as:
Fe^ = X,(t X,M. (9)
Fen+i(<) = GDP (t -ft)'— GDP (t). (10)
In an input-output model, the construction of differential equations is related to the sources of investment. In general, the generators of investment include GFCF, III, government expenditures imports, and exports. However, for investment, we usually consider only two factors, GFCF and III, because they are sources of investment based on internal demand.
Based on the formula (7), we regard the investment as part of the net profit, which is written as:
Cp, = rnj ■ Prhj, j = 1,
(11)
where rnj denotes the share of net profit used for investment. The formula (11) can be written as follows:
Op, = GFCFj + IIIj, Cp
rn = -,
GDP1
(12) (13)
where rn denotes the volume of investment Cp as a share of GDP.
Based on (3), (4), (5), (7), we can obtain the set of dynamic differential equations describing the economic development of each sector considered.
X
(1 - apj — tdj ) (1 - rwj — rtj )
Fe,
YJak
\j=1
■ Xj + Yrj ■ GDP
(14)
GDP = Jn I V; (1 - apj ) Xj + td ■ GDP
Fe.
n+l
\j=1
In vector form, the system (14),(15) is shown below:
where,
X = DX, D = MQ
(15)
(16)
Q
Xll
X21
xn1
X12 X22
I
xn1
1 — api 1 — ap2
xin Yri X2n Yr2
xnn Yrn
1 - apn td
M
ai Fe\ 0
a2 Fe2
an Fen 0
0
td
a.j = rnj (1 — rwj — rtj) (1 — tdj — apj), j = 1, ■ ■ ■ , n.
Considering the direct consumption coefficient apj can better estimate the production cost and resource utilization, considering the wage rate rwj can assist in analyzing the relationship between production cost and output, and considering the net production tax rtj can more accurately calculate the production cost and profit, the model can help us to more accurately predict the total output as well as the GDP by considering all of these factors, which is helpful to the government and the enterprises to formulate the economic policy and strategic planning.
j
0
e
n+1
2.3. Input-output Analysis in the United States
For the U.S. input-output tables 2000-2018, the current paper utilized R Programming Language and Python to calculate capital intensity and control for errors between total output and actual total output, as detailed in (Dan, 2023).
We aggregated 45 sectors into three major categories, namely advanced manufacturing, abbreviated as AMI (comprising 12 sectors), modern services, abbreviated as MSI (comprising 12 sectors), and other industries (23 sectors). Based on Eqs. (14), (15), we used linear regression to predict capital intensity under inflation, denoted as Fa, Fen+1. Additionally, we employed time series analysis to forecast GDP and the total output of the three aggregated sectors from 2019 to 2024.
According to Fig. 2, the year 2007 required the highest capital expenditure per unit of total output for both advanced manufacturing and modern services, while the years 2001 and 2008 saw the lowest capital expenditure per unit of total output for these sectors. The minimal capital required for advanced manufacturing in 2001 reflects the sector's sluggish development during that period. The financial crisis of 2007-2008 is also one of the key reasons for the notably low capital expenditure in modern services in 2008. From the projections, considering inflation, the capital intensity for advanced manufacturing is consistently higher than that for modern services and other industries.
According to Fig. 3, from the historical data, the unit capital expenditure required to increase the unit of GDP was the highest in 2007, and the unit capital expenditure required to increase the unit of GDP was the lowest in 2008. For the 2019-2024 forecast, the lowest unit capital expenditure is required to increase unit GDP in 2022.
2000 2005 2010 2015 2020
Fig. 2. Predicted values of capital density Fei under real GDP
According to Fig. 4, nominal GDP is consistent with forecast GDP from 2000 to 2018. At this stage, the official U.S. GDP was released only until 2021. For our forecasts of 215,906,93 , 224, 878, 63, 233,095,49 and 240,621, 01 million for the years 2019, 2020, 2021, 2022, respectively, the variances are: 0.97%, 6.34%, 0.04% and 1.84%. The projected nominal GDP for 2023 and 2024 will be $ 247, 513, 36 and $ 253, 825, 81 million respectively.
According to Fig. 5, historical data from 2000 to 2018 shows that the actual values are consistent with the predicted values. Overall, the total output of the modern services industry has always been higher than that of advanced Manufacturing. Notably, the year 2008 serves as an inflection point for both sectors, with the financial
Fig. 3. Predicted values of capital density Fen+i under real GDP
Fig. 4. Nominal GDP and predicted GDP
crisis being an evident and significant cause. In the forecast for 2019-2024, the total output for modern services industry is predicted to rise, while the total output for advanced manufacturing industry is expected to decline. This trend further confirms that, in developed countries, the modern services industry contributes more significantly to the economy.
Fig. 5. Total outputs and predicted total inputs of the 3 sectors
3. Equilibrium Economic Program Control
The control system:
X(t) = DX (t) + Q u + f (t),
(17)
where X = (xi
,Xn)T e u = (ui,..., ur)T e
- state vector, T- transposition, control vector,
P(t), Q(t) are (n x n) and (n x r) matrices with continuous components, f(t) - n-dimensional continuous vector function - a disturbance.
This is a particular case of continuous one when D = constant, Q = constant. The general solution of the Cauchy form of system (17) is:
x (t, 0,x0) = Y(t)
xo + f Y-1(t)(Q(t)u(t) + f(t))dT
0
(18)
Let two points x0, x1 and interval t e [0, T] are given. It is necessary to find an admissible control u(t) such that
x (T, xo, u(^)) = xi.
(19)
The pair x0, xi is said to be controllable on interval t e [0, T] if there exists an admissible control that is a solution of equation (18).
3.1. Program control algorithm
© Check the complete controllable of the system.
Stationary systems: the Kalman criterion. rang [Q, DQ,..., Dn-1Q] = n. Nonstationary systems: Sufficient conditions for complete controllability. rangS (t*) = n.
© Calculate the fundamental matrix Y(t) of the homogeneous system X = D(t)x.
© Construct the matrix B(t) = Y-1(t)Q(t). © Calculate the Gramian A(T) = /0T B(t)Bt(t)dr.
© Calculate the vector n = Y-1(T)xi - xo - /QT Y-1(r)f (r)dr. © Solve the system A(T)c = r/. © Solve the integral equation JT B(r)v(r)dr = 0 or take v(t) = 0.
© Form the program control u(t) = BT(t)c + v(t).
3.2. The United States input-output program control
Informed by the computational analysis presented in Section 2.3, we have acquired data (n = 4) for the United States during 2017-2018. We set control for the advanced manufacturing industry to obtain the following non-homogeneous equation.
"0.0195 0.0013 0.0031 0.0047 1
0.0081 0.0155 0.0098 0.0151 X + 0
0.0240 0.0067 0.0114 0.0336 0
0.0375 0.0379 0.0315 0.0011 0
3823518.3 4164850
11922242.2 12656899.5
xo = 18063252.8 , xi = 19074639.5
19541523.1 20611103.4
To assess the system's complete controllability, we implement program control targeting a specific coordinate pair {x0,x1} within the interval [0,1]. Here, x0 and
x1 represent the total output of the advanced manufacturing industry, the modern service industry, others, as well as the GDP for the years 2017 and 2018, respectively.
D
0.0195 0.0013 0.0031 0.0047 0.0081 0.0155 0.0098 0.0151 0.0240 0.0067 0.0114 0.0336 0.0375 0.0379 0.0315 0.0011
Q
, f(t)
T = 1. xo
3823518.3 11922242.2 18063252.8 19541523.1
xi
4164850 12656899.5 19074639.5 20611103.4
Initially, the system's complete controllability is verified using the Kalman criterion for stationary systems, which necessitates a non-zero determinant of the controllability matrix (det D = 0 ).
Rank(Q, DQ, D2Q, D3Q)
1 0.0195 0.00064143 0.00002 0 0.0081 0.00108495 0.00006 0 0.024 0.00108495 0.00010 0 0.0375 0.00182829 0.00013
Therefore, it is completely controllable.
Due to the intricate nature of the computational formulae involved, intermediate steps are omitted for brevity. Consequently, the control input u(t) is presented directly. For practical purposes, the constant e is approximated by 2.7183, facilitating a numerical solution to the control problem.
u(t) = [[1.6118 x 1018 x ((-0.7006 x 2.71830 003t + 0.0113 x 2.71830 015t +0.3318 x 2.7183°'059i +0.3574/2.718°'03i)x(-0.2937x2.718°'003i) + - • -+(0.6445/2.7183°'03i)x (0.4216 x 2.7180'003i+0.0018 x 2.71830'015i +0.3076 x2.71830'059i +0.269/2.7180'03i) x (• • • + 0.218 x 2.7183a015i + 0.217 x 2.7183a059i + 0.076/2.7183a03i)))]j.
To evaluate the effectiveness of the proposed control strategy, Table 1 delineates the discrepancies between the aggregated outputs of the advanced manufacturing industry, modern service industry, and other sectors, as influenced by the control variables, in comparison to the actual GDP of the United States for the year 2018.
Figure 6 reveals significant amplitude of fluctuations in the trajectory of key economic indicators under the proposed control scheme. The control mechanism we have constructed here is merely the initial step towards achieving the ultimate objective. In future work, we will consider a discrete form of control functions in order to smooth out fluctuations in the program mode.
Table 1. The error rate between the control total output and the actual total output
Year 2018 actual total output 2018 controlled total output Error rate
AMI 4164850 4165722.87 0.02%
MSI 12656899.5 12667558.54 0.08%
Others 19074639.5 19085804.38 0.05%
GDP 20611103.4 20620479.56 0.04%
4
*10s
Fig. 6. The control process of total output and GDP among 3 sectors in 2017-2018
4. Input-output Model of Saddle Point Equilibrium Strategy
Drawing upon the Leontief input-output relational framework, we can delineate a general dynamic input-output model. The input-output processes of individual sectors are conceptualized as dynamic economic activities, which determine that the output levels of each sector, as well as the final consumption products, are functions of time, denoted by t. This temporal dimension introduces a dynamic perspective to the traditionally static input-output analysis, allowing for a more nuanced understanding of the economic interdependencies and their evolution over time.
X (t) = AX (t)+ BX(t) + Y (t), (20)
where, t G [t*,t*](product planning period), n dimension vectors X(t) and Y(t) respectively is the output vector and the final consumption product vector (excluding the investment part). The n time-invariant matrix A is the direct consumption coefficient matrix, whose elements satisfy aj < 1 and 0 < aj < 1, i, j = 1, 2, • • • , n. X(t) is the rate of change in investments in each sector. The n time-invariant matrix B is the investment coefficient matrix.
Building on the foundational work of Mao, 1992, our study takes a holistic approach to addressing the stochastic uncertainties inherent in the national economy. We delve into the formulation of a continuous dynamic input-output model that operates under the most adverse conditions of ambiguity. By incorporating concepts from game theory, we frame the model as a linear quadratic differential game, where the control strategy of the first player, denoted by X(t), is juxtaposed against the uncertain factor-represented by the stochastic variable z(t) — which serves as the control strategy for the second player.
Our analysis is geared towards identifying the optimal control strategy for the dynamic input-output system. This is achieved by rigorously deriving the saddle point equilibrium strategy, which encapsulates the most favorable course of action
for both players under the given conditions of uncertainty. The resolution of this equilibrium strategy provides valuable insights into the optimization of the dynamic input-output system in the face of stochastic disturbances.
4.1. The game model of dynamic input-output system
Owing to the intricate nature of real economic activities, formula (20) is insufficient to accurately capture genuine economic laws. Consequently, we propose an amendment to the model by introducing a random factor, denoted as variable z(t), to account for uncertain factors in social and economic activities. This results in the formulation of a new model represented by formula (21).
X (t) = AX (t)+ BX(t) + Y (t) + z(t). (21)
The investment growth rates X (t) of each sector determine the changes in their output capabilities, and by adjusting them, the final consumption product vector Y(t) can be controlled. Therefore, the continuous dynamic input-output model (21) can be transformed into the following state-space representation.
r X(t) = u(t) (22)
\Y(t) = (I - A)X(t) - Bu(t) - z(t) (22)
Among them, t e [t*,t*]. I is the identity matrix of order n. u(t) e Rn is the decision-making control variable of player 1, z(t) e Rn is the "natural" decision-making control variable of player 2, output level vector X(t) is the state variable of the system, and final consumption product vector Y(t) is the output variable of the system.
We use the n dimensional column vector G(t) to represent the social demand vector for products, assuming it is a known continuous function vector. When the national economy is in a dynamic equilibrium, the social demand vector G(t) is equal to the system output vector Y(t). However, it is difficult to eliminate the supply-demand imbalance. At this time, on the one hand, we hope to make the performance index J(u,z) obtain a minimum value by adjusting the control variable u(t). On the other hand, it is hoped that under the worst interference of the random factor z(t), the performance index J(u, z) takes the maximum value, where,
J(u, z) = 2 ¡[ [Y(t) - G(t)]TP[Y(t) - G(t)] + uT (t)Ru(t)dt, (23)
where, P and R are positive definite matrices of order n respectively.
Their practical significance is to distinguish between the difference between the number of final consumer products provided by each sector and the social demand, and the difference in the primary and secondary degrees required by the changes in the output capacity of each sector. At this point, a complete dynamic input-output system game model is formed. In the following chapters, we will use the method of solving the saddle point equilibrium strategy to solve this problem.
4.2. Optimal control construction
The above model is a normal linear quadratic differential game model. Next, we use the result of Deissenberg to solve it, and construct the Hamilton function of the system (23) as follows:
H(X, u, Z, 0,t) = 1 [Y(t) - G(t)]TP[Y(t) - G(t)] + 1 uTRu + 0Tu. (24)
Then, d H
— = (R + BTPB) u - BTP(I - A)X + A + BTPG(t) + BTPz, du
d H
— = Pz - P(I - A)X + PG(t) + PBu, dz
d H
— = (I - A)TP(I - A)x - (I - A)TPBu - (I - A)TPz - (I - A)TPG(t). dX
According to the minimum principle, the optimal control of the system (23) u*(t), z*(t) and the optimal trajectory X*(t) and the corresponding co-state variable 0*(i) satisfy:
XT = u,X (to) = Xo,
9 = -(I - A)TP(I - A)X + (I - A)TPBu + (I - A)TPz + (I - A)TPG(t), 0 (t*) = 0, u = -R-1(9,
z = (I - A)X - BR-1(9 - G(t).
(25)
Let:
0(t) = S(t)X (t) + v(t). (26)
Then after proper calculation, S(t) and v(t) respectively satisfy the following matrix Riccati differential equation:
S(t) + (I - A)TPBR-1S(t) + S(t)R-1BTPT(1 - A) - S(t)R-1S(t) = 0, S (t*) = 0.
(27)
v(t) + (1 - A)TPBR-1r(t) + v(t)R-1BTPT(1 - A) - S(t)R-1v(t) = 0, v (t*) = 0.
Substituting equation (26) into equation (25), the optimal control law of the system is:
u*(t) = - R-1S(t)X(t)- R-1v(t), (28)
z*(t)=[l - BR-1S(t)) X(t)-BR-1v(t) - G(t), ( )
where, S(t) and v(t) are uniquely determined by equation (27).
Hence, this section uses the saddle point equilibrium theory in the differential game to study the multi-sector dynamic input-output problem in macroeconomic decision-making and designs the realization method of the control strategy.
5. Conclusion
The present study commences with the formulation of a dynamic input-output model, building upon the framework established by the OECD. Using the United States as a case study, the model, in conjunction with the least squares method, is employed to forecast capital intensity amidst inflationary conditions. For the U.S.
input-output tables, time series analysis is utilized to project the trajectory of the maximum aggregate total output as well as the GDP for the period spanning 2019 to 2024. Our findings indicate that the contribution of modern services to the U.S. GDP surpasses that of advanced manufacturing. In light of this, we recommend that the United States should relax trade restrictions and proactively engage in international economic and technological collaborations to foster robust economic growth.
Furthermore, the study integrates classical program control theory with a linear non chi-square differential equation that encompasses advanced manufacturing, modern services, other sectors, and GDP. By applying control to the aggregated advanced manufacturing sector in the United States for the years 2017 to 2018, we derive equations that account for the control of multiple influencing factors and validate the effectiveness of program control through numerical analysis.
As a final point, we explore the optimal strategy design for the continuous variant based on the dynamic input-output model proposed by Leontief. The dynamic input-output system is conceptualized as a saddle-point equilibrium game model, and the saddle-point equilibrium strategy is adopted to devise an innovative approach for resolving dynamic input-output challenges.
References
Dan, W. (2023). https://github.com/WangDan199810/Master-thesis-Wang-Dan.git. Jin, M. Z., Tang, R. Z., Ji, YJ, Liu, F., Gao, L., Huisingh, D. (2017). Impact of advanced manufacturing on sustainability: An overview of the special volume on advanced manufacturing for sustainability and low fossil carbon emissions. Journal of Cleaner Production, 161, 69-74.
Kang, J.-L., Mao, Y. Y., and Yang, Z. F. (2002). An optimal control model for dynamic
input and output with constraints. Journal of Tianjin University, 35(3), 384-388. Leontief, W. (1956). Factor Proportions and the Structure of American Trade: Further Theoretical and Empirical Analysis. Review of Economics and Statistics, 38(4), 386407.
Mao, Y., Fuping, B. (1992). Dynamic input-output optimal control model. Mathematical
Practice and Understanding, 1992(1), 58-64. Miller, R. E. Blair, P. D. (2009). Dynamic Input-Output Models: Basic Concepts and Solution Methods. In: Input-Output Analysis: Foundations and Extensions Cambridge University Press, 3, 331-369. Mourtzis, D., Vlachou, E., Zogopoulos, V., Gupta, R., Belkadi, F., Debbache, A., Bernard, A. (2018). Customer feedback gathering and management tools for product-service system design. Procedia CIRP, 67, 577-582. Qu, B., Huang, X. (1999). Optional production strategy for consumption tracking dynamic
input-output systems. Journal of Systems Engineering, 12(4), 325-329. Singhry, H.B., Abd, Rahman, A., Imm, N. S. (2016). Effect of advanced manufacturing technology, concurrent engineering of product design, and supply chain performance of manufacturing companies. The International Journal of Advanced Manufacturing Technology, 86, 663-669. Smirnov, N. V., Peresada, V. P., Smirnova, T. E. (2021). Dynamic models of inter-industrial balance. In: Intersectoral balance: analysis of dynamics and management of macroeco-nomic trends (N. V. Smirnov., ed), Vol. 2021.2, pp. 58-80. Lan Press (in Russian). Tamasyan, G. Sh. (2008). Program management and observability. Software controls, 2008, 28-32.