实验11Lingo中数据的导入与导出

1、实验报告(11)学生姓名学号指导老师郑成勇实验时间地点综合实验大楼实验题目Lingo中数据的导入与导出实验目的l 熟悉Lingo 中数据的导入与导出实验内容财政收入预测问题：财政收入与国民收入、工业总产值、农业总产值、总人口、就业人口、固定资产投资等因素有关。下表列出了1952-1981年的原始数据，试构造预测模型。年份国民收入（亿元）x1工业总产值(亿元)x2农业总产值（亿元）x3总人口（万人）x4就业人口（万人）x5固定资产投资（亿元）x6财政收入(亿元)y1952598349461574822072944184195358645547558796213648921619547075204

2、916026621832972481955737558529614652232898254195682571555662828230181502681957837798575646532371113928619581028123559865994266002563571959111416815096720726173338444196010791870444662072588038050619617571156434658592559013827119626779644616729525110662301963779104651469172266408526619649431250584704

3、992773612932319651152158163272538286701753931966132219116877454229805212466196712491647697763683081415635219681187156568078534319151273031969137221016888067133225207447197016382747767829923443231256419711780315679085229356203556381972183333657898717735854354658197319783684855892113665237469119741993

4、36968919085937369393655197521214254932924213816846269219762052430995593717388344436571977218949259719497439377454723197824755590105896259398565509221979270260651150975424058156489019802791659211949870541896568826198129276862127310007273280496810设财政收入y与国民收入x1、工业总产值x2、农业总产值x3、总人口x4、就业人口x5、固定资产投资x6之间的的

5、函数关系为。问题转化为根据已知数据估计出，根据最小二乘法，作离差平方和 选择使Q达到最小，即可求得财政收入的预测模型。请用lingo求解。要求：分别使用file函数和ole函数实现数据的导入，并将结果用text导出到result.txt文件中。实验指导1、文件输入输出函数输入：file(filename) 输出：text(filename)2、Lingo中的变量默认都是非负的，使用free(x)可取消x的非负限制。3、从Excel文件中导入数据ole函数只能用在模型的集合定义段、数据段和初始段，使用格式可以分成以下几种类型：（1）变量名1,变量名2=ole(文件名,数据块名称1,数据块名称2)

6、;从指定的Excel文件读取数据，文件名可以包括扩展名（.xls），还可以包含完整的路径目录名称，如果没有指定路径，则默认路径是LINGO的当前工作目录。该文件中定义了两个数据块，其中的数据分别用来对变量1和变量2初始化。如果变量名是集合，则对应数据块应该是文本格式表示的集合成员名，如果变量名是集合的属性，则对应数据块应该是一系列数字，并且，若变量是初始集合的属性，则对应的数据块应当是一列数据，若变量是二维衍生集合的属性，则对应数据块应当是二维矩形数据区域。ole函数无法读取三维数据区域。（2）变量名1,变量名2=ole(文件名,数据块名称);左边的两个变量必须定义在同一个集合中，ole的参数

7、仅指定一个数据块名称，该数据块应当包含类型相同的两列数据，第1列赋值给变量1，第2列赋值给变量2。（3）变量名1,变量名2=ole(文件名);没有指定数据块名称，默认使用Excel文件中与变量名同名的数据块。4将计算结果导出到Excel文件中使用ole函数能把计算结果写入Excel文件，使用格式也有以下三种：（1）ole(文件名,数据块名称1,数据块名称2)=变量名1,变量名2;将两个变量的内容分别写入指定文件的两个预先已经定义了名称的数据块，数据块的长度（大小）不应小于变量所包含的数据，如果数据块原来有数据，则ole写入语句运行后原来的数据将被新的数据覆盖。（2）ole(文件名,数据块名称)

8、=变量名1,变量名2;两个变量的数据写入同一数据块（不止1列），先写变量1，变量2写入另外1列。（3）ole(文件名)=变量名1,变量名2;不指定数据块的名称，默认使用Excel文件中与变量名同名的数据块。实验结果与分析首先将上面的原始数据复制到表格中并命名为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.运行结果：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