项目设计论文微弱信号检测

1、 006A013微弱信号检测装置参赛学生：联系方式：联系方式：联系方式：指导教师：理工大学自动化专业摘要本设计是检测在强噪声背景下已知频率的微弱正弦波信号的幅度值，并数字显示出该幅度值。以锁相放大器作为微弱信号检测装置的核心，利用模块化设计，分别实现系统各部分功能。加法器模块，将输入信号淹没在强噪声里；衰减模块是将混合信号变成小信号；跟随器与放大模块，利用OP27实现了电压跟随并将小信号进行放大至原来的幅值；锁相放大器是一种基于AD633芯片与低通滤波的电路，利用AD633乘法器对同频同相参考信号与探测信号进行处理，再通过低通滤波电路，使输出为直流量；放大与显示模块把输出信号放大到适合STM3

2、2A/D采样的电压值并经过单片机处理后用数码管显示。通过硬件和软件的调试，结果与预期的一样。关键词：微弱信号 锁相放大 强噪声 STM32 AD633Abstract：This design is in the strong noise background of weak sinusoidal signal of known frequency amplitude value and digital display of amplitude value. The core of weak signal detection device is lock-in amplifier, using

3、of modular design, respectively realizing the function of each part in the system. Adder module submerges the input signal in the noise. Pressure module declines the mixed signal into a small signal. The follower and amplification module use OP27 to achieve a voltage follower andenlarge the small si

4、gnal to the original amplitude. The lock-in amplifier is based on the AD633 chip and the low pass filter circuit, using the AD633 multiplier to handle reference signal and the detection signal on the same frequency and phase, and then through the low pass filter circuit, the output becomes a straigh

5、t flow. Amplification and display module is to amplify the output to the STM32A/D for sampling the voltage value and processed by single chip microcomputer with digital delay. Through the hardware and software testing, the result is identical to the expected.Keywords: weaksignal； lock-in amplifier ；

6、 strong noise ； STM32 ； AD633目录一、 设计要求 1二、 总体方案设计 22.1总体方案概述 2 2.2总体方案分析 3 三、 各模块方案分析与论证 33.1噪音信号处理 33.2加法器设计 33.3纯电阻网络衰减设计 33.4微弱信号检测设计 43.4.1前置放大电路设计 43.4.2微弱信号检测电路设计 3.4.3低通滤波设计 73.5软件设计 7四、 实验分析与测试结果五、 总结六、 参考文献七、 附录一、设计要求1任务正弦波信号源微弱信号检测电路噪声源VS VN VC Vo 纯电阻压网络显示电路加法器Vi A B C D E 设计并制作一套微弱信号检测装置，

7、用以检测在强噪声背景下已知频率的微弱正弦波信号的幅度值，并数字显示出该幅度值。为便于测评比较，统一规定显示峰值。整个系统的示意图如图1所示。正弦波信号源可以由函数信号发生器来代替。噪声源采用给定的标准噪声（wav文件）来产生，通过PC机的音频播放器或MP3播放噪声文件，从音频输出端口获得噪声源，噪声幅度通过调节播放器的音量来进行控制。图中A、B、C、D和E分别为五个测试端点。图1 微弱信号检测装置示意图 2基本要求（1） 噪声源输出VN的均方根电压值固定为1V±0.1V；加法器的输出VC=VS+VN，带宽大于1MHz；纯电阻分压网络的衰减系数不低于100。 （2）微弱信号检测电路的输

8、入阻抗Ri1 MW。 （3）当输入正弦波信号VS 的频率为1 kHz、幅度峰峰值在200mV 2V围时，检测并显示正弦波信号的幅度值，要求误差不超过5%。 3发挥部分 （1）提高正弦波信号的识别能力，当输入正弦波信号VS 的频率在100Hz10kHz围、幅度峰峰值在20mV 200mV围时，检测并显示正弦波信号的幅度值，误差不超过5%。 （2）在发挥部分（1）的条件下，要求检测误差不超过2%。 （3）当输入正弦波信号VS 的频率在100Hz10kHz围时，进一步降低VS的幅度，检测并显示正弦波信号的幅度值，误差不超过2%。二、总体方案设计2.1总体方案概述 本设计以锁相放大器为微弱信号检测装置

9、的核心。输入信号与强噪声经过加法器叠加，经由衰减电路、电压跟随器电路、前置放大电路，然后和参考信号输入到乘法器电路中，之后通过低通滤波滤除交流分量，再经过放大电路，最后经STM32A/D采样、处理，输出检测信号幅值并由数码管显示。2.2总体结构设计三、各模块方案分析与论证3.1硬件设计 3.1噪音信号处理题目要求噪声源输出VN的均方根电压值固定为1V±0.1V，噪声源采用给定的标准噪声（wav文件）来产生，通过PC机的音频播放器或MP3播放噪声文件，从音频输出端口获得噪声源，噪声幅度通过调节播放器的音量来进行控制。而现有设备播放噪声信号均方根幅值不能达到要求，因此设置一前置放大电路对

10、噪声信号进行处理，电路设置放大倍数为十倍，通过调节音量大小来控制噪声均方根幅值大小。放大电路原理图见附图3.1.1。3.2加法器设计 按照题目要求需要设计加法器电路将噪声信号与标准正弦信号叠加实现此要求电路有以下两种方案：方案一：设计一个同向加法器电路，R1，R2，R3的阻值相等，R4的阻值应为R1/R2/R3的值，试验中选用R的阻值为10K，电路原理图见附图3.2（a）。但实际电路中对于微弱信号同向加法电路幅值波动较大，故而舍弃此方案方案二：设计一个反相加法器电路，由于改变了输入信号的相位故此再加以反相器电路，由此两部分组成加法器电路，电路原理图见3.2（b）。在实际调试时，反向加法器电路幅

11、值波动小，波形更为稳定，因此我们选择方案二。3.3纯电阻网络衰减设计题目要求采用纯电阻分压网络实现分压网络的衰减系数不低于100，通过查阅资料，现有以下几种方案，可以供我们选择。方案一: 由一个电位器和一个固定的电阻组成的分压网络。通过改变电位器的阻值就可以改变其衰减系数。这样就可以避免衰减系数达不到或者更换元器件的情况。但是由于控制衰减系数要用手调节，这样就造成可能调节不准确,造成产生误差，使测量不准确。因此舍弃此方案。方案三: T型网络能实现衰减并保持输入阻抗不变, 从输入端看进去的输入电阻应该等于R0，所以可得此电路具有良好的衰减效果，而且还不影响后面电路的输入，仿真效果非常好。综上所述

12、：选择方案三。3.4.1微弱信号检测前置放大电路设计方案一：采用两级放大电路，电路图原理图见3.4.1(a),第一级电路为同相比例运算电路，第二级为反相运算放大电路，利用叠加原理放大100倍，然后接一个反相电路以保证相位不变。此方案满足题目对输入阻抗必须大于1 MW的要求，但是在电阻选择上有困难，造成的误差比较大，而且用反相电路来保证相位不变结构繁琐。故舍弃。方案二：综合以上优缺点，设计出电路图见3.4.1(b)先用电压跟随器保证输入阻抗大于1 MW，然后用两级放大电路以方便电阻的选择并减小误差，同时保证相位不变。此方案满足题目要求，而且结构简单，误差较小，切实可行，故采用此方案。3.4.2微

13、弱信号检测电路设计方案一：带通滤波：但在检测噪声波形与频率围分析时发现，噪声频率大概在0到20kHz，而要求检测的信号频率在100Hz10kHz之间，完全被噪声信号淹没，因此不能采用带通滤波提取微弱信号。方案二：锁相放大：使用模拟乘法器AD633JN，该方案是应用互相干原理，可以提取同频同相的正弦信号，并且输出值为含有直流分量的正弦波，与输入的正弦波的幅值有一定的对应关系。因此采用此方案锁相放大原理：锁相放大器是由乘法器和积分器组成，其中乘法器采用模拟乘法器，积分器通常由低通滤波器组成，下图给出锁相放大器的构成原理图。设待测信号，其中为待测信号中的有效信号，n为噪声。参考信号，则经乘法器的输出

14、信号为，积分器的输出信号为：=上式中、分别是待测有效信号与参考信号与参考信号与噪声之间的相关函数，对于项，由于噪声的频率和相位都是随机量，可以认为参考信号与噪声互相独立，它的相关函数为零：即经过长时间的积分，项为零，实际上积分时间不可能无限长，也就是说实际计算中不可能为零，它表现为剩余噪声，但噪声对测量结果的影响已经大大降低，对于项为：=由上式可知，积分的两项是待测有效信号与参考信号的和频与差频相关，当时，且噪声相关函数为零时，积分器输出的信号可表示为：式中K为低通滤波器的传递系数，为待测信号与参考信号的相位差，当待测信号与参考信号同频率时，锁相放大器输出的信号与待测有效信号的幅度有关，也与待

15、测信号与参考信号的相位差有关，调整参考信号的相位，当=时，锁相放大器的输出信号与待测信号的幅度成正比。总之，由于输入信号中的噪声与参考信号不相关，经积分器积分后，噪声得到了很好的抑制，而有效信号得到了放大，从而提高了测量的精度和信噪比。锁相放大电路：AD633乘法器输出与输入的关系：Vo=，故输出为含有直流分量的交流量，经过低通滤波电路除去交流量后为直流量。2号和4号管脚接地，8号管脚和5号管脚分别串联一个0.1uf的电容接+15v与-15v电源，电容的作用是去噪。3.4.3低通滤波器的选择方案1：无源低通滤波器此电路无放大作用带负载后，通带放大倍数的数值减小，通带截止频率升高，不符合信号处理

16、的要求。方案2：有源低通滤波器要求：Aup=-1，fo尽量小。我们选择： C2=C1=1F R1=R2=R3=51k为了增强滤波效果，我们将两个二阶低通滤波器串联起来，得到四阶低通滤波器。3.5软件设计 3.5.1主控芯片方案方案一: 宏晶公司的STC89C52。STC89C52为8位单片机，价格低，程序易写，但I/O口少,存小，且需接外部AD，运算速度低，功耗高。故舍弃。方案二: ST公司的STM32f103。STM32为32位单片机，频率高，速度快，功耗低，性价比高，I/O口数量多达51个，存为128k,同时部自带12位AD，AD转换精度高，可简化电路，易于电路板设计和程序编写。故选此方案

17、。3.2.2显示电路方案设计方案一: 液晶显示。液晶可以显示数字、字符和汉字，但本项目只需显示数字便可完成要，同时液晶的程序控制相对复杂。方案二: 数码管显示。数码管只能显示数字，但程序控制相对简单易于实现。综上所述: 由于本项目只需显示峰峰值，只显示数字即可，所以采用方案二。因此采用此方案系统通过多次AD采样，将采集后的数据进行处理除去最大值、最小值，然后将剩下的数据取平均值，以减小由于信号波动带来的影响。再将平均值数据进行相应处理转换成对应的峰峰值，将峰峰值通过数码管显示出来。程序流程图如下：四、实验分析与测试结果先分别改变输入正弦波的峰峰值和频率，观察数码管显示的峰峰值，再同时改峰峰值和

18、频率观察数码管显示的峰峰值。结果分析通过反复采集数据和调试，对于基本部分要求200mv2v、误差0.5%基本完成，而对于发挥部分20mv200mv、误差0.2%只有部分几个段落的数据能够满足要求，其他部分的数据则不能达到误差要求。同时频率在该项目中的影响不大。 发挥部分不能达到要求的原因可能如下： 1、该项目为微信号检测，所以当测试小信号是稍微有一点信号的波动和噪音的影响都会对显示数据产生一定影响。而在实际的硬件电路中，波形发生器的稍微波动和线路中的噪音以与环境中的噪音都会对最后的结果产生一定的影响。 2、由于知识的欠缺，在软件方面，还有许多值得提高。在该函数程序中只采用了计算平均值、除去最大

19、值、最小值等最为简单方法来消除干扰。而对于其他高级算法像FFT、PID等，由于知识上的不足没有被应用到程序中。五、总结本次设计以锁相放大为检测微弱信号装置的核心，在参考信号与输入信号同频同相的前提下，使输出信号为直流量，然后经过单片机采样与数据处理通过数码管实现微弱信号的幅值显示。经过软件与硬件的调试，基本实现了题目设计的要求。但是本次设计还存在很多缺憾，比如发挥部分没有很好的完成，希望在以后的学习中能逐步完善。此次设计是我们团队共同努力的成果，渗透着每个队员的心血，从一开始的迷茫到讨论慢慢熟悉，我们在其中一点点共同进步，一步步完成了既定的目标。我相信只要我们不畏艰难定会取得更大的进步。六、参

20、考文献1 蒙博宇 STM32自学笔记 ：航空航天大学，2012.22 童诗白 华成英 模拟电子技术基础 ：高等教育，2006.53 AD633使用手册附件附件一：各模块电路原理图图3.1 图3. 2（a）图3.2（b） 图3.4.1（a） 图3.4.1（b）开始102次AD采样？100次取平均值处理数据数码管显示新的AD采样数据？否是是部分测试数据f=500Hz输入正弦波峰峰值数码管显示峰峰值相对误差0.100v0.145v45%0.150v0.165v10%0.200v0.216v8%0.300v0.302v0.70%0.400v0.405v1.25%0.500v0.502v0.40%0.6

21、00v0.625v4.17%0.700v0.721v3.00%0.800v0.815v1.88%0.900v0.912v1.33%1.000v1.011v1.10%1.500v1.513v0.86%2.000v1.999v0.05%f=5kHz 输入正弦波峰峰值数码管显示峰峰值相对误差0.100v0.140v40%0.150v0.175v16.7%0.200v0.213v6.5%0.300v0.309v3%0.400v0.407v1.75%0.500v0.512v2.4%0.600v0.617v2.83%0.700v0.715v2.14%0.800v0.812v1.50%0.900v0.912

22、v1.30%1.000v1.013v1.30%1.500v1.505v0.33%2.000v2.002v0.10%附件二:软件程序#include "stm32f10x_lib.h"#include "math.h"#define ADC1_DR_Address (u32)0x4001244C)/*数码管选通端为PB0，PB1,PB2，PA1*/#define PB0_H GPIO_SetBits(GPIOB,GPIO_Pin_0)#define PB0_L GPIO_ResetBits(GPIOB,GPIO_Pin_0)#define PB1_H GP

23、IO_SetBits(GPIOB,GPIO_Pin_1)#define PB1_L GPIO_ResetBits(GPIOB,GPIO_Pin_1)#define PB2_H GPIO_SetBits(GPIOB,GPIO_Pin_2)#define PB2_L GPIO_ResetBits(GPIOB,GPIO_Pin_2)#define PA1_H GPIO_SetBits(GPIOA,GPIO_Pin_1)#define PA1_L GPIO_ResetBits(GPIOA,GPIO_Pin_1)ErrorStatus HSEStartUpStatus;/* 函数申明 -*/void R

24、CC_Configuration(void);void GPIO_Configuration(void);void TIM1_Configuration(void);void NVIC_Configuration(void);void ADC1_Configuration(void);/void WWDG_Configuration(void);/void Delay(vu32 nCount); void Delayms (vu32 m); /* 变量定义 -*/vu16 ADC_ConvertedValue;double AD_temp,max,min;double AD_value;dou

25、ble temp102;int i,a=0,b=0,c=0,d=0;double data=0;double t=0;unsigned char table=0;u16 table010=0xBF,0x86,0xDB,0xCF,0xE6,0xED,0xFD,0x87,0xFF,0xEF;u16 table110=0x3F,0x06,0x5B,0x4F,0x66,0x6D,0x7D,0x07,0x7F,0x6F;/* 函数名称 ： main(void)* 函数功能 ： 主函数* 输 入： 无* 输 出： 无* 返 回： 无*/int main(void) RCC_Configuration();

26、 GPIO_Configuration();/ WWDG_Configuration(); ADC1_Configuration(); TIM1_Configuration(); NVIC_Configuration(); while(1) a=(int)AD_temp;PB0_L; PB1_H; PB2_H; PA1_H; GPIO_Write(GPIOC,table0a); Delayms(1);b=(int)(AD_temp*10);b=b%10;PB0_H; PB1_L; PB2_H; PA1_H; GPIO_Write(GPIOC,table1b); Delayms(1);c=(in

27、t)(AD_temp*100);c=c%10; PB0_H; PB1_H; PB2_L; PA1_H; GPIO_Write(GPIOC,table1c); Delayms(1);d=(int)(AD_temp*1000);d=d%10;PB0_H; PB1_H; PB2_H; PA1_L; GPIO_Write(GPIOC,table1d); Delayms(1); /* 函数名称 ： RCC_Configuration(void)* 函数功能 ： 时钟初始化* 输 入： 无* 输 出： 无* 返 回： 无*/void RCC_Configuration(void) /* RCC syste

28、m reset(for debug purpose) */ RCC_DeInit(); /* Enable HSE */ RCC_HSEConfig(RCC_HSE_ON); /* Wait till HSE is ready */ HSEStartUpStatus = RCC_WaitForHSEStartUp(); if(HSEStartUpStatus = SUCCESS) /* HCLK = SYSCLK */ RCC_HCLKConfig(RCC_SYSCLK_Div1); /* PCLK2 = HCLK */ RCC_PCLK2Config(RCC_HCLK_Div1); /* P

29、CLK1 = HCLK/2 */ RCC_PCLK1Config(RCC_HCLK_Div2); /* Flash 2 wait state */ FLASH_SetLatency(FLASH_Latency_2); /* Enable Prefetch Buffer */ FLASH_PrefetchBufferCmd(FLASH_PrefetchBuffer_Enable); /* PLLCLK = 8MHz * 9 = 72 MHz */ RCC_PLLConfig(RCC_PLLSource_HSE_Div1, RCC_PLLMul_9); /* Enable PLL */ RCC_P

30、LLCmd(ENABLE); /* Wait till PLL is ready */ while(RCC_GetFlagStatus(RCC_FLAG_PLLRDY) = RESET) /* Select PLL as system clock source */ RCC_SYSCLKConfig(RCC_SYSCLKSource_PLLCLK); /* Wait till PLL is used as system clock source */ while(RCC_GetSYSCLKSource() != 0x08) /* 函数名称 ： GPIO_Configuration(void)*

31、 函数功能 ： 端口初始化* 输 入： 无* 输 出： 无* 返 回： 无*/void GPIO_Configuration(void) GPIO_InitTypeDef GPIO_InitStructure; RCC_APB2PeriphClockCmd( RCC_APB2Periph_GPIOA | RCC_APB2Periph_GPIOB | RCC_APB2Periph_GPIOC | RCC_APB2Periph_GPIOD, ENABLE ); /* Configure USART1 Tx (PA.09) as alternate function push-pull */ GPI

32、O_InitStructure.GPIO_Pin = GPIO_Pin_0| /PC0到PC7管脚数码管输出 GPIO_Pin_1|GPIO_Pin_2|GPIO_Pin_3|GPIO_Pin_4|GPIO_Pin_5|GPIO_Pin_6|GPIO_Pin_7; GPIO_InitStructure.GPIO_Mode = GPIO_Mode_Out_PP; / 推挽输出 GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz; / 最高输出速率50MHz GPIO_Init(GPIOC, &GPIO_InitStructure); / 选择

33、C端口 GPIO_InitStructure.GPIO_Pin = GPIO_Pin_0| GPIO_Pin_1|GPIO_Pin_2; /数码管片选端 GPIO_InitStructure.GPIO_Mode = GPIO_Mode_Out_PP; /浮空输入 GPIO_Init(GPIOB, &GPIO_InitStructure); /选择B端口 GPIO_InitStructure.GPIO_Pin = GPIO_Pin_1; /数码管片选端 GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz; / 最高输出速率50MHz GPIO

34、_InitStructure.GPIO_Mode = GPIO_Mode_Out_PP; /推挽输出 GPIO_Init(GPIOA, &GPIO_InitStructure); /* Configure PC.04 (ADC Channel14) as analog input*/ GPIO_InitStructure.GPIO_Pin = GPIO_Pin_0;/ 选中管脚0 GPIO_InitStructure.GPIO_Mode = GPIO_Mode_AIN;/ 模拟输入 GPIO_Init(GPIOA, &GPIO_InitStructure);/选择A端口/* 函

35、数名称 ： ADC1_Configuration(void)* 函数功能 ： ADC1初始化* 输 入： 无* 输 出： 无* 返 回： 无*/void ADC1_Configuration(void) ADC_InitTypeDef ADC_InitStructure; RCC_APB2PeriphClockCmd(RCC_APB2Periph_ADC1, ENABLE ); /* Time Base configuration */ ADC_InitStructure.ADC_Mode = ADC_Mode_Independent;/独立工作模式 ADC_InitStructure.ADC

36、_ScanConvMode = ENABLE;/扫描方式 ADC_InitStructure.ADC_ContinuousConvMode = ENABLE;/连续转换 ADC_InitStructure.ADC_ExternalTrigConv = ADC_ExternalTrigConv_None;/外部触发禁止 ADC_InitStructure.ADC_DataAlign = ADC_DataAlign_Right;/数据右对齐 ADC_InitStructure.ADC_NbrOfChannel = 1;/用于转换的通道数 ADC_Init(ADC1, &ADC_InitSt

37、ructure); ADC_RegularChannelConfig(ADC1, ADC_Channel_0, 1, ADC_SampleTime_55Cycles5);/采样时间为55.5周期 ADC_DMACmd(ADC1, ENABLE); / 使能 ADC1 DMA ADC_Cmd(ADC1, ENABLE); / 使能 ADC1 ADC_ResetCalibration(ADC1); / 重置指定的ADC的校准寄存器 while(ADC_GetResetCalibrationStatus(ADC1); /被选中的ADC重置校准寄存器状态 Check the end of ADC1 r

38、eset calibration register ADC_StartCalibration(ADC1); / Start ADC1 calibration 开始ADC1校准 while(ADC_GetCalibrationStatus(ADC1); /检查ADC1校准的结束Check the end of ADC1 calibration ADC_SoftwareStartConvCmd(ADC1, ENABLE);/开始ADC1软件转换Start ADC1 Software Conversion /* 函数名称 ： TIM1_Configuration(void)* 函数功能 ： 定时器1

39、初始化* 输 入： 无* 输 出： 无* 返 回： 无*/void TIM1_Configuration(void) TIM_TimeBaseInitTypeDef TIM_TimeBaseStructure; RCC_APB2PeriphClockCmd(RCC_APB2Periph_TIM1, ENABLE ); /* Time Base configuration */ TIM_DeInit(TIM1); TIM_TimeBaseStructure.TIM_Prescaler = 7200-1; /设置预分频器分频系数71，即APB2=72M, TIM1_CLK=72/72=1MHz ，

40、 /它的取值必须在0x0000和0xFFFF之间 TIM_TimeBaseStructure.TIM_CounterMode = TIM_CounterMode_Up; / 向上计数模式 TIM_TimeBaseStructure.TIM_Period = 5000-1; / 500ms定时 它的取值必须在0x0000和0xFFFF之间设置了在下一个更新事件装入活动的自动重装载寄存器周期的值 TIM_TimeBaseStructure.TIM_ClockDivision = 0x0; /设置了定时器时钟分割， 时钟分频 TIM_TimeBaseStructure.TIM_RepetitionC

41、ounter = 0x0; / 多少次更新一次 周期计数 重复计数 设置了周期计数器值，它的取值必须在0x00和0xFF之间。 TIM_TimeBaseInit(TIM1,&TIM_TimeBaseStructure); / 根据TIM_TimeBaseInitStruct中指定的参数初始化TIMx的时间基数单位 TIM_ClearFlag(TIM1, TIM_FLAG_Update); /清中断，以免一启用中断后立即产生中断 TIM_ITConfig(TIM1, TIM_IT_Update, ENABLE); /使能TIM1中断源 TIM_Cmd(TIM1, ENABLE); /TI

42、M1总开关：开启 /* 函数名称 ： TIM1_UP_IRQHandler(void)* 函数功能 ： 定时器1中断* 输 入： 无* 输 出： 无* 返 回： 无*/ void TIM1_UP_IRQHandler(void) double AD_value=0;for(i=0;i<102;i+) tempi=ADC_GetConversionValue(ADC1);/ AD_value=AD_value+ADC_GetConversionValue(ADC1);/ AD_value=AD_value/100.0; max=temp0; min=temp0; for(i=0;i<

43、102;i+) if(max<tempi) max=tempi;else if(min>tempi) min=tempi; AD_value=AD_value+tempi; AD_value=(AD_value-max-min)/100.0;AD_temp=AD_value*3.3/4096;if(AD_temp<0.256890) /600mv if(AD_temp<0.095681) /300mv AD_temp=-19789*pow(AD_temp,4)+7315.6*pow(AD_temp,3)-969.03*pow(AD_temp,2)+58.221*AD_temp-1.1463; else AD_temp=-79839*pow(AD_temp,6)+76318*pow(AD_temp,5)-29490*pow(AD_temp,4)+5909.5*pow(AD_temp,3)-655.45*pow(AD_temp,2)+41.088*AD_temp-0.8917;else AD_temp=-0.004*pow(AD_temp,6)+0.0503*pow(AD_temp,5)-0.2647*pow(AD_temp,4)+0.7533*pow(AD_temp,3)-1.2979*pow(AD_temp,2)+1.8342*AD_temp+0.2063