• Nie Znaleziono Wyników

,Q WU R G X F WLR Q DQG SUHOLP LQDULHV /HW XV FRQVLGHU D V\VWHP RI LQ

N/A
N/A
Protected

Academic year: 2021

Share " ,Q WU R G X F WLR Q DQG SUHOLP LQDULHV /HW XV FRQVLGHU D V\VWHP RI LQ"

Copied!
5
0
0

Pełen tekst

(1)

52&=1,.O 32/6.,(*2 72:$5=<67:$ 0$7(0$7<&=1(*2

6HULD ,,, 0$7(0$7<.$ 67262:$1$ ;;;; 

+H Q U \ N %U ] H V N Z L Q L H Z L F ]D Q G .U ] \ V ] W R I .O D F ] \ Q V N L

3R]QDQ

&RQVLVWHQF\ H[DPLQDWLRQ RIOLQHDU LQHTXDOLW\ V\VWHP

3UDFD ZSW\QWID GR 5HGDNFML      

6 X P P D U \  $Q LQYHVWLJDWLRQ RI WKH FRQVLVWHQF\ RI OLQHDU LQHTXDOLW\ V\VWHP LV FRQ

VLGHUHG ,W KDV EHHQ SURYHQ WKDW WKH V\VWHP RI OLQHDU LQHTXDOLWLHV$[!   LV FRQVLV

WHQW LI DQG RQO\ LI IRU DQ\ JHQHUDOL]HG LQYHUVH$7RI PDWUL[$WKH V\VWHP RI HTXDWLRQV

^,³$$a Y ³ ,³$$a EKDV D QRQQHJDWLYH VROXWLRQ GXH WR YHFWRUY&RQVLVWHQF\ RI

WKH DERYH V\VWHP GRHV QRW GHSHQG RQ WKH FKRLFH RIPDWUL[$7

7KH SDSHU SUHVHQWV DOVR PHWKRGV RI WKH H[DPLQDWLRQ LQ WKH H[LVWHQFH RI QRQQHJDWLYH

VROXWLRQV RI OLQHDU HTXDWLRQ V\VWHP

 ,Q WU R G X F WLR Q  DQG SUHOLP LQDULHV /HW XV FRQVLGHU D V\VWHP RI LQ

HTXDOLWLHV

 

$ [  ! IW

ZKHUH $  LV D P [ S PDWUL[ IW LV D P [  YHFWRU DQG WKH QRWDWLRQ $ [  ! IW

GHQRWHV WKH IDFW WKDW YHFWRU $ [ ³IWKDV DOO QRQQHJDWLYHHOHPHQWV 6\VWHPRI   IRUP LV HQFRXQWHUHG LQ PDQ\ SUDFWLFDO SUREOHPV 2QH RIWKHP LV WKH

HVWLPDWLRQ RISDUDPHWHUV LQ WKH IL[HG OLQHDU VWDWLVWLFDO PRGHO LQ ZKLFK V\V

WHP $ [ ! IW SUHVHQWV DQ DGGLWLRQDO QRQVWRFKDVWLF LQIRUPDWLRQ DERXW YHFWRU

RIXQNQRZQ SDUDPHWHUV 6XFK PRGHO ILQGV DSSOLFDWLRQ ERWK LQ HFRQRPHWULF

SUREOHPV VHH HJ ,WR  6FKPLGW   DQG LQELRPHWULFSUREOHPV VHH

HJ $UPVWURQJ DQG )URPH   $Q DQDO\VLV RIWKLV PRGHOLVSUHVHQWHG E\

-XGJH DQG 7DND\DPD   /RZHOO DQG 3UHVFRWW   /HYLV DQG 2GHOO

  /LHZ D DQG E  (VFREDU DQG 6NDUSHQHHV   :HUQHU

  DQG .ODF]\LLVNL   7KH V\VWHP RIOLQHDU LQHTXDOLWLHV LQWKH IRUP

RI   LV VRPHWLPHV DGGHG WR D PL[HG VWDWLVWLFDO OLQHDU PRGHO VHH 6HDO\

. H\ ZRUGV DQG SKUDVHV FRQVLVWHQF\ RI OLQHDU LQHTXDOLW\ V\VWHP QRQQHJDWLYH VHPL

SRVLWLYH  VROXWLRQV RIOLQHDU HTXDWLRQ V\VWHPV OLQHDU PRGHOV JHQHUDOL]HG LQYHUVH RI PDWUL[

(2)

1971), in which system (1.1) can indicate, as in Brown (1978) an additional nonstochastic information about the vector of variance components.

System of inequality (1.1) occurs also in linear programming problems as a system of linear constraints (see Bazaraa and Shetty, 1979; Chap. I and II).

In this paper for the research of inconsistency system of linear inequalities and their solution the notion of generalized inverse of matrix is used. It will give an idea of linear inequalities from the matrix point of view. Furthermore, in case where the generalized inverse of matrix is known (it frequently occurs in linear models) it may have a practical importance.

The consistency of the linear inequality system (1.1) will be solved on the basis of the existence or nonexistence of nonnegative solutions of the corresponding linear equation system. Hence, it contains some information about the linear equations, mainly in the context of the existence of their nonnegative solutions.

2. R esults. Let A~ denote a generalized inverse of matrix A, i.e. any matrix satisfying condition A = AA~A.

First, we consider a system of equations

(2.1) Rx = a,

where R is q x p matrix and 8 is q x 1 vector. In linear - equation theory well known is

LEMMA

1. (a) The equation system (2.1) is consistent if and only if there occurs one of the following conditions:

(i) r(R. 8) = r(R ) ( see known Kronecker-Capelly Theorem), (ii) 8 E 7Z(R) or, equivalently,

(iii) R R ~s = 8, for some R (see Rao and Mitra, 1971, p. 23),

where 'IZ(-) andr(-) denotes the column space and rank of a matrix argument, respectively.

(b) Consistent, system (2.1) has a general solution in the form

(2.2) x = R ~s + (I - R~R)v,

where v is any vector from space TV’ (see Rao and Mitra, 1971, p. 23).

Considering the consistency of system A x > b , defined in (1.1), we start with a simple statement, which is a direct consequence of linear-equation theory.

LEMMA 2.

The following conditions are equivalent:

(i) A x > h for some vector x,

(ii) 3v > 0 : A x = b + v for some vector x,

(3)

Consistency examination of linear inequality system 85

(iii) 3v > 0 : b + v e 7Z(A) (see Lemma 1(a) (ii)),

(iv) 3v > 0 :r(A ,6 4- v) — r (A) (see Lemma 1(a) (i)),

(v) 3v > 0 : AA (b + v) — b + v for some A~ (see Lemma 1(a) (iii)).

As conclusions t,o the Lemma, let us formulate the sufficient, conditions for consistency of system (1.1).

Lemma 2 (iii) with v = 0 implies that if

(2.2) b e 71(A),

then system (1.1) is consistent.

Note relation (2.2) is not a necessary condition of the consistency of system (1.1).

For example the system of inequalities - x i - x 2 > -4

6.'zq T

2x2

^ 8

X\ + hx2 > 8

is consistent (e.g. vector x — (1,2)' is its solution), but condition (2.2) is not satisfied.

It is easily visible that relation (2.2) is particularly satisfied if matrix A is of a full row rank (i.e. if r(A) = m).

Now wo are going to show that about the consistency of system (1.1) decides the existence of nonnegative solution of some linear equation system.

Th e o r e m

2.1. The system of inequalities (1.1) is consistent if and only if the system of equations

(2.3) (I - A A ~)v = - ( I - AA~)b

has a nonneyative solution due to vector v for any generalized inverse A of matrix A. Consistency of system (2.3) does not depend on the choice of matrix A ~ .

P ro o f. From Lemma 2(v) we have (I — AA~)(b + v) = 0 , which leads to (2.3).

To complete the proof consider two following systems of equations (2.4) (I - A A f)v = - ( I - A A f)b

(2.5) (I - A A f)v = - ( I - A A f)b,

where A f and A f are two different generalized inverses of matrix A. Since (I - A, A ~ )(I - A tA~) = I - A tA ~ , i = 1,2, it follows from Lemma 1(b) that for any vector Z\ and z2 from space 7ZP

Vi = - ( I - A A f)b + A A fz i and

v2

=:

—(I

A A 2)b

T ^4^4.2 ^2

(4)

are general solutions of (2.4) and (2.5), respectively. If we assume that V\ > 0 then substituting Z

2

= A A f

Z\

— AA^b-^ A A fb we get V

2

= V\ > 0

,

which completes the proof.

The Theorem 2.1 and Lemma 1(b) yields

C O R O L L A R Y

2.1. Consistent system (1.1) has a general solution in the form

(2.6) x 0 = A~{b + v0) + (I - A~ A)z,

where Vo is nonnegative solution of (2.3) and z is any vector from space 1ZP.

P r o o f . Substituting (2.6) into (1.1) we obtain

Ax

0 =

AA~ (b

+

V

q

) + A(I

A~A)z = b

+ vq

> 6

,

what completes the proof.

Note, that in practical situations nonnegative solution of R x = s (or of (2.3)) can be obtained basing on the solution of the following linear pro- gramming problem:

(2.7.) minimization l'tt

subject to the constraints

(2.8) R x + u — s x > 0 and u > 0

.

It is known (see Gale, 1960, Chap. IV.5 and Bazaraa and Shetty, 1979, ]). 70) that nonnegative solution of R x = s exists if and only if the optimal solution of the problem (2.7)-(2.8) is attained at the point, in which the minimized function reaches the value zero.

References

A r m s tr o n g , R. D. and F roin e, E. L .‘(1976): A branch-and-bound solution of a restricted least squares problem, Technometrics 18, 447-450.

B row n , K. G. (1978): Estimation of variance components using residuals, JASA 73, 141 146.

B a z a r a a , M. S. and S h e tt y , C. M. (1979): Nonlinear Programming. Theory and Algo- rithms, .John Wiley & Sons., New York.

E sc o b a r , L. A. and S k a r p e n e e s, B. (1984): A closed form solution for the least squares regression problem with linear inequality constraints, Commnn. Statist. - Theor. Meth.

13, 1127 1134.

G a le , D. (1960): The theory of linear economic models, McGraw-Hill Book Company, Inc.

New York.

Ito , T. (1980): Methods of estimation for multi-market disequilrium models, Econometrica 48, 97 125.

J u d g e, G. G. and T a k a y a m a , T. (1966): Inequality restrictions in regression analysis, JASA 61, 166 -181.

K la c z y n s k i, K. (1994): On estimation of parameter functions in a weakly singular linear model with linear inequality restrictions (in Polish), Matematyka Stosowana 37, 49-66.

(5)

Consistency examination of linear inequality system 87

L ew is, T. O. and O d e ll, P.L. (1971): Estimation in Linear Models, Prentice Hall, Inc.

Englewood Cliffs, New Jersey.

Liew , Cli. K. (1976 a): Inequality constrained least-squares estimation, JASA 71, 746-751.

Liew , Cli. Iv. (1976 b): A tuio-stage least-squares estimation with inequality restrictions on para,meters, The Review of Economics and Statistics LVIII No. 2, 234-238.

L o w ell, M. C. and P r e s c o tt, E. (1970): Multiple regression with inequality constraints:

Pretesting bias, hypothesis testing and efficiency, JASA 65, 913-925.

lla o . C. R. and M itra , K . M. (1971): Generalized Inverse of Matrix and its Applications, Wiley, New York.

Schm idt., P. (1981): Constraints on the parameters in simultaneous tobit and probit mod- els. in Structural Analysis of Discrete Data with Econometric Applications (C.F. Man- ski et al., Eds.), M IT Press, Mas., Cambridge 1981.

S c a ly . .1. (1971): Linear spaces and unbiased estimation application to the mixed linear model, The Annals of Mathematical Methods 41, 1735-1748.

T h e il, H. (1971): Principles of Econometrics, Wiley, New York.

W a t e r m a n . M.S. (1974): A restricted least squares problem, Technometrics 16, 135-136.

W e r n e r . H. .1. (1990): On Inequality Constrained Generalized Least Squares Estimation, Linear Algebra Appl. 127, 379-392.

DKPAHTMKNT OF MATHEMATICAL AND STATISTICAL METHODS AGRICULTURAL UNIVERSITY OF POZNAN

UL. WO.ISKA POLSKIEGO 28 POLAND

Cytaty

Powiązane dokumenty

STUDENTS OF THE FACULTY OF LETTERS: ENGLISH-LANGUAGE MAJORS Second-year MA students of English-language majors will attend the Entrepreneurial Skills: Work, Business, Career course

In Section 4 we consider the linearized problem (3.1); first in 4(a) we prove the existence of solutions in a half-space, in 4(b) we obtain the regu- larity of solutions and in the

Key words and phrases : elliptic system, nonnegative radial solution, existence, fixed point theorem in cones.. The work was supported by NNSF

Di Blasio, Differentiability of spatially homogeneous solution of the Boltzmann equation in the non Maxwellian case, Comm.. Ehlers, Survey of general relativity theory,

Applying process II to the ν-summation, we have to consider the prime term and the almost-prime term... Then we have an

We first notice that if the condition (1.7) is satisfied then the a priori estimates for u − ε 1 (x) given in Corollary 3.3 can be modified so as to be independent of ε... Below

Let G be Green’s function of C\J (that trivially coincides with Siciak’s function defined before (3)).. The implication 3)⇒2) is proved.. This was the most

Our version of the proof does not use the Poisson integral representation of harmonic functions in the unit disk D2. In order to make our method easily understandable, we have