刘训臣燃烧诊断实验室主页
  • Alumni
  • People
  • Post
  • Research

Post

科研笔记和一些想法

十一月 8, 2025

HITRAN 数据库用法

在上一篇 吸收光谱简介 ,我们介绍了分子吸收光谱的基本知识。 常见的分子光谱的数据库有:

  • HITRAN, USA Harvard-Smithsonian Center for Astrophysics Molecular Data
  • NASA JPL Molecular Spectroscopy, USA
  • The Cologne Database for Molecular Spectroscopy, CDMS, Germany
  • GEISA Spectroscopic Database, France
  • Pacific Northwest National Labs Spectral Database, USA
  • Ames Research Center Infrared Spectral Database of PAH, USA
  • MPI-Mainz-UV-VIS Spectral Atlas of Gaseous Molecules, Germany
  • surface types: US Geological Society Digital Spectral Laboratory, USA
  • NIST Chemistry WebBook, USA
  • AIST Integrated Spectral Data Base System for Organic Compounds, Japan
  • Bremen, Institute of Environmental Physics, Germany
  • LISA database (SCOOP), France
  • Institute of Atmospheric Optics (CO2 O3), Tomsk, Russia
  • Instituto Superior Tecnico SPATRAN, Gas and Plasma Radiation Database, Italy

HITRAN数据库是我们最常用的分子光谱谱线的数据库。NASA JPL 数据库中有很多星际分子的谱线,另外JPL数据库用的 SPCAT/SPFIT 软件也是功能非常强大。未来我们会尝试用 SPCAT 模拟高温反应流体的分子光谱。Tomsk 有很好的CO2的光谱数据库。

继续阅读
十月 3, 2025

吸收光谱简介

简介

在燃烧诊断中,吸收光谱的测量方法从光源上分非激光和激光两大类。 传统的非激光的吸收光谱测量,最早从 Hertzberg 那时侯开始,多采用宽频的非相干光源如闪光灯和红外热幅射光源,结合光栅、滤片或者干涉仪进行分光,通过测量不同波长电磁波的透射光强被火焰吸收减弱得到吸收光谱定量信息,就可以获得燃烧体系中分子浓度、温度、压力和速度的信息,能得到很多对火焰结构的测量结果和深刻理解。 比如Gaydon的书《Flames,Their Structure, Radiation, and Temperature》

激光在 20 世纪下半叶的发明和应用,极大地促进了吸收光谱的发展。激光吸收光谱用高亮度、窄线宽的可调激光光源,当扫描激光波长时,光电探测器测量不同波长激光在通过燃烧体系后的光强损失,从得到所扫描的激光波长范围内的高分辨率的激光吸收光谱。 现代的激光吸收光谱多指采用近红外和中红外半导体激光光源,通过调制方法扫描激光的波长,光电探测器直接测量不同波长激光在通过燃烧体系后的光强损失,从得到所扫描的激光波长范围内的可调谐半导体吸收光谱(Tunable Diode Laser Absorption Spectroscopy, TDLAS)。

吸收光谱研究的是光与物质作用过程中最基本的吸收过程,是最基本的燃烧诊断方法,有易于操作、选择性强、定量准确、灵敏度高的优点。 吸收光谱所测量的光谱有位置、强度、和线宽三个基本信息,可以定量地研究反应流中分子浓度、温度、压力和流速等各种信息。 谱线的位置由所测量分子的能级结构决定,与燃烧体系的温度压力等热力学状态与流场结构无关。(在高压燃烧环境中,分子碰撞导致的能级降低和多普勒效应还会引起的微小的吸收峰位置的移动。) 因此可以通过测量燃烧体系的吸收谱准确地区分所测量的不同分子,具有极高的选择性(selectivity)。 在谱线强度测量方面,结合光腔共振等方法增加吸收路径的长度,吸收光谱可以测量ppt量级的分子浓度,具有很高的灵敏度(sensitivity)。 研究吸收光谱所测量的谱线吸收强度只需考虑体系状态从低能级到高能级的光子吸收这一简单过程,可以忽略上能级过程,实验和理论计算均可较准确的定量研究。因此吸收光谱的最大特点为定量准确(absolute)。

本文主要介绍吸收光谱中常用的一些定量关系。包括 Einstein A B 系数的介绍和HITRAN 数据库的使用。 可以参考Peter Bernath 的课本《Spectra of Atoms and Molecules》和这篇 HITRAN 官方的介绍文献[1]。

Einstein’s A and B constant

这里从量子力学能级的角度对光子吸收过程稍做描述。 分子对光子的吸收模型是上下两个能级之间的跃迁,因此将分子看作包含$|\phi_0\rangle|\phi_1\rangle$两个态的两能级系统(two-level system),将光用经典的电磁波理论看作一个正弦波描述的波动电场。由于吸收光谱中所用激光的瞬时电场强度远远小于分子内原子和电子之间的作用强度,可以将激光的作用做含时微扰处理。当电磁波达到共振频率时,有较大的跃迁速率。(电磁波的频率以 [$\mathrm{Hz}$]为单位。)可以得到下能级到上能级的吸收(absorption)与上能级到下能级的受激辐射(stimulated emission)的跃迁速率$B$ [$\mathrm{s^{-1} /(J \cdot m^{-3}\cdot Hz^{-1}}$][2]:

$$ \begin{aligned} \frac{d N_1}{d t} &= -\frac{d N_0}{d t} = B_{1\leftarrow0} \cdot \rho(\nu) \cdot N_0 \\ -\frac{d N_1}{d t} &= \frac{d N_0}{d t} = B_{1\rightarrow0} \cdot \rho(\nu) \cdot N_1 \\ B &= \frac{2 \pi^2}{3\epsilon_0h^2} |\vec{p}|^2 \end{aligned} $$

其中, $\rho_\nu$ 为普朗克定理中光场辐射强度,单位[ $\mathrm{J \cdot m^{-3} \cdot Hz^{-1}}$]。

继续阅读
九月 16, 2025

用Witoszynski曲线设计收缩段

今天阅读了Witoszynski提出witoszynski曲线的文章,witoszynski曲线最初的设计是用于优化扩压器设计,目的是减小气体或液体的流速而没有显著能量损失。

Witoszynski工作的核心成果是描述扩压器内流线径向坐标r的表达式: $r = \frac{r_0}{\sqrt{1 - \left(1 - \frac{r_0^2}{r_1^2}\right) \frac{\left(1 - 3\frac{z^2}{a^2}\right)^2}{\left(1 + \frac{z^2}{a^2}\right)^3}}}$

其中,

  • r:径向坐标
  • r₀:入口半径,
  • r₁:出口半径,
  • z:表示沿流向的轴向坐标,范围是0到$\frac{a}{\sqrt{3}}$
  • a:几何常数参数, 流线示意图如下:

根据文章,流线图v对应的为公式中的r,α对应公式中的a。 流线径向坐标r表达式的推导始于轴对称理想流体在柱坐标下的微分方程: $\frac{\partial^2 \psi}{\partial r^2} - \frac{1}{r} \frac{\partial \psi}{\partial r} + \frac{\partial^2 \psi}{\partial z^2} = 0$

其解为流函数的级数形式: $\psi = r^2 f(z) - \frac{r^4}{2 \cdot 4} \frac{d^2 f(z)}{dz^2} + \frac{r^6}{2 \cdot 4 ^2\cdot 6} \frac{d^4 f(z)}{dz^4} + \cdots$

$f(z)$是一个光滑的函数,满足$f_{\infty}=1$.

则假设一个符合条件的$f(z)$函数,将其带入流函数得

A由边界条件而定 把A带入 得出r的表达式

$$ R = \frac{R_0}{\sqrt{1 - \left[1 - \left( \frac{R_0}{R_1} \right)^2 \right] \frac{\left( 1 - \frac{z^2}{z_0^2} \right)^2}{\left( 1 + \frac{z^2}{3z_0^2} \right)^3}}} $$

与原式子的不同点在于,公式轴向坐标z的范围简化,变成0到$z_0$,于是公式随之变化,但其实两式是相同的。 那么如何用wito曲线设计收缩段呢

继续阅读
八月 21, 2025

解决流量计跳数问题总结和感悟

在进行第二次轴向分级燃烧第一级点火燃烧实验时,主流空气更换了两个大量程流量计2000slpm,在实验过程中,返回数据出现了位数错乱现象,如下图。 流量的读取是通过截取固定位置及位数的返回码而后经过换算实现的,因此流量显示出现了频繁的掉数。刚开始考虑可能是硬件问题,例如信号线接线不牢固,拓展坞插口异常,usb转485设备传输不稳定,流量计信号线之间存在电磁干扰等问题。于是急忙重新接了一次线,更换了另一个拓展坞,发现问题依然存在,但是好在,在实验中发现了,只有流量计流量在二百多到三百的时候才会出现掉数,大流量则不会出现这种情况,这起码证明掉数的原因并非是设备密集而导致的干扰问题,也非接线的问题。至此,除usb转485接口外,其他硬件因素基本已经排查完毕。

当时因为只有两个大量程流量计出现了跳数,考虑可能是用大流量计测较小流量会出现错误,询问流量计客服后,用串口助手测试,发现收发功能码都是正常的,但是用labview则会出现返回码不完整和冗余的现象,考虑可能是labview的编写问题。后来,购买的高质量usb转485设备到了,迫不及待的换上试了一下,确实质量好,也确实好用了很多,但是问题并没有解决,于是便开始一心解决软件及程序问题。

经过多途径的大量搜索,怀疑可能是大量程流量计在收发小流量时,写入速度大于读取速度,累计数据量增大导致缓冲区溢出,进而导致跳数。于是基于这一猜想,想出了几种解决方案:

  1. 在读取前,增加一个清除读取缓冲区,保证每次的读数不会被前面影响
  2. 在写入和读取之间增加一个延时,避免读取过快过慢导致的错位
  3. 使用条件结构,在只有返回数为7时,才进行数据的读取和写入
  4. 使用条件结构,判断起始位是否是01,如果是,再进行读数和存储

实践后发现无论哪种方法,都无法避免其在特定的流量范围内出现返回码不完整的问题,总是会频繁的出现0103020A这个不完整的返回码。还有一个比较麻烦的方法就是在visa读取后,增加一个队列,数据存入队列,长度等于7则处理,不足7则继续等待,这个学习了一会,还没学会,后续可以再学习一下。

在验证这些方案是否可行的过程中,突然发现,每次返回数的位数出现偏差,都是以0A为结尾,于是搜索了一下发现,labview的默认终止符是0A!在前面的blog中提到,读取流量的位数是第4,5位,如图。

在十六进制中,0A的取值范围为0A00-0AFF,对应十进制的值为2560-2815,再除以流量计的放大倍数为256-281.5,即在这个流量范围内都会因触发终止符而中断,导致后续的错位,致使无法正常显示流量。解决办法即是在visa配置串口将启用终止符创建常量,并将其关闭,如图。

完成修改后,再次测试程序,返回码中断,错位问题得到解决,如图。

经过此次经历让我发现,labview的学习任重而道远。一步就能解决的问题,竟然走了这么多弯路。

继续阅读
八月 20, 2025

Blow-out limits of nonpremixed turbulent jet flames in a cross flow at atmospheric and sub-atmospheric pressures

文章重点研究了非预混湍流射流火焰的吹熄极限。

实验采用的设备如图: 风洞的横截面高度1.5m,宽度1.2m,长度10m,高压机械风扇产生横向气流,蜂窝结构用以稳定气流,另一端是开放式,燃料喷嘴的内径范围 3 - 8mm,壁厚 1.5mm,喷嘴放在风洞出口0.5m 处,距离底部 0.3m,燃料是丙烷,流量由气体流量计调节。燃料温度与环境温度相同。

在合肥和拉萨分别进行了实验,火焰用 CCD 相机记录。将实验结果进行了比较,以评估压力对横流火焰吹熄极限的影响,所有实验重复了三次。

实验结果如下: a,b 图分别为合肥和拉萨实验火焰形状变化图 (a) 图工况:射流速度39.3m/s,喷管内径3mm,环境压力 100KPa 保持喷嘴形状和射流速度不变,随横流速度的增加,火焰有三种状态, 0 - 2m/s,火焰高度长度减少,弯曲程度变得更明显,蓝色火焰区扩大。 这归因于横流增大,增强了喷嘴和抬升火焰边缘区域空气的卷吸和混合,使燃料核心区域局部预混更强,整体火焰亮度减低。 2 - 3m/s,火焰整体形态变化不大,火焰主体与横流几乎平行,火焰前缘近乎垂直,火焰变得更蓝,亮度显著降低,在水平方向略有拉长。 结果同 Moore 等人在共流条件下观察到的过程相似,当横流超过一个阈值,火焰从稳定火焰过渡到熄灭和消失。火焰前缘向下游移动,尺寸变小;扩散火焰尺寸显著减少;整个火焰以火焰前缘的消失而结束,此为吹熄。 临界吹熄横流速度定义为火焰不能持续存在,发生熄灭,所以吹熄极限是 3.37m/s (b) 图工况:射流速度 16.9m/s,喷管内径 3mm,环境压力 64KPa 火焰出现与 a 相似的三种状态 (c) 图是火焰边缘特写,工况为:射流速度 19.7m/s,环境压力 100KPa 由于迎风侧与背风侧的混合程度不同,导致火焰抬升高度不同。 (由射流与横流的相互作用,横流中射流火焰的流场比静止空气中自由射流火焰更复杂。迎风侧的火焰边缘由于横流,会遭受更高的应变效应,横向正应变相较于背风测更大。)

整理实验数据绘制了横流速度与燃料射流速度的吹熄极限图 由图可知 100KPa,吹灭极限处的气流速度随燃料喷射速度增大而增大(位于吹熄曲线的横流主导区域) 64KPa,吹灭极限处的气流速度随燃料喷射速度增大而增大,然后又减少。呈非单调行为,与 kalghatgi 观察的一样

将速度与喷嘴直径的比值与 kalghatgi 的数据(100KPa)作比较,如下图: 结果显示,在标准大气压下两者数据相吻合。 由此图可知,在标准大气压和亚大气压两种情况下,数据都显示出非线性行为。 有两个区域: 1.横流主导:在 “横流主导” 区域,熄灭主要由横流产生的应变引起。在这种情况下,首先火焰被吹向平行于横流的方向,然后,火焰沿水平方向被横流吹灭。 在点 A,燃料射流速度接近静止空气中自熄灭极限的值(对应点 C),自 A 点以后向燃料射流主导转变 2.燃料射流主导:熄灭主要由高速燃料射流本身产生的应变引起。 B 点是可持续火焰的最大值,射流速度超过点 B 时,即使在横流中,火焰也无法维持。若要维持火焰,射流速度必须降低 C 点,吹熄完全由燃料射流本身产生的应变效应引起的。 吹熄极限曲线包围的是火焰存在区域。 环境压力降低到 64kpa,吹熄极限曲线明显缩小,这说明随着环境压力的降低,火焰对于燃料射流速度和横流速度的吹熄更敏感,更容易熄灭。吹熄时(横流速度 / 喷管内径)局部最大值约为 400s⁻¹ 约是标准大气压的三分之一,射流速度 / 直径的值也有所下降。

继续阅读
八月 15, 2025

利用JANAF表估算氢氧火焰温度

一、扫描求解:氢氧反应平衡与温度搜索

$$ \mathrm{H_2 + \tfrac{1}{2}O_2 \rightleftharpoons H_2O} $$$$ n_{\mathrm{H_2O}} = \xi,\quad n_{\mathrm{H_2}} = 2 - \xi,\quad n_{\mathrm{O_2}} = 1 - \frac{\xi}{2} $$$$ K_p(T) = \frac{n_{\mathrm{H_2O}}}{n_{\mathrm{H_2}} \cdot (n_{\mathrm{O_2}})^{1/2}} \left(\frac{P}{P_0\, n_{\text{total}}}\right)^{-1/2} = \frac{\xi}{(2-\xi)\sqrt{1-\xi/2}} \cdot \sqrt{3 - \frac{\xi}{2}} $$

(式中压力比 $P/P_0$ 取 1)。对每一个猜测温度 $T$,由 JANAF 数据计算 $\Delta_r G^\circ(T) = \Delta_f G^\circ_{\mathrm{H_2O}}(T)$,进而得 $K_p(T) = \exp[-\Delta_r G^\circ(T)/(RT)]$。求解上述方程即可得到该温度下的平衡反应进度 $\xi_{\mathrm{eq}}(T)$,同时也就确定了平衡组成。

$$ H_{\mathrm{prod}}(T) = \xi_{\mathrm{eq}} \cdot H_{\mathrm{H_2O}}(T) + (2-\xi_{\mathrm{eq}}) \cdot H_{\mathrm{H_2}}(T) + \left(1-\frac{\xi_{\mathrm{eq}}}{2}\right) \cdot H_{\mathrm{O_2}}(T) $$

其中 $H_i(T) = \Delta H_i(T) + \Delta_f H_i^\circ$,单质的 $\Delta_f H^\circ$ 为零。反应物的初始总焓 $H_{\mathrm{react}} = 0$(因为反应物就是 298.15 K 的单质)。因此,需要寻找使 $H_{\mathrm{prod}}(T) = 0$ 的温度。

继续阅读
六月 30, 2025

激光车

一般纳秒固体激光就是指Nd:YAG激光。 频率多为10Hz,工作模式是1个脉冲pump,1个脉冲发光。

就像高速相机拍摄一串图像一样,串脉冲激光在1个长时间的宏脉冲中可以发出多个脉冲,可以实现高能高频的脉冲,特别适合湍流激光成像,如PIV PLIF等。非常适合用于燃烧台、风洞等 test rig 。

我的看法是应当发展这种先进的测试设备,用于现有的实验台。 于是我们研发了激光器,并把激光器需要的洁净室功能集成到一个自走的底盘上做成了一个激光车,可以很方便地用集装箱运输。到了实验场地插上一个电缆就可以工作。非常方便。 激光就像是灯泡一样。可以把实验台照亮。

1 2

而不能反过来,不能因为灯泡去盖一个房子。 但这样往往比较容易申请到“大项目”,可以下一步争取到“大帽子”,名利双收。 这种科研上的盲目扩张和恒大爆雷没有本质区别。 趋之若鹜,其实是浪费资源,也浪费了自己宝贵的科研时间。

继续阅读
六月 29, 2025

低空经济会议

2025年6月28日我们学院组织了一个低空经济为主题的会。并且有一个低空经济现代产业学院的揭牌仪式。欢迎大家来参会!感谢来做报告的专家!

2025低空技术与装备研讨会在青岛召开

机电学院举办低空技术与装备专业建设论坛

1 1 1 1 1 1 1 1


未来我们也会做一些无人机方面的工作。最近要做的就是有一个风墙吹风的无人机飞行测试笼子。 侧面采用很多个风扇吹风产生风速和分布可变可控的强风环境,模拟无人机在强风环境下飞行控制实验。

初步可以通过少量几个相机进行初步粗略的轨迹追踪,后续可以增加风墙面积、无刷电机数量、相机数量,逐步提高实验能力和实验精度。

和图像有关的主要做两个事情:

  1. 无人机周围多个摄像头, 通过图像重构方法实现无人机空间定位和三维姿态动作捕捉。重构无人机轨迹姿态,在实验室可以作为复杂控制算法的测试平台,提供精确的位置和姿态信息。

  2. 无人机上一个或者两个摄像头,拍外部画面,计算自己的规矩姿态。这个和无人驾驶就很像了。现在无人驾驶用的很多的是端到端的算法,通过这个项目可以学学纯视觉的算法。有大规模应用的前景。当然无人机的话要有地图去匹配。

继续阅读
六月 23, 2025

横向射流轨迹的拟合

横向射流(crossflow)是很常见的一种流体混合结构,当一股小尺寸的射流垂直进入大尺寸的横流时,垂直的射流就会逐渐转到横向过来,形成横向射流。 在这个过程中形成的涡系主要包括剪切涡、马蹄涡、尾迹涡和对旋涡,其中剪切涡和对旋涡的生成发展、传播和最终破碎是横向射流混合增强的主要因素。通过控制射流和主流的动量比,可以控制轨迹,改变穿透深度和形成的火焰长度。

远场的横向射流的轨迹一般可以用幂函数表示[1,2]。如图1,横向射流轨迹的可以由速度比$r$和射流直径$d$进行无量纲化。

$$ \left(\frac{y}{rd}\right) = A \left(\frac{x}{rd}\right)^{B} $$

典型的轨迹方程中,$B\approx \frac{1}{3}$。

以我们实验测量的横向射流的轨迹为例。

图1 实验原始图片、拟合轨迹、和以轨迹为坐标转化后的图像


拟合得到轨迹方程包含如下几步:

  1. 在每个高度,也就是$y$坐标上,找到浓度或速度的最大值对应的$x$值。这里找到的最大点坐标还比较粗糙,最好再和相邻的几个点通过抛物线拟合的方法,得到更精确的坐标。 如图2 ,得到的最大点基本描绘出了轨迹,但还需要修剪一下数据。

图2 找最大点

  1. 去掉头尾部分。去掉轨迹后面几个点比较好理解,因为到轨迹最后信噪比很低了;去除开头几个点是为什么呢? 以实验测量的双喷嘴横向射流为例子,图3中给出了初步找到的中心线和修正后的中心线。

图3 修正后的中心线

  1. 下一步就是拟合了。首先要把找到的轨迹起始点归零。然后求log10 ,得到拟合直线,此时会发现需要去掉开头的非常接近零的几个点,也就是近场不符合幂函数规律的点,便于拟合。 从图4中可以看出,拟合的直线需要确定两个拟合变量,一个是拟合长度范围,一个是拟合起始位置。 不一定拟合长度越长越好,需要确定恰当的范围。 我们将这两个变量都在一定范围内变化,看得到的结果情况。

图4 变量拟合

图4中,1)为待拟合的直线,2)横坐标为拟合位置、纵坐标为拟合长度,当拟合长度大时,拟合位置变化的就小。3)为2)图的梯度,可以看到有趣的梯度变化的图案,说明了拟合参数的影响。 可见当拟合长度和位置恰当时,得到的参数有一个稳定的区域,数值基本不变。表现为梯度接近零,可以用这个指标来确定恰当的拟合长度和位置。

在这个例子中,得到直线的方程为:

$$ \ln(y) = 0.30\times \ln(x) + 1.36 $$

而 $r = 7$, $d = 10mm$,得到

$$ \left(\frac{y}{rd}\right) = 0.2 \left(\frac{x}{rd}\right)^{0.3} $$

参考文献

  1. BROADWELL J E, BREIDENTHAL R E. Structure and mixing of a transverse jet in incompressible flow[J]. Journal of Fluid Mechanics, 1984, 148: 405-412.
  2. HASSELBRINK E F, MUNGAL M G. Transverse jets and jet flames. Part 1. Scaling laws for strong transverse jets[J]. Journal of Fluid Mechanics, 2001, 443: 1-25.
  3. SMITH S H, MUNGAL M G. Mixing, structure and scaling of the jet in crossflow[J]. Journal of Fluid Mechanics, 1998, 357: 83-122.
  4. SU L K, MUNGAL M G. Simultaneous measurements of scalar and velocity field evolution in turbulent crossflowing jets[J]. Journal of Fluid Mechanics, 2004, 513: 1-45.
继续阅读
六月 22, 2025

利用HDM3065 Hantek 6位半万用表对5V电源输出电压结果的测试

在电子设备的使用中,电源输出电压的稳定性至关重要。今天,我将通过实际测量和数据分析,使用最小二乘法拟合它的输出电压值,结合得到的方差来评判它的输出电压效果。

测试方法

首先,介绍一下测试方法。我使用万用表作为测量工具,将电源输出电压从 0V 开始,以 0.1V 为步长逐步调节,直至达到 5V。我这里用到的是 Hantek 的 HDM3065 系列的万用表来测量小电源的输出电压,我们将万用表连接到测电压端,红黑表笔分别接到小电源的正负极,用来测量小电源的输出电压。连接到其他的地方以后也可以测试电流、电阻等。我们将红色线接到电源正极,黑色线连接到电源负极。

该万用表有专用的 IO 控制软件,将万用表连接到电脑以后,我们可以通过电脑的 IO 软件发送指令,万用表便会执行相应的操作。为了保证数据的可靠性,针对小电源的每个电压值,我都进行了五十次测量,并将测量的数据都记录下来。

测量连接示意图

数据处理

数据记录下来之后,我们需要用到 Python 来处理这些数据,运用最小二乘法对数据进行处理。最小二乘法的核心思想是通过最小化误差的平方和寻找数据的最佳函数匹配。在这里,我们假设电源输出电压与设定电压之间存在线性关系,通过最小二乘法拟合出一条直线,来描述这种关系。经过计算,可以得到拟合直线的方程。

通过计算方差来评估电源输出电压的稳定性。方差是用来衡量一组数据离散程度的统计量,方差越小,说明数据越集中,电源输出电压越稳定。分别计算每个设定电压点下十个测量数值的方差,发现大部分电压点的方差都较小,说明该电源在不同输出电压下的稳定性表现良好。

测试结果

测试结果如下图所示:

最小二乘法拟合图

拟合值相对于拟合直线理论值残差图

拟合值相对于设定值的残差图

我们可以看到该小电源的线性效果还是很好的,方差也很小,说明电压输出比较稳定。从残差图中我们可以看出,残差图上的点代表对应设定电压下,输出电压实际值与拟合模型预测值的差值,可以看出这些点的值都比较小,残差绝对值大多在 0.006V 以内,数值很小,表明实际值和拟合预测值偏差不大,模型对该电源电压数据的拟合精度较高。这张残差图显示模型拟合效果不错,误差随机且整体很小,能较好描述设定电压和输出电压的关系。并且我们可以在图中看出,拟合值相对于设定值的插值都在-0.01内,数据很小,也没有某些值与设定值差距较大,说明该小电源的输出效果还是很不错的。

另外,我还对每个电压值下测量的电压输出值查看其分布规律,根据测试结果生成的数据大致都有如下图的趋势。这里由于数据过多,我就不列出来了,拿出几个图来看一下分布规律即可:

0.9V电压分布

1.1V电压分布

2.9V电压分布

4.3V电压分布

图中每个蓝色柱子代表“输出电压落在某一小区间内的测量次数密度”,且直方图与曲线贴合度较高,故可以看出,电压的波动符合正态分布规律。

结论

通过最小二乘法得到了输出电压的估算值,我们大致了解到了该小电源输出电压的大致水平;结合方差的计算结果,我们对其稳定性也有了清晰的认知。从本次测量结果来看,该小电源的输出电压效果整体表现良好,在每个电压值下都能提供接近于小电源显示电压值的输出电压;但会存在一定的波动,所以在对于电压稳定性较高要求的条件下不适用。在今后的电子项目的实验中,掌握这种评估电源输出电压效果的方法,有助于我们更好地选择合适的电源,保证项目的顺利进行。

继续阅读
  • ««
  • «
  • 1
  • 2
  • 3
  • 4
  • »
  • »»
© 刘训臣燃烧诊断实验室主页 2026