基于matlab的相关器仿真程序编写.docx

PAGE1

微弱信号检测技术

实验一基于matlab的相关器仿真程序编写

一、目的要求

了解相关器的原理

掌握基本的matlab信号处理编程方法

测量相关器的相位响应曲线

二、基本原理

相关器由乘法器与积分器组成，是锁定放大器的核心部件，其中，积分器也可以看作是一个低通滤波器。锁定放大器中的相关器的结构如图1。

图1 相关器结构

图中，Vin为信号输入端，其可以包含被测信号，也可以包含外来噪声信号。通常情况下，我们以傅里叶分解的观点，以其包含的一个频率分量对其进行描述

V

VR为参考信号输入端，一般情况下为方波，其基频的角频率记为ωR。由于相关器中使用的乘法器一般不用模拟乘法器，而是采用开关电路实现，参考信号

V

当输入信号Vin的角频率ω=ωR时，表达被测信号，而当ω≠ωR时则为噪声或干扰。Vin和VR之间的相位差在锁定放大器中可以通过参考通道的相移电路调节。图1

V1=Vin?VR

V

(1-2)

式中，Q

各符号R1、R0、C0等如图1所示。其中，

从公式（1-2）可见，对相关器进行直接的解析计算，即使输入信号是给定频率的单频正弦，也会非常繁琐。当输入信号包含复杂情况噪声时，解析计算变的不可实现。在本实验中，采取将相关器分为乘法器和积分器两个独立的模块，并通过求取积分器的冲激响应的方式，利用卷积的方法得出其时域响应信号。具体方法参考实验步骤与内容。

实验步骤与内容

1.基于matlab的相关器功能仿真程序编写

(1)根据图1的相关器结构可知，需要的考虑的信号有：

a.信号输入端Vin与参考输入端VR，涉及到两者的频率f和fR，以及两者之间的相位差?

b.与输入信号可能同时存在的噪声Vn，以及输入信噪比

c.电阻、电容R1、R0、

d.乘法器输出V1

e.相关器输出Vout

除此之外，由于程序需要用离散信号的方式处理，需要设定采样频率fs

由如上信号，可以编写相关器功能仿真程序的顶层代码，参考main.m。

(2)生成输入信号：根据采样频率生成时间序列t，总持续时长建议设定为10tc，生成被测信号Vs=sin2πft+?

(3)乘法器模块：以函数形式编写，Vin和VR进行逐点相乘（matlab中的.

乘法器模块编写完成后，将main.m中乘法器之后的所有语句屏蔽掉，设定f=fr，根据要求调整?的取值，并对每一个相位差绘制输入信号、参考信号、以及乘法器输出信号，并将所绘信号粘贴入表1（请调整横坐标范围，使曲线清晰

表1.不同相位差下的乘法器输出信号波形

?=0?：

?=π

?=π

?=π

?=π

?=π?：

(4)积分器模块：根据图1所示的相关器结构，积分器模块的输入信号为V1，关键影响参数为Tc

C

将其改写为s域方程：

C

可得积分器的s域传输函数为：

H

所以，积分器的冲激响应为：

h

设R0=

h

将h(t)在程序中生成为与t等长的数组，则相关器输出为：

V

Matlab中的卷积运算为y=conv(x1,x2)，输出结果长度为输入信号长度的2倍，仅需保留前一半的结果。

程序编写完成后，设定f=fr，tc=0.1s，fs=8000Hz，?=0，分别绘制fr=100Hz与fr=10Hz时的相关器输出波形，包括0到1秒之间的整体波形，以及0.7至0.8秒间的细节波形，贴入表2，对波形进行对比，并对其差异进行分析。

表2.相关器输出波形

fr=100Hz，整体波形

fr=10Hz，整体波形

fr=100Hz，细节波形

fr=10Hz，细节波形

波形主要区别：

答：频率较低的时候，波形的幅度比较大。当频率较高的时候，波形的幅度较小。

原因分析：

答：输入频率增高，积分时间缩短，在积分时间常数不变的情况下，积分所得的电压下降，输出的幅度就减小了。

2.相关器的相位响应曲线测量

（1）利用前一步编写的相关器仿真程序，设定?=0，观察输入信号频率为参考信号频率的1、2、3、4、5倍时的输出情况，并设法计算其稳态输出幅度，设计表格，记录输出波形与稳态输出幅度

?

频率倍数

输出波形

输出幅度

1

-0.6337

2

0.0164

3

-0.215

4

0.0081

5

-0.125

（2）观察输入信号频率为参考信号频率的1.1倍、3.2倍、4.3倍时的输出情况，并计算其输出幅度（交流幅度），设计表格，记录输出波形与输出幅度；

频率倍数

输出波形

输出幅度

1.1

0.105

3.2

0.021

4.3

0.0085

（3）总结输入信号频率与参考信号频率之间呈谐波关系及无谐波关系情况下的输出幅度的规律，从f=0.1fr为起点，频率步进值为0.1fr，编写程序，绘制[0.1fr,10fr]范围内的相关器频率响应曲线。

对所生成的频率响应曲线的描述与分析解释：

答：随着频率的升高，幅度开始逐渐降低