#/垃圾桶/桶#这是一个shell存档。去掉这行之前的任何东西,然后打开包装#将其保存到文件中并键入“sh文件”。覆盖现有#文件,键入“sh file-c”。您也可以通过以下方式将其作为标准输入#取消shar,或键入“sh“MANIFEST”<<“END_OF_FILE”X文件名存档编号说明X(X)-----------------------------------------------------------X MANIFEST 1此发货清单X生成文件1X自述文件1X通道1 1X ch1eg1/amstat.dat 1号机组X信道1eg1/detail.bak 3X ch1eg1/结束.bak 1X ch1eg1/型号.cc 1X ch1eg1/型号.h 1X ch1eg1/启动日期1X通道1 1X ch5eg1/detail.bak 2X ch5eg1/electa.dat 4号X ch5eg1/electb.dat 1X ch5eg1/electrc1.dat 5号X ch5eg1/electrc2.dat 5X通道5eg1/end3l.bak 2X ch5eg1/结束.bak 2X ch5eg1/型号.cc 3X ch5eg1/型号.h 1X ch5eg1/启动.dat 1X通道1 1X ch6eg1/自述文件1X ch6eg1/detail.bak 6X ch6eg1/结束.bak 2X ch6eg1/启动.dat 1X通道2 1X ch6eg2/细节.bak 3X ch6eg2/详图3l.bak 2X通道2/端3l.bak 1X ch6eg2/ending.bak 1号X ch6eg2/hansen.dat 3号文件X ch6eg2/型号.cc 1X ch6eg2/型号.h 1X ch6eg2/启动3l.dat 1X ch6eg2/启动.dat 1X dcond.cc 1号X显示.cc 2X数据挖掘.cc 1X nlmdl.cc 1号X nlmdl.h 2X nlmdl.opr 3号机组X补丁级别.h 1X realmat.cc 4X realmat.h 2X星盒.cc 1X状态.cc 4X状态.h 1X工具.h 1X通常.h 1文件结束如果测试1664-ne `wc-c<'MANIFEST'`;然后echo shar:\“'MANIFEST'”以错误的尺寸解包!fi(菲涅耳)#“MANIFEST”结束fi(菲涅耳)如果测试-f“生成文件”-a“${1}”!=“-c”;那么echo shar:不会破坏现有文件“'Makefile'”其他的echo shar:解压缩“生成文件”(1692个字符)sed“s/^X//”>“生成文件”<<'END_OF_FILE'XCC=克++XCFLAGS=-OXRMDIR=。。XNLDIR=。。XMDLDIR=。XLDFLAGS=-lmX(X)Xnlmdl:nlmdl.o状态.o显示.o模型.o\X realmat.o dcond.o dsweep.o starbox.oX$(CC)-o nlmdl nlmtl.o状态.o显示.o型号.o\X realmat.o dcond.o dsweep.o starbox.o$(LDFLAGS)X(X)Xnlmdl.o:$(NLDIR)/nlmdl.cc$(NLDIR)/nlmdl.opr\X$(MDLDIR)/model.h$(NLDIR)/nlmdl.h$(NLDIR)/status.h\X$(RMDIR)/usual.h$(RMDIR)/realmat.h$(RM DIR)/tools.hX$(CC)$(CFLAGS)-c-I$(MDLDIR)-I$(NLDIR)-I$X(X)Xmodel.o:$(MDLDIR)/model.cc\X$(MDLDIR)/model.h$(NLDIR)/nlmdl.h$(NLDIR)/status.h\X$(RMDIR)/usual.h$(RMDIR)/realmat.h$(RM DIR)/tools.hX$(CC)$(CFLAGS)-c-I$(MDLDIR)-I$(NLDIR)-I$X(X)Xstatus.o:$(NLDIR)/status.cc\X$(NLDIR)/状态.h\X$(RMDIR)/usual.h$(RMDIR)/realmat.h$(RM DIR)/tools.hX$(CC)$(CFLAGS)-c-I$(MDLDIR)-I$(NLDIR)-I$X(X)Xdisplay.o:$(NLDIR)/display.cc\X$(NLDIR)/状态.h\X$(RMDIR)/usual.h$(RMDIR)/realmat.h$(RM DIR)/tools.hX$(CC)$(CFLAGS)-c-I$(MDLDIR)-I$(NLDIR)-I$(RMDIR)$(NLDIR)/diisplay.CCX(X)Xrealmat.o:$(RMDIR)/realmat.cc\X$(RMDIR)/usual.h$(RMDIR)/realmat.h$(RM DIR)/tools.hX$(CC)$(CFLAGS)-c-I$(RMDIR)$(RMDIR)/realmat.CCX(X)Xstarbox.o:$(RMDIR)/starbox.cc\X$(RMDIR)/usual.h$(RMDIR)/tools.hX$(CC)$(CFLAGS)-c-I$(RMDIR)$(RMDIR)/starbox.CCX(X)Xdcond.o:$(RMDIR)/dcond.cc\X$(RMDIR)/usual.h$(RMDIR)/tools.hX$(CC)$(CFLAGS)-c-I$(RMDIR)$(RMDIR)/dcond.CCX(X)Xdsweep.o:$(RMDIR)/dsweep.cc\X$(RMDIR)/usual.h$(RMDIR)/tools.hX$(CC)$(CFLAGS)-c-I$(RMDIR)$(RMDIR)/dsweep.CCX(X)Xclean公司:X rm-f*.o内核X(X)Xveryclean公司:X rm-f*.o核心nlmdl详细信息.dat ending.dat tmp.dat文件结束如果测试1692-ne `wc-c<'Makefile'`;然后echo shar:\“'Makefile'\”以错误的大小解包!fi(菲涅耳)#“生成文件”结束fi(菲涅耳)如果测试-f“自述”-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“README”其他的echo shar:提取“README”(2856个字符)sed“s/^X//”>“自述”<<'END_OF_FILE'X(X)Xnlmdl是a.Ronald Gallant中统计方法的C++实现,X“非线性统计模型”,纽约:John Wiley and Sons,1987年,使用Xmatrix类,realmat,与之一起分发。矩阵类可以是一些独立的利益。X(X)X该程序计算单变量非线性的最小二乘估计多元的X回归模型和广义最小二乘估计X非线性回归模型。对两者都使用SUR选项。其他选项可以Xbe集用于校正异方差和/或序列的方差估计X相关性。对于可能是隐式形式的同时非线性系统,Xit计算三阶段最小二乘估计、TSLS选项和非线性Xdynamic系统,可能以隐式形式,广义矩量法Xestimates,GMM选项。X(X)X该文档由该文件、两个头文件nlmdl.h和Xrealmat.h,以及目录ch1eg1、ch5eg1、ch6eg1和Xch6eg2。X(X)XA是分布式的,代码设置为在XUnix盒。使用其中一个示例测试代码,例如ch1eg1,cd到Xch1eg1并键入make-f/生成文件。每个目录都有一个文件detail.bakXcreated using nlmdl>detail.bak,可用于比较。对于版本g++的X1.35,通常.h中I/O样式的#undef和#define必须是X设置为GNU。该代码尚未与其他版本一起测试。X(X)X该代码也可以与Borland的Turbo C++编译器一起使用。要使用Turbo C++,Xpermute#undef和#define,以便在通常的.h中选择编译器,复制例如,将ch1eg1的Xcontents发送到工作目录,并使用中的Xa.cc扩展。PRJ文件。使用大内存选项编译。这个Xcode还没有在其他C++编译器上进行过测试。X(X)XUpdates和bug修复程序可通过ftp匿名获取,网址为ccvr1.cc.ncsu.edu目录pub/arg/nlmdl中的X(128.109.21.20)。X(X)X版权所有(C)1990,作者:A.Ronald Gallant。X(X)X使用、复制、修改和分发本软件及其特此免费授予用于任何目的的Xdocumentation,前提是X上述版权声明出现在所有副本中,并且Xnotice和此许可通知出现在支持文档中。X(X)XT本软件按“原样”提供,无任何明示或暗示担保。X(X)X(X)-------------------------------------------------------------------------------X教授A.Ronald Gallant电话:919-737-2531X统计部传真:919-737-7591XNorth Carolina州立大学电子邮件:arg@ccvr1.cc.ncsu.eduXCampus Box 8203 ftp:ccvr1.cc.ncsu.edu(128.109.21.20),XRaleigh NC 27695-8203美国用户:匿名,目录:pub/argX(X)-------------------------------------------------------------------------------文件结束如果测试2856-ne `wc-c<'README'`;然后echo shar:\“'README'\”以错误的大小解包!fi(菲涅耳)#“自述”结束fi(菲涅耳)if测试-d“ch1eg1”;然后echo shar:创建目录“'ch1eg1'\”mkdir“ch1eg1”fi(菲涅耳)如果测试-f'ch1eg1/amstat.dat'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch1eg1/amstat.dat'\”其他的echo shar:提取“'ch1eg1/amstat.dat'\”\(1020个字符\)sed“s/^X//”>'ch1eg1/amstat.dat'<<'END_OF_FILE'X 1 0.98610 1 1 6.28X 2 1.03848 0 1 9.86号X 3 0.95482 1 1 9.11X 4 1.04184 0 1 8.43号X 5 1.02324 1 1 8.11号X 6 0.90475 0 1 1.82X 7 0.96263 1 1 6.58X 8 1.05026 0 1 5.02号X 9 0.98861 1 1 6.52X10 1.03437 0 1 3.75X11 0.98982 1 1 9.86X12 1.01214 0 1 7.31X13 0.66768 1 1 0.47X14 0.55107 0 1 0.07X15 0.96822 1 1 4.07X16 0.98823 0 1 4.61X17 0.59759 1 1 0.17X18 0.99418 0 1 6.99X19 1.01962 1 1 4.39X20 0.69163 0 1 0.39X21 1.04255 1 1 4.73X22 1.04343 0 1 9.42X23 0.97526 1 1 8.90X24 1.04969 0 1 3.02X25 0.80219 1 1 0.77X26 1.01046 0 1 3.31X27 0.95196 1 1 4.51X28 0.97658 0 1 2.65X29 0.50811 1 1 0.08X30 0.91840 0 1 6.11文件结束如果测试1020-ne `wc-c<'ch1eg1/amstat.dat'';然后echo shar:\“'ch1eg1/amstat.dat'”以错误的大小解包!fi(菲涅耳)#“ch1eg1/amstat.dat”结束fi(菲涅耳)如果测试-f'ch1eg1/ending.bak'-a“${1}”!=“-c”;然后echo-shar:不会阻塞现有文件“'ch1eg1/ending.bak'”其他的echo shar:提取“'ch1eg1/ending.bak'”\(2523个字符\)sed“s/^X//”>'ch1eg1/ending.bak'XGallant,“非线性统计模型”,第1章,图5,第35-36页。XSUR什么估算方法?代码SUR、TSLS或GMM。X30观察次数,t=1。。。,n.(名词)。X1方程式数量,即e的尺寸。X0仪器数量,即Z的尺寸。X4参数的数量,即θ的维数。X50高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。同型或异型编码。X0用于var估计的移动平均项MA的数量。Xnone代码none或Parzen,当MA>0时为none是不明智的。X1.000000e-13收敛公差,tol=1.0e-8是合理的。X1.000000e-10反转公差,eps=1.0e-13是合理的Xfull输出多少?无代码、最小代码或完整代码。X4 V的计算秩。变量的X修正除数,修正或未修正。X空行。X空行。X空行。X-2.5889698083278231e-02-0.02588970θ(1)X 1.0156796693196859e+00 1.01567967θ(2)X-1.1156971398981310e+00-1.11569714θ(3)X-5.0490286151154495e-01-0.50490286θ(4)X 1.1729052766553795e-03 0.00117291变量(1,1)X 1.5936130118638265e-04 0.00015936伏(1,1)X-7.8715725112262180e-05-0.00007872伏(2.1)X-1.7710704773277823e-04-0.00017711伏(3.1)X-4.4094836963940918e-05-0.00004409伏(4.1)X-7.8715725112262180e-05-0.00007872伏(1,2)X 9.8762401664581004e-05 0.00009876伏(2,2)X 6.0702486057973514e-04 0.00060702伏(3.2)X-1.8513810538785802e-06-0.00000185伏(4.2)X-1.7710704773277823e-04-0.00017711伏(1,3)X 6.0702486057973514e-04 0.00060702伏(2,3)X 2.6745983237711428e-02 0.02674598伏(3,3)X 2.35621195000099977e-03 0.00235621伏(4.3)X-4.4094836963940918e-05-0.00004409伏(1.4)X-1.8513810538785802e-06-0.00000185伏(2,4)X 2.35621195000099977e-03 0.00235621伏(3,4)X 6.5829241588721008e-04 0.00065829伏(4,4)X-3.1122653130181253e-13-0.00000000 D(1)X 1.0666063107489048e-12 0.00000000 D(2)X 4.6999396801498302e-11 0.00000000 D(3)X 4.1405085967972343e-12 0.00000000 D(4)X 2.5999999999999996e+01 26.00000000对象文件结束如果测试2523-ne `wc-c<'ch1eg1/ending.bak'`;然后echo shar:\“'ch1eg1/ending.bak'”以错误的大小解包!fi(菲涅耳)#“ch1eg1/ending.bak”的末尾fi(菲涅耳)如果测试-f'ch1eg1/model.cc'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch1eg1/model.cc'”其他的echo-shar:提取“'ch1eg1/model.cc'”(1660个字符)sed“s/^X//”>'ch1eg1/model.cc'<<'END_OF_FILE'X#包含“model.h”X(X)Xmodel::model(){}X(X)Xmodel::~model(){}X(X)X(X)Xvoid模型::e(整数t)X(X){X _元素(1,1)X=数据元素(t,2)X-数据elem(t,3)*sθelem(1,1)X-数据元素(t,4)*sθ元素(2,1)X-sθ元素(4,1)*exp(数据元素(t,5)*sθ元(3,1));X}(X)X(X)Xvoid模型::dele(整数t)X(X){X REAL tmp;X(X)X tmp=exp(数据元素(t,5)*sθ元素(3,1));X(X)X dele_t.elem(1,1)=-数据元素(t,3);X dele_t.elem(1,2)=-data.elem(t,4);X dele_t.elem(1,3)=-s.theta.elem;X代表元素(1,4)=-tmp;X}(X)X(X)X(X)Xint模型::initialize()X(X){X(X)X#ifdef GNU_GPP_COMPILERX(X)X#ifdef USE_ATT_STYLE_IO_WITH_GNU(带GNU)X文件buf数据buf;X if(databuf.open(“amstat.dat”,输入)==0){X cerr<<“错误,模型::model,无法打开amstat.dat\n”;X出口(1);X}(X)X istream数据获取(&databuf);X#结尾X(X)X#ifdef USE_GNU_STYLE_IO_WITH_GNU(带GNU)X istream数据获取(“amstat.dat”,io_readonly,a_useonly);X if(!dataget){X cerr<<“错误,模型::model\n”;X出口(1);X}X#结尾X(X)X#结尾X(X)X#ifdef涡轮_CPP_COMPILERX ifstream数据获取(“amstat.dat”);X if(!dataget){X cerr<<“错误,模型::model\n”;X出口(1);X}(X)X#结束X(X)X数据。调整大小(30.5);X(X)X双X;X(X)X代表(整数i=1;i<=30;i++){X代表(整数j=1;j<=5;j++){X数据获取>>X;X数据元素(i,j)=(真实)X;X}(X)X}(X)X(X)X如果(strcmp(s.detail,“full”)==0)cout<“ch1eg1/model.h”<<“END_OF_FILE”X#ifndef __FILE_MODEL_H_SEEN__X#pragma一次X#定义__FILE_MODEL_H_SEEN__1X(X)X#包含“status.h”X(X)Xextern状态;X(X)Xclass模型X(X){X(X)X保护:X realmat数据;X(X)Xpublic公司:X模型();X~模型();X(X)X realmat(_t);X realmat dele_t;X实数Z_t;X(X)X int initialize();X int terminate();X(X)X空位e(整数t);X无效代表(整数t);X空隙Z(整数t);X(X)十} ;X(X)X#结尾文件结束如果试验455-ne‘wc-c<‘ch1eg1/h型’;然后echo shar:\“'ch1eg1/model.h'”打开时尺寸错误!fi(菲涅耳)#“ch1eg1/model.h”结束fi(菲涅耳)如果test-f'ch1eg1/starting.dat'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch1eg1/starting.dat'\”其他的echo shar:提取“'ch1eg1/starting.dat'”\(1097个字符\)sed“s/^X//”>'ch1eg1/starting.dat'XGallant,“非线性统计模型”,第1章,图5,第35-36页。XSUR什么估算方法?代码SUR、TSLS或GMM。X30观测次数,t=1。。。,n.(名词)。X1方程式数量,即e的尺寸。X0仪器数量,即Z的尺寸。X4参数的数量,即θ的维数。X50高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。同型或异型编码。X0用于var估计的移动平均项MA的数量。Xnone代码none或Parzen,当MA>0时为none是不明智的。X1.000000e-13收敛公差,tol=1.0e-8是合理的。X1.000000e-10反转公差,eps=1.0e-13是合理的Xfull输出多少?无代码、最小代码或完整代码。X空行。X空行。X空行。X空行。X空行。X-0.04866号X 1.03884号X-0.73792毫米X-0.51362号文件结束如果测试1097-ne `wc-c<'ch1eg1/starting.dat'';然后echo shar:\“'ch1eg1/starting.dat'”以错误的大小解包!fi(菲涅耳)#“ch1eg1/starting.dat”结束fi(菲涅耳)如果测试-d“ch5eg1”;然后echo shar:创建目录“'ch5eg1'\”mkdir“ch5eg1”fi(菲涅耳)如果测试-f'ch5eg1/electb.dat'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch5eg1/electb.dat'\”其他的echo shar:提取“'ch5eg1/electb.dat'\”\(612个字符\)sed“s/^X//”>'ch5eg1/electb.dat'<<'END_OF_FILE'X 1 1.06 2.86 3.90X 2 1.78 2.86 3.90X 3 1.06 3.90 3.90X 4 1.78 3.90 3.90X 5 1.37 3.34 5.06X 6 1.06 2.86 6.56号X 7 1.78 2.86 6.56X 8 1.06 3.90 6.56X 9 1.78 3.90 6.56文件结束如果测试612-ne“wc-c<‘ch5eg1/electb.dat’”;然后echo shar:\“'ch5eg1/electb.dat'”以错误的大小解包!fi(菲涅耳)#“ch5eg1/electb.dat”结束fi(菲涅耳)如果测试-f'ch5eg1/model.h'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch5eg1/model.h'”其他的echo shar:提取“'ch5eg1/model.h'”\(465个字符\)sed“s/^X//”>'ch5eg1/model.h'<<'END_OF_FILE'X#ifndef __FILE_MODEL_H_SEEN__X#pragma一次X#定义__FILE_MODEL_H_SEEN__1X(X)X#包含“status.h”X(X)Xextern状态;X(X)Xclass模型X(X){X(X)X保护:X realmat数据;X(X)Xpublic公司:X模型();X~模型();X(X)X实际材料(_t);//e_t是2乘1X realmat dele_t;//dele_是2乘8X realmat Z_t;//Z_t是17乘1X(X)X int initialize();X int terminate();X(X)X空隙e(整数t);X无效删除(整数t);X空隙Z(整数t);X(X)十} ;X(X)X#结尾文件结束如果试验465-ne‘wc-c<‘ch5eg1/h型’;然后echo shar:\“'ch5eg1/model.h'”打开时尺寸错误!fi(菲涅耳)#“ch5eg1/型号.h”结束fi(菲涅耳)如果测试-f'ch5eg1/starting.dat'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch5eg1/starting.dat'\”其他的echo shar:提取“'ch5eg1/starting.dat'”\(1113个字符\)sed“s/^X//”>'ch5eg1/starting.dat'<<'END_OF_FILE'XGallant,“非线性统计模型”,第5章,图3,第324-326页XSUR什么估算方法?代码SUR、TSLS或GMM。X224观察次数,t=1。。。,n.(名词)。X2方程式数量,即e的尺寸。X0仪器数量,即Z的尺寸。X8参数的数量,即θ的尺寸。X50高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。同型或异型编码。X0用于var估计的移动平均项MA的数量。Xnone代码none或Parzen,当MA>0时为none是不明智的。X1.000000e-13收敛公差,tol=1.0e-8是合理的。X1.000000e-10反转公差,eps=1.0e-13是合理的Xminimal输出多少?无代码、最小代码或完整代码。X空行。X空行。X空行。X空行。X空行。X-2.90倍X-1.30号X 0.82X 0.36X-1.50倍X-1.00号X-0.03号X-0.47倍文件结束如果测试1113-ne `wc-c<'ch5eg1/starting.dat'';然后echo shar:\“'ch5eg1/starting.dat'”以错误的大小解包!fi(菲涅耳)#“ch5eg1/starting.dat”结束fi(菲涅耳)如果测试-d“ch6eg1”;然后echo shar:创建目录“'ch6eg1'\”mkdir“ch6eg1”fi(菲涅耳)如果测试-f'ch6eg1/README'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg1/README'\”其他的echo shar:提取“'ch6eg1/README'”(90个字符)sed“s/^X//”>'ch6eg1/README'<<'END_OF_FILE'X(X)X执行此操作以获取丢失的文件:X(X)X cp/ch5eg1/el*.日期。X坐标/ch5eg1/型号*。X(X)文件结束如果测试90-ne `wc-c<'ch6eg1/README'`;然后echo shar:\“'ch6eg1/README'\”打开时尺寸错误!fi(菲涅耳)#“ch6eg1/README”结束fi(菲涅耳)如果test-f'ch6eg1/starting.dat'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg1/starting.dat'\”其他的echo shar:提取“'ch6eg1/starting.dat'”\(1113个字符\)sed“s/^X//”>'ch6eg1/starting.dat'<<'END_OF_FILE'XGallant,“非线性统计模型”,第6章,图1,第436页。XTSLS什么估算方法?代码SUR、TSLS或GMM。X224观测次数,t=1。。。,n.(名词)。X2方程式数量,即e的尺寸。X17仪器数量,即Z的尺寸。X8参数的数量,即θ的维数。X20高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。同型或异型编码。X0用于var估计的移动平均项MA的数量。Xnone代码none或Parzen,当MA>0时为none是不明智的。X1.000000e-08收敛公差,tol=1.0e-8是合理的。X1.000000e-10反转公差,eps=1.0e-13是合理的Xfull输出多少?无代码、最小代码或完整代码。X空行。X空行。X空行。X空行。X空行。X-2.98年X-1.16年X 0.787X 0.353X-1.51号X-1.0X 0.054X-0.474倍文件结束如果测试1113-ne `wc-c<'ch6eg1/starting.dat'';然后echo shar:\“'ch6eg1/starting.dat'”以错误的大小解包!fi(菲涅耳)#“ch6eg1/starting.dat”结束fi(菲涅耳)如果测试-d“ch6eg2”;然后echo shar:创建目录“'ch6eg2'\”mkdir“ch6eg2”fi(菲涅耳)如果测试-f'ch6eg2/end3l.bak'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg2/end3l.bak'\”其他的echo-shar:提取“'ch6eg2/end3l.bak'”(2101个字符)sed“s/^X//”>'ch6eg2/end3l.bak'XGallant,“非线性统计模型”,第6章,图2,第449-450页。XGMM什么估算方法?代码SUR、TSLS或GMM。X238观察次数,t=1。。。,n.(名词)。X1方程式数量,即e的尺寸。X3仪器数量,即Z的尺寸。X2参数的数量,即θ的维数。X20高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。异方差代码同方差或异方差。X3用于var估计的移动平均项MA的数量。XParzen Code none或Parzen,当MA>0时为none是不明智的。X1.000000e-05收敛公差,tol=1.0e-8是合理的。X1.000000e-13反转公差,eps=1.0e-13合理Xfull输出多少?无代码、最小代码或完整代码。X2 V的计算秩。Xuncorrected var的除数,修正或未修正。X空行。X空行。X空行。X-8.5525378703973232e-01-0.85525379θ(1)X 9.980456031317236e-01 0.99804560θ(2)X 4.397820854904504e-01 0.43978221变量(1,1)X 4.4050039895600307e-01 0.44050040变量(2,1)X 4.3420956263910121e-01 0.43420956 var(3.1)X 4.4050039895600307e-01 0.44050040变量(1,2)X 4.4122700259798669e-01 0.44122700 var(2,2)X 4.349363977747488e-01 0.43493664 var(3,2)X 4.3420956263910121e-01 0.43420956 var(1,3)X 4.349363977747488e-01 0.43493664变量(2,3)X 4.2967081964723475e-01 0.42967082 var(3,3)X 3.6066325678816527e+00 3.60663257伏(1,1)X-7.1576798102674808e-03-0.00715768伏(2,1)X-7.1576798102674808e-03-0.00715768伏(1,2)X 2.1170843618445885e-05 0.00002117伏(2,2)X 2.6747016119385693e-06 0.00000267 D(1)X-1.7073275785211286e-08-0.00000002 D(2)X 1.0861341198092420e+00 1.08613412对象文件结束如果测试2101-ne `wc-c<'ch6eg2/end3l.bak'';然后echo shar:\“'ch6eg2/end3l.bak'”打开时大小不正确!fi(菲涅耳)#“ch6eg2/end3l.bak”的末尾fi(菲涅耳)如果测试-f'ch6eg2/ending.bak'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg2/ending.bak'\”其他的echo shar:提取“'ch6eg2/ending.bak'”\(2101个字符\)sed“s/^X//”>'ch6eg2/ending.bak'XGallant,“非线性统计模型”,第6章,图2,第449-450页。XGMM什么估算方法?代码SUR、TSLS或GMM。X238观察次数,t=1。。。,n.(名词)。X1方程式数量,即e的尺寸。X3仪器数量,即Z的尺寸。X2参数的数量,即θ的维数。X20高斯-牛顿迭代的上限。X1数值var迭代,ivar=0表示无。同型或异型编码。X0用于var估计的移动平均项MA的数量。XParzen Code none或Parzen,当MA>0时为none是不明智的。X1.000000e-05收敛公差,tol=1.0e-8是合理的。X1.000000e-13反转公差,eps=1.0e-13合理Xfull输出多少?无代码、最小代码或完整代码。X2 V的计算秩。Xuncorrected var的除数,修正或未修正。X空行。X空行。X空行。X-1.0335217392035390e+00-1.03352174θ(1)X 9.9825572364240200e-01 0.99825572θ(2)X 4.0582168346665826e-01 0.40582168变量(1,1)X 4.0643410965960231e-01 0.40643411变量(2,1)X 3.9873658267692480e-01 0.39873658变量(3,1)X 4.0643410965960231e-01 0.40643411变量(1,2)X 4.0705487414860658e-01 0.40705487变量(2,2)X 3.9936326231932223e-01 0.39936326变量(3,2)X 3.9873658267692480e-01 0.39873658 var(1,3)X 3.9936326231932223e-01 0.39936326变量(2,3)X 3.9272327136336027e-01 0.39272327变量(3,3)X 3.5800932072477121e+00 3.58009321伏(1,1)X-7.2126817691390913e-03-0.00721268伏(2.1)X-7.2126817691390913e-03-0.00721268伏(1,2)X 2.0603180900457831e-05 0.00002060伏(2,2)X-9.1931859179544446e-08-0.00000009 D(1)X 1.8525999753608607e-10 0.00000000 D(2)X 1.0569217148505041e+00 1.056921171对象文件结束如果测试2101-ne `wc-c<'ch6eg2/ending.bak'';然后echo shar:\“'ch6eg2/ending.bak'”以错误的大小解包!fi(菲涅耳)#“ch6eg2/ending.bak”的末尾fi(菲涅耳)如果测试-f'ch6eg2/model.cc'-a“${1}”!=“-c”;那么echo shar:不会破坏现有文件“'ch6eg2/model.cc'”其他的echo shar:提取“'ch6eg2/model.cc'”\(1679个字符\)sed“s/^X//”>'ch6eg2/model.cc'<<'END_OF_FILE'X#包含“model.h”X(X)Xmodel::model(){}X(X)Xmodel::~model(){}X(X)X(X)Xvoid模型::e(整数t)X(X){X实数a,b,yta,yt,xt;X(X)X整数i=t+2;X(X)X yt=(nds[i]/pop[i])/(nds[1]/pop[1]);X xt=(1.0+nyse[i])*(dflat[i-1]/dflat[i];X(X)X a=s.θ[1];X b=s.θ[2];X(X)X yta=功率(yt,a);X(X)X e_t[1]=b*yta*xt-1.0;X}X(X)Xvoid模型::dele(整数t)X(X){X真实a,b,yta,yt,xt;X(X)X整数i=t+2;X(X)X yt=(nds[i]/pop[i])/(nds[1]/pop[1]);X xt=(1.0+nyse[i])*(dflat[i-1]/dflat[i];X(X)X a=s.θ[1];X b=s.θ[2];X(X)X yta=功率(yt,a);X(X)X dele_t.elem(1,1)=b*log(yt)*yta*xt;X dele_t.elem(1,2)=yta*xt;X}(X)X(X)X(X)Xint模型::initialize()X(X){X整数i;X双tmp;X(X)X#ifdef GNU_GPP_COMPILERX#ifdef USE_ATT_STYLE_IO_WITH_GNU(带GNU)X(X)X文件缓冲区cdata_buf;X if(cdata_buf.open(“hansen.dat”,输入)==0){X cerr<<“错误,模型::model,无法打开hansen.dat.\n”;X出口(1);X}(X)X istream cdata(&cdata_buf);X(X)X#结尾X(X)X#ifdef USE_GNU_STYLE_IO_WITH_GNU(带GNU)X(X)X istream cdata(“hansen.dat”,io_readonly,a_useonly);X如果(!cdata){X cerr<<“错误,模型::model,无法打开hansen.dat.\n”;X出口(1);X}(X)X(X)X#结尾X(X)X#结尾X(X)X#ifdef涡轮_CPP_COMPILERX(X)X ifstream cdata(“hansen.dat”);X if(!cdata){X cerr<<“错误,model::model,无法打开hansen.dat。\n”;X出口(1);X}(X)X(X)X#结尾X(X)X代表(i=1;i<=240;i++){X cdata>>tmp;nds[i]=(真实)tmp;X cdata>>tmp;pop[i]=(真实)tmp;X cdata>>tmp;nyse[i]=(真实)tmp;X cdata>>tmp;dflat[i]=(真实)tmp;X}(X)X(X)X返回1;X}(X)X(X)Xvoid模型::Z(整数t)X(X){X整数i=t+2;X Z_t[1]=1.0;X Z_t[2]=(nds[i-1]/pop[i-1])/(nds[i-2]/pop[i-2]);X Z_t[3]=(1.0+nyse[i-1])*(dflat[i-2]/dflat[i-1]]);X}(X)X(X)Xint模型::terminate(){return 0;}X(X)X(X)X(X)X(X)X(X)文件结束如果测试1679-ne `wc-c<'ch6eg2/model.cc'';然后echo shar:\“'ch6eg2/model.cc'”打开时尺寸错误!fi(菲涅耳)#“ch6eg2/model.cc”结束fi(菲涅耳)如果测试-f'ch6eg2/model.h'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg2/model.h'”其他的echo shar:提取“'ch6eg2/model.h'”\(536个字符\)sed“s/^X//”>'ch6eg2/model.h'<<'END_OF_FILE'X#ifndef __FILE_MODEL_H_SEEN__X#pragma一次X#定义__FILE_MODEL_H_SEEN__1X(X)X#包含“status.h”X(X)Xextern状态s;X(X)Xclass模型X(X){X(X)X保护:X实数[256];X REAL流行音乐[256];X真实纽约证券交易所[256];X实dflat[256];X(X)Xpublic公司:X realmat e_t;//e_t是1乘1X realmat dele_t;//dele_t是1乘2X realmat Z_t;//Z_t是3乘1X(X)X模型();X~模型();X(X)X int initialize();X int terminate();X(X)X空隙e(整数t);X无效删除(整数t);X空隙Z(整数t);X(X)十} ;X(X)X#结尾文件结束如果试验536-ne‘wc-c<‘ch6eg2/h型’;然后echo shar:\“'ch6eg2/model.h'”打开时尺寸错误!fi(菲涅耳)#“ch6eg2/model.h”结束fi(菲涅耳)如果测试-f'ch6eg2/start3l.dat'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg2/start3l.dat'\”其他的echo shar:提取“'ch6eg2/start3l.dat'\”\(1068个字符\)sed“s/^X//”>'ch6eg2/start3l.dat'XGallant,“非线性统计模型”,第6章,图2,第449-450页。XGMM什么估算方法?代码SUR、TSLS或GMM。X238观察次数,t=1。。。,n.(名词)。X1方程的数量,即e的维数。X3仪器数量,即Z的尺寸。X2参数的数量,即θ的维数。X20高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。异方差代码同方差或异方差。X3用于var估计的移动平均项MA的数量。XParzen Code none或Parzen,当MA>0时为none是不明智的。X1.000000e-05收敛公差,tol=1.0e-8是合理的。X1.000000e-13反转公差,eps=1.0e-13合理Xfull输出多少?无代码、最小代码或完整代码。X空行。X空行。X空行。X空行。X空行。X-0.4倍X 0.9倍文件结束如果测试1068-ne“wc-c<”ch6eg2/start3l.dat“;然后echo shar:\“'ch6eg2/start3l.dat'”以错误的大小解包!fi(菲涅耳)#“ch6eg2/start3l.dat”结束fi(菲涅耳)如果测试-f'ch6eg2/starting.dat'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg2/starting.dat'\”其他的echo shar:提取“'ch6eg2/starting.dat'”\(1068个字符\)sed“s/^X//”>'ch6eg2/starting.dat'<<'END_OF_FILE'XGallant,“非线性统计模型”,第6章,图2,第449-450页。XGMM什么估算方法?代码SUR、TSLS或GMM。X238观察次数,t=1。。。,n.(名词)。X1方程式数量,即e的尺寸。X3仪器数量,即Z的尺寸。X2参数的数量,即θ的维数。X20高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。Xhomoskedastic代码同源或异基。X0用于var估计的移动平均项MA的数量。XParzen Code none或Parzen,当MA>0时为none是不明智的。X1.000000e-05收敛公差,tol=1.0e-8是合理的。X1.000000e-13反转公差,eps=1.0e-13合理Xfull输出多少?无代码、最小代码或完整代码。X空行。X空行。X空行。X空行。X空行。X-0.4倍X 0.9倍文件结束如果测试1068-ne `wc-c<'ch6eg2/starting.dat'';然后echo shar:\“'ch6eg2/starting.dat'”以错误的大小解包!fi(菲涅耳)#“ch6eg2/starting.dat”结束fi(菲涅耳)如果测试-f“dcond.cc”-a“${1}”!=“-c”;那么echo shar:不会破坏现有文件“'dcond.cc'\”其他的echo shar:提取“'dcond.cc'”(1599个字符)sed“s/^X//”>'第二个cc'<<'END_OF_FILE'X(X)/*XName dcond-提高了当对角线为a的X元素的数量级不同。X(X)XUsage#包括“usual.h”X#包括“tools.h”X无效数据端(实数*a,整数n,实数*s,整数isw)X(X)X刀具类型.hX(X)X描述a是一个存储的对称正定n乘n矩阵X列式,没有未使用的空格和[0]中的第一个元素;即。,X代表(j=1;j<=n;j++)代表(i=1;i<=n,i++)aij=a[n*(j-1)+(i-1)];X将遍历矩阵,aij是第X行和第j列。s是长度为n的功向量。X预期的调用顺序是:X数据端(a,n,s,0);X ier=dsewe(a,n,eps);X d键(a,n,s,1);X用法X dcond(a,n,s,0);X将返回中的相关矩阵和标准错误当a是方差-方差矩阵时,X为s。X(X)XRemark dscond.c是将dcond.f转换为c++的。X(X)X返回值无。X(X)XFunctions库:(无)调用的子库:(无)X(X)*/X(X)X#包括“usual.h”X#包括“tools.h”X(X)Xvoid dcond(实数*a,整数n,实数*s,整数isw)X(X){X整数i;X REAL*si,*sj,*top,*aij;X(X)如果(isw==1)转到L20,则为X;X(X)X代表(i=0;i<n;++i)X s[i]=sqrt(a[i+i*n]);X(X)X L20:X(X)X sj=秒;X顶部=s+n;X aij=a;X(X)X while(sj<顶部){X如果(*sj!=0.){X si=s;X while(si<顶部)X如果(*si!=0.)*aij++/=(*si++)*(*sj);X}(X)X sj++;X}(X)X(X)X返回;X(X)X}(X)文件结束如果测试1599-ne `wc-c<'cond.cc'';然后echo shar:“'dcond.cc'”解压缩时大小错误!fi(菲涅耳)#“dcond.cc”结尾fi(菲涅耳)如果测试-f'dsweep.cc'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'dspew.cc'\”其他的echo shar:提取“'dspew.cc'”(4560个字符)sed“s/^X//”>'dsweep.cc'<<'END_OF_FILE'X(X)/*XName dsweep-通过X对角线扫描。X(X)XUsage#包括“usual.h”X#包括“tools.h”X整数数据扫描(REAL*a,INTEGER n,REAL eps)X(X)X刀具类型.hX(X)X描述a是一个存储的对称正定n乘n矩阵X列式,没有未使用的空格和[0]中的第一个元素;即。,X代表(j=1;j<=n;j++)代表(i=1;i<=n,i++)aij=a[n*(j-1)+(i-1)];X将遍历矩阵,aij是第X行和第j列。返回时,包含X按列存储。eps是用作相对值的输入常量退化秩检验中的X公差。合理的价值eps的X为1.e-13。X(X)XRemark dsweep.cc是dsweep.f翻译成C++的。X(X)XReference Schatzoff,M.、R.Tsao和S.Fienberg(1968)“高效X所有可能回归的计算,“技术计量学10,X 769-779。X(X)X返回值返回值ier是一个错误参数,编码如下:X ier=0,无误,a的秩为n;ier>0,a是单数,秩a的X是n-ier。如果ier>0,则dsweet返回g-inverse。X(X)XFunctions库:(无)调用的子库:(无)X(X)*/X(X)X#包括“usual.h”X#包括“工具.h”X(X)X(X)XINTEGER数据流(REAL*a,INTEGER n,REAL eps)X(X){X整数i,j,k,kk,ier;X REAL测试,d,阴性,akj,akk,tol;X REAL*p_aij、*p_aji、*p_aik、*p _a1j、*p _akk、*top;X(X)X a--;//调整偏移量以允许Fortran样式索引。X(X)X公差=0.0;X(X)X代表(i=1;i<=n;i++){X检验=a[n*(i-1)+i];X如果(test>tol)tol=测试;X}(X)X(X)X tol=tol*eps;X ier=0;X(X)X(X)X(X)/*十、目的是:X(X)X代表(k=1;k<=n;k++){X(X)X kk=n*(k-1)+k;X(X)X akk=a[kk];X(X)X如果(akk<tol){X(X)X ier=ier+1;X(X)X代表(j=k;j<=n;j++)//上三角形的第k行零点X a[n*(j-1)+k]=0.0;X(X)X代表(i=1;i<=k-1;i++)//上三角的k列零点X a[kk-i]=0.0;X(X)X}(X)X else{//扫描X(X)X d=1.0/akk;X(X)X表示(j=1;j<k;j++){X akj=-a[n*(k-1)+j];X代表(i=1;i<=j;i++){X aik=a[n*(k-1)+i];X a[n*(j-1)+i]-=aik*akj*d;X}(X)X}(X)X(X)X代表(j=k+1;j<=n;j++){X akj=a[n*(j-1)+k];X表示(i=1;i<k;i++){X aik=a[n*(k-1)+i];X a[n*(j-1)+i]-=aik*akj*d;X}(X)X代表(i=k+1;i<=j;i++){X aik=a[n*(i-1)+k];X a[n*(j-1)+i]-=aik*akj*d;X}(X)X}(X)X(X)X代表(j=k;j<=n;j++)Xα[n*(j-1)+k]*=d;X(X)X代表(i=1;i“nlmdl.cc”<<“END_OF_FILE”X/*----------------------------------------------------------------------------X(X)Xnlmdl:nlmdl.ccX(X)Xnlmdl是a.Ronald中统计方法的C++实现XGallant,“非线性统计模型”,纽约:John Wiley and Sons,X1987,ISBN 0-471-80260-3,使用分布式的矩阵类realmat头文件nlmdl.h和realmat.h描述了Xprogram和矩阵类。X(X)X版权(C)1990。X(X)XA公司。罗纳德·格兰特XP。邮政信箱5513XRaleigh北卡罗来纳州27650-5513XUSA公司X(X)X使用、复制、修改和分发本软件及其特此免费授予用于任何目的的Xdocumentation,前提是X上述版权声明出现在所有副本中,并且X版权声明和本许可声明出现在支持X文档。X(X)XT本软件按“原样”提供,无任何明示或暗示担保。X(X)X----------------------------------------------------------------------------*/X#包含“nlmdl.h”X(X)Xtypedef无效(*ZERO_ARGUMENT_NL_OPR)();Xtypedef void(*ONE_ARGUMENT_NL_OPR)(int);X(X)X状态;X模型m;X(X)XZERO_ARGUMENT_NL_OPR OPR_mgn;XZERO_ARGUMENT_NL_OPR OPR_obj;XONE_ARGUMENT_NL_OPR操作变量;XZERO_ARGUMENT_NL_OPR OPR_V;X(X)X#包括“nlmdl.opr”X(X)Xint main()X(X){X字符消息[MAX_STATUS_LINE];X(X)X s.from(s.starting);X(X)X重复:X(X)X如果(m.initialize()==0)退出(0);X(X)X s.from(s.starting);X(X)X if(strcmp(s.detail,“none”)!=0) {X冲刺(消息,“/nlmdl%s//_”,版本);X cout<<星空(msg);X s.显示(启动);X cout<<starbox(“/起始θ//_”)<<s.theta;X}(X)X(X)X如果(strcmp(s.method,“SUR”)==0){X opr_mgn=&SUR_mgn;X opr_obj=&SUR_obj;X opr_var=&SUR_var;X opr_V=&SUR_V;X}(X)X else if(strcmp(s.method,“TSLS”)==0){X(X)X if(strcmp(s.vartype,“异方差”)==0||s.MA>0){X cerr<<“\n错误,nlmdl,在这种情况下使用GMM,而不是TSLS。\n\n”;X出口(1);X}(X)X opr_mgn=&TSLS_mgn;X opr_obj=&TSLS_obj;X opr_var=&SUR_var;X opr_V=&TSLS_V;X(X)X qZ.调整大小(s.M*s.K,1);X QZ.调整大小(s.M*s.K,s.p);X}(X)X else如果(strcmp(s.method,“GMM”)==0){X opr_mgn=&GMM_mgn;X opr_obj=&GMM_obj;X opr_var=&GMM_var;X opr_V=&GMM_V;X(X)X qZ.调整大小(s.M*s.K,1);X QZ.调整大小(s.M*s.K,s.p);X}(X)X其他{X cerr<<“错误,nlmdl,s.method设置错误。\n”;X出口(1);X}(X)X(X)X m.e_t.调整大小(s.m,1);X m.dele_t.调整大小(s.m,s.p);X,如果(s.K>0)m.Z_t.调整大小(s.K,1);X(X)X s.V.resize(s.p、s.p、REAL_ZERO);X表示(整数i=1;i<=s.p;i++)s.V.elem(i,i)=1.0;X s.D.调整大小(s.p,1,REAL_ZERO);X(X)X代表(整数var_loop=0;var_loop<=s.ivar;var_loop++){X(X)X(*opr_var)(var_loop);X(X)X如果(strcmp(s.detail,“full”)==0){X cout<<“\nvar_loop”<<var_loop<<“\n”;X s.显示(VAR_ITERATE);X}(X)X(X)X代表(整数theta_loop=0;theta_loop<=s.itheta;theta_ loop++){X(X)X(*opr_mgn)();X(X)X if(strcmp(s.detail,“full”)==0){X cout<<“\ntheta_loop”<<theta_loop<<“\n”;X s.显示(THETA_ITERATE);X}X(X)X int rv=线路搜索(&msg[0]);X(X)X if(strcmp(s.detail,“full”)==0)cout<<msg;X,如果(rv!=0)断裂;X}(X)X}(X)X(X)X(*opr_mgn)();X(*opr_V)();X(X)X s.to(s.ending);X(X)X if(strcmp(s.detail,“none”)!=0)s.display(终止);X(X)X如果(m.terminate()!=0)转到重复;X(X)X出口(0);X}(X)文件结束如果测试3428-ne `wc-c<'nlmdl.cc'`;然后echo shar:“'nlmdl.cc'”未打包,大小错误!fi(菲涅耳)#“nlmdl.cc”的结尾fi(菲涅耳)如果测试-f'pathlevel.h'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“patchlevel.h”其他的echo shar:提取“patchlevel.h”(1006个字符)sed“s/^X//”>'patchlevel.h'<<'END_OF_FILE'X/*--------------------------------------------------------------------------X(X)Xnlmdl:patchlevel.hX(X)Xnlmdl是a.Ronald中统计方法的C++实现XGallant,“非线性统计模型”,纽约:John Wiley and Sons,X1987,ISBN 0-471-80260-3,使用分布式的矩阵类realmat头文件nlmdl.h和realmat.h描述了Xprogram和矩阵类。X(X)X版权所有(C)1990。X(X)XA公司。罗纳德·格兰特XP。邮政信箱5513XRaleigh北卡罗来纳州27650-5513XUSA公司X(X)X使用、复制、修改和分发本软件及其特此免费授予用于任何目的的Xdocumentation,前提是X上述版权声明出现在所有副本中,并且X版权声明和本许可声明出现在支持X文档。X(X)XT本软件按“原样”提供,无任何明示或暗示担保。X(X)X-------------------------------------------------------------------------*/X(X)X#定义PATCHLEVEL 1X(X)X(X)文件结束如果测试1006-ne“wc-c<”atchlevel.h“”;然后echo shar:“'patchlevel.h'”打开时大小错误!fi(菲涅耳)#“pathlevel.h”结束fi(菲涅耳)如果测试-f“starbox.cc”-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'starbox.cc'\”其他的echo shar:提取“'starbox.cc'”(2486个字符)sed“s/^X//”>'starbox.cc'<<'END_OF_FILE'X(X)/*XName星形框-将标题放在一个中心的星形框中。X(X)XUsage#包括“usual.h”X#包括“tools.h”X char*starbox(char*title,char终止符);X(X)X刀具类型.hX(X)XDescription title是一个字符串。该字符串由X条连接线。一行由0-68个字符组成X后跟斜线(默认值)。例如:X cout<<starbox(“//第1行/第2行/第3行///”);X cout<<starbox(“\n\nline 1\nline 2\nline 3\n\n\n”,'\n');X(X)X返回值字符串。X(X)XFunctions库:strcat、strchr、strlen调用的工具:(无)X子库:(无)X(X)*/X(X)X#包括“usual.h”X#包括“工具.h”X(X)Xchar*starbox(char*title,char终止符)X(X){X字符边框[81],制表符[81]、工件[256]、中间[81];X char*line,*begin,*end,*t;X int i,pad,mpad,linesize;X long int length,cum_length;X(X)X char starbox_str_pointer[1024];X long int starbox_str_length=1024;X(X)X累计长度=0;X*starbox_str_pointer=“\0”;X(X)X代表(i=0;i<70;i++)X边框[i]=“*”;X边框[70]=“\n”;border[71]=“\0”;X(X)X((行数<72)||(行数>133))?(线条尺寸=133):(线条尺寸=线条尺寸);X(X)X焊盘=(linesize-72)/2+1;X代表(i=0;i68)return“错误,starbox,行太长。”;X中间[0]=“*”;X mpad=(68-长度)/2+1;X代表(i=1;i<=mpad;i++)中间[i]='';X t=中间[mpad];X while(begin<end)*t++=*begin++;X开始++;X代表(i=mpad+长度;i<=68;i++)中间[i]='';X中间[69]=“*”;middle[70]='\n';中间[71]=“\0”;X功[0]=“\0”;X线=strcat(工作,制表符);X线=strcat(工作,中间);X(X)X(X)X如果(cum_length<starbox_str_length)X strcat(starbox_str_pointer,line);X其他X return“错误,starbox,标题过长。”;X(X)X}X(X)X功[0]=“\0”;X线=strcat(工作,制表符);X线=strcat(工作,边界);X(X)X如果(cum_length<starbox_str_length)X strcat(starbox_str_pointer,line);X其他X return“错误,starbox,标题过长。”;X(X)X返回starbox_str_pointer;X(X)X}X(X)X(X)X(X)文件结束如果测试2486-ne `wc-c<'starbox.cc'`;然后echo shar:\“'starbox.cc'\”未打包,大小不正确!fi(菲涅耳)#“starbox.cc”结束fi(菲涅耳)如果测试-f“status.h”-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“status.h”其他的echo shar:提取“status.h”(2046个字符)sed“s/^X//”>“状态.h”<<'END_OF_FILE'X/*----------------------------------------------------------------------------X(X)Xnlmdl:状态.hX(X)Xnlmdl是a.Ronald中统计方法的C++实现XGallant,“非线性统计模型”,纽约:John Wiley and Sons,X1987,ISBN 0-471-80260-3,使用分布式的矩阵类realmat头文件nlmdl.h和realmat.h描述了Xprogram和矩阵类。X(X)X版权所有(C)1990。X(X)XA公司。罗纳德·加兰特XP。邮政信箱5513XRaleigh北卡罗来纳州27650-5513XUSA公司X(X)X使用、复制、修改和分发本软件及其特此免费授予用于任何目的的Xdocumentation,前提是X上述版权声明出现在所有副本中,并且X版权声明和本许可声明出现在支持X文档。X(X)XT本软件按“原样”提供,无任何明示或暗示担保。X(X)X----------------------------------------------------------------------------*/X/*status是nlmdl使用的类*/X(X)X#ifndef __FILE_STATUS_H_SEEN__X#pragma一次X#定义__FILE_STATUS_H_SEEN__1X(X)X#包括“usual.h”X#包括“realmat.h”X#包括“工具.h”X(X)X(X)X(X)X#定义MAX_STATUS_LINE 256X#define STARTING_FILENAME“开始.dat”X#定义ENDING_FILENAME“ENDING.dat”X(X)氙气显示模式{START_UP,VAR_ITERATE,THETA_ITERATE,TERMINATION};X(X)Xclass状态X(X){Xpublic公司:X状态();X~状态();X(X)X字符开关[MAX_STATUS_LINE];X字符方法[MAX_STATUS_LINE];X整数n;X整数M;X整数K;X整数p;X整数itheta;X整数ivar;X字符变量类型[MAX_STATUS_LINE];X整数MA;X字符权重[MAX_STATUS_LINE];X真实公差;X REAL每股收益;X字符细节[MAX_STATUS_LINE];X整数秩;X字符df[MAX_STATUS_LINE];X realmatθ;X realmat变量;X realmat V;X实际材料D;X真实对象;X(X)X字符*开始;X字符*结尾;X(X)X int from(char*filename);X int to(char*filename);X int显示(display_mode);X(X)十} ;X(X)X#结尾X(X)文件结束如果测试2046-ne‘wc-c<‘status.h’’;然后echo shar:“'status.h'”大小错误!fi(菲涅耳)#“status.h”结束fi(菲涅耳)如果测试-f'tools.h'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'tools.h'\”其他的echo shar:提取“tools.h”(290个字符)sed“s/^X//”>'tools.h'<<'END_OF_FILE'X#ifndef __FILE_TOOLS_H_SEEN__X#pragma一次X#定义__FILE_TOOLS_H_SEEN__1X(X)X#包括“usual.h”X(X)Xchar*starbox(char*title,char终止符='/');X(X)Xvoid dcond(实数*a,整数n,实数*s,整数isw);XINTEGER dswew(REAL*a,INTEGER n,REAL eps);X(X)X#结尾/*__FILE_TOOLS_H_SEEN__*/文件结束如果测试290-ne“wc-c<”tools.h“”;然后echo-shar:\“'tools.h'\”解压缩时大小错误!fi(菲涅耳)#“tools.h”结束fi(菲涅耳)如果测试-f'usual.h'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'usual.h'”其他的echo shar:提取“'usual.h'”(4267个字符)sed“s/^X//”>'通常.h'<<'END_OF_FILE'X#ifndef __FILE_USUAL_H_SEEN__X#杂注一次X#定义__FILE_USUAL_H_SEEN__1X(X)X(X)X#未定义TURBO_CPP_COMPILERX#定义GNU_GPP_COMPILER 1X(X)X#定义USE_ATT_STYLE_IO_WITH_GNU 1X#取消USE_GNU_STYLE_IO_WITH_GNUX(X)X(X)X/*-------------------------参考---------------------------------X(X)XBjarne Stroustrup(1986),《C++编程语言》。阅读马萨诸塞州:XAddison Wesley出版公司。X(X)X头文件:X(X)X(X)第123、173页X(X)第77、237页X(X)第22、115页X(X)第126页X(X)第124页X(X)第12、226页X(X)第111页X(X)第35页X(X)X(X)X(X)XKernigham、Brian W.和Dennis M.Ritchie(1988年)。C编程语言,X第二版。新泽西州恩格尔伍德悬崖:普伦蒂斯·霍尔。X(X)X(X)X附录B标准库X(X)X B1输入和输出:第241页X B2字符等级测试:第248页X B3字符串函数:第249页X B4数学函数:第250页X B5实用程序功能:第251页X B6诊断第253页X B7变量参数列表:第254页X B8非本地跳跃第254页X B9信号:第255页X B10日期和时间功能第255页X B11实施定义的限制:第257页X(X)X(X)---------------------------------------------------------------------------*/X(X)X/*-------------------------g++-------------------------------------------X(X)XGNU的g++遵循Stroustrup。这些是g++等价物以及X它们的依赖性。X(X)X(X)X(X) X(X) X(X)X(X)X(X)X(X)X(X)X(X)X(X) X(X)X(X)X(X)X(X)X(X)X(X) X(X)X(X)X(X) X(X) X(X)X(X)X(X) X(X)X(X)&X(X)X(X)XTUSE_ATT_STYLE_IO_WITH_GNU选项使用符合以下条件的filebuf类斯特劳斯特鲁普第8章。打开文件的语法:X文件buf f1;X if(f1.open(“filename”,input)==0)错误(“无法打开文件”);X来自(&f1);X//选项:输入、输出、附加X(X)X USE_GNU_STYLE_IO_WITH_GNU选项使用File类。打开的语法X文件:X来自(“文件名”,io_readonly,a_useonly);X if(!from)错误(“无法打开文件”);X//选项:io_readonly、io_writeonly、io_readwrite、io_appendonly、io_appendX//选项:a_createonly、a_create、a_useonly、a_useX(X)----------------------------------------------------------------------------*/X(X)X#ifdef GNU_GPP_COMPILERX#包括X#结尾X(X)X(X)X/*----------------------------Turbo C++------------------------------------X(X)XTurbo C++使用了令人讨厌的ANSI C头文件和Stroustrup头文件的混合。有关说明,请参阅《Turbo C++程序员指南》第155页。X(X)XSyntax打开文件:X ifstream来自(“文件名”);X if(!from)错误(“无法打开文件”);流的X到(“文件名”)X if(!to)error(“无法打开文件”);X--------------------------------------------------------------------------------------------*/X(X)X#ifdef涡轮_CPP_COMPILERX#包含X#包含X#包含X#包含X#包含X#包含X#结束X(X)X#define LINESIZE 80//定义输出例程的行大小X(X)Xtypedef long int INTEGER;//定义整数算术的精度Xtypedef双REAL;//定义浮点运算的精度X(X)X#define REAL_ZERO 0.0//为REAL、0.0F、0.0或0.0L更正类型0X(X)Xtypedef void(*ONE_ARG_ERROR_HANDLER_T)(常量字符*);Xtypedef void(*TWO_ARG_ERROR_HANDLER_T)(常量字符*,常量字符*);X(X)X#结尾/*__FILE_USUAL_H_SEEN__*/文件结束如果测试4267-ne‘wc-c<‘通常.h’’;然后echo shar:“'usual.h'”打开时大小错误!fi(菲涅耳)#“usual.h”的结尾fi(菲涅耳)echo shar:存档1结束(共6个)。cp/dev/null ark1完成缺少=“”对于1 2 3 4 5 6中的I;如果测试-方舟${I} isdone公司 ; 然后MISSING=“${MISSING}${I}”fi(菲涅耳)完成如果测试“${MISSING}”=“”;然后echo您已打开所有6个档案。rm-f方格[1-9]完成其他的echo您仍然需要解压缩以下存档文件:echo“”${MISSING}fi(菲涅耳)##shell存档结束。退出0