东南大学计算机科学试验中心,受迫振动实验数据处理的研究-东南大学物理实验中心.doc...

受迫振动实验数据处理的研究-东南大学物理实验中心.doc

受迫振动实验数据处理的研究

牛星

(东南大学 计算机科学与工程学院,南京 210096)

通过计算机软件Mathematica 6对实验记录数据进行处理，得到拟合曲线和阻尼系数．指出不确定的固有角频率造成的系统误差很难避免．对实验处理方法的选取做出建议．

波尔共振仪；受迫振动；数据处理；拟合曲线；系统误差

Researches in Experimental Data Processing of Forced Vibration

NiuXing

(Department of Computer Science and Engineering, Southeast University, Nanjing 210096)

Through using computer software Mathematica6 to process experimental data, we can get fitting curve and damping coefficient. Pointing out the system error caused by uncertain natural angular frequency is very difficult to avoid. Making a suggestion to the selection of experiment data processing method.

Pohl resonator; Forced vibration; Data processing; Fitting curve; System error

大学物理实验“受迫振动的研究”在数据处理方面较为复杂，所得结果极大程度上取决于作图的精度．我考虑使用Mathematica 6软件进行数据辅助处理，以期得到较为理想的结果．

计算机拟合

幅频特性曲线拟合法

受迫振动的振幅和频率的关系如下：

作者简介：牛星,1988年，男，本科生, oxstar.student@.

(1)

经过上述变换即可得与的关系方程．为做出拟合曲线，可以将和设为待定系数．令=， =，根据最小二乘法原理，容易得到、的值．

幅频特性曲线拟合法

受迫振动的相位和频率的关系如下：

(2)

根据与幅频特性曲线相同的分析方法，得到系数=．

拟合曲线

选取实验测得的阻尼1和阻尼2两组数据，使用Mathematica 6软件，按上述方法得到两种阻尼各两条拟合曲线(图1--图4)．

图1 阻尼1幅频特性拟合曲线

图2 阻尼1相频特性拟合曲线

图3 阻尼2幅频特性拟合曲线

图4 阻尼2相频特性拟合曲线

数据处理

“幅频特性曲线拟合法”和“相频特性曲线拟合法”所能求得的分别是和，这里应该如何取值是个棘手的问题，现暂取共振峰值处的值，并在下文详细讨论．

数据处理结果见表1：(这里的峰值法即为教材中利用峰值的作图计算方法)

表1 各种数据处理方法的比较

处理

方法的值与阻尼法所得阻尼

系数的相对误差阻尼1峰值法0.0549841.81%=0.0560幅频曲线拟合0.0553281.20%相频曲线拟合0.0563530.630%阻尼2峰值法0.0612510.242%=0.0614幅频曲线拟合0.0615910.311%相频曲线拟合0.0631142.79%从表中误差栏可以看到，这3种处理方法所得到的结果误差都不大，而且在这两次不同的实验中不能判断那种方法更精确．

峰值法的详细分析

峰值法核心公式

为什么用拟合方法遇到了取值问题，而峰值法却可以通过考虑随的变化关系修正系统误差？要弄清这个问题需要详细分析峰值法的原理．

首先考察计算公式：

(3)

(4)

(5)

可见，满足(5)的基础是，即需要在共振峰附近．

实验数据处理后画出的幅频特性曲线以为唯一变量，那么和必须被当成定值系数，也即这种处理方法默认是一个定值．实际上，这种处理方法只用到了很小的一段(图5)，取值范围很小，确实可以做这样的近似．不过，在(1)式中同时将看成定量和可变修正量依然是一种矛盾的做法．

不确定度的分析

我们不妨再从这种修正对不确定度的影响来考察．设的不确定度为，根据间接测量量不确定度公式(这里只考虑的影响)：

(6)

(1)’

按(１)＇直接运算得

(7)

将视为与(１)＇中地位相同的参数，并设，只对其它运算得

(8)

为方便比较(7)(8)两式大小，先求出各无关参数的取值范围(阻尼1实验)得：，0.0560，当，［0.6321,0.6487］，［0.9867,1.013］．代入计算得：［0,377］，［