模拟互相关滤波笔记

实验学习：

1、了解互相关滤波原理和物理实现，xcorr(x)函数；

2、相关计算、低通滤波及其检测；

参考文献：

信号处理趣学D1——相关函数的意义&利用自相关函数消除噪声-CSDN博客

实验原理：

        使用信号发生器产生，Vpp=5V，f=50HZ，相位0°、28°、45°、 69°、90°、118°、135°的正弦波，接入滤波信号端 IN1。参考信号端接入 Vpp=5V，f=50HZ，的方波信号。依次测得不同相角下的输出 X，Y，并利用 X、Y 计算A0，φ0，与输入做比较。当正弦信号与方波信号相位差为0°、90°时， 用示波器观察两信号相乘后的全波、半波检波结果。

% 本代码用于研究自相关函数
clc; clear; close all;
% %% 生成函数序列
% fs = 1000;
% T = 1/fs;
% N = 6000; % 信号长度
% t = (0:N-1)*T; %时间序列
% f0 = 8;
% x = sin(2*pi*f0*t); %原信号
% z=x+randn(size(x)); %加入噪音之后的信号
% % x1=circshift(x,100); 
% figure;subplot(211); plot(t,z);title('时域信号');
% 
% %% 自相关R（0）
% % 时延为0
% [c,lags]=xcorr(z,x,600,'coeff'); %600是指定时延
% subplot(212); plot(lags,c); title('无差自相关');

%% 自己实现自相关
% y=zeros(1,2*N-1);
% for k=1:2*N-1
%     for n=1:N
%         if k-n>=1 && k-n<=N % 保持自相关和自己有交点
%             y(k)=y(k)+x(n).*x(k-n);
%         end
%     end
% end
% t1=1:2*N-1; t1=t1-500;
% figure; plot(t1,y);

%% 参考信号时方波
fs = 1000;
T = 1/fs;
N = 6000; % 信号长度
t = (0:N-1)*T; %时间序列
f0 = 50;
aa = 0*pi/180; %设置相位弧度制
x = 5*cos(2*pi*f0*t+aa); %原信号 相位是弧度
a = 100*pi; %周期参数
y = 5*square(a*t); %方波周期 T=2pi/a
y1 = circshift(y,25); %产生相位差90°的方波
subplot(311);
plot(t,x,t,y,t,y1); title('Sin(x) & Square'); axis([0 0.1 -15 15]);

%% 相关计算
c = x.*y;
c1 = x.*y1;
subplot(312);
plot(t,c,t,c1); title('全波整流');axis([0 0.1 -40 40]);
% [c,lags]=xcorr(y,x,200);
% [c1,lags1]=xcorr(y1,x,200);
% subplot(312);
% plot(lags,abs(c),lags1,abs(c1)); title('全波整流');

%% 低通滤波
% d=lowpass(c,0.01,fs);
% d1=lowpass(c1,1,fs);

fc_low = 1.5; % 低截止频率
fc_high = 30; % 高截止频率
NN = 7; % 滤波器阶数
[b1, a1] = butter(NN, fc_low/(fs/2), 'low'); % 计算低通滤波器系数
% 使用中通滤波器对信号进行滤波
d = filter(b1, a1, c);
d1 = filter(b1,a1,c1);

%% 低通滤波器检测
F_y=fft(d); F_x=fft(d1);
P2=abs(F_x/N); %双边频谱
P1=P2(1:N/2+1); %单边频谱
P1(2:end-1)=2*P1(2:end-1); %单边频谱翻倍
f=fs*(0:(N/2))/N;
% subplot(313);
% plot(f,P1);

subplot(313);
plot(t,d,t,d1); title('低通滤波');
disp('通道Y的值为'); mean(abs(d))
disp('通道X的值为'); mean(abs(d1))

    

实验中有三个问题没解决：

1.为什么改变相位（参数aa）为0的时候，并不是理论上Y通道的最小值，而是在aa=9附近有最小值?

2.通过低通滤波之后为什么值是从0开始的呢？

3.改变fs（增大）也会明显改变结果，但是理论上采样频率更高应该不会降低实验结果？