试验11Lingo中数据的导入与导出课案

1、学生姓名实验报告 （11）学号指导老师 郑成勇实验时间 地点 综合实验大楼实验题目 Lingo 中数据的导入与导出实 验 目 的熟悉 Lingo 中数据的导入与导出实验内年份国民收入（亿元） x1工业总产 值（亿元） x2农业总产 值（亿元） x3总人口 （万人） x4就业人口（万人） x5固定资产投 资（亿元） x6财政收入（亿元 ） y195259834946157482207294418419535864554755879621364892161954707520491602662183297248195573755852961465223289825419568257155566282

2、823018150268195783779857564653237111392861958102812355986599426600256357195911141681509672072617333844419601079187044466207258803805061961757115643465859255901382711962677964461672952511066230196377910465146917226640852661964943125058470499277361293231965115215816327253828670175393196613221911687745

3、422980521246619671249164769776368308141563521968118715656807853431915127303196913722101688806713322520744719701638274776782992344323125641971178031567908522935620355638197218333365789871773585435465819731978368485589211366523746911974199336968919085937369393655197521214254932924213816846269219762052

4、430995593717388344436571977218949259719497439377454723197824755590105896259398565509221979270260651150975424058156489019802791659211949870541896568826198129276862127310007273280496810财政收入预测问题：财政收入与国民收入、工业总产值、农业总产值、总人口、就 业人口、固定资产投资等因素有关。下表列出了1952-1981 年的原始数据，试构造预测模型。设财政收入 y 与国民收入 x1、工业总产值 x2、农业总产值 x3

5、、总人口 x4、就业人 口 x5、固定资产投资 x6 之间的的函数关系为 y 0 1x1 . 6x6 。问题转化为根 据已知数据估计出 0,., 6 ，根据最小二乘法，作离差平方和30 2Qyi 01xi1 . k xi6i1选择 0,., 6使 Q达到最小，即可求得财政收入的预测模型。请用 lingo 求解 0 ,., 6 。要求：分别使用 file 函数和 ole 函数实现数据的导入，并将结果 0 ,., 6 用text 导出到 result.txt 文件中。实 验 指 导1、文件输入输出函数输入： file（filename）输出： text（filename）2、 Lingo 中的变量

6、默认都是非负的，使用 free（x） 可取消 x 的非负限制。3、 从 Excel 文件中导入数据ole函数只能用在模型的集合定义段、 数据段和初始段， 使用格式可以 分成以下几种类型：（1）变量名 1, 变量名 2=ole （ 文件名 , 数据块名称 1, 数据块 名称 2）;从指定的 Excel 文件读取数据，文件名可以包括扩展名（ .xls），还可以 包含完整的路径目录名称，如果没有指定路径，则默认路径是 LINGO 的当 前工作目录。该文件中定义了两个数据块， 其中的数据分别用来对变量 1和变量 2 初 始化。如果变量名是集合， 则对应数据块应该是文本格式表示的集合成员名， 如果变量名

7、是集合的属性，则对应数据块应该是一系列数字，并且，若变量 是初始集合的属性，则对应的数据块应当是一列数据，若变量是二维衍生集 合的属性， 则对应数据块应当是二维矩形数据区域。 ole函数无法读取三维 数据区域。（2）变量名 1, 变量名 2= ole （ 文件名, 数据块名称 ）; 左边的两个变量必须定义在同一个集合中， ole 的参数仅指定一个数据 块名称，该数据块应当包含类型相同的两列数据，第 1 列赋值给变量 1，第 2 列赋值给变量 2。（3）变量名 1, 变量名 2= ole （ 文件名 ）;没有指定数据块名称， 默认使用 Excel 文件中与变量名同名的数据块。 4将计算结果导出到

8、 Excel 文件中使用ole函数能把计算结果写入 Excel 文件，使用格式也有以下三种： （1）ole （ 文件名 , 数据块名称 1, 数据块名称 2）= 变量名 1, 变量名 2;将两个变量的内容分别写入指定文件的两个预先已经定义了名称的数 据块，数据块的长度（大小）不应小于变量所包含的数据，如果数据块原来 有数据，则 ole 写入语句运行后原来的数据将被新的数据覆盖。（2）ole （ 文件名, 数据块名称）= 变量名 1,变量名 2; 两个变量的数据写入同一数据块 （不止 1列），先写变量 1，变量 2写入 另外 1 列。（3）ole （ 文件名 ）= 变量名 1, 变量名 2; 不

9、指定数据块的名称，默认使用 Excel 文件中与变量名同名的数据块。实 验 结 果 与 分 析首先将上面的原始数据复制到表格中并命名为data存放至 c盘中代码： 1.sets :yearSet/ ole (c:data.xls,years)/:y;ItemSet/x1.x6/:b; link(yearset,Itemset):X;endsetsdata :y= ole (c:data.xls,incoming);X= ole (c:data.xls,data);enddatamin = sum(yearset(i):(y(i)-b0-sum(ItemSet(k):b(k)*X(i,k)2);

10、free (b0);for (ItemSet: free (b);2.运行结果： 1.Variable Value Reduced Cost B0 160.2137 0.000000Y( 1952)184.00000.000000Y( 1953)216.00000.000000Y( 1954)248.00000.000000Y( 1955)254.00000.000000Y( 1956)268.00000.000000Y( 1957)286.00000.000000Y( 1958)357.00000.000000Y( 1959)444.00000.000000Y( 1960)506.00000

11、.000000Y( 1961)271.00000.000000Y( 1962)230.00000.000000Y( 1963)266.00000.000000Y( 1964)323.00000.000000Y( 1965)393.00000.000000Y( 1966)466.00000.000000Y( 1967)352.00000.000000Y( 1968)303.00000.000000Y( 1969)447.00000.000000Y( 1970)564.00000.000000Y( 1971)638.00000.000000Y( 1972)658.00000.000000Y( 19

12、73)691.00000.000000Y( 1974)655.00000.000000Y( 1975)692.00000.000000Y( 1976)657.00000.000000Y( 1977)723.00000.000000Y( 1978)922.00000.000000Y( 1979)890.00000.000000Y( 1980)826.00000.000000Y( 1981)810.00000.000000B( X1)0.45923720.000000B( X2)-0.1121112E-010.000000B( X3)-0.51307430.000000B( X4)0.788502

13、5E-030.000000B( X5)-0.2779719E-020.000000B( X6)0.31570470.000000X( 1952, X1)598.00000.000000X( 1952, X2)349.00000.000000X( 1952, X3)461.00000.000000X( 1952, X4)57482.000.000000X( 1952, X5)20729.000.000000X( 1952, X6)44.000000.000000X( 1953, X1)586.00000.000000X( 1953, X2)455.00000.000000X( 1953, X3)

14、475.00000.000000X( 1953, X4)58796.000.000000X( 1953, X5)21364.000.000000X( 1953, X6)89.000000.000000X( 1954, X1)707.00000.000000X( 1954, X2)520.00000.000000X( 1954, X3)491.00000.000000X( 1954, X4)60266.000.000000X( 1954, X5)21832.000.000000X( 1954, X6)97.000000.000000X( 1955, X1)737.00000.000000X( 1

15、955, X2)558.00000.000000X( 1955, X3)529.00000.000000X( 1955, X4)61465.000.000000X( 1955, X5)22328.000.000000X( 1955, X6)98.000000.000000X( 1956, X1)825.00000.000000X( 1956, X2)715.00000.000000X( 1956, X3)556.00000.000000X( 1956, X4)62828.000.000000X( 1956, X5)23018.000.000000X( 1956, X6)150.00000.00

16、0000X( 1957, X1)837.00000.000000X( 1957, X2)798.00000.000000X( 1957, X3)575.00000.000000X( 1957, X4)64653.000.000000X( 1957, X5)23711.000.000000X( 1957, X6)139.00000.000000X( 1958, X1)1028.0000.000000X( 1958, X2)1235.0000.000000X( 1958, X3)598.00000.000000X( 1958, X4)65994.000.000000X( 1958, X5)2660

17、0.000.000000X( 1958, X6)256.00000.000000X( 1959, X1)1114.0000.000000X( 1959, X2)1681.0000.000000X( 1959, X3)509.00000.000000X( 1959, X4)67207.000.000000X( 1959, X5)26173.000.000000X( 1959, X6)338.00000.000000X( 1960, X1)1079.0000.000000X( 1960, X2)1870.0000.000000X( 1960, X3)444.00000.000000X( 1960,

18、 X4)66207.000.000000X( 1960, X5)25880.000.000000X( 1960, X6)380.00000.000000X( 1961, X1)757.00000.000000X( 1961, X2)1156.0000.000000X( 1961, X3)434.00000.000000X( 1961, X4)65859.000.000000X( 1961, X5)25590.000.000000X( 1961, X6)138.00000.000000X( 1962, X1)677.00000.000000X( 1962, X2)964.00000.000000

19、X( 1962, X3)461.00000.000000X( 1962, X4)67295.000.000000X( 1962, X5)25110.000.000000X( 1962, X6)66.000000.000000X( 1963, X1)779.00000.000000X( 1963, X2)1046.0000.000000X( 1963, X3)514.00000.000000X( 1963, X4)69172.000.000000X( 1963, X5)26640.000.000000X( 1963, X6)85.000000.000000X( 1964, X1)943.0000

20、0.000000X( 1964, X2)1250.0000.000000X( 1964, X3)584.00000.000000X( 1964, X4)70499.000.000000X( 1964, X5)27736.000.000000X( 1964, X6)129.00000.000000X( 1965, X1)1152.0000.000000X( 1965, X2)1581.0000.000000X( 1965, X3)632.00000.000000X( 1965, X4)72538.000.000000X( 1965, X5)28670.000.000000X( 1965, X6)

21、175.00000.000000X( 1966, X1)1322.0000.000000X( 1966, X2)1911.0000.000000X( 1966, X3)687.00000.000000X( 1966, X4)74542.000.000000X( 1966, X5)29805.000.000000X( 1966, X6)212.00000.000000X( 1967, X1)1249.0000.000000X( 1967, X2)1647.0000.000000X( 1967, X3)697.00000.000000X( 1967, X4)76368.000.000000X( 1

22、967, X5)30814.000.000000X( 1967, X6)156.00000.000000X( 1968, X1)1187.0000.000000X( 1968, X2)1565.0000.000000X( 1968, X3)680.00000.000000X( 1968, X4)78534.000.000000X( 1968, X5)31915.000.000000X( 1968, X6)127.00000.000000X( 1969, X1)1372.0000.000000X( 1969, X2)2101.0000.000000X( 1969, X3)688.00000.00

23、0000X( 1969, X4)80671.000.000000X( 1969, X5)33225.000.000000X( 1969, X6)207.00000.000000X( 1970, X1)1638.0000.000000X( 1970, X2)2747.0000.000000X( 1970, X3)767.00000.000000X( 1970, X4)82992.000.000000X( 1970, X5)34432.000.000000X( 1970, X6)312.00000.000000X( 1971, X1)1780.0000.000000X( 1971, X2)3156

24、.0000.000000X( 1971, X3)790.00000.000000X( 1971, X4)85229.000.000000X( 1971, X5)35620.000.000000X( 1971, X6)355.00000.000000X( 1972, X1)1833.0000.000000X( 1972, X2)3365.0000.000000X( 1972, X3)789.00000.000000X( 1972, X4)87177.000.000000X( 1972, X5)35854.000.000000X( 1972, X6)354.00000.000000X( 1973,

25、 X1)1978.0000.000000X( 1973, X2)3684.0000.000000X( 1973, X3)855.00000.000000X( 1973, X4)89211.000.000000X( 1973, X5)36652.000.000000X( 1973, X6)374.00000.000000X( 1974, X1)1993.0000.000000X( 1974, X2)3696.0000.000000X( 1974, X3)891.00000.000000X( 1974, X4)90859.000.000000X( 1974, X5)37369.000.000000

26、X( 1974, X6)393.00000.000000X( 1975, X1)2121.0000.000000X( 1975, X2)4254.0000.000000X( 1975, X3)932.00000.000000X( 1975, X4)92421.000.000000X( 1975, X5)38168.000.000000X( 1975, X6)462.00000.000000X( 1976, X1)2052.0000.000000X( 1976, X2)4309.0000.000000X( 1976, X3)955.00000.000000X( 1976, X4)93717.00

27、0.000000X( 1976, X5)38834.000.000000X( 1976, X6)443.00000.000000X( 1977, X1)2189.0000.000000X( 1977, X2)4925.0000.000000X( 1977, X3)971.00000.000000X( 1977, X4)94974.000.000000X( 1977, X5)39377.000.000000X( 1977, X6)454.00000.000000X( 1978, X1)2475.0000.000000X( 1978, X2)5590.0000.000000X( 1978, X3)1058.0000.000000X( 1978, X4)96259.000.000000X( 1978, X5)39856.000.000000X( 1978, X6)550.00000.000000X( 1979, X1)2702.0000.000000X( 1979, X2)6065.00