實(shí)驗(yàn)11Lingo中數(shù)據(jù)的導(dǎo)入與導(dǎo)出

1、實(shí)驗(yàn)報告(11)學(xué)生姓名學(xué)號指導(dǎo)老師鄭成勇實(shí)驗(yàn)時間地點(diǎn)綜合實(shí)驗(yàn)大樓實(shí)驗(yàn)題目Lingo中數(shù)據(jù)的導(dǎo)入與導(dǎo)出實(shí)驗(yàn)?zāi)康膌 熟悉Lingo 中數(shù)據(jù)的導(dǎo)入與導(dǎo)出實(shí)驗(yàn)內(nèi)容財政收入預(yù)測問題：財政收入與國民收入、工業(yè)總產(chǎn)值、農(nóng)業(yè)總產(chǎn)值、總?cè)丝凇⒕蜆I(yè)人口、固定資產(chǎn)投資等因素有關(guān)。下表列出了1952-1981年的原始數(shù)據(jù)，試構(gòu)造預(yù)測模型。年份國民收入（億元）x1工業(yè)總產(chǎn)值(億元)x2農(nóng)業(yè)總產(chǎn)值（億元）x3總?cè)丝冢ㄈf人）x4就業(yè)人口（萬人）x5固定資產(chǎn)投資（億元）x6財政收入(億元)y1952598349461574822072944184195358645547558796213648921619547075204

2、916026621832972481955737558529614652232898254195682571555662828230181502681957837798575646532371113928619581028123559865994266002563571959111416815096720726173338444196010791870444662072588038050619617571156434658592559013827119626779644616729525110662301963779104651469172266408526619649431250584704

3、992773612932319651152158163272538286701753931966132219116877454229805212466196712491647697763683081415635219681187156568078534319151273031969137221016888067133225207447197016382747767829923443231256419711780315679085229356203556381972183333657898717735854354658197319783684855892113665237469119741993

4、36968919085937369393655197521214254932924213816846269219762052430995593717388344436571977218949259719497439377454723197824755590105896259398565509221979270260651150975424058156489019802791659211949870541896568826198129276862127310007273280496810設(shè)財政收入y與國民收入x1、工業(yè)總產(chǎn)值x2、農(nóng)業(yè)總產(chǎn)值x3、總?cè)丝趚4、就業(yè)人口x5、固定資產(chǎn)投資x6之間的的

5、函數(shù)關(guān)系為。問題轉(zhuǎn)化為根據(jù)已知數(shù)據(jù)估計(jì)出，根據(jù)最小二乘法，作離差平方和 選擇使Q達(dá)到最小，即可求得財政收入的預(yù)測模型。請用lingo求解。要求：分別使用file函數(shù)和ole函數(shù)實(shí)現(xiàn)數(shù)據(jù)的導(dǎo)入，并將結(jié)果用text導(dǎo)出到result.txt文件中。實(shí)驗(yàn)指導(dǎo)1、文件輸入輸出函數(shù)輸入：file(filename) 輸出：text(filename)2、Lingo中的變量默認(rèn)都是非負(fù)的，使用free(x)可取消x的非負(fù)限制。3、從Excel文件中導(dǎo)入數(shù)據(jù)ole函數(shù)只能用在模型的集合定義段、數(shù)據(jù)段和初始段，使用格式可以分成以下幾種類型：（1）變量名1,變量名2=ole(文件名,數(shù)據(jù)塊名稱1,數(shù)據(jù)塊名稱2)

6、;從指定的Excel文件讀取數(shù)據(jù)，文件名可以包括擴(kuò)展名（.xls），還可以包含完整的路徑目錄名稱，如果沒有指定路徑，則默認(rèn)路徑是LINGO的當(dāng)前工作目錄。該文件中定義了兩個數(shù)據(jù)塊，其中的數(shù)據(jù)分別用來對變量1和變量2初始化。如果變量名是集合，則對應(yīng)數(shù)據(jù)塊應(yīng)該是文本格式表示的集合成員名，如果變量名是集合的屬性，則對應(yīng)數(shù)據(jù)塊應(yīng)該是一系列數(shù)字，并且，若變量是初始集合的屬性，則對應(yīng)的數(shù)據(jù)塊應(yīng)當(dāng)是一列數(shù)據(jù)，若變量是二維衍生集合的屬性，則對應(yīng)數(shù)據(jù)塊應(yīng)當(dāng)是二維矩形數(shù)據(jù)區(qū)域。ole函數(shù)無法讀取三維數(shù)據(jù)區(qū)域。（2）變量名1,變量名2=ole(文件名,數(shù)據(jù)塊名稱);左邊的兩個變量必須定義在同一個集合中，ole的參數(shù)

7、僅指定一個數(shù)據(jù)塊名稱，該數(shù)據(jù)塊應(yīng)當(dāng)包含類型相同的兩列數(shù)據(jù)，第1列賦值給變量1，第2列賦值給變量2。（3）變量名1,變量名2=ole(文件名);沒有指定數(shù)據(jù)塊名稱，默認(rèn)使用Excel文件中與變量名同名的數(shù)據(jù)塊。4將計(jì)算結(jié)果導(dǎo)出到Excel文件中使用ole函數(shù)能把計(jì)算結(jié)果寫入Excel文件，使用格式也有以下三種：（1）ole(文件名,數(shù)據(jù)塊名稱1,數(shù)據(jù)塊名稱2)=變量名1,變量名2;將兩個變量的內(nèi)容分別寫入指定文件的兩個預(yù)先已經(jīng)定義了名稱的數(shù)據(jù)塊，數(shù)據(jù)塊的長度（大?。┎粦?yīng)小于變量所包含的數(shù)據(jù)，如果數(shù)據(jù)塊原來有數(shù)據(jù)，則ole寫入語句運(yùn)行后原來的數(shù)據(jù)將被新的數(shù)據(jù)覆蓋。（2）ole(文件名,數(shù)據(jù)塊名稱)

8、=變量名1,變量名2;兩個變量的數(shù)據(jù)寫入同一數(shù)據(jù)塊（不止1列），先寫變量1，變量2寫入另外1列。（3）ole(文件名)=變量名1,變量名2;不指定數(shù)據(jù)塊的名稱，默認(rèn)使用Excel文件中與變量名同名的數(shù)據(jù)塊。實(shí)驗(yàn)結(jié)果與分析首先將上面的原始數(shù)據(jù)復(fù)制到表格中并命名為data存放至c盤中代碼：1.sets: yearSet/ole(c:data.xls,years)/:y; ItemSet/x1.x6/:b; link(yearset,Itemset):X; endsets data: y=ole(c:data.xls,incoming); X=ole(c:data.xls,data); enddat

9、a min=sum(yearset(i):(y(i)-b0-sum(ItemSet(k):b(k)*X(i,k)2); free(b0);for(ItemSet:free(b);2.運(yùn)行結(jié)果：1. Variable Value Reduced Cost B0 160.2137 0.000000 Y( 1952) 184.0000 0.000000 Y( 1953) 216.0000 0.000000 Y( 1954) 248.0000 0.000000 Y( 1955) 254.0000 0.000000 Y( 1956) 268.0000 0.000000 Y( 1957) 286.0000

10、 0.000000 Y( 1958) 357.0000 0.000000 Y( 1959) 444.0000 0.000000 Y( 1960) 506.0000 0.000000 Y( 1961) 271.0000 0.000000 Y( 1962) 230.0000 0.000000 Y( 1963) 266.0000 0.000000 Y( 1964) 323.0000 0.000000 Y( 1965) 393.0000 0.000000 Y( 1966) 466.0000 0.000000 Y( 1967) 352.0000 0.000000 Y( 1968) 303.0000 0.

11、000000 Y( 1969) 447.0000 0.000000 Y( 1970) 564.0000 0.000000 Y( 1971) 638.0000 0.000000 Y( 1972) 658.0000 0.000000 Y( 1973) 691.0000 0.000000 Y( 1974) 655.0000 0.000000 Y( 1975) 692.0000 0.000000 Y( 1976) 657.0000 0.000000 Y( 1977) 723.0000 0.000000 Y( 1978) 922.0000 0.000000 Y( 1979) 890.0000 0.000

12、000 Y( 1980) 826.0000 0.000000 Y( 1981) 810.0000 0.000000 B( X1) 0.4592372 0.000000 B( X2) -0.1121112E-01 0.000000 B( X3) -0.5130743 0.000000 B( X4) 0.7885025E-03 0.000000 B( X5) -0.2779719E-02 0.000000 B( X6) 0.3157047 0.000000 X( 1952, X1) 598.0000 0.000000 X( 1952, X2) 349.0000 0.000000 X( 1952,

13、X3) 461.0000 0.000000 X( 1952, X4) 57482.00 0.000000 X( 1952, X5) 20729.00 0.000000 X( 1952, X6) 44.00000 0.000000 X( 1953, X1) 586.0000 0.000000 X( 1953, X2) 455.0000 0.000000 X( 1953, X3) 475.0000 0.000000 X( 1953, X4) 58796.00 0.000000 X( 1953, X5) 21364.00 0.000000 X( 1953, X6) 89.00000 0.000000

14、 X( 1954, X1) 707.0000 0.000000 X( 1954, X2) 520.0000 0.000000 X( 1954, X3) 491.0000 0.000000 X( 1954, X4) 60266.00 0.000000 X( 1954, X5) 21832.00 0.000000 X( 1954, X6) 97.00000 0.000000 X( 1955, X1) 737.0000 0.000000 X( 1955, X2) 558.0000 0.000000 X( 1955, X3) 529.0000 0.000000 X( 1955, X4) 61465.0

15、0 0.000000 X( 1955, X5) 22328.00 0.000000 X( 1955, X6) 98.00000 0.000000 X( 1956, X1) 825.0000 0.000000 X( 1956, X2) 715.0000 0.000000 X( 1956, X3) 556.0000 0.000000 X( 1956, X4) 62828.00 0.000000 X( 1956, X5) 23018.00 0.000000 X( 1956, X6) 150.0000 0.000000 X( 1957, X1) 837.0000 0.000000 X( 1957, X

16、2) 798.0000 0.000000 X( 1957, X3) 575.0000 0.000000 X( 1957, X4) 64653.00 0.000000 X( 1957, X5) 23711.00 0.000000 X( 1957, X6) 139.0000 0.000000 X( 1958, X1) 1028.000 0.000000 X( 1958, X2) 1235.000 0.000000 X( 1958, X3) 598.0000 0.000000 X( 1958, X4) 65994.00 0.000000 X( 1958, X5) 26600.00 0.000000

17、X( 1958, X6) 256.0000 0.000000 X( 1959, X1) 1114.000 0.000000 X( 1959, X2) 1681.000 0.000000 X( 1959, X3) 509.0000 0.000000 X( 1959, X4) 67207.00 0.000000 X( 1959, X5) 26173.00 0.000000 X( 1959, X6) 338.0000 0.000000 X( 1960, X1) 1079.000 0.000000 X( 1960, X2) 1870.000 0.000000 X( 1960, X3) 444.0000

18、 0.000000 X( 1960, X4) 66207.00 0.000000 X( 1960, X5) 25880.00 0.000000 X( 1960, X6) 380.0000 0.000000 X( 1961, X1) 757.0000 0.000000 X( 1961, X2) 1156.000 0.000000 X( 1961, X3) 434.0000 0.000000 X( 1961, X4) 65859.00 0.000000 X( 1961, X5) 25590.00 0.000000 X( 1961, X6) 138.0000 0.000000 X( 1962, X1

19、) 677.0000 0.000000 X( 1962, X2) 964.0000 0.000000 X( 1962, X3) 461.0000 0.000000 X( 1962, X4) 67295.00 0.000000 X( 1962, X5) 25110.00 0.000000 X( 1962, X6) 66.00000 0.000000 X( 1963, X1) 779.0000 0.000000 X( 1963, X2) 1046.000 0.000000 X( 1963, X3) 514.0000 0.000000 X( 1963, X4) 69172.00 0.000000 X

20、( 1963, X5) 26640.00 0.000000 X( 1963, X6) 85.00000 0.000000 X( 1964, X1) 943.0000 0.000000 X( 1964, X2) 1250.000 0.000000 X( 1964, X3) 584.0000 0.000000 X( 1964, X4) 70499.00 0.000000 X( 1964, X5) 27736.00 0.000000 X( 1964, X6) 129.0000 0.000000 X( 1965, X1) 1152.000 0.000000 X( 1965, X2) 1581.000

21、0.000000 X( 1965, X3) 632.0000 0.000000 X( 1965, X4) 72538.00 0.000000 X( 1965, X5) 28670.00 0.000000 X( 1965, X6) 175.0000 0.000000 X( 1966, X1) 1322.000 0.000000 X( 1966, X2) 1911.000 0.000000 X( 1966, X3) 687.0000 0.000000 X( 1966, X4) 74542.00 0.000000 X( 1966, X5) 29805.00 0.000000 X( 1966, X6)

22、 212.0000 0.000000 X( 1967, X1) 1249.000 0.000000 X( 1967, X2) 1647.000 0.000000 X( 1967, X3) 697.0000 0.000000 X( 1967, X4) 76368.00 0.000000 X( 1967, X5) 30814.00 0.000000 X( 1967, X6) 156.0000 0.000000 X( 1968, X1) 1187.000 0.000000 X( 1968, X2) 1565.000 0.000000 X( 1968, X3) 680.0000 0.000000 X(

23、 1968, X4) 78534.00 0.000000 X( 1968, X5) 31915.00 0.000000 X( 1968, X6) 127.0000 0.000000 X( 1969, X1) 1372.000 0.000000 X( 1969, X2) 2101.000 0.000000 X( 1969, X3) 688.0000 0.000000 X( 1969, X4) 80671.00 0.000000 X( 1969, X5) 33225.00 0.000000 X( 1969, X6) 207.0000 0.000000 X( 1970, X1) 1638.000 0

24、.000000 X( 1970, X2) 2747.000 0.000000 X( 1970, X3) 767.0000 0.000000 X( 1970, X4) 82992.00 0.000000 X( 1970, X5) 34432.00 0.000000 X( 1970, X6) 312.0000 0.000000 X( 1971, X1) 1780.000 0.000000 X( 1971, X2) 3156.000 0.000000 X( 1971, X3) 790.0000 0.000000 X( 1971, X4) 85229.00 0.000000 X( 1971, X5)

25、35620.00 0.000000 X( 1971, X6) 355.0000 0.000000 X( 1972, X1) 1833.000 0.000000 X( 1972, X2) 3365.000 0.000000 X( 1972, X3) 789.0000 0.000000 X( 1972, X4) 87177.00 0.000000 X( 1972, X5) 35854.00 0.000000 X( 1972, X6) 354.0000 0.000000 X( 1973, X1) 1978.000 0.000000 X( 1973, X2) 3684.000 0.000000 X(

26、1973, X3) 855.0000 0.000000 X( 1973, X4) 89211.00 0.000000 X( 1973, X5) 36652.00 0.000000 X( 1973, X6) 374.0000 0.000000 X( 1974, X1) 1993.000 0.000000 X( 1974, X2) 3696.000 0.000000 X( 1974, X3) 891.0000 0.000000 X( 1974, X4) 90859.00 0.000000 X( 1974, X5) 37369.00 0.000000 X( 1974, X6) 393.0000 0.

27、000000 X( 1975, X1) 2121.000 0.000000 X( 1975, X2) 4254.000 0.000000 X( 1975, X3) 932.0000 0.000000 X( 1975, X4) 92421.00 0.000000 X( 1975, X5) 38168.00 0.000000 X( 1975, X6) 462.0000 0.000000 X( 1976, X1) 2052.000 0.000000 X( 1976, X2) 4309.000 0.000000 X( 1976, X3) 955.0000 0.000000 X( 1976, X4) 9

28、3717.00 0.000000 X( 1976, X5) 38834.00 0.000000 X( 1976, X6) 443.0000 0.000000 X( 1977, X1) 2189.000 0.000000 X( 1977, X2) 4925.000 0.000000 X( 1977, X3) 971.0000 0.000000 X( 1977, X4) 94974.00 0.000000 X( 1977, X5) 39377.00 0.000000 X( 1977, X6) 454.0000 0.000000 X( 1978, X1) 2475.000 0.000000 X( 1978, X2) 5590.000 0.000000 X( 1978, X3) 1058.000 0.000000 X( 1978, X4) 96259.00 0.000000 X( 1978, X5) 39856.00 0.000000 X( 1978, X6) 550.0000 0.000000 X( 1979, X1) 2702.000 0.000000 X( 1979, X2) 6065.000 0.0