分形维数算法
分形维数算法
分形包括规则分形和无规则分形两种。规则分形是指可以由简单的迭代或者是按一定规律所生成的分形，如Cantor集，Koch曲线，Sierpinski海绵等。这些分形图形具有严格的自相似性。无规则分形是指不光滑的，随机生成的分形，如蜿蜒曲折的海岸线，变换无穷的布朗运动轨迹等。这类曲线的自相似性是近似的或统计意义上的，这种自相似性只存于标度不变区域。
对于规则分形，其自相似性、标度不变性理论上是无限的（观测尺度可以趋于无限小）。不管我们怎样缩小（或放大）尺度（标度）去观察图形，其组成部分和原来的图形没有区别，也就是说它具有无限的膨胀和收缩对称性。因些对于这类分形，其计算方法比较简单，可以用缩小测量尺度的或者不断放大图形而得到。分形维数
D=lnN(λ)/ln(1/λ) （2-20）
如Cantor集，分数维D=ln2/ln3=0.631；Koch曲线分数维D=ln4/ln3=1.262; Sierpinski海绵分数维D=ln20/ln3=2.777。
对于不规则分形，它只具有统计意义下的自相似性。不规则分形种类繁多，它可以是离散的点集、粗糙曲线、多枝权的二维图形、粗糙曲面、以至三维的点集和多枝权的三维图形，下面介绍一些常用的测定方法[26]。
（1）尺码法
用某个选定尺码沿曲线以分规方式测量，保持尺码分规两端的落点始终在曲线上。不断改变尺码λ，得到一系列长度N（λ），λ越小、N越大。如果作lnN～lnλ图后得到斜率为负的直线，这表明存在如下的幂函数关系
N～λ-D（2-21）
上式也就是Mandelbrot在《分形：形状、机遇与维数》专著中引用的Richardson公式。Richardson是根据挪威、澳大利亚、南非、德国、不列颠西部、葡萄牙的海岸线丈量结果得出此公式的，使用的测量长度单位一般在1公里到4公里之间。海岸线绝对长度L被表示为：
L=Nλ～λ1-D（2-22）
他得到挪威东南部海岸线的分维D≈1.52，而不列颠西部海岸线的分维D≈1.3。这说明挪威的海岸线更曲折一些[27]。
（2）小岛法
如果粗糙曲线都是封闭的，例如海洋中的许多小岛，就可以利用周长-面积关系求分维，因此这个方法又被称为小岛法。
对于规则图形的周长与测量单位尺寸λ的一次方成正比，而面积A 则与λ的二次方成正比。通常我们可以把它们写成一个简单的比例关系：
P ∝A 1/2 (2-23)
对于二维空间内的不规则分形的周长和面积的关系显然更复杂一些，Mandelbrot 提出，应该用分形周长曲线来代替原来的光滑周长，从而给出了下述关系式：
1/1/2D 1/200[()](1)/[()][()]D P a D D A a A 1/-1λ=λ-λ=λλλ （2-24）
这里的分维D 大于1（周长光滑时D=1，上式转化成为（2.23）式），使P 的变化减缓，a 0是和岛的形状有关的常数，λ是测量尺寸，一般取λ为小于1的数
值（如取岛的最大直径为1），使因子λ（1-D ）/D 随测量尺寸λ减小而增大。作log[P(λ)/λ]～log[A(λ)1/2/λ]图，从其中直线部分的斜率的倒数，可以得到分维D 。
这个方法也可以推广到粗糙曲线（表面积-体积法）。
（3）计盒维数法[28]
这是一种常用的计算分形图形分维数的实用方法。取边长为r 的小盒子，把分形曲线覆盖起来。则有些小盒子是空的，有些小盒子覆盖了曲线的一部分。计数多少小盒子不是空的，所得的非空盒子数记为N （r ）。然后缩小盒子的尺寸，所得N （r ）自然要增大，当r→0时，得到分形维数：
0log ()lim log r N r D r
→=- （2-25） 实际计算中只能取有限的r ，通常的做法与尺码法类似，求一系列r 和N （r ），然后在双对数坐标中用最小二乘法拟合直线，所得直线的斜率即所求分形维数。
（４）结构函数法[29]
具有分形特征的时间序列能使其采样数据的结构函数满足：
242()[()()]D S z x z x C τττ-=+-= （2-26）
式中：
2[()()]z x z x τ+-表示差方的算术平均值。
τ是数据间隔的任意选择值。
针对若干尺度τ对分形曲线的离散信号计算出相应的Ｓ（τ），然后在对数坐标中得logS(τ)～log τ直线的斜率Ｗ，则分形维数：
42
W D -= （2-27） 2.2.4系统所采用的二种计算维数的方法
以上介绍的各种测量不规则分形的分维方法，在原理上都是利用了它们的自相似性和被测量是随测量尺度的改变而改变的特性。因此选择哪一种方法来测定和计算分维只能从实际问题出发，没有统一的标准。但在计算分维时存在的共同点是在计算原则上要求图形象素尽量多以及相似的层次尽量多。但实际图形往往达不到这样的要求，计算机模拟结果原则上可以有大得多的线性范围，但限于计算时，一般双对数图上的线性范围是2～3个量级。因此我们在实际的研究工作中，对研究对象使用分形或分维等概念时一定要注意它的适用范围。
下面介绍在系统中所使用的二种求分形的方法。
a 、 半方差法
半方差法用于复杂的分形曲线的计算，适用于对随机过程数据的处理。该方法简单易行，适合于计算机处理，是一种较实用的计算方法。
设在某一测量距离或测量时间序列上得到一族z （t ），且随机变量的平均差表示为：
1()[()()]m a z t z t t n
=-+?∑ （2-28） 其中：
m(a)为平均差；
z(t)为在t 位置函数曲线的测量值；
z(t+Δt)为在t+Δt 位置函数曲线的测量值；
Δt 为一对数据的间据
n 为数据对数。
方差表示为：
21()[()()]s a z t z t t n
=-+?∑ （2-29） 半方差表示为：
211()()[()()]22r a s a z t z t t n
==-+?∑ （2-30） 式中数据的对数n 的确定方法是：若以等间距Δt 连续测量某一距离的各点数值时，得到一随机数据z(1),z(2),…,z(k),如图2-6所示
当一对数据的间距t 1=Δt 时，数据的对数n=k-1,如图2-6 (a)所示。
当一对数据的间距t 2=2Δt 时，计算相应的半方差时，数据的对数n 2=k-2,如图2-6 (b)所示。
F ig 2-6 the definition of n in semi-variance
method
当试验数据较多时，往下依次类推。
每当改变一对数据的间距时，由式（2-30）可以得到相应的半方差r(a)。对于分形曲线，a 与r(a)存在如下的幂型关系：
r(a)∝h W (2-31)
t 1=Δt n 1=k-1 z(1) z(2) z(3) z(4) z(5) z(6) z(7) z(8) z(9) z(k) (a) t 2=2Δt n 2=k-2 z(1) z(2) z(3) z(4) z(5) z(6) z(7) z(8) z(9) z(k) (b) t 3=3Δt
其中，W 是幂指数，是分形维数D 的一种逼近，把h 和r(h)绘到双对数坐标图上，并进行线性回归，得到回归方程，其斜率即为W 。而斜率W 与分形维数D 有如下关系[23]：
W=4-2D （2-32）
则 42
W D -= （2-33） b 、 变换法
这是Dubuc 等[29]介绍的方法，在本质上它与计盒维数法相似，但对已知分形曲线运用此法得到的结果比计盒维数法准确，。后来Spanos 和Irene [25]把此方法推广应用于粗糙曲面，也得到很好的结果。
此法设置宽为Ｒ的矩形（盒子）覆盖到分形曲线上，矩形的高度由分形曲线在框内的最高点和最低点决定（图2-7），一步一步移动矩形遍及所有象素点，将所有矩形的高和宽相乘并且相加起来得到总面积S （R ），系列改变Ｒ的大小重复以上操作，得到一系列S （R ）。注意上述操作过程中矩形经过的范围应远远大于矩形的宽度。将
图2-7 变换法求分维
Fig 2-7 dimension calculating using
variation
S （R ）除以R ２得到N （R ）＝S （R ）／R ２，作lnN(R)～ln(1/R)曲线，取其中线性部分的斜率为分维Ｄ，因为在线性范围内存在N （R ）～R －Ｄ的关系。或者直接作lnS(R)～lnR 曲线，其中线性部分斜率为W ，并且由此斜率得到分维D 。
R R1
D＝２－W （2-34）
变换法也可以推广到粗糙曲面的分维计算。此时测量用的矩形被正方柱代替。变换法和计盒维数法在本质上是相同的，它们都是用不断改变尺寸的盒子去覆盖图形。其较为准确的原因在于它允许二维或三维的盒子数N（R）为非整数，同时N（R）也是遍及所有象素点得到的数值。
C语言的学习要从基础，100个经典的算法 题目：古典问题：有一对兔子，从出生后第3个月起每个月都生一对兔子，小兔子长到第三个月后每个月又生一对兔子，假如兔子都不死，问每个月的兔子总数为多少？ ________________________________________________________________ 程序分析：兔子的规律为数列1,1,2,3,5,8,13,21.... __________________________________________________________________ 程序源代码： main() { long f1,f2; int i; f1=f2=1; for(i=1;i<=20;i++) { printf("%12ld %12ld",f1,f2); if(i%2==0) printf("\n");/*控制输出，每行四个*/ f1=f1+f2;/*前两个月加起来赋值给第三个月*/ f2=f1+f2;/*前两个月加起来赋值给第三个月*/ } } 上题还可用一维数组处理，you try! 题目：判断101-200之间有多少个素数，并输出所有素数。 _________________________________________________________________ 程序分析：判断素数的方法：用一个数分别去除2到sqrt(这个数)，如果能被整除，则表明此数不是素数，反之是素数。 ___________________________________________________________________ 程序源代码： #include "math.h" main() { int m,i,k,h=0,leap=1; printf("\n"); for(m=101;m<=200;m++) { k=sqrt(m+1); for(i=2;i<=k;i++) if(m%i==0) {leap=0;break;} if(leap) {printf("%-4d",m);h++; if(h%10==0) printf("\n"); } leap=1; } printf("\nThe total is %d",h); } 题目：打印出所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字立方和等于该数本身。例如：153是一个“水仙花数”，因为153=1的三次方＋5的三次方＋3的三次方。 __________________________________________________________________ 程序分析：利用for循环控制100-999个数，每个数分解出个位，十位，百位。 ___________________________________________________________________ 程序源代码：
分形维数算法． 分形维数算法 分形包括规则分形和无规则分形两种。规则分形是指可以由简单的迭代或者是按一定规律所生成的分形，如Cantor集，Koch曲线，Sierpinski海绵等。这些分形图形具有严格的自相似性。无规则分形是指不光滑的，随机生成的分形，
如蜿蜒曲折的海岸线，变换无穷的布朗运动轨迹等。这类曲线的自相似性是近 似的或统计意义上的，这种自相似性只存于标度不变区域。 对于规则分形，其自相似性、标度不变性理论上是无限的（观测尺度可以趋于无限小）。不管我们怎样缩小（或放大）尺度（标度）去观察图形，其组成部分和原来的图形没有区别，也就是说它具有无限的膨胀和收缩对称性。因些对于这类分形，其计算方法比较简单，可以用缩小测量尺度的或者不断放大图形而得到。分形维数 D=lnN(λ)/ln(1/λ) （2-20） 如Cantor集，分数维D=ln2/ln3=0.631；Koch曲线分数维 D=ln4/ln3=1.262; Sierpinski海绵分数维D=ln20/ln3=2.777。 对于不规则分形，它只具有统计意义下的自相似性。不规则分形种类繁多，它可以是离散的点集、粗糙曲线、多枝权的二维图形、粗糙曲面、以至三维的[26]。点 集和多枝权的三维图形，下面介绍一些常用的测定方法（1）尺码法 用某个选定尺码沿曲线以分规方式测量，保持尺码分规两端的落点始终在曲线上。不断改变尺码λ，得到一系列长度N（λ），λ越小、N越大。如果作lnN～lnλ图后得到斜率为负的直线，这表明存在如下的幂函数关系
-D（2-21） N～λ上式也就是Mandelbrot在《分形：形状、机遇与维数》专著中引用的Richardson公式。Richardson是根据挪威、澳大利亚、南非、德国、不列颠西部、葡萄牙的海岸线丈量结果得出此公式的，使用的测量长度单位一般在1公里到4公里之间。海岸线绝对长度L被表示为： 1-D（2-22）L=Nλ～λ 他得到挪威东南部海岸线的分维D≈1.52，而不列颠西部海岸线的分维D≈[27]。。这说明挪威的海岸线更曲折一些1.3． ）小岛法（2面积如果粗糙曲线都是封闭的，例如海洋中的许多小岛，就可以利用周长-关系求分维，因此这个方法又被称为小岛法。则与λ的而面积A对于规则图形的周长与测量单位尺寸λ的一次方成正比， 二次方成正比。通常我们可以把它们写成一个简单的比例关系：1/2 (2-23) AP∝对于二维空间内的不规则分形的周长和面积的关系显然更复杂一些，提出，应该用分形周长曲线来代替原来的光滑周长，从而给出了下Mandelbrot 述关系式：21/??D??1/1/D2）（2-24)]?(?)]?[a?AP[(?)][??a(1?D)/DA(?00的P）式），使1（周长光滑时D=1，上式转化成为（2.23这里的分维D大于??的数1变化减缓，a是和岛的形状有关的常数，为小于是测量尺寸，一般取0/D）（1-D??减小而增大。作随测
看懂一个程序，分三步：1、流程；2、每个语句的功能；3、试数； 小程序：1、尝试编程去解决他；2、看答案；3、修改程序，不同的输出结果； 4、照答案去敲； 5、调试错误； 6、不看答案，自己把答案敲出来； 7、实在不会就背会。。。。。周而复始，反复的敲。。。。。 【程序1】 题目：有1、2、3、4个数字，能组成多少个互不相同且无重复数字的三位数？都是多少？ ============================================================== 【程序2】 题目：企业发放的奖金根据利润提成。利润(I)低于或等于10万元时，奖金可提10%；利润高 于10万元，低于20万元时，低于10万元的部分按10%提成，高于10万元的部分，可可提 成7.5%；20万到40万之间时，高于20万元的部分，可提成5%；40万到60万之间时高于 40万元的部分，可提成3%；60万到100万之间时，高于60万元的部分，可提成1.5%，高于 100万元时，超过100万元的部分按1%提成，从键盘输入当月利润I，求应发放奖金总数？ ============================================================== 【程序3】 题目：一个整数，它加上100后是一个完全平方数，再加上168又是一个完全平方数，请问该数是多少？============================================================== 【程序4】 题目：输入某年某月某日，判断这一天是这一年的第几天？ ============================================================== 【程序5】 题目：输入三个整数x,y,z，请把这三个数由小到大输出。 ============================================================== 【程序6】 题目：用*号输出字母C的图案。 ============================================================== 【程序7】 题目：输出特殊图案，请在c环境中运行，看一看，Very Beautiful! ============================================================== 【程序8】 题目：输出9*9口诀。 ============================================================== 【程序9】 题目：要求输出国际象棋棋盘。 ============================================================== 【程序10】 题目：打印楼梯，同时在楼梯上方打印两个笑脸。 -------------------------------------------------------------------------------- 【程序11】 题目：古典问题：有一对兔子，从出生后第3个月起每个月都生一对兔子，小兔子长到第三个月 后每个月又生一对兔子，假如兔子都不死，问每个月的兔子总数为多少？ ==============================================================
分形维数算法
分形维数算法 分形包括规则分形和无规则分形两种。规则分形是指可以由简单的迭代或者是按一定规律所生成的分形，如Cantor集，Koch曲线，Sierpinski海绵等。这些分形图形具有严格的自相似性。无规则分形是指不光滑的，随机生成的分形，如蜿蜒曲折的海岸线，变换无穷的布朗运动轨迹等。这类曲线的自相似性是近似的或统计意义上的，这种自相似性只存于标度不变区域。 对于规则分形，其自相似性、标度不变性理论上是无限的（观测尺度可以趋于无限小）。不管我们怎样缩小（或放大）尺度（标度）去观察图形，其组成部分和原来的图形没有区别，也就是说它具有无限的膨胀和收缩对称性。因些对于这类分形，其计算方法比较简单，可以用缩小测量尺度的或者不断放大图形而得到。分形维数 D=lnN(λ)/ln(1/λ) （2-20） 如Cantor集，分数维D=ln2/ln3=0.631；Koch曲线分数维D=ln4/ln3=1.262; Sierpinski海绵分数维D=ln20/ln3=2.777。 对于不规则分形，它只具有统计意义下的自相似性。不规则分形种类繁多，它可以是离散的点集、粗糙曲线、多枝权的二维图形、粗糙曲面、以至三维的点集和多枝权的三维图形，下面介绍一些常用的测定方法[26]。 （1）尺码法 用某个选定尺码沿曲线以分规方式测量，保持尺码分规两端的落点始终在曲线上。不断改变尺码λ，得到一系列长度N（λ），λ越小、N越大。如果作lnN～lnλ图后得到斜率为负的直线，这表明存在如下的幂函数关系 N～λ-D（2-21） 上式也就是Mandelbrot在《分形：形状、机遇与维数》专著中引用的Richardson公式。Richardson是根据挪威、澳大利亚、南非、德国、不列颠西部、葡萄牙的海岸线丈量结果得出此公式的，使用的测量长度单位一般在1公里到4公里之间。海岸线绝对长度L被表示为： L=Nλ～λ1-D（2-22） 他得到挪威东南部海岸线的分维D≈1.52，而不列颠西部海岸线的分维D≈1.3。这说明挪威的海岸线更曲折一些[27]。
2008-02-18 18:48 【程序1】 题目：有1、2、3、4个数字，能组成多少个互不相同且无重复数字的三位数？都是多少？ 1.程序分析：可填在百位、十位、个位的数字都是1、2、3、4。组成所有的排列后再去 掉不满足条件的排列。 2.程序源代码： main() { int i,j,k; printf("\n"); for(i=1;i<5;i++) ／*以下为三重循环*/ for(j=1;j<5;j++) for (k=1;k<5;k++) { if (i!=k&&i!=j&&j!=k) /*确保i、j、k三位互不相同*/ printf("%d,%d,%d\n",i,j,k); } } ============================================================== 【程序2】 题目：企业发放的奖金根据利润提成。利润(I)低于或等于10万元时，奖金可提10%；利润高 于10万元，低于20万元时，低于10万元的部分按10%提成，高于10万元的部分，可可提 成7.5%；20万到40万之间时，高于20万元的部分，可提成5%；40万到60万之间时高于 40万元的部分，可提成3%；60万到100万之间时，高于60万元的部分，可提成1.5%，高于 100万元时，超过100万元的部分按1%提成，从键盘输入当月利润I，求应发放奖金总数？ 1.程序分析：请利用数轴来分界，定位。注意定义时需把奖金定义成长整型。 2.程序源代码： main() { long int i; int bonus1,bonus2,bonus4,bonus6,bonus10,bonus; scanf("%ld",&i); bonus1=100000*0.1;bonus2=bonus1+100000*0.75; bonus4=bonus2+200000*0.5; bonus6=bonus4+200000*0.3; bonus10=bonus6+400000*0.15; if(i<=100000)
计算题 1．已知某地铁线路车辆定员每节240人，列车为6节编组，高峰小时满载率为120％，且单向最大断面旅客数量为29376人，试求该小时内单向应开行的列车数。 2、已知某地铁线路采用三显示带防护区段的固定闭塞列车运行控制方式，假设各闭塞分区长度相等，均为1000米，已知列车长 度为420米，列车制动距离为100米，列车运行速度为70km/h，制动减速度为2米/秒2，列车启动加速度为1.8米/秒2，列车最大停站时间为40秒。试求该线路的通过能力是多少？ 若该线路改成四显示自动闭塞，每个闭塞分区长度为600米，则此时线路的通过能力是多少？ 3．已知某地铁线路采用移动闭塞列车运行控制方式，已知列车长度为420米，车站闭塞分区为750米，安全防护距离为 200米，列车进站规定速度为60km/h，制动空驶时间为1.6秒，制动减速度为2米/秒2，列车启动加速度为1.8米/秒2，列车最大停站时间为40秒。试求该线路的通过能力是多少？ 4．已知某地铁线路为双线线路，列车采用非自动闭塞的连发方式运行，已知列车在各区间的运行时分和停站时分如下表，线路的连发间隔时间为12秒。试求该线路的通过能力是多少？
5．已知地铁列车在某车站采用站后折返，相关时间如下：前一列车离去时间1.5分钟，办理进路作业时间0.5分钟，确认信号时间0.5分钟，列车出折返线时间1.5分钟，停站时间1分钟。试计算该折返站通过能力。 6．已知某终点折返站采用站前交替折返，已知列车直到时间 为40秒，列车侧到时间为1分10秒，列车直发时间为40秒，列车侧发时间为1分20秒，列车反应时间为10秒， 办理接车进路的时间为15秒，办理发车进路的时间为15秒。试分别计算考虑发车时间均衡时和不考虑发车时间均衡时，该折返站的折返能力是多少？ 7．已知线路上有大小交路两种列车，小交路列车在某中间折返 站采用站前折返（直到侧发），已知小交路列车侧发时间为1分20秒，办理接车进路的时间为15秒，办理发车进路的时间为15秒，列车反应时间为10秒，列车直到时间为25 秒，列车停站时间为40秒；长交路列车进站时间为25秒。试分别计算该中间折返站的最小折返能力和最大折返能力分别是多少？ 8．已知线路上有大小交路两种列车，小交路列车在某中间折返站采用站后折返，已知小交路列车的相关时分为：列车驶出车站 闭塞分区时间为1分15秒，办理出折返线调车进路的时间 为20秒，列车从折返线至车站出发正线时间为40秒，列车反应时间为10秒，列车停站时间为40秒。
遥感图象分形维数的几种估计算法研究1 张凯选1，郭嗣琮2 1辽宁工程技术大学测绘与地理科学学院，辽宁阜新（123000） 2辽宁工程技术大学理学院，辽宁阜新（123000） E-mail：zhangkaixuan@https://www.doczj.com/doc/0412894294.html, 摘要：美籍法国数学家曼德布罗特（B.Mandelbrot）首次引入分形这个新术语，今天分形理论已经成为一门描述自然界中许多不规则事物规律性的科学，在遥感影象学中也有很大的用途。在研究遥感图像的分形维数时,通常把图像看作一个由许多像素点的灰度值构成的曲面来进行估算和分析，本文给出了遥感图象分形维数的几种估算方法，并作了相关实验。关键词：分形，分形维数，遥感图象 中图分类号：TP7 1．引言 分形理论始创立于20世纪70年代中期[1],创立伊始就引起人们极大的兴趣,与耗散结构、混沌并称为70年代科学史上的三大发现。作为一门独立的学科,该理论只有大约30多年的历史。 基于对复杂景物自相似性的描述,Mandelbrot创立了分形几何学理论,提出用分形维数( fractal dimension)D来度量自然现象的不规则程度。分形理论借助相似性原理洞察隐藏于混乱现象中的精细结构,为人们从局部认识整体、从有限认识无限提供新的方法论,为不同的学科发现的规律提供了崭新的语言和定量的描述,为现代科学技术提供了新的思想方法。近年来,分形理论在自然科学、社会科学以及遥感的许多领域中得到了广泛的应用,并逐步成为连结现代各学科的纬线。 2．分形与分形维数的定义 美籍法国数学家曼德布罗特(B.Mandelbrot) 于1967 年在《科学》杂志上发表了一篇题为“英国的海岸线有多长? 统计自相似性与分数维数” 的论文[2], 通常被认为是“分形”学科诞生的标志。自然界的许多物体在某一范围内都具有统计的自相似性,即每一部分都被认为是整体的一个缩小图像。曼德布罗特在随后两本著作《自然界的分形几何学》和《分形、形状、机遇与维数》中第一次提出了fractal这个英文词,其原意是“不规则的”、“分数的”、“支离破碎的”物体,并阐述分形理论的基本思想,即分形研究的对象是具有自相似性的无序系统,其维数的变化是连续的。 关于分形,目前还没有严格的数学定义,只能给出描述性的定义。粗略地说,分形是对没有特征长度但具有一定意义下的自相似性图形和结构的总称。它具有两个基本性质:自相似性和标度不变性。自相似性是指局部是整体成比例缩小的性质。形象地说,就是当用不同倍数的照相机拍摄研究对象时,无论放大倍数如何改变,看到的照片都是相似的(统计意义) ,而从相片上也无法断定所用的相机的倍数,即标度不变性或全息性。严格按一定的数学方法生成的许多经典的分形(如图1) 具有严格的自相似性,称之为有规分形。而一般情况下的分形都是无规分形,即自相似性并不是严格的,只是统计意义下的自相似性,其局部经放大或缩小操作可能得到与整体完全不同的表现形式,但表征自相似结构或系统的定量参数如分形维数,并 本课题得到辽宁工程技术大学青年基金（05－124），辽宁省教育厅基金项目（05L181）,辽宁省高等学校重点实验室项目基金（20060370）的资助。
线性方程组AX=B 的数值计算方法实验 一、 实验描述： 随着科学技术的发展，线性代数作为高等数学的一个重要组成部分， 在科学实践中得到广泛的应用。本实验的通过C 语言的算法设计以及编程，来实现高斯消元法、三角分解法和解线性方程组的迭代法（雅可比迭代法和高斯-赛德尔迭代法），对指定方程组进行求解。 二、 实验原理： 1、高斯消去法： 运用高斯消去法解方程组，通常会用到初等变换，以此来得到与原系数矩阵等价的系数矩阵，达到消元的目的。初等变换有三种：(a)、(交换变换)对调方程组两行；(b)、用非零常数乘以方程组的某一行；(c)、将方程组的某一行乘以一个非零常数，再加到另一行。 通常利用(c)，即用一个方程乘以一个常数，再减去另一个方程来置换另一个方程。在方程组的增广矩阵中用类似的变换，可以化简系数矩阵，求出其中一个解，然后利用回代法，就可以解出所有的解。 2、选主元： 若在解方程组过程中，系数矩阵上的对角元素为零的话，会导致解出的结果不正确。所以在解方程组过程中要避免此种情况的出现，这就需要选择行的判定条件。经过行变换，使矩阵对角元素均不为零。这个过程称为选主元。选主元分平凡选主元和偏序选主元两种。平凡选主元： 如果()0p pp a ≠，不交换行；如果()0p pp a =，寻找第p 行下满足() 0p pp a ≠的第一 行，设行数为k ，然后交换第k 行和第p 行。这样新主元就是非零主元。偏序选主元：为了减小误差的传播，偏序选主元策略首先检查位于主对角线或主对角线下方第p 列的所有元素，确定行k ，它的元素绝对值最大。然后如果k p >，则交换第k 行和第p 行。通常用偏序选主元，可以减小计算误差。 3、三角分解法： 由于求解上三角或下三角线性方程组很容易所以在解线性方程组时，可将系数矩阵分解为下三角矩阵和上三角矩阵。其中下三角矩阵的主对角线为1，上三角矩阵的对角线元素非零。有如下定理： 如果非奇异矩阵A 可表示为下三角矩阵L 和上三角矩阵U 的乘积： A LU = (1) 则A 存在一个三角分解。而且，L 的对角线元素为1，U 的对角线元素非零。得到L 和U 后，可通过以下步骤得到X ： (1)、利用前向替换法对方程组LY B =求解Y 。 (2)、利用回代法对方程组UX Y =求解X 。 4、雅可比迭代：
车站通过能力 车站通过能力是在车站现有设备条件下，采用合理的技术作业过程，一昼夜能接发和方向的货物（旅客）列车数和运行图规定的旅客（货物）列车数。 车站通过能力包括咽喉通过能力和到发线通过能力。 咽喉通过能力是指车站某咽喉区各衔接方向接、发车进路咽喉道岔组通过能力之和，咽喉道岔通过能力是指在合理固定到发线使用方案及作业进路条件下，某衔接方向接、发车进路上最繁忙的道岔组一昼夜能够接、发该方向的货物（旅客）列车数和运行图规定的旅客（货物）列车数。 到发线通过能力是指到达场、出发场、通过场或到发场内办理列车到发作业的线路，采用合理的技术作业过程和线路固定使用方案，一昼夜能够接、发各衔接方向的货物（旅客）列车数和运行图规定的旅客（货物）列车数。 车站咽喉通过能力计算 咽喉占用时间标准 表咽喉道岔占用时间表 顺序作业名称时间标准 （min） 顺序作业名称 时间标准 （min） 1 货物列车接车占用6~8 4 旅客列车出发占用4~6 2 旅客列车接车占用5~7 5 单机占用2~4 3 货物列车出发占用5~7 6 调车作业占用4~6 道岔组占用时间计算 表到发线固定使用方案 线路编号固定用途 一昼夜 接发列车数 线路 编号 固定用途 一昼夜 接发列车数 1 接甲到乙、丙旅客列车8 7 接乙到甲直通、区段货物列车9 4 接乙到甲旅客列车 5 8 接甲、乙到丙直通、区段货物列车10 接丙到甲旅客列车 3 9 接丙到甲、乙直通、区段货物列车10 5 接甲到乙直通、区段货物列车11 10 接发甲、乙、丙摘挂货物列车10 表甲端咽喉区占用时间计算表 编号作业进路名称 占用 次数 每次 占用时间 总占用 时间 咽喉区道岔组占用时间 1 3 5 7 9 固定作业 1 1道接甲-乙,丙旅客列车8 7 56 56 2 4道发乙-甲旅客列车 5 6 30 30 30 3 4道发丙-甲旅客列车 3 6 18 30 30 5 往机务段送车 3 6 18 18 6 从机务段取车 2 6 12 12
分形维数算法 分形包括规则分形和无规则分形两种。规则分形是指可以由简单的迭代或者是按一定规律所生成的分形，如Cantor集，Koch曲线，Sierpinski海绵等。这些分形图形具有严格的自相似性。无规则分形是指不光滑的，随机生成的分形，如蜿蜒曲折的海岸线，变换无穷的布朗运动轨迹等。这类曲线的自相似性是近似的或统计意义上的，这种自相似性只存于标度不变区域。 对于规则分形，其自相似性、标度不变性理论上是无限的（观测尺度可以趋于无限小）。不管我们怎样缩小（或放大）尺度（标度）去观察图形，其组成部分和原来的图形没有区别，也就是说它具有无限的膨胀和收缩对称性。因些对于这类分形，其计算方法比较简单，可以用缩小测量尺度的或者不断放大图形而得到。分形维数 D=lnN(λ)/ln(1/λ) （2-20） 如Cantor集，分数维D=ln2/ln3=0.631；Koch曲线分数维D=ln4/ln3=1.262; Sierpinski海绵分数维D=ln20/ln3=2.777。 对于不规则分形，它只具有统计意义下的自相似性。不规则分形种类繁多，它可以是离散的点集、粗糙曲线、多枝权的二维图形、粗糙曲面、以至三维的点集和多枝权的三维图形，下面介绍一些常用的测定方法[26]。 （1）尺码法 用某个选定尺码沿曲线以分规方式测量，保持尺码分规两端的落点始终在曲线上。不断改变尺码λ，得到一系列长度N（λ），λ越小、N越大。如果作lnN～lnλ图后得到斜率为负的直线，这表明存在如下的幂函数关系 N～λ-D（2-21） 上式也就是Mandelbrot在《分形：形状、机遇与维数》专著中引用的Richardson公式。Richardson是根据挪威、澳大利亚、南非、德国、不列颠西部、葡萄牙的海岸线丈量结果得出此公式的，使用的测量长度单位一般在1公里到4公里之间。海岸线绝对长度L被表示为： L=Nλ～λ1-D（2-22） 他得到挪威东南部海岸线的分维D≈1.52，而不列颠西部海岸线的分维D≈1.3。这说明挪威的海岸线更曲折一些[27]。
根据交叉口的现场交通调查数据，通过各流向流量的构成关系，可推得各路段流量，从而得到饱和度V/C 比。路段通行能力的确定采用建设部《城市道路设计规范》（CJJ 37-90）的方法，该方法的计算公式为：单条机动车道设计通行能力n C N N a ????=ηγ0，其中N a 为车道可能通行能力，该值由设计车速来确定，如表2.2所示。 表2.13 一条车道的理论通行能力 其中γ为自行车修正系数，有机非隔离时取1，无机非隔离时取0.8。η为车道宽度影响系数，C 为交叉口影响修正系数，取决于交叉口控制方式及交叉口间距。修正系数由下式计算： s 为交叉口间距(m)，C 0为交叉口有效通行时间比。 车道修正系数采用表 2.3所示 表2.3 车道数修正系数采用值 路段服务水平评价标准采用美国《道路通行能力手册》，如表2.4所示 表2.4 路段服务水平评价标准
由路段流量的调查结果，并且根据交叉口的间距、路段等级、车道数等对路段的通行能力进行了修正。在此基础上对路段的交通负荷进行了分析。 路段机动车车道设计通行能力的计算如下： δ m c p m k a N N = （1） 式中: m N —— 路段机动车单向车道的设计通行能力（pcu/h ） p N —— 一条机动车车道的路段可能通行能力（pcu/h ） c a —— 机动车通行能力的分类系数，快速路分类系数为0.75；主干道分类 系数为0.80；次干路分类系数为0.85；支路分类系数为0.90。 m k —— 车道折减系数，第一条车道折减系数为 1.0；第二条车道折减系数 为0.85；第三条车道折减系数为0.75；第四条车道折减系数为0.65.经过累加，可取单向二车道 m k ＝1.85；单向三车道 m k ＝2.6；单向四车道 m k ＝3.25； δ—— 交叉口影响通行能力的折减系数，不受交叉口影响的道路（如高架 道路和地面快速路）δ＝1；该系数与两交叉口之间的距离、行车速度、绿信比和车辆起动、制动时的平均加、减速度有关，其计算公式如下： ?+++= b v a v v l v l 2/2///δ （2） l —— 两交叉口之间的距离（m ）； a —— 车辆起动时的平均加速度，此处取为小汽车0.82/s m ； b —— 车辆制动时的平均加速度，此处取为小汽车1.662/s m ； ?—— 车辆在交叉口处平均停车时间，取红灯时间的一半。 Np 为车道可能通行能力，其值由路段车速来确定： 表4.1 Np 的确定
C语言经典算法大全 老掉牙 河内塔 费式数列 巴斯卡三角形 三色棋 老鼠走迷官（一） 老鼠走迷官（二） 骑士走棋盘 八个皇后 八枚银币 生命游戏 字串核对 双色、三色河内塔 背包问题（Knapsack Problem） 数、运算 蒙地卡罗法求PI Eratosthenes筛选求质数 超长整数运算（大数运算） 长PI 最大公因数、最小公倍数、因式分解 完美数 阿姆斯壮数 最大访客数 中序式转后序式（前序式） 后序式的运算 关于赌博 洗扑克牌（乱数排列） Craps赌博游戏 约瑟夫问题（Josephus Problem） 集合问题 排列组合 格雷码（Gray Code） 产生可能的集合
m元素集合的n个元素子集 数字拆解 排序 得分排行 选择、插入、气泡排序 Shell 排序法- 改良的插入排序Shaker 排序法- 改良的气泡排序Heap 排序法- 改良的选择排序快速排序法（一） 快速排序法（二） 快速排序法（三） 合并排序法 基数排序法 搜寻 循序搜寻法（使用卫兵） 二分搜寻法（搜寻原则的代表）插补搜寻法 费氏搜寻法 矩阵 稀疏矩阵 多维矩阵转一维矩阵 上三角、下三角、对称矩阵 奇数魔方阵 4N 魔方阵 2(2N+1) 魔方阵
1.河内之塔 说明河内之塔(Towers of Hanoi)是法国人M.Claus(Lucas)于1883年从泰国带至法国的，河内为越 战时北越的首都，即现在的胡志明市；1883年法国数学家Edouard Lucas曾提及这个故事，据说创世纪时Benares有一座波罗教塔，是由三支钻石棒（Pag）所支撑，开始时神在第一根棒上放置64个由上至下依由小至大排列的金盘（Disc），并命令僧侣将所有的金盘从第一根石棒移至第三根石棒，且搬运过程中遵守大盘子在小盘子之下的原则，若每日仅搬一个盘子，则当盘子全数搬运完毕之时，此塔将毁损，而也就是世界末日来临之时。 解法如果柱子标为ABC，要由A搬至C，在只有一个盘子时，就将它直接搬至C，当有两个盘 子，就将B当作辅助柱。如果盘数超过2个，将第三个以下的盘子遮起来，就很简单了，每次处理两个盘子，也就是：A->B、A ->C、B->C这三个步骤，而被遮住的部份，其实就是进入程式的递回处理。事实上，若有n个盘子，则移动完毕所需之次数为2^n - 1，所以当盘数为64时，则所需次数为：264- 1 = 18446744073709551615为5.05390248594782e+16年，也就是约5000世纪，如果对这数字没什幺概念，就假设每秒钟搬一个盘子好了，也要约5850亿年左右。 #include
简单分形及维数的研究 （河南大学，物理与电子学院，物理学，河南开封，475004）摘要：本文介绍了分形、维数的相关知识，并以简单分形做例子进行了演示，又求得了Sierpinski三角分形及埃侬映射的维数。 关键词：分形，维数，程序设计。 一、分形 分形（fractal）是指由各部分组成的形态，每个部分以某种方式与整体相似。对这一描述加以引伸，它可以包括以下含义： 分形可以是几何图形，也可以是由“功能”或“信息”架起的数理模型；分形可以同时具有形态、功能和信息三方面的自相似性，也可以只有其中某一方面的自相似性。 分形的创建历史： (1)曼德勃罗在美国《科学》杂志上发表论文《英国的海岸线有多长》震惊学术界(1967 年)。 (2)法兰西学院讲演报(1973年)。 (3)“病态”“数学怪物”命名——分形(Fractal)(1975年)。 (4)法文版《分形对象：形、机遇和维数》出版(1975年)。 (5)英文版《分形：形、机遇和维数》出版(1977年)。 (6)英文版《大自然的几何学》出版(1982年) 。 分形是由Mandelbrot在20世纪70年代为了表征复杂图形和复杂过程而引入自然领域的。原意是破碎的、不规则的物体。分形分为两类，规则分形，又称决定类的分形，它是按一定的规则构造出的具有严格自相思的分形；另一类是无规则的分形，它是在生长现象中和许多物理问题中产生的分形，其特点是不具备严格意义上的自相似，只是在统计意义上是自相似的。本文研究的是规则分形。 有以上可知，自相似性是分形最大的几何特征。下面我们就科赫曲线和Sierpinski对此进行讨论。 1、科赫曲线 科赫曲线的生成方法：把一条曲线三等分，中间的一段用夹角为60的折线替代，得到第一个生成元；把第一个生成元中的每一条直线都用生成元迭代，得到第二个生成元；经过无数次迭代，即可得到科赫曲线。 实现程序如下： s=[0,1];t=[0,0];n=8; for j=1:n
60.题目:古典问题:有一对兔子,从出生后第3个月 起每个月都生一对兔子,小兔 子长到第三个月后每个月又生一对兔子,假如兔子都不死,问每个月的兔子总 数 为多少? _________________________________________________________________ _ 程序分析:兔子的规律为数列1,1,2,3,5,8,13,21.... _________________________________________________________________ __ 程序源代码: main() { long f1,f2; int i; f1=f2=1; for(i=1;i<=20;i++) { printf("%12ld %12ld",f1,f2); if(i%2==0) printf("\n");/*控制输出,每行四个*/
f1=f1+f2;/*前两个月加起来赋值给第三个月*/ f2=f1+f2;/*前两个月加起来赋值给第三个月*/ } } 上题还可用一维数组处理,you try! 61.题目:判断101-200之间有多少个素数,并输出所有素数。 _________________________________________________________________ _ 1 程序分析:判断素数的方法:用一个数分别去除2到sqrt(这个数),如果能被 整 除,则表明此数不是素数,反之是素数。 _________________________________________________________________ __ 程序源代码: #include "math.h" main() { int m,i,k,h=0,leap=1;
题目：有1、2、3、4个数字，能组成多少个互不相同且无重复数字的三位数？都是多少？ 1.程序分析：可填在百位、十位、个位的数字都是1、2、3、4。组成所有的排列后再去 掉不满足条件的排列。 2.程序源代码： main() { int i,j,k; printf("\n"); for(i=1;i<5;i++) ／*以下为三重循环*/ for(j=1;j<5;j++) for (k=1;k<5;k++) { if (i!=k&&i!=j&&j!=k) printf("%d,%d,%d\n",i,j,k); } } ============================================================== 【程序2】 题目：企业发放的奖金根据利润提成。利润(I)低于或等于10万元时，奖金可提10%；利润高 于10万元，低于20万元时，低于10万元的部分按10%提成，高于10万元的部分，可可提成7.5%；20万到40万之间时，高于20万元的部分，可提成5%；40万到60万之间时高于40万元的部分，可提成3%；60万到100万之间时，高于60万元的部分，可提成1.5%，高于 100万元时，超过100万元的部分按1%提成，从键盘输入当月利润I，求应发放奖金总数？ 1.程序分析：请利用数轴来分界，定位。注意定义时需把奖金定义成长整型。 2.程序源代码： main() { long int i; int bonus1,bonus2,bonus4,bonus6,bonus10,bonus; scanf("%ld",&i); bonus1=100000*0.1;bonus2=bonus1+100000*0.75; bonus4=bonus2+200000*0.5; bonus6=bonus4+200000*0.3; bonus10=bonus6+400000*0.15; if(i<=100000) bonus=i*0.1; else if(i<=200000) bonus=bonus1+(i-100000)*0.075; else if(i<=400000) bonus=bonus2+(i-200000)*0.05; else if(i<=600000) bonus=bonus4+(i-400000)*0.03;
路区间通过能力计算办法 1984年10月1日，铁道部 第一章总则 第1条为了保证铁路完成和超额完成不断增长的运输任务，以适应国民经济发展和国防建设对铁路运输的需要，铁路必须大力加强运输组织工作，采取有效措施，积极提高铁路线路通过能力。 铁路线路通过能力，是根据现有技术设备、行车组织方法及规定的技术作业过程确定的在一昼夜内所能通过的最大列车对数或列数。 铁路线路通过能力，系按区间、车站、机务段设备和整备设备、车站给水设备、电气化铁路的供电设备分别确定，以其中最小的通过能力，作为该区段的限制通过能力。 为了计算铁路区间通过能力，本办法规定了铁路区间通过能力的计算办法。 第2条铁路区间通过能力，是指每一区间在一昼夜内所能通过的列车数量（列数或对数）。 区间通过能力的大小，在一定的行车组织条件下，主要取决于正线数目、区间长度、线路纵断面、信联闭设备、牵引机车类型和列车运行速度等因素。 第3条计算区间通过能力时，应先计算平行运行图通过能力，再计算非平行运行图通过能力。 平行运行图通过能力，一般应按货物列车对数或列数计算；非平行运行
图通过能力，系在规定旅客列车数量的基础上，以扣除系数的方法计算出旅客列车和货物列车的对数或列数。 第4条铁路区间通过能力，由各铁路局或分局负责计算，并填制区间通过能力计算表及区间通过能力汇总表，经铁路局审核后报铁道部运输局。 第5条本办法系根据我国铁路现有技术设备条件及多年来编制和执行列车运行图的经验，规定了铁路区间通过能力的一般计算方法。个别特殊情况，由铁路局根据具体情况和特点，进行图解和计算。 第二章平行运行图区间通过能力 第6条平行运行图区间通过能力，应分别对区段内每一区间计算。运行图周期最大的区间通过能力，即为该区段的限制区间通过能力。 运行图周期，是指一定类型运行图的一组列车占用区间的总时间。其组成因素，在非自动闭塞区段包括：列车区间运行时分，起停车附加时分及列车在车站的间隔时间。在自动闭塞区段为追踪列车间隔时间。 平行运行图区间通过能力的基本关系式如下： 1440 Ｎ＝―――― (1) Ｔ周 式中：Ｎ――平行运行图通过能力（对数或列数）； 1440――一昼夜时分； Ｔ周――运行图周期。 电力牵引区段，由于每日须进行接触网检修，因此，其计算公式为：
1.河床稳定计算及河相分析 1.1.河床稳定计算 河床稳定指标可采用横向稳定指标、纵向稳定指标及综合稳定指标 3 种形式分析，以确定河道特性。 1.1.1.河道横向稳定分析 河道横向稳定系数按下式计算: 式中： 横向稳定系数； Q造床流量， m3/s ； J河床比降； B 相当于造床流量的平摊河宽，m。 1.1. 2.河道纵向稳定分析 水流对河床泥沙的拖曳力与床面泥沙抵抗运动的摩阻力之间的相互作用，决定河床的纵向稳定性。根据黄河水利出版社出版《治河及泥沙工程》中河道纵向稳定系数采用爱因斯坦水流强度函数按下式计算： 式中： 纵向稳定系数； D床沙平均粒径，mm； J河床纵比降； H河流平摊水深，m。
1.1.3.综合稳定指标 综合稳定指标是综合考虑河床的纵、横向稳定性。建议采用的公式为 2 （b）*h 1.2.河床演变分析与河相关系 调查工程区河道历史主流及河道变迁，分析工程区河道形态。共分为蜿蜒型河道、游荡型河道两种形式。 蜿蜒型河段一般凹岸崩退，凸岸淤长，凹岸深槽和过渡段浅滩在年内发生互相交替的冲淤变化。 游荡型河道的河岸及河床抗冲性较差，从长距离来看河道往往呈藕节状，其中窄段水流 归顺，有控制河势的作用，宽段则河床宽浅，洲滩密布，汊道交织，水流散乱，主流迁徙不 定。河道的平面状态可用“宽、浅、散、乱”四个字概括。 在水流长期作用下形成的河床，其形态有一定的规律，大量资料表明，表征河床形态的 水深、河宽、比降等，与来水来沙条件及河床地质条件之间，有一定函数关系，这种关系便 称为河相关系。 根据俄罗斯国立水文所提出公式，河道横断面河相关系公式为: B H 式中 : ξ 河相相关系数； B 造床流量下的水面宽（m）； H造床流量下的平均水深（m）； （蜿蜒型河道ζ 约为2~4，较为顺直的过渡性河段约为8~12，游荡型河道ζ 约为20~30）2.护岸结构设计 2.1.护岸顶高程确定 根据《堤防工程设计规范》（GB50286-2013）（以下简称《堤防规范》）要求，堤顶高程为设计洪水位加超高值确定。堤顶超高按下式计算：
经典的滤波算法(转) 1、限幅滤波法（又称程序判断滤波法） A、方法： 根据经验判断，确定两次采样允许的最大偏差值（设为A） 每次检测到新值时判断： 如果本次值与上次值之差<=A,则本次值有效 如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值 B、优点： 能有效克服因偶然因素引起的脉冲干扰 C、缺点 无法抑制那种周期性的干扰 平滑度差 2、中位值滤波法 A、方法： 连续采样N次（N取奇数） 把N次采样值按大小排列 取中间值为本次有效值 B、优点： 能有效克服因偶然因素引起的波动干扰 对温度、液位的变化缓慢的被测参数有良好的滤波效果 C、缺点： 对流量、速度等快速变化的参数不宜 3、算术平均滤波法 A、方法： 连续取N个采样值进行算术平均运算 N值较大时：信号平滑度较高，但灵敏度较低 N值较小时：信号平滑度较低，但灵敏度较高 N值的选取：一般流量，N=12；压力：N=4 B、优点： 适用于对一般具有随机干扰的信号进行滤波 这样信号的特点是有一个平均值，信号在某一数值范围附近上下波动 C、缺点： 对于测量速度较慢或要求数据计算速度较快的实时控制不适用 比较浪费RAM
递推平均滤波法对偶然出现的脉冲性干扰的抑制作用较差 4、递推平均滤波法（又称滑动平均滤波法） A、方法： 把连续取N个采样值看成一个队列 队列的长度固定为N 每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进先出原则) 把队列中的N个数据进行算术平均运算,就可获得新的滤波结果 N值的选取：流量，N=12；压力：N=4；液面，N=4~12；温度，N=1~4 B、优点： 对周期性干扰有良好的抑制作用，平滑度高 适用于高频振荡的系统 C、缺点： 灵敏度低 对偶然出现的脉冲性干扰的抑制作用较差 不易消除由于脉冲干扰所引起的采样值偏差 不适用于脉冲干扰比较严重的场合 比较浪费RAM 5、中位值平均滤波法（又称防脉冲干扰平均滤波法） A、方法： 相当于“中位值滤波法”+“算术平均滤波法” 连续采样N个数据，去掉一个最大值和一个最小值 然后计算N-2个数据的算术平均值 N值的选取：3~14 B、优点： 融合了两种滤波法的优点 对于偶然出现的脉冲性干扰，可消除由于脉冲干扰所引起的采样值偏差 C、缺点： 测量速度较慢，和算术平均滤波法一样 比较浪费RAM 6、限幅平均滤波法 A、方法： 相当于“限幅滤波法”+“递推平均滤波法” 每次采样到的新数据先进行限幅处理， 再送入队列进行递推平均滤波处理 B、优点： 融合了两种滤波法的优点 对于偶然出现的脉冲性干扰，可消除由于脉冲干扰所引起的采样值偏差 C、缺点： 比较浪费RAM