标准击实试验的曲线拟合方法
曲线拟合的方法及过程
一、课程设计题目： 对于函数 xex x f --=)(从00=x 开始，取步长1.0=h 的20个数据点，求五次最小二乘拟合多项式5522105)()()()(x x a x x a x x a a x P -++-+-+= 其中 ∑===1995.020i ix x 二、原理分析 （1）最小二乘法的提法当数据量大且由实验提供时，不宜要求近似曲线）（x y φ=严格地经过所有数据点),(i i y x ，亦即不应要求拟合函数）（x ϕ在i x 处的偏差（又称残差）i i i y x -=)(φδ(i=1,2,…,m)都严格的等于零，但是，为了使近似曲线能尽量反应所给数据点的变化趋势，要求偏差i δ适当的小还是必要的，达到这一目标的途径很多，例如，可以通过使最大偏差i δmax 最小来实现，也可以通过使偏差绝对值之和∑ii δ最小来实现……，考虑到计算方便等因素，通常用使得偏差平方和∑ii 2δ最小（成为最小二乘原则）来实现。
按最小二乘原则选择近似函数的方法称为最小二乘法。
用最小二乘法求近似函数的问题可以归结为：对于给定数据),(i i y x(i=1,2,…,m)，要求在某个函数类Φ中寻求一个函数）（x *ϕ，使[][]21)(21*)()(min ∑∑=Φ∈=-=-mi iix mi iiy x y x ϕϕϕ（1-1） 其中）（x ϕ为函数类Φ中任意函数。
（1）确定函数类Φ，即确定）（x ϕ的形式。
这不是一个单纯的数学问题，还与其他领域的一些专业知识有关。
在数学上，通常的做法是将数据点),(i i y x 描绘在坐标纸上，然后根据这些点的分布情况来选择的）（x ϕ形式。
（2）球最小二乘法的解，即求满足条件（1-1）的近似函数）（x *ϕ。
（3）最小二乘法的实验原理 设）（x ϕ具有如下形式）（x ϕ=F),,,,10x a a a n ⋅⋅⋅（ （1-2） 其中n<m, k a (k=0,1,…,n)是待定参数，求具有这种形式的最小二乘解的实质，就是要适当选择k a =*k a (k=0,1,…,n)，使相应的函数 ),,,,()(**1*0*x a a a x n ⋅⋅⋅=ϕ （1-3） 满足条件（1-1），也就是说，点),,,(**1*0n a a a ⋅⋅⋅是多元函数 []211010),,,,),,,∑=-⋅⋅⋅=⋅⋅⋅mi i i n n y x a a a F a a a s （（的极小点，从而使**1*0,,,n a a a ⋅⋅⋅满足方程组0S=∂∂ka ，(k=0,1,…,n) （1-4） 因此，可以通过解上述方程组（称为法方程组）来求取*k a (k=0,1,…,n)，以便获得最小二乘解。
关于击实曲线的理解
关于土体击实曲线的理解清华大学李广信老师版的《土力学（第2版）》中，对土体的击实曲线有所叙述，实际上击实曲线是一般土力学教材中都会涉及到的问题，一般教材在叙述此部分时通常也是分为细粒土和粗粒土来讲的，今天笔者对这两种类型的击实曲线的理解做一点阐述，以供后面同学参考。
对细粒土：我们通过试验能够得到的击实曲线如下我们通过试验的曲线可以直观的看出一下几点：①在某一含水量、某一击实功的条件下，土体的干密度能达到最大值②随着击实功的增大，土体所能到达的最大干密度越来越大，相应的最优含水量会逐渐减小③击实曲线始终位于饱和曲线的下方对于上述的三点是我们直接从试验曲线中直观的得到的结论，下面我们对这三点结论给出一些阐述： 首先，我们先来理解饱和曲线，关于含水量=w s m m ω和干密度s d m Vρ=中，对任一土体，在s m 不变的情况下，水含量的多少会影响在不限制压实条件下最终的V 值，即我们假定（注意是假定！），在压实过程中，我们有手段或者措施对土体实现100%的压实，那么随着水的增多，压实后土体的剩余体积就在增大，最终如下图，也就是饱和曲线但是！但是！现实的情况并非如此，也就是说我们不可能有充分的手段对土体进行100%的击实，事实上在试验中我们反而要控制变量，即我们只能付出一定量的功来对土体进行压实，我们还是来定量分析下，下图的土颗粒，我们假设只用100J 进行压实（这里的100焦耳只是一个定性的概念，即指100个单位的能量）当无水加入时，我们假设土体在100J 的作用下体积压实了5个单位，此时压实功的效率为55%/100J =个,此时的干密度500.5395d ρ==。
当有10%的水加入时，体积压实了15个单位，此时功的效率为1515%/100J =个，相应的干密度500.58885d ρ== 再加入10%的水时，体积压实了20个单位，此时功的效率为2020%/100J =个，相应的干密度500.6380d ρ== 当再加入10%水时，由于水的继续增多，在土颗粒内形成封闭的气泡，会降低功的利用率，体积的压实会下降！此时体积由于封闭气泡的存在，体积只被压实了15个单位，功的利用率1515%/100J =个，相应的干密度500.58885d ρ==，以此类推，当水的含量继续增加时，压实功的利用率会进一步的降低，从而不能得到更大的干密度。
击实曲线与理论饱和线HDW20200915
击实曲线与理论饱和线在学习粘性土压实特性时，会遇到黏性土击实曲线和理论饱和线两条曲线。
学生在理解理论饱和线的物理含义，以及理论饱和线和击实曲线的关系时往往会遇到一定的困难。
因此，本文将对这一知识点进行扩展讲解。
图1即为常见的黏性土击实曲线和理论饱和线的关系图，黏性土的击实曲线反映的是：在一定击实功作用下，黏性土的干密度随含水率变化关系。
曲线中存在一个最优含水率，与之对应的是最大干密度。
黏性土的击实曲线是一条试验曲线，其数据点来自击实试验，通过数据的拟合可以得到一条光滑曲线。
粘性土的击实曲线反映出黏性土的击实特性。
在一定击实功作用下，当含水率较少时，黏土颗粒之间水膜的润滑作用不足，颗粒之间不易错动，土体不易被压实；当含水率较高时，颗粒间水膜润滑作用充足，但土体孔隙中水体不易被挤出而占据较多的孔隙空间，土体仍然不易被压实。
只有当土体含水率达到该击实功作用下的最优含水率时，黏土颗粒之间既有一定的润滑作用，水分也没有挤占太多的孔隙空间，土体可以达到该击实功作用下的最大干密度。
图1典型击实曲线与理论饱和线位置关系图当击实功变化时，黏性土的最优含水率也会相应变化。
因为较大的击实功可以抵抗较大的粒间错动阻力，使土颗粒相互错动形成更紧密的排列，从而达到压实的目的。
因此，当击实功增大，最优含水率会相应减小；击实功减小，最优含水率会相应增大。
并且，同一种土，击实功越大，所能达到的最大干密度也越大。
理论饱和线描述的是饱和含水率和干密度的函数关系，是一条理论曲线，而非试验曲线，与最大干密度没有必然联系。
理论饱和线的函数关系可以由三相草图进行推导。
理论饱和线中的干密度必大于击实曲线中相同含水率对应的干密度。
因为理论饱和线中对应的点表征的是土体处于饱和状态，而击实曲线中对应的点表征的是土体处于非饱和状态。
同一土体在相同含水率下，饱和状态肯定更密实，具有更大的干密度。
理论饱和线的含义：1，对应于纵坐标任一干密度做一水平线，与理论饱和线相交的点对应的横坐标含水率就是该土在该干密度下的饱和含水率。
曲线拟合的数值计算方法实验教材
曲线拟合的数值计算方法实验【摘要】实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。
曲线拟合（curve fitting）是指选择适当的曲线类型来拟合观测数据，并用拟合的曲线方程分析两变量间的关系。
曲线直线化是曲线拟合的重要手段之一。
对于某些非线性的资料可以通过简单的变量变换使之直线化，这样就可以按最小二乘法原理求出变换后变量的直线方程，在实际工作中常利用此直线方程绘制资料的标准工作曲线，同时根据需要可将此直线方程还原为曲线方程，实现对资料的曲线拟合。
常用的曲线拟合有最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束。
关键词曲线拟合、最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束一、实验目的1.掌握曲线拟合方式及其常用函数指数函数、幂函数、对数函数的拟合。
2.掌握最小二乘法、线性插值、三次样条插值、端点约束等。
3.掌握实现曲线拟合的编程技巧。
二、实验原理1.曲线拟合曲线拟合是平面上离散点组所表示的坐标之间的函数关系的一种数据处理方法。
用解析表达式逼近离散数据的一种方法。
在科学实验或社会活动中，通过实验或观测得到量x与y的一组数据对（X i，Y i)（i=1，2，...m），其中各X i 是彼此不同的。
人们希望用一类与数据的背景材料规律相适应的解析表达式，y=f(x，c）来反映量x与y之间的依赖关系，即在一定意义下“最佳”地逼近或拟合已知数据。
f(x，c)常称作拟合模型，式中c=(c1，c2，…c n)是一些待定参数。
当c在f中线性出现时，称为线性模型，否则称为非线性模型。
有许多衡量拟合优度的标准，最常用的一种做法是选择参数c使得拟合模型与实际观测值在各点的残差(或离差)，c）-f （f y e k k k =的加权平方和达到最小，此时所求曲线称作在加权最小二乘意义下对数据的拟合曲线。
击实实验报告标准击实试验报告
击实实验报告标准击实试验报告贵XX大学学生实验报告xxx五——xxx六学年度第一学期材料与建筑工程学院三年级专业土木工程班级课程名称土力学试验组别组员成绩 XX大学学生实验报告实验项目名称_ _击实试验_________ 组长 _______ 组别同组学生实验时间_20XX_年_X月_X日指导教师 _______ 实验报告需包含以下内容：①实验目的与原理；②主要实验仪器及设备；③实验步骤；④实验数据及处理；⑤实验结论；⑥其它；⑦思考题（一）实验目的与原理1.实验目的击实实验的目的是测定试样在一定击实次数下或某种压实功能下的干密度与含水率之间的关系，从而确定土的最大干密度和最优含水率。
2.试验原理土在一定的压实效应下，若含水率不同，则密度也会不同。
当压实功能和压实方法不变时，土的密度随含水率的增加而增加，但当含水率增大到一定程度后，土的密度反而减小，这是因为细粒土在含水率较低时，颗粒表面形成薄膜水，摩擦力大，不易压实；当含水率继续增加时，颗粒表面结合水膜渐渐加厚，其润滑作用也增大。
在外力作用下，容易移动，易于压实；而继续增加水量，只会增加土的孔隙体积，从而使干密度降低。
能使土体达到最大干密度的含水率称为最优含水率。
（二）主要实验仪器及设备1.击实仪：由击实筒、击锤及导筒组成。
本次实验选用轻型击实仪。
2.天平：称量200g，最小分度值0.01g；称量500g，最小分度值1.0g。
3.台秤：称量10kg，最小分度值5g。
4.圆孔筛：孔径为40mm、20mm和5mm标准筛各一个。
5.试样推出器：本实验用修土刀和刮土刀从击实筒中取出试样。
6.烘箱。
7.其他：喷水设备、碾土器、盛土盘、修土刀、平直尺、量筒和保湿设备等。
（三）实验步骤1.试样制备（湿法制备）① 称取天然含水率的代表性土体碾散，按要求过筛，将筛下的土样拌匀，并测定其含水率。
② 根据土样的塑限预估最佳含水率，按干法制备的原则选择至少5个含水率试样，分别将天然含水率土样风干或加水进行制备，并让制备好的土样喷洒的水分能均匀分布。
曲线拟合和插值运算原理和方法
实验10 曲线拟合和插值运算一． 实验目的学会MATLAB 软件中软件拟合与插值运算的方法。
二． 实验内容与要求在生产和科学实验中，自变量x 与因变量y=f(x)的关系式有时不能直接写出表达式，而只能得到函数在若干个点的函数值或导数值。
当要求知道观测点之外的函数值时，需要估计函数值在该点的值。
要根据观测点的值，构造一个比较简单的函数y=t (x),使函数在观测点的值等于已知的数值或导数值，寻找这样的函数t(x),办法是很多的。
根据测量数据的类型有如下两种处理观测数据的方法。
（1） 测量值是准确的，没有误差，一般用插值。
（2） 测量值与真实值有误差，一般用曲线拟合。
MATLAB 中提供了众多的数据处理命令，有插值命令，拟合命令。
1.曲线拟合已知离散点上的数据集[(1x ,1y ),………(n x ,n y )]，求得一解析函数y=f (x),使f(x)在原离散点i x 上尽可能接近给定i y 的值，之一过程叫曲线拟合。
最常用的的曲线拟合是最小二乘法曲线拟合，拟合结果可使误差的平方和最小，即使求使21|()|n i ii f x y =-∑ 最小的f(x).格式：p=polyfit(x,Y ,n).说明：求出已知数据x,Y 的n 阶拟合多项式f(x)的系数p ，x 必须是单调的。
[例 1.9]>>x=[0.5,1.0,1.5,2.0,2.5,3.0]; %给出数据点的x 值>>y=[1.75,2.45,3.81,4.80,7.00,8.60]; %给出数据点的y 值>>p=polyfit (x,y,2); %求出二阶拟合多项式f(x)的系数>>x1=0.5:0.05:3.0; %给出x 在0.5~3.0之间的离散值>>y1=polyval(p,1x ); %求出f(x)在1x 的值>>plot(x,y,‟*r ‟, 11,x y ‟-b ‟) %比较拟合曲线效果计算结果为：p=0.5614 0.8287 1.1560即用f(x)=0.56142x +0.8287x+1.1560拟合已知数据，拟合曲线效果如图所示。
巧用Excel软件处理击实试验数据并绘制曲线效果图
巧用Excel软件处理击实试验数据并绘制曲线效果图永利(新疆生产建设兵团农二师勘测新疆库尔勒市841000)【摘要】应用Excel对击实试验数据成果进行简化处理，绘制击实曲线图，并利用Excel强大的分析计算功能，求解出曲线的趋势方程及相关系数，进一步求解出最优含水率和最大干密度，简化了绘制曲线的烦琐工作，同时减少手工绘图的误差。
本文对此作了详细介绍。
【关键词】Excel 击实试验曲线饱和趋势方程相关系数最大干密度最优含水率饱和度1 前言土的击实试验是衡量土的压实性的一项很重要指标，在铁路、公路、桥梁、机场、水利（水库坝体）等工程建设的填方地基设计中都是必不可少的。
此外，土的渗透、压缩、剪切等试验在一定程度上又是由击实试验来控制的。
所以击实试验的成果对设计、施工以至于工程造价都会带来影响。
土在外力作用下（碾压或夯实），孔隙度变小，密度增加，强度提高，压缩性及渗透性降低，使土的工程性质得到改善。
而土的压实性就是指在一定的含水量条件下，以人工或机械的办法，使土能够达到某种密实程度的性能。
土的压实强度与含水量、压实功能和压实方法有着密切的关系。
当压实功能和压实方法不变时，则土的干密度随含水量的增加而增加；当干密度达到某一最大值后，含水量的继续增加反而使干密度减小。
此时干密度的最大值称为最大干密度，其相应的含水量称为最优含水量。
击实试验的目的，就是模拟工地压实条件，用标准击实的方法，测定在某种压实功能下土的含水量和干密度的关系，确定土的最优含水量和相应的最大干密度，以求用最小的压实功能，得到符合工程要求的密实度。
图1 击实试验计算和曲线效果图长期以来，工作在试验场所的技术人员面对大量击实试验数据，采用手工描点，曲线尺绘图的办法，不仅工作量大烦琐，并且极易受个人因素影响，得出的试验结果可信度比较差，不同程度地困绕着工程技术人员。
目前计算机已经普及，但专业绘制击实效果图的软件并不多见，且有的软件在输入数据时相对较为复杂。
标准击实试验的曲线拟合方法
[8 ]
4 结语
解决最大干密度和最佳含水量的方法时用三角函 数拟合是一种尝试 。本文是对于无机结合料的击实数 据进行的拟合 , 但对于不同材料的击实数据更适合于 哪种数据的拟合还有待进一步研究 。 参考文献 :
[ 1 ] 宫凤强 , 李夕兵 , 邓建 . 填土最大干密度的拟合正交多项式法推断 . 金
y ( x) =
6
n
cφ j j ( x)
( 1)
j =0
中寻找一曲线按照某种原则去拟合数据 ( 1 ) , 用所得的拟合曲线去代替数据 ( 1 ) 所反映的函数关 系 。若曲线
y ( x) =
m n
3
6
3
n
cj φj ( x )
m n
3
( 2)
j =0
使得 ( 3 ) 式成立
6 [6
i =0
j =0
(p - y, p
3
- y) = m ψ in ( p - y, p - y )
其中 p
A Ac
T
=
6
T
n
j =0
cj < j 。故满足条件 ( 3 ) 的拟合曲
3
3 线 y ( x ) 存在且唯一 , 并且从法方程 3
中解出 cj
3
( 4) =A y 3 3 3 T = ( c0 , c1 , Λ, cn ) , 就得到拟合曲线 ( 2 )
出的问题和风险 , 隧道施工应对可燃天然气的安全措 施绝不能减弱 。
5 建议与结语
[ 1 ] 合武铁路片 (麻 ) 岩隧道可燃气评价及处理措施 [M ]. 铁道第四勘察
可燃天然气隧道要安全施工 , 地质工作应该先 行 , 准确的地质数据是设计与施工措施的前提与基
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
拟合曲线 y ( x ) 对数据的拟合精度 , 可用均方差 σ 即二范数的平方来描述 σ =
3 实例
ρ w % = 41727 m ax = 21257 -4 σ = 2173 × 误差为 10 当 n = 2时 3 π π x × x × y ( x ) = c0 + c1 sin + c2 cos + 91842 91842 π π 2x × 2x × c3 sin + c4 cos 91842 91842 根据 ( 4 ) 可得 c0 = 3177537 c1 = - 21198849 c2 = - 0160240 c3 = 0132563 c4 = - 0168303 故所求方程为 3 y ( x ) = 2176118 - 0180536 sin ( 013192 x ) 0137426cos ( 013192 x ) + 0120053 sin ( 016382 x ) 0130054cos ( 016382 x ) 对方程求导 , 得 ρ m ax = 21257, w % = 41727 -4 σ = 1129 × 误差为 10 为了方便计算 , 取 a = 2 时已经完全满足工程的 需要了 。 而用幂函数作为插值基函数时 , 当 n = 2 时 σ = -4 -4 2129 × 10 ; 当 n = 3 时 σ = 1188 × 10 现将三角函数插值法和图解法列图如图 1。
p = [ y ( x0 ) , y ( x1 ) , Λ, y ( xm ) ]
T m +1
=
6
2
n
j =0
cj < j
由于
( p - y, p - y ) = =
3 3
6
m
m
[ y ( xi ) - yi ]
n
i =0
6 6
[
i =0 p∈
cφ j j ( xi ) - yi ]
2
j =0
所以 , 条件 ( 3 ) 可写成如下形式
冯 林等 : 标准击实试验的曲线拟合方法 = 0, 1, Λ, n ) 。由于基函数组 {φj ( x ) , j = 0, 1, Λ, n}
・73・
在点集 { x i , j = 0, 1, Λ, m } 上线性无关等价于向量组 T < j = (φj ( x0 ) , (φj ( x i ) , Λ, (φj ( xm ) ) ) , j = 0, 1, Λ, n 线性无关 , 故可用此向量组作基底构成 R 的一个子 空间 。记 ψ = Sp an { <0 , <1 , Λ, < n } 。 T A = [ <0 , <1 , Λ, < n ] y = ( y0 , y1 , Λ, ym )
[ 6 ] 王永良 . 曲线拟合法在土工数据处理中的应用 . 西部探矿工程 , 2004, 95 (4) : 38 - 39. [ 7 ] 颜庆 津 . 数 值 分 析 2000: 207 - 209. [ 8 ] 公路工程无机结合料稳定材料实验规程 (JTJ 057). 北京 : 人民交通出 (修 订 版 ). 北 京 : 北 京 航 空 航 天 大 学 出 版 社 ,
版社 , 1994.
收稿日期 : 2006 - 09 - 18
x = a + b - a
[8 ]
4 结语
解决最大干密度和最佳含水量的方法时用三角函 数拟合是一种尝试 。本文是对于无机结合料的击实数 据进行的拟合 , 但对于不同材料的击实数据更适合于 哪种数据的拟合还有待进一步研究 。 参考文献 :
[ 1 ] 宫凤强 , 李夕兵 , 邓建 . 填土最大干密度的拟合正交多项式法推断 . 金
・72・
全国中文核心期刊 路基工程 2007 年第 6 期 (总第 135 期 )
标准击实试验的曲线拟合方法
冯 林 刘 根
(长沙理工大学 湖南长沙 410076 ) (浙江核工业井巷建设公司 )
摘 要 对击实试验中矿料和含水量之间的函数关系 , 基于数值逼近原理 , 构造三角函数进行曲 线拟合 , 而后对其求导得到最大干密度和最佳含水量 。其误差用均方差描述 , 计算很符合矿料击实曲 线的分布规律 。结合实例为处理击实试验数据提供了一种方法 。 关键词 曲线拟合 三角函数 最小二乘法 最大干密度 最佳含水量
出的问题和风险 , 隧道施工应对可燃天然气的安全措 施绝不能减弱 。
5 建议与结语
[ 1 ] 合武铁路片 (麻 ) 岩隧道可燃气评价及处理措施 [M ]. 铁道第四勘察
可燃天然气隧道要安全施工 , 地质工作应该先 行 , 准确的地质数据是设计与施工措施的前提与基
设计院 , 2006.
收稿日期 : 2006 - 08 - 21
π
t 故 t =
x - a π b - a
现分别取 a = 1, 2 时的情况对表 1 拟合 。
ρ 表 1 d - ω标准击实试验数据表
含水量 w / %
21291 21690 31413 41921 61158 61497
属矿山 , 2005, 352 (10) : 51 - 53.
[ 2 ] 卿启湘 , 王永和 , 李光耀 . 标准击实试验的一种数据处理方法 . 湖南
1 引言
目前国内击实试验数据处理常用图解法 , 即先通 过室内标准试验 , 得到关于干密度和含水量的系列数 据 , 然后绘制曲线 , 曲线的峰值点即为所求的最大干 密度和最佳含水量 。但由于人为因素的影响或试验设 备精度不足 , 其数据离散性比较大 , 过所有点要绘制 一条光滑曲线几乎不可能 , 图解法所得到的最大干密 度也要小于或等于该样品的实际最大干密度 。为了得 到更精 确 的 数 据 , 国 内 同 行 大 多 用 幂 函 数 拟 合 曲 [ 1 ]～[ 7 ] 线 , 作者通过用三角函数多项式拟合曲线 , 发 现它很符合击实曲线的分布规律 , 其中室内试验数据 [8 ] 按有关规程 所得 。
2 最小二乘法
标 ) : ( x i , y i ) , 其中 a = x0 < x1 <Λ xm = b 。又 选定 n + 1 个在区间 [ a, b ] 上连续且在点 { x i , i = 0, 1, Λ, m } 上线性无关的基函数 φj ( x ) ( j = 0, 1, ∧, n ) , 其中 n ≤m 。问题是要在曲线族
(p - y, p
3
- y) = m ψ in ( p - y, p - y )
其中 p
A Ac
T
=
6
T
n
j =0
cj < j 。故满足条件 ( 3 ) 的拟合曲
3
3 线 y ( x ) 存在且唯一 , 并且从法方程 3
中解出 cj
3
( 4) =A y 3 3 3 T = ( c0 , c1 , Λ, cn ) , 就得到拟合曲线 ( 2 )
91842
根据 ( 4 ) 可得 c0 = 1187107 c1 = 0138771 c2 = - 01021901 故所求方程为 3 y ( x ) = 1187107 + 0138771 sin ( 013192 x ) 01021901cos ( 013192 x ) 对方程求导 , 得
学报 , 2002, 24 (4) : 538 - 540.
y ( x) =
( 1)
j =0
中寻找一曲线按照某种原则去拟合数据 ( 1 ) , 用所得的拟合曲线去代替数据 ( 1 ) 所反映的函数关 系 。若曲线
y ( x) =
m n
3
6
3
n
cj φj ( x )
m n
3
( 2)
j =0
使得 ( 3 ) 式成立
6 [6
i =0
j =0
长安大学学报 , 2002, 22 (2) : 10 - 13.
[ 4 ] 吕鹏 , 邓海 , 张昀青 . 填料最大干密度的曲线拟合解 . 岩土工程技术 , 2003, (3) : 230 - 232. [ 5 ] 黄英 , 符必昌 . 确定土的最大干密度和最优含水率的数解法 . 岩土工程
91842
6
m
[ y ( xi ) - yi ]
3
2
i =0
设一组击实试验数据如表 1 , 取 ρ d - ω函数关 系曲线的三角函数拟合 。 取三角函数组 ψ = { 1, sin x, co sx, sin2 x, cos2 x, Λ, sin ax, cosax} , a = 1, 2, Λ, n 作为插值基函数 。取含水 量为自变量 , 取干密度为因变量 。大量实验表明 , ρ d - ω曲线形状成抛物线 , 而无机结合料的试验所取含 水量为 [ a, b ] , a = 0, b为所测最大干密度对应含水 量的 2 倍 , 故 b = 91842, 其中选取函数区间时需将其 变换成三 角 函 数 的 半 个 正 弦 区 间 即 [ 0, π ] , 变 换 公式
cj φj ( x i ) - y i ]
2
= m in
{ cj }
6 [6
i =0
j =0
cφ j j ( xi ) -
设在 xoy直角坐标系中给定 m + 1 对数据 (即坐
yi ]
2
( 3)
3
冯 林 , 男 , 硕士研究生 。
则称 y ( x ) 为在曲线族 ( 1 ) 中按最小二乘原则 确定的对于数据 ( 1 ) 的拟合曲线 。所谓最小二乘法 3 3 求拟合曲线 y ( x ) , 就是按条件 ( 3 ) 求出系数 cj ( j 础 , 是建设方拨付措施费用数量和监理方采取相关监 控手段和处置方案的依据 。隧道地质工作可以从以下 几个方面开展 。 ( 1 ) 针对该区地质条件特点有效地开展超前地 质预测预报工作 。技术手段上主要立足于在隧道施工 过程中采用 TSP、超前水平钻 、隧道开挖风钻钻进情 况记录和掌子面地质编录 , 结合使用超前 15 m 的声 波探测 、地质雷达 、 HSP 等进行 。 ( 2 ) 补充对已出现燃烧事件的地段进行岩石学 与岩石组构学研究 , 重点是 : ① 有效孔隙度 ; ② 渗透 性; ③ 孔隙结构和微构造裂隙的定性定量分析 。 ( 3 ) 补充对埋深较大的隧道进行区域性的地表 破裂系统的调查与分析 , 以便找出裂隙发育规律和为 评价近地表可能的圈闭规模提供数据 。 ( 4 ) 天 然 气 (瓦 斯 ) 隧 道 施 工 , 应 严 格 执 行 《 铁路瓦斯隧道技术规范 》安全标准 , 加强瓦斯监测 与预测预报 , 加强隧道通风管理 , 制定并认真落实瓦 斯隧道安全施工预案 , 确保施工安全 参考文献 :