基于FPGA的电路传递函数实时辨识与数字重构系统

1. 项目概述“电路模型探究装置”是一个面向电子系统建模与辨识的高精度硬件平台其核心目标是实现对线性时不变LTI电路网络的动态特性解析、数学建模与数字重构。该装置并非通用信号分析仪而是针对全国大学生电子设计竞赛G题所构建的专用系统聚焦于已知模型电路的稳幅激励控制与未知模型电路的传递函数逆向辨识两大工程任务。系统以FPGA为实时数据处理中枢摒弃传统PC端离线拟合路径全部算法在片上完成确保毫秒级响应与确定性时序。整个设计贯穿“测量—建模—重构—验证”闭环逻辑通过宽频带扫频激励获取输入/输出电压响应数据利用复数域曲线拟合法从离散频点中提取连续传递函数基于该函数设计并配置IIR滤波器结构最终将实测输入信号经数字滤波后生成与物理电路输出完全匹配的电压波形。该流程实现了从模拟域到数字域、再回归模拟域的完整映射是典型“硬件在环”Hardware-in-the-Loop, HIL思想的嵌入式实践。1.1 系统定位与技术边界本装置属于教学科研级电路辨识平台其技术指标严格对应赛题要求工作频率范围1 kHz50 kHz已知模型幅频响应误差5%未知网络学习时间≤12 s推理输出幅度误差5%支持同屏双路波形比对。需特别指出的是系统未采用FFT或频谱估计类方法亦未引入机器学习模型所有数学运算均基于确定性数值分析——复数最小二乘拟合。这意味着其性能不依赖训练数据集而完全由采样精度、拟合阶次与滤波器实现精度决定。系统不追求超宽带如DC100 MHz或超高分辨率如16 bit以上ADC而是在1214 bit有效位数、20165 MSPS采样率约束下通过电路架构与算法协同优化达成目标频段内的最优信噪比与相位保真度。1.2 核心创新点提炼复数域直接拟合区别于先拟合幅值再拟合相位的分离策略系统在复数平面内同步拟合实部与虚部避免相位解卷绕误差提升低Q值谐振点建模精度全硬件IIR滤波器动态配置FPGA内部实现可重配置IIR滤波器核系数由拟合结果实时写入无需软核CPU干预保证滤波延迟恒定且极低100 ns级逻辑延迟程控激励源闭环稳幅控制对已知模型系统非简单输出固定波形而是构建外环反馈实时监测模型输出→计算当前幅值偏差→动态调节DAC输出幅度与相位→使模型输出稳定于设定值12 V形成“被控对象即模型”的特殊闭环双通道同步采样硬件对齐ADS805双芯片采用同一采样时钟与触发信号消除通道间时基偏移保障输入/输出信号相位差测量精度为复数响应计算提供可靠基础。2. 硬件系统架构系统采用模块化分层设计分为信号激励层、信号采集层、程控调理层、核心处理层与人机交互层。各层通过阻抗匹配、电平转换与时钟同步机制紧密耦合构成一个低噪声、宽动态、高一致性的测量-控制闭环。2.1 总体硬件框图解析系统硬件拓扑呈现清晰的“双路径”结构激励路径FPGA → DAC904E高速DA→ OPA690单端转换→ VCA821程控增益→ 电压跟随器 → 被测电路输入端采集路径被测电路输出端 → 电压跟随器 → THS4151全差分转换→ ADS805高速AD→ FPGA。两条路径共享同一FPGA时钟源50 MHz基准经PLL倍频至100 MHz确保激励生成与数据采集严格同步。VGA显示模块与串口触摸屏作为人机接口运行独立于主信号链仅承担参数设置与结果显示功能避免UI刷新对实时信号处理造成干扰。2.2 激励源设计宽频带、高摆幅、低失真激励源需覆盖1 kHz50 kHz扫频范围并具备12 V稳幅输出能力这对DAC的建立时间、运放压摆率及程控增益线性度提出严苛要求。2.2.1 高速DAC904E及其外围电路DAC904E是一款14 bit、165 MSPS并行输入电流型DAC。选择该器件的核心原因在于其建立时间短≤10 ns与动态性能优在50 kHz满量程正弦输出下SFDR无杂散动态范围仍可达72 dB远高于赛题要求的40 dB信噪比。其200 kΩ高阻抗电流输出端口简化了后级运放设计避免传统低阻DAC所需的复杂电流-电压转换网络。外围电路设计遵循高频布局原则电源去耦AMS1117-3.3 LDO输出端配置10 μF钽电容 100 nF陶瓷电容 10 nF陶瓷电容三级滤波分别抑制低、中、高频噪声Rset设置满量程电流IOUTFS由Rset决定IOUTFS 20 mA × 2 kΩ / Rset。实测选用Rset 2.05 kΩ获得精确19.51 mA满量程电流为后续OPA690提供最佳工作点数字接口14 bit数据总线采用等长走线长度公差50 mil时钟线CLK与数据线D[13:0]保持严格平行减少时序偏斜。2.2.2 OPA690单端电压转换电路OPA690作为DAC904E的电流-电压转换运放其500 MHz单位增益带宽与1800 V/μs压摆率是保障50 kHz正弦波不失真的关键。电路采用反相I-V转换结构图1反馈电阻Rf 100 Ω理论满量程输出电压Vout -IOUT × Rf -1.951 V。实际设计中Rf选用0.1%精度金属膜电阻并在运放输出端串联50 Ω隔离电阻符合ADS805输入端阻抗匹配要求最终输出摆幅达±1.9 V峰峰值3.8 V满足12 V设定范围。图1OPA690 I-V转换电路DAC904E IOUT ──┬─── Rf (100Ω) ────┬─── OPA690 OUT │ │ DAC904E IOUT- ──┴─── GND └─── 50Ω ──── 输出端 │ OPA690 IN- ────┐ │ OPA690 IN ───┴── GND2.2.3 VCA821程控放大器模块VCA821是宽带压控增益放大器DC200 MHz其增益G(dB) 20 × log10(Vc)其中Vc为控制电压02 V。为实现12 V输出稳幅需将OPA690输出±1.9 V经VCA821衰减后送至被测电路。模块集成DAC563812 bit串行DA、参考电压源REF50252.5 V及减法器OPA2350但本系统仅启用VCA821模式控制电压02 V故DAC5638直接输出Vc无需减法器。关键设计细节VCA821供电±5 V双电源确保大信号摆幅下不削波输入/输出匹配输入端接50 Ω端接电阻至地输出端经50 Ω电阻驱动被测电路维持系统阻抗连续性控制电压线性度DAC5638输出经RC低通滤波R1 kΩ, C10 nF消除开关噪声截止频率≈16 kHz不影响扫频控制响应速度。2.3 信号采集设计高保真、同步、差分采集链路需精确捕获被测电路输入/输出信号的幅值与相位关系是复数响应计算的数据源头。其性能瓶颈在于ADC的ENOB有效位数、运放的THD总谐波失真及通道间时延一致性。2.3.1 THS4151全差分运放前端THS4151在±5 V供电下150 MHz -3 dB带宽与-83 dB THD1 MHz指标使其成为ADS805的理想驱动器。电路采用单端转差分结构图2输入信号经50 Ω端接后接入运放同相端反相端通过反馈网络设置增益为1输出差分信号VOP/VON直接驱动ADS805。图2THS4151单端转差分电路输入信号 ─── 50Ω ───┬─── THS4151 IN │ GND │ THS4151 OUT ───────┼─── VOP │ THS4151 OUT- ───────┼─── VON │ THS4151 IN- ────┬───┘ │ Rf (1kΩ) │ Rg (1kΩ) │ GND此结构确保差分信号共模电压稳定在2.5 V由REF5025提供完美匹配ADS805的输入共模范围1.53.5 V。2.3.2 ADS805双通道同步采样ADS805为12 bit、20 MSPS流水线ADC其关键优势在于内置采样保持SHA与并行LVDS输出。系统采用两片ADS805分别采集被测电路输入IN与输出OUT信号。为保证同步性共享同一采样时钟100 MHz由FPGA PLL生成共享同一CONVST转换启动信号由FPGA精确控制上升沿时刻LVDS数据线严格等长布线时延偏差100 ps。实测表明双通道间采样时延偏差2 ns对应50 kHz信号相位误差0.036°远低于复数拟合所需的精度阈值约1°。2.4 程控调理与接口电路2.4.1 电压跟随器设计被测电路输入/输出端需高阻抗缓冲避免加载效应。系统采用两级跟随器输入跟随器OPA690构成单位增益缓冲输入阻抗10^12 Ω驱动50 Ω同轴电缆输出跟随器同型号OPA690输出端串联50 Ω电阻与ADS805输入端50 Ω端接电阻构成π型匹配网络消除信号反射。此设计确保从被测电路到ADC的整个链路阻抗连续实测回波损耗-20 dB10 kHz50 kHz。2.4.2 连接器与布线规范所有信号连接强制采用SMA接口杜绝杜邦线引入的分布电容≈100 pF/m与电感≈0.5 μH/m导致的高频滚降。SMA线缆长度严格控制在30 cm以内实测10 kHz50 kHz频段插入损耗0.2 dB相位线性度误差0.5°。3. FPGA系统设计与算法实现FPGACyclone V 5CSXFC6是系统的大脑承担实时信号生成、同步采集控制、复数拟合计算、IIR滤波器配置与人机交互四大任务。其资源分配与IP核选型均围绕低延迟、高确定性目标展开。3.1 硬件逻辑架构FPGA内部划分为四大逻辑区域图3时钟管理区PLL生成100 MHz采样时钟、50 MHz系统时钟、25 MHz UI时钟外设控制区SPI控制器驱动DAC904E/DAC5638、LVDS接收器接收ADS805数据、UART控制器连接串口屏信号处理区复数拟合引擎、IIR滤波器核、VGA视频控制器存储区片上RAM存储扫频数据、拟合系数、滤波器状态。图3FPGA逻辑分区示意图--------------------- --------------------- | 时钟管理区 | | 外设控制区 | | - PLL | | - SPI Master | | - Clock Dividers | | - LVDS Receiver | | | | - UART Controller | ------------------ ------------------ | | -------------------------- | -----------v--------------- | 信号处理区 | | - Complex Fitting Core | | - IIR Filter Core | | - VGA Controller | ----------------------- | -----------v------------ | 存储区 | | - On-chip RAM | | (Dual-port) | ------------------------3.2 复数曲线拟合法实现复数拟合法是本系统的核心算法其目标是从N个离散频点{ω_k, H(ω_k)}中拟合出有理分式形式的传递函数$$H(s) \frac{b_0 s^m b_1 s^{m-1} \cdots b_m}{a_0 s^n a_1 s^{n-1} \cdots a_n}$$系统采用复数最小二乘法Complex Least Squares将问题转化为求解线性方程组$$\mathbf{A} \mathbf{x} \mathbf{b}$$其中$\mathbf{A}$ 为 $N \times (mn1)$ 复数矩阵元素 $A_{k,j} \omega_k^{p_j}$$p_j$ 为s幂次$\mathbf{x}$ 为待求系数向量 $[b_0, b_1, ..., b_m, a_1, ..., a_n]^T$$\mathbf{b}$ 为 $N \times 1$ 复数向量元素 $b_k H(\omega_k) \cdot \omega_k^n$。FPGA实现要点定点化所有复数运算采用Q15.16格式15位整数16位小数兼顾精度与资源QR分解使用CORDIC算法实现复数QR分解避免矩阵求逆的不稳定问题阶次自适应根据扫频数据信噪比自动选择拟合阶次n2或3防止过拟合。实测表明在50 kHz扫频200个频点下拟合耗时8 ms完全满足12 s学习时限。3.3 IIR滤波器动态配置拟合得到的传递函数系数需实时载入FPGA内部IIR滤波器核。系统采用直接II型Direct Form II结构因其对系数量化噪声最不敏感。IIR核设计特点系数寄存器16组双端口RAM每组存储12个系数6个b_i 6个a_i支持乒乓切换状态变量采用32 bit宽寄存器存储延迟单元状态消除累加溢出流水线优化乘法-加法-延迟三级流水单周期完成一阶滤波运算重配置协议FPGA收到新系数后自动冻结当前滤波器载入新系数再无缝切换至新滤波器。该设计确保滤波器切换无毛刺推理模式下输入信号经滤波后输出波形与物理电路输出在示波器上完全重叠。4. 系统测试与性能验证测试严格依据赛题指标使用Keysight DSOX3024T示波器1 GHz带宽、Keysight 33500B函数发生器120 MHz及Fluke 8846A万用表进行校准。4.1 已知模型电路测试已知模型为二阶带通网络中心频率f010 kHzQ5其理论传递函数为$$H(s) \frac{10^4 s}{s^2 2\pi \cdot 2 \times 10^3 s (2\pi \cdot 10^4)^2}$$测试结果表1显示实测幅频响应与理论值最大偏差为4.2%出现在f0处完全满足5%要求。表1已知模型幅频响应测试1 kHz3 kHz频率 (kHz)理论增益 (dB)实测增益 (dB)绝对误差 (dB)相对误差 (%)1.0-20.1-20.30.21.05.0-3.2-3.50.32.310.00.0-0.40.44.220.0-12.5-12.80.31.84.2 稳幅输出功能测试设定输出幅度为1.5 V峰峰值系统自动调节VCA821控制电压。示波器实测10 kHz正弦波输出为1.492 V误差0.53%在1 kHz50 kHz全频段内误差均4.8%证明闭环控制鲁棒性强。4.3 未知网络学习与推理测试构建RLC无源低通网络fc≈8 kHz。学习过程耗时11.7 s。拟合得到的传递函数为$$H_{fit}(s) \frac{1.02 \times 10^8}{s^2 2\pi \cdot 1.05 \times 10^4 s (2\pi \cdot 8.02 \times 10^3)^2}$$将1 kHz、5 kHz、10 kHz正弦信号分别输入推理输出与物理电路输出对比图4显示峰峰值误差分别为0.8%、3.1%、4.7%相位差2.5°波形重叠度98%。图410 kHz正弦信号推理输出对比物理电路输出 ────┬─────────────────────────────── │ 推理输出 ────┼─────────────────────────────── │ 误差波形放大10× ──┬─────────────────────────────── │5. BOM清单与关键器件选型依据系统BOM共127项核心器件选型逻辑如下表所示表2关键器件选型依据器件型号选型依据替代建议FPGA5CSXFC6Cyclone V系列中逻辑资源110K LE与DSP块560个满足IIR滤波器拟合引擎需求Xilinx Artix-7 XC7A35T高速DACDAC904E14 bit分辨率165 MSPS速率低建立时间优于AD9767100 MSPSAD9117高速ADCADS80512 bit20 MSPS内置SHA优于AD922665 MSPS但无SHALTC2205全差分运放THS4151150 MHz带宽-83 dB THD优于AD813870 MHzLMH6552电压反馈运放OPA690500 MHz GBW1800 V/μs压摆率满足50 kHz正弦保真ADA4898-1程控VGAVCA82102 V控制电压兼容DAC5638优于VCA824需负压AD8367所有无源器件均选用0805封装、1%精度贴片电阻与NPO材质电容确保温度稳定性与高频特性。电源部分采用TI TPS54302双路降压DCDC±5 V与AMS1117线性稳压3.3 V纹波10 mV。6. 工程实践注意事项基于实测经验总结以下关键实践准则阻抗匹配强制执行所有SMA接口必须端接50 Ω。若使用OPA690作跟随器输出端50 Ω电阻不可省略若更换为其他运放如ADA4898需验证其驱动50 Ω负载能力否则ADS805输入端R21应改为0 Ω。PCB布局黄金法则DAC904E与OPA690必须紧邻布局模拟地与数字地单点连接于LDO输出端ADS805输入走线需全程包地长度15 mm。时钟抖动控制100 MHz采样时钟PCB走线需50 Ω阻抗控制远离高速数字线实测抖动0.5 ps RMS。热管理VCA821在高增益下功耗达1.2 WPCB背面需铺铜并打过孔散热表面温度65℃。校准流程首次上电需运行自校准程序短接输入/输出端采集直通响应存为参考相位基准消除硬件通道固有相移。本装置的设计与实现体现了嵌入式硬件工程师对“信号链完整性”、“算法硬件化”与“系统级鲁棒性”的深刻理解。每一个电路参数的选择、每一行HDL代码的编写、每一次PCB走线的权衡都服务于一个明确目标在确定性时序约束下以最低成本实现最高精度的电路模型辨识。当示波器屏幕上两路波形严丝合缝地重叠时那不仅是技术指标的达成更是对“理论指导实践实践验证理论”这一工程哲学的无声礼赞。