基于多元回归的炉龄问题的实现 2.0

1、题目：基于多元回归的炉龄问题的实现成员信息：姓名班级学号张队长队员&答辩人日期：2019年12月30日摘要工业生产在生产出产品的同时，也生产出大量关于生产过程的信息。应充分利用这些信息并在此基础上寻求如何才能使生产得到进一步改善的方法。因为钢铁生产过程特别复杂，出现的因素特别多，各因素间又相互牵连1。因此找出某一时期影响某一实际问题的最主要原因，提供优化决策，不仅对生产技术管理有很大的参考价值，而且也有助于积累经验，还对机理性问题的研究具有启发和促进作用。现有某钢铁公司炼钢转炉的炉龄按30天炉/天炼钢规模，大约一个月就需等炉一次进行检修。为了减少消耗，厂家通过实际测定，得到33组数据。要求对炉

2、龄进行预测和分析。本文针对此问题，进行分析并建立合适的数学模型，为解决炉龄问题提供方案。针对问题一，分析各个因素横向之间以及各个因素纵向内部数据之间的关系采用灰度预测模型，将缺失的数据进行补全。针对问题二，通过主成分分析法找出与影响炉龄的主要因素，剔除对炉龄影响较小的因素。建立以多个影响因素为自变量和以炉龄为因变量的多元线性回归模型，并且建立了支持向量机的回归模型。针对问题三，对于多元线性回归模型，主要考虑其预测的残差来分析模型的可靠性，模型预测残差的分布越没有趋势，表示回归的结果越是可靠，所以多元线性回归模型对本问题有一定的适用性。对于支持向量机回归模型，将其与BP神经网络进行对比，发现其性

3、能明显优于BP神经网络。针对问题四，分析找出对转炉寿命影响较大的因素为喷补料量、炼钢时间和渣中含铁量，并将设计延长炉龄方案时重点放在这三个因素。关键词：转炉炉龄灰度预测多元线性型回归支持向量机目录TOC o 1-5 h z HYPERLINK l bookmark8 一、问题重述4 HYPERLINK l bookmark10 问题背景4 HYPERLINK l bookmark12 问题的提出4二、问题的分析4 HYPERLINK l bookmark14 三、基本假设5四、符号说明5五、问题一的模型建立与求解5问题的分析5 HYPERLINK l bookmark50 模型的建立与求解6六

4、、问题二的模型建立与求解9 HYPERLINK l bookmark110 问题的分析9模型的建立与求解9七、问题三的求解19 HYPERLINK l bookmark114 问题的分析19 HYPERLINK l bookmark116 问题的求解19八、问题四的求解21问题的分析21问题的求解21 HYPERLINK l bookmark118 九、参考文献22、问题重述问题背景自1952和1953年氧气顶吹转炉炼钢在奥地利的钢铁股份公司林茨钢厂与砂冶公司多纳维茨钢厂先后建成投产。转炉炼钢以其生产率高、品种多、质量好、热效率高、原材料适应性好、消耗低、成本低、基建投资少、建设速度快等优点，

5、被国内外钢铁冶炼行业广泛采用，成为现代炼钢的主要方法之一。转炉炉龄是炼钢车间的一项综合性技术经济指标。炉龄的离低直接影响转炉钢产及耐火材料消耗和成本等指标。因此，炉龄的提高对于技经指标的改善和炼钢成本的降低具有十分重要的意义2。问题的提出某钢铁公司炼钢转炉的炉龄按30天炉/天炼钢规模，大约一个月就需等炉一次进行检修。为了减少消耗，厂家通过实际测定，得到下表所示的数据，其中q为喷补料量、X2为吹炉时间、%3为炼钢时间、X4为钢水中含锰量、X5为渣中含铁量、%6为作业率、目标变量y为炉龄（炼钢炉次/炉）。要求完成如下四问：问题1：由于某种原因，造成个别数据缺失，试对这些缺失数据（表中用表示）进行补

6、全。问题2：试根据附表数据建立炉龄的预测模型。问题3：采用适当的指标和方法对第2问建立的模型，和第1问修正数据之后的改进模型进行可靠性分析，说明模型对实际问题的适用性。问题4：为钢铁公司提出延长炉龄的方案。二、问题的分析此问题属于多元回归问题。回归问题是建立因变量Y与自变量X之间关系的模型。我们需要利用数据统计原理，对大量统计数据进行数学处理，并确定因变量与某些自变量的相关关系，建立一个相关性较好的回归方程（函数表达式），并加以外推，用于预测今后的因变量的变化。根据因变量和自变量的个数分为：一元回归分析和多元回归分析。此问题属于多元回归。三、基本假设1.所得数据可靠准确，不包含人为误差。2.炉

7、龄的影响因素只有喷补料量、吹炉时间、炼钢时间、钢水中含锰量、渣中含铁量、作业率六个因素。四、符号说明符号说明单位%1喷补料量无2吹炉时间3炼钢时间%4钢水中含锰量无5渣中含铁量作业率y炉龄炼钢炉次/炉A(k)序列的级比b.j信息贡献率Q残差平方和K(x.,x.)=(x冲*)核函数表1：符号说明五、问题一的模型建立与求解5.1问题的分析问题一首先要求我们根据一项已知的数据，对个别数据缺失，试这些缺失数据（表中用表示）进行补全。由附表可以看出，给出的数据是一些离散的，规律性不强的数据，此时我们采用灰色预测的方法对其进行补全。灰色系统理论认为，尽管客观表象复杂，但总是有整体功能的，因此必然蕴含某种内

8、在规律。关键在于如何选择适当的方式去挖掘和利用它。灰色系统是通过对原始数据的整理来寻求其变化规律的，这是一种就数据寻求数据的现实规律的途径，即为灰色序列的生成。一切灰色序列都能通过某种生成弱化其随机性，显现其规律性。灰色预测的主要特点是模型使用的不是原始数据序列，而是生成的数据序列。其核心体系是灰色模型(GreyModel，简称GM)，即对原始数据作累加生成(或其它方法生成)得到近似的指数规律再进行建模的方法。5.2模型的建立与求解数据的检验与处理首先，为了保证建模方法的可行性，需要对已知数据列作必要的检验处理。设参考数据为x(o)=(x(o)(1),x(o)(2),，尤仇)，计算序列的级比：

9、x(o)(k1)x(o)(k)k=2,3,5.1)如果所有的级比A(k)都落在可容覆盖0=(e-訂,e才)内，贝U序列x(o)可以作为模型GM(1,1)的数据进行灰色预测。如果所有的级比A(k)都落在可容覆盖0=(e-訂,e才)内，则序列兀(。)可以作为模型GM(1,1)的数据进行灰色预测。否则，需要对序列x(o)做必要的变换处理，使其落入可容覆盖内。即取适当的常数c,作平移变换:y(o)(k)=x(o)(k)c，k=1,2,小5.2)使序列y(o)=(y(o)(1),y(o)(2),y(o)(n)的级比:=2,3,n5.3)5.2.2建立模型1.已知参考数据列5.4)x(o)=(x(o)(1

10、),x(o)(2),x(o)(n)进行1次累加生成序列：%(1)=(%(1)(1),%(1)(2),%(1)(n)=(x(o)(1),x(o)(1)X(o)(2),x(o)(1)X(o)(n)5.5)其中x(1)(k)=Sf=1(0)(i)Ck=1,2,-,n)o尤的均值生成序列：z(i)=(z(i)(2),z(i)(3),，z(i)(n)(5.6)其中z(i)(k)=0.5x(i)(k)+0.5x(i)(k1),k=2,3,兀。2.建立灰微分方程x(o)(k)+az(i)(k)=bk=2,3,皿(5.7)相应的白化微分方程为：dx(i)吐+ax(i)(t)=b(5.8)TOC o 1-5 h

11、 zz(1)(2)1记=a,bT，Y=x(0)(2)x(0)(3),-,x(0)(n)T，B=-z(3)1，则由最小-z(i)(n)1二乘法，求得使J(u)=(Y-Bu)r(YBu)达到最小值的的估计值u=a,bT=(BtB)1BtY于是求解方程(5-8)得：(1)(k+1)=(x(0)(1)e(ik+JLk=0,1,n1,(5.9)CLCL5.2.3模型求解1.对于无4(钢水中含锰量)，先选取前19个数据对第20个数据进行的级比图1:残差A(k)，发现不满足级比要求。多次测试后，选择：x(o)=(50，44，46，46，48，45，42，40)此时级比满足要求。由此补全所缺失的值：502.对

12、于无6(作业率)，先选取前19个数据对第27个数据进行的级比A(k)，发现不满足级比要求。多次测试后，选择：%(0)=(35，33.3，37.9，42.9)0.21.21415US2222.4Z62.83图2：残差此时级比、残差满足要求。由此补全所缺失的值：37.777.3对于y(炉龄)，先选取前25各数据对第27个数据进行的级比A(k)，发现不满足级比要求。多次测试后，选择%(0)=(832，1076，1376，914)此时级比满足要求。由此补全所缺失的值：1366.5。补全后的数据如下：序号x1x2x3x4x5x6y200.215520.240.25016.841.71098280.243

13、617.737.24516.237.771105300.19717.335.94613.857.41366.5表2：应补全的数据六、问题二的模型建立与求解6.1问题的分析问题二要求根据附表数据建立炉龄的预测模型。由表中数据可以看出，我们需要找出多个变量间是否相关、相关方向与强度，并建立数学模型以便观察特定变量。对此我们采用回归分析的方法。回归分析（RegressionAnalysis是一种统计学上分析数据的方法，回归分析是建立因变量Y（或称依变量，反因变量）与自变量X（或称独变量，解释变量）之间关系的模型。在建立模型前，由于自变量过多，可先对自变量进行主成分分析。然后可先建立多元线性回归模型来

14、观察。模型的建立与求解6.2.1数据的处理主成分分析的主要目的是希望用较少的变量去解释原来资料中的大部分变异，将我们手中许多相关性很高的变量转化成彼此相互独立或不相关的变量。通常是选出比原始变量个数少，能解释大部分资料中的变异的几个新变量，即所谓主成分，并用以解释资料的综合性指标。由此可见，主成分分析实际上是一种降维方法。主成分分析法的步骤如下对原始数据进行标准化处理假设进行主成分分析的指标变量有m个，分别为xlfx2，,xm，共有个评价对象，第i个评价对象的第丿个指标的取值为。将各指标值%转换成标准化指标值孔，I)a仏/、&.=L(i=1,2,，n；j=1,2,，m)iJs.j(6.1)其中

15、出=1Lf_ia，s.=丄&(%-比)2，(丿=1,2,m),即岀，Sj为第丿九t丄丿丿n1i丄U丿丿丿丿个指标的样本均值和样本标准差。对应地有，X.U.x.L(j1,2，m)丿s.丿(6.2)计算相关系数矩阵R相关系数矩阵R_(r.)mxm九a.r.k1且乩，Ci,j1,2,m)，ijn1(6.3)式中q1，冷r.i，r.是第i个指标与第丿个指标的相关系数。计算特征值和特征向量计算相关系数矩阵R的特征值0，及对应的特征向量u1fu2,，其中u.u1j,u2j,umjT，由特征向量组成m个新的指标变量y1un1+u212+um1my2u121+u222+um2myux+ux+uxm1m12m2

16、mmm式中歹是第1主成分，y2是第2主成分，ym是第m主成分。选择p(pSm)个主成分，计算综合评价值i)计算特征值A.(j=1,2,m)的信息贡献率和累积贡献率。称A=丄()=1,2,，m)丿电人6.4)为主成分丹的信息贡献率,k=1k6.5)为主成分片y2,，儿的累积贡献率，当接近于1(一般取丄2ppP0.85,0.90,0.95)时，则选择前p个指标变量yi，y2，-,yp作为p个主成分，代替原来m个指标变量，从而可对p个主成分进行综合分析。ii)计算综合得分pZ=Yb.y.j=i6.6)其中吗为第丿个主成分的信息贡献率，根据综合得分值就可进行评价。主成分分析法的结果分析如下：各于I:戒

17、分的黃献率菩T主成分图4：各成分贡献率由结果可知，坷（喷补料量）、X2（吹炉时间）、%3（炼钢时间）、X4（钢水中含锰量）、x5（渣中含铁量）和耳（作业率）的贡献率分别为30.3143,25.7442,5619.7510，14.1698，7.6058，2.4148。因此前五个主成分的贡献率已达97.5852%，实际应用中可只取前五个成分。6.2.2建立模型1.多元线性回归预测模型（1）建立的模型为卩=00+耳卩+氐5+&、W（O,02）,（6.1）式中00心,鶴心都是与X1,X2,Xm无关的未知参数，其中W,B血称为回归系数。现得到个独立观测数据切Qj，Q庆,其中q为y的观察值，%分别为的观察

18、值，i=1,皿，nm，由（7.1）得b.=Po+Pa.i+-+Pma.m+se耳N（0,o2）,i=1,n.6.2)记：111X=1an1b1Y=:bn=1，疋BB0“1”血Tairnanm6.3)（6.1）即可表示为a=xp+&、WN（0,O2E）（6.4）其中琶为九阶单位矩阵。2）参数估计模型（6.1）中的参数仔0，仔肌用最小二乘法估计，即应选取估计值Bj使当仔丿.=片，）=0,1,m时，误差平方和i=1i=1i=1（6.5）达到最小。为此，令：西=0，j=0,1,2,，九，曙=一2(-0。一叭鶴g=0dQ一仔0_仔1唱1maim)aij=0j=1,2，m.6.6)经整理化为以下正规方程组

19、：恥+0茂勺1+02$勺2+血乞咕工TOC o 1-5 h zHHHH和n和n和.0。aH+久a2+02anai2+血anatmZaHbi，i.1i.1i.1i.1i.1、n、n、n、n、力B/a+B/aa+B/aa+B/a2ab,H。im1imH2imi2mimimi,i.1i.1i.1i.1i.1（6.7）正规方程组的矩阵形式为XtXBXtY（6.8）当矩阵X列满秩时，XTX为可逆方阵，式（7.8）的解为B（XtX）1XtY.（6.9）将“代回原模型得到y的估计值TOC o 1-5 h zy仁+Rx+Rx011m（6.10）而这组数据的拟合值为b.R+pa+Ba.（I1,，Q,i01i1i

20、mbr-,bvT,拟合误差YP称为残差，可作为随机误差的估计，HH6.11）称Q为残差平方和。（3）模型求解通过计算（MATLAB程序参考了文献7）,得到：y=3527.1+1005.8*xl-24.335*x2-58.476x3-0.081657x4-6.7537x5+0.5902x6。得到真实值和拟合值的结果如图：图5：真实值和拟合值结果残差的置信区间如下图所示：复相关系数R为0.6相关关系不成立。同时，系数95%置信区间bint包含零点，可认为在95%置信度下，相应变量对y影响不显著，寻找异常点，考虑改进模型。600400200总0-200-600ResidualCaseOrderPlo

21、t51015202530CaseNumber图7：去除异常点发现3个异常点：10、12、29。从原始数据中将其剔除后，重新做多元线性图8：去除异常点，真实值和拟合值结果回归，拟合结果如图：此时计算得：R=0.803，相关关系成立。同时F值有了较大的提升。此时，得到y=仔。+臥+B6x6+&中B的值为:2387.0891，1508.4431,11.9415,-42.6046,1.9757,-25.5838，0.4669。同时B6为0.4，验证了主成分分析的结果。除此之外，我们还建立了支持向量机的回归预测模型，下面介绍如下。2.支持向量机的回归预测模型3,4,8,9支持向量机(sup-portve

22、ctormachines,SVM)是-一种基于结构风险最小化原理的机器学习技术，能较好地解决小样本、非线性、高维数和局部极小点等问题5。模型建立设含有i个数据的集合为(叫,号川=1,2.2，其中叫eRd是第i个训练样本的输入列向量，叫二尤严2，厂此eR为对应的输出值。设在高维特征空间中建立的线性回归函数为/(%)=30(叫)+b。其中，0(%)为非线性映射函数,e表示权重，b表示偏差。引入松弛因子g.,?，原问题转化为：I*引入松弛因子召，原问题转化为：m加訓3112+CZ(+*)y.-w0(x.)-b+勺，30(叫)y.+g.,i=1,2I(6.12)式中：表示不敏感损失系数；C为惩罚因子，

23、C越大表示对训练误差大于的样本惩罚越大，规定了回归函数的误差要求，越小表示回归函数的误差越小。J=1i=1上述模型对应的对偶形式：Q*)q-0*)x(%,亏)一zs+。*)+ZS+即儿J=1即=，j=1i=1a.C0a.a.C(6.13)其中：K(x.,x.)=e(x)e(x为核函数，本文选取应用广泛、精度较高的i丿ijRBF核函数，其表达式为一比一%.|2K(Xj,Xj)=exp(/)(6.14)通过(3)得最优解a=aia2勺,a-=aia2叩，并根据此求出权重e*和b*，最后得出支持向量机得回归函数f(x)=w*(x)+b*=a*)(x.)(x)+b*=/(a*)K(x.,x)+b*1=

24、11=1(6.15)模型求解先将附表中的数据去除缺失行，利用余下的数据建立回归模型。通过多次计算及参数选择(MATLAB程序参考了文献5-6,得出计算过程中的的最优训练集和测试集样本的预测仿真效果见下图。训练礙预測斜里对比m昭二0.00010183R2二0999811600I.III3t?400111051015202530样本跚号40-1O20OCJ18DO60图9：训练集预测结果对比从两图可以看出，训练集和测试集的均方误差分别为0.00010183和0.00063372,决定系数分别达到0.99981和0.92449。这表明利用六个自变量和一个因变量所建立的SVR回归模型具有非常好的泛化能

25、力。此时选择的训练集为27个,测试集为3个。七、问题三的求解7.1问题的分析问题三有以下两个要求：a.采用适当的指标和方法对第2问建立的模型进行评价和分析；b.对利用第1问修正数据之后的改进模型进行可靠性分析，说明模型对实际问题的适用性。7.2问题的求解7.2.1多元线性回归预测模型根据问题1，补全后的数据见附表，利用补全后的数据进行回归分析，步骤同6.2.2。此时，得到y=仔。+臥+B6x6+&中B的值为：2458.792,1647.034,40.609，-56.875，2.282，-31.774，1.086。Q*6005002530700140D13001100二.1OOD*900BOO1

26、01520样本躺号缶.-0拟命價hrnr图11：真实值和拟合值结果回归模型的可靠性分析主要考虑其预测的残差，由其原理可知，模型预测残差的分布越没有趋势，就表示回归的结果越是可靠。图12：残差分布图由残差的分布图可知，每个点代表一个残差，在途中可以看到残差点的分布是非常分散的。所以该模型的可靠性具有一定的显著性，对本问题有一定的适用性。支持向量机的回归预测模型为了对比SVR回归模型的性能，这里将之与BP神经网络对比。与图1和图2对应的训练集和测试集相同，建立的BP神经网络6对测试集的预测结果如图所不难发现，SVR回归模型的性能要明显优于BP神经网络。八、问题四的求解问题的分析问题四要求为钢铁公司

27、提出延长炉龄的方案。对于此问题我们需要分析六个自变量中哪个对因变量的影响更大。问题的求解根据问题作以下考虑：自变量变动一个单位，因变量变动多少个单位。单从此方面来说，如果系数很小，比如0.00001这种，我们说虽然具有统计上的显著性，但是没有经济的显著性。由本文“六、问题二的模型建立与求解”知，xl，x2，x3,x4，x5,x6的系数分别为1508.4431,11.9415,-42.6046,1.9757,-25.5838,0.4669。可知x1,x3,x5对炉龄影响较大。据此建议钢铁公司设计延长炉龄方案时应该把重点放在关注喷补料量、炼钢时间和渣中含铁量三方面提出改进措施。九、参考文献刘苏杭提

28、高卧式转炉炉龄的方法和措施探究J.湖北农机化,2018(05):52.文灏，区勇铭转炉炉龄影响因素的探讨J.南方钢铁，1996(04):13-15.薛磊，吴炜一种基于支持向量回归的多层原油分层产能贡献预测模型J.石油化工高等学校学报,2006,19(4):88-92.徐飞，徐卫亚.基于支持向量机-马尔可夫链的位移时序预测J.岩土力学,2010，31(3):944-948.蔡煜东遗传程序设计用于工业调优一一炼钢转炉炉龄预测J.自动化学报，1997，23(1):51-54 HYPERLINK /lyxleft/article/details/82880860 /lyxleft/article/de

29、tails/82880860司守奎，孙玺菁.数学建模算法与应用M.2011.陈剑勇，苏浩益.结合支持向量机和马尔可夫链算法的中长期电力负荷预测模型J.南方电网技术，2012，6(1):54-58.李丹玲，陈平雁，周凤麒.基于线性v-支持向量回归机的异常数据检测J.数理统计与管理，2011，30(01):59-63.附录：1、附表：序号x1x2x3x4x5x6y10.292218.541.4581883.3103020.267218.441511891.7100630.268517.738.65217.378.9100040.183518.941.81812.847.270250.2348183

30、9.45117.457.4108760.138618.940.53912.822.590070.208318.339.86417.152.670880.41818.8416416.426.7122390.10318.439.22012.335803100.489319.341.44919.131.3715110.205819404018.841.2784120.092517.938.75014.366.7535130.18541940.8442128.6949140.196318.137.24615.3631012150.100818.2374616.833.9716160.270218.93

31、9.54820.231.3858170.146519.138.64517.828.1826180.13531938.64216.739.71015190.224418.837.74017.449861200.215520.240.24616.841.71098210.031620.941.24817.452.6580220.049120.340.65619.735573230.148719.439.54218.333.3832240.244518.236.64115.237.91076250.222218.4374013.742.91376260.129818.437.24517.244.39

32、14270.2318.437.14722.921.6861280.243617.737.24516.237.7771105290.280418.337.54617.320.31013300.19717.335.94613.857.41366.5310.18416.235.34316.644.81039320.167917.134.64320.337.31502330.152417.6365114.236.711282、源码1）灰度预测clc,clearx0=load(mydatax.txt);%注意这里为列向量n=length(x0);lamda=x0(1:n-1)./x0(2:n)%计算级比

33、range二minmax(lamda)%计算级比的范围x1=cumsum(x0);%累加运算B=-0.5*(x1(1:n-1)+x1(2:n),ones(nT,1);&Y=x0(2:n);u=BY%拟合参数u(1)=a,u(2)=bx=dsolve(Dx+a*x=b,x(0)=x0);%求微分方程的符号解x=subs(x,a,b,x0,u(1),u(2),x0(1);%代入估计参数值和初始值yuce1=subs(x,t,0:nT);%求已知数据的预测值y=vpa(x,6)%其中的6表示显示6位数字yuce=x0(1),diff(yuce1)%差分运算，还原数据epsilon=x0-yuce%计

34、算残差delta二abs(epsilon./xO)%计算相对误差rho=l-(l-0.5*u(l)/(l+0.5*u(l)*lamda%计算级比偏差值,u(1)=a主成分分析clc,cleargj=load(mydatax.txt);%把原始数据保存在纯文本文件gj=zscore(gj);%数据标准化r=corrcoef(gj);%计算相关系数矩阵%下面利用相关系数矩阵进行主成分分析，x的列为r的特征向量，即主成分的系数x,y,z二pcacov(r)%y%r的特征值，z为各个主成分的贡献率f=repmat(sign(sum(x),size(x,1),1);%构造与x同维数的元素为1的矩阵&x=

35、x.*f%修改特征向量的正负号，每个特征向量乘以所有分量和的符号函数值num=6;%num为选取的主成分的个数df=gj*x(:,1:num);%计算各个主成分的得分tf=df*z(1:num)/100;%计算综合得分stf,ind=sort(tf,descend);%把得分按照从高到低的次序排列stf二stf,ind二indbar(z)title(各个主成分的贡献率)xlabel(各个主成分)ylabel(贡献率);多元线性回归clc,clearx=load(reallx.txt);%提取x的观察值y=load(really.txt);%提取因变量y的观察值Y=nonzeros(y);%去掉

36、y后面的0,并变成列向量n=length(Y);%样本点的个数N=1:n;%标记样本点X=ones(n,1),x;%构造多元线性回归分析的数据矩阵X&beta,betaint,r,rint,st=regress(Y,X)；%计算回归系数和统计量等，st的第2个分量就是F统计量，下面根据统计量的表达式重新计算的结果和这里是一样的。q=sum(r.2);%计算残差平方和ybar二mean(Y);%计算y的观察值的平均值yhat=X*beta;%计算y的估计值u=sum(yhat-ybar).2);%计算回归平方和MM=size(x);m=MM(2);%变量的个数，拟合参数的个数为m+1F=u/m/

37、(q/(n-m-l);%计算F统计量的值，自由度为样本点的个数减拟合参数的个数fwl二finv(0.025,m,n-m-l);%计算上alpha/2分位数fw2二finv(0.975,m,n-m-l);%计算上l-alpha/2分位数c=diag(inv(X*X);%计算c(j,j)的值t二beta./sqrt(c)/sqrt(q/(n-m-l);%计算t统计量的值tfw二tinv(0.975,n-m-1);%计算t分布的上alpha/2分位数savexydataYx;%把Y和x123保存到mat文件xydata中，供问题(3)的二次模型使用PY=beta(1)+beta(2)*x(:,1)+

38、beta(3)*x(:,2)+beta(4)*x(:,3)+beta(5)*x(:,+beta(6)*x(:,5)+beta(7)*x(:,6);%计算拟合值%二次模型%loadxydata%rstool(x,Y)%模型相关数据R=sqrt(u/(u+q);%复相关系数figure(1)plot(N,rint(:,1),N,rint(:,2);xlabelC残差的置信区间)figure(2)%rcoplot(r,rint);plot(r,o)legend(残差值)xlabel(样本编号)%拟合值真实值对比figure(3)plot(N,y,r-*,N,PY,b:o)legend(真实值,拟合值

39、)xlabel(样本编号)ylabel(炉龄)4）支持向量机的回归%I.清空环境变量clcclearall%II.导入数据%loadconcrete_data.matattributes=load(attributes2.txt);strength=load(strength2.txt);&%1.随机产生训练集和测试集n=randperm(size(attributes,2);%2.训练集xl=input(输入训练集的个数(小于30);p_train=attributes(:,n(1:xl);t_train=strength(:,n(l:xl);%3.测试集cs=xl+1;p_test二att

40、ributes(:,n(cs:end);t_test二strength(:,n(cs:end);%III.数据归一化%1.训练集pn_train,inputps=mapminmax(p_train);%归一化(mapminmax)pn_train=pn_train;pn_test二mapminmax(apply,p_test,inputps);pn_test二pn_test;%2.测试集tn_train,outputps=mapminmax(t_train);tn_train=tn_train;tn_test二mapminmax(apply,t_test,outputps);tn_test二tn_test;%IV.SVM模型创建/训练%1.寻找最佳c参数/g参数%c,g=meshgrid(0:3:120,500:3:620);(T0:0.5:10,T0:0.5:10)c,