Homework 61. Analysis of a power control algorithm. In this problem we consider again the powercontrol method described in homework 1 problem 1. Please refer to this problem for the setupand background. In that problem, you expressed the power control method as a discrete-timelinear dynamical system, and simulated it for a specific set of parameters, with several valuesof initial power levels, and two target SINRs. You found that for the target SINR value γ = 3,the powers converged to values for which each SINR exceeded γ, no matter what the initialpower was, whereas for the larger target SINR value γ = 5, the powers appeared to diverge,and the SINRs did not appear to converge. You are going to analyze this, now that you knowalot more about linear systems.a) Explain the simulations. Explain your simulation results from the problem 1(b) for thegiven values of G, α, σ, and the two SINR threshold levels γ = 3 and γ = 5.b) Critical SINR threshold level. Let us consider fixed values of G, α, and σ. It turns outthat the power control algorithm works provided the SINR threshold γ is less than somecritical value γcrit (which might depend on G, α, σ), and doesn’t work for γ > γcrit.(‘Works’ means that no matter what the initial powers are, they converge to values forwhich each SINR exceeds γ.) Find an expression for γcrit in terms of G ∈ Rn×n, α, andσ. Give the simplest expression you can. Of course you must explain how you came upwith your expression.2. Real modal form. We learned about the modal form of a system in class. Show that whensome of eigenvalues of the dynamics matrix A are complex, the system can be put in realmodal form (Assuming the eigenvectors of A are independent):S1AS = diag (Λr, Mr+1, Mr+3, . . . , Mn?1)where Λr = diag(λ1, . . . , λr) are the real eigenvalues, andMj =�σj ωjωj σj�, λj = σj + iωj , j = r + 1, r + 3, . . . , n ? 1where λj are the complex eigenvalues (one from each conjugate pair). Clearly explain whatthe matrix S is.Generate a matrix A in R10×10 using A=randn(10). (The entries of A will be drawn froma unit normal distribution.) Find the eigenvalues of A. If by chance they are all real, generatea new instance of A. Find the real modal form of A, i.e., a matrix S such that S1AS hasthe real modal form. Your solution should include the source code that you use to find S, andsome code that checks the results (i.e., computes S1AS to verify it has the required form).13. Output response envelope for linear system with uncertain initial condition. Weconsider the autonomous linear dynamical system x˙ = Ax, y(t) = Cx(t), where x(t) ∈ Rn andy(t) ∈ R. We do not know the initial condition exactly; we only know that it lies in a ball ofradius r centered at the point x0:kx(0) x0k ≤ r.We call x0 the nominal initial condition, and the resulting output, ynom(t) = CetAx0, thenominal output. We define the maximum output or upper output envelope asy(t) = max{y(t) | kx(0) x0k ≤ r},i.e., the maximum possible value of the output at time t, over all possible initial conditions.(Here you can choose a different initial condition for each t; you are not required to find asingle initial condition.) In a similar way, we define the minimum output or lower outputenvelope asy(t) = min{y(t) | kx(0) x0k ≤ r},i.e., the minimum possible value of the output at time t, over all possible initial conditions.a) Explain how to find y(t) and y(t), given the problem data A, C, x0, and r.b) Carry out your method on the problem data in uie_data.m. On the same axes, plotynom, y, and y, versus t, over the range 0 ≤ t ≤ 10.4. Spectral mapping theorem. Suppose f : R → R is analytic, i.e., given by a power seriesexpansionf(u) = a0 + a1u + a2u2 + · · ·(where ai = f(i)(0)/(i!)). (You can assume that we only consider values of u for which thisseries converges.) For A ∈ Rn×n, we define f(A) asf(A) = a0I + a1A + a2A2 + · · ·(again, we’ll just assume that this converges).Suppose that Av = λv, where v 6= 0, and λ ∈ C. Show that f(A)v = f(λ)v (ignoring theissue of convergence of series). We conclude that if λ is an eigenvalue of A, then f(λ) is aneigenvalue of f(A). This is called the spectral mapping theorem.To illustrate this with an example, generate a random 3 × 3 matrix, for example usingA=randn(3). Find the eigenvalues of (I + A)(I ? A)1 by first computing this matrix, thenfinding its eigenvalues, and also by using the spectral mapping theorem. (You should get veryclose agreement; any difference is due to numerical round-off errors in the various compuations.)5. Interconnection of linear systems. Often a linear system is described in terms of ablock diagram showing the interconnections between components or subsystems, which are2themselves linear systems. In this problem you consider the specific interconnection s代写SINR作业、代做Matlab实验作业、代写Matlab语言程序作业、代做linear systems作业 代做留学hownbelow:S Tu v yw1w2Here, there are two subsystems S and T. Subsystem S is characterized byx˙ = Ax + B1u + B2w1, w2 = Cx + D1u + D2w1,and subsystem T is characterized byz˙ = F z + G1v + G2w2, w1 = H1z, y = H2z + Jw2.We don’t specify the dimensions of the signals (which can be vectors) or matrices here. You canassume all the matrices are the correct (i.e., compatible) dimensions. Note that the subscriptsin the matrices above, as in B1 and B2, refer to different matrices. Now the problem. Expressthe overall system as a single linear dynamical system with input, state, and output given by�respectively. Be sure to explicitly give the input, dynamics, output, and feedthrough matricesof the overall system. If you need to make any assumptions about the rank or invertibilityof any matrix you encounter in your derivations, go ahead. But be sure to let us know whatassumptions you are making.6. Analysis of investment allocation strategies. Each year or period (denoted t = 0, 1, . . .)an investor buys certain amounts of one-, two-, and three-year certificates of deposit (CDs)with interest rates 5%, 6%, and 7%, respectively. (We ignore minimum purchase requirements,and assume they can be bought in any amount.) B1(t) denotes the amount of one-year CDs bought at period t. B2(t) denotes the amount of two-year CDs bought at period t. B3(t) denotes the amount of three-year CDs bought at period t.We assume that B1(0) + B2(0) + B3(0) = 1, i.e., a total of 1 is to be invested at t = 0. (Youcan take Bj (t) to be zero for t sum of six terms: 1.05B1(t1), i.e., principle plus 5% interest on the amount of one-year CDs bought oneyear ago. 1.06B2(t2), i.e., principle plus 6% interest on the amount of two-year CDs bought twoyears ago.3 1.07B3(t 3), i.e., principle plus 7% interest on the amount of three-year CDs boughtthree years ago. 0.06B2(t 1), i.e., 6% interest on the amount of two-year CDs bought one year ago. 0.07B3(t 1), i.e., 7% interest on the amount of three-year CDs bought one year ago. 0.07B3(t 2), i.e., 7% interest on the amount of three-year CDs bought two years ago.The total wealth held by the investor at period t is given byw(t) = B1(t) + B2(t) + B2(t 1) + B3(t) + B3(t1) + B3(t2).Two re-investment allocation strategies are suggested. The 35-35-30 strategy. The total payout is re-invested 35% in one-year CDs, 35% in twoyearCDs, and 30% in three-year CDs. The initial investment allocation is the same:B1(0) = 0.35, B2(0) = 0.35, and B3(0) = 0.30. The 60-20-20 strategy. The total payout is re-invested 60% in one-year CDs, 20% in twoyearCDs, and 20% in three-year CDs. The initial investment allocation is B1(0) = 0.60,B2(0) = 0.20, and B3(0) = 0.20.a) Describe the investments over time as a linear dynamical system x(t + 1) = Ax(t),y(t) = Cx(t) with y(t) equal to the total wealth at time t. Be very clear about what thestate x(t) is, and what the matrices A and C are. You will have two such linear systems:one for the 35-35-30 strategy and one for the 60-20-20 strategy.b) Asymptotic wealth growth rate. For each of the two strategies described above, determinethe asymptotic growth rate, defined as limt→∞ w(t+ 1)/w(t). (If this limit doesn’t exist,say so.) Note: simple numerical simulation of the strategies (e.g., plotting w(t+ 1)/w(t)versus t to guess its limit) is not acceptable. (You can, of course, simulate the strategiesto check your answer.)c) Asymptotic liquidity. The total wealth at time t can be divided into three parts: B1(t) + B2(t�1) + B3(t2) is the amount that matures in one year (i.e., theamount of principle we could get back next year) B2(t) + B3(t1) is the amount that matures in two years B3(t) is the amount that matures in three years (i.e., is least liquid)We define liquidity ratios as the ratio of these amounts to the total wealth:L1(t) = (B1(t) + B2(t�1) + B3(t2))/w(t),L2(t) = (B2(t) + B3(t�1))/w(t),L3(t) = B3(t)/w(t).For the two strategies above, do the liquidity ratios converge as t → ∞ If so, to whatvalues? Note: as above, simple numerical simulation alone is not acceptable.4d) Suppose you could change the initial investment allocation for the 35-35-30 strategy,i.e., choose some other nonnegative values for B1(0), B2(0), and B3(0) that satisfyB1(0) + B2(0) + B3(0) = 1. What allocation would you pick, and how would it bebetter than the (0.35, 0.35, 0.30) initial allocation? (For example, would the asymptoticgrowth rate be larger?) How much better is your choice of initial investment allocations?Hint for part d: think very carefully about this one. Hint for whole problem: watchout for nondiagonalizable, or nearly nondiagonalizable, matrices. Don’t just blindlytype in matlab commands; check to make sure you’re computing what you think you’recomputing.转自：http://ass.3daixie.com/2018111723601656.html
