根据示波器存储的波形数据得到两列信号的相位差（MATLAB源码）

最新推荐文章于 2023-11-20 15:12:02 发布

liuz_notes

最新推荐文章于 2023-11-20 15:12:02 发布

阅读量3.5k

点赞数 5

分类专栏：杂项文章标签： matlab 信号处理

本文链接：https://blog.csdn.net/liuz_notes/article/details/103801378

版权

杂项专栏收录该内容

11 篇文章 9 订阅

订阅专栏

目的

示波器显示的波形如下，需要通过MATLAB处理得到这两列信号的相位差。
在这里插入图片描述

步骤

将U盘插入到示波器上，存储信号序列。在MATLAB中读取这个.CSV 信号序列文件，处理得到相位差。代码如下：

clear all
close all
%读取数据
t=xlsread('A2.CSV','A2','A1:A500');%读取时间序列
x=xlsread('A2.CSV','A2','B1:B500');%读取参考信号
x1=xlsread('A2.CSV','A2','C1:C500');%读取待测信号
subplot 211;
plot(t,x,'b-');hold on;
plot(t,x1,'r--');grid;
legend('Reference','Signal');title('Origin');

%数据归一化
m1=(max(x)+min(x))/2;
m2=(max(x1)+min(x1))/2;
x=(x-m1)/max(x);   
x1=(x1-m2)/max(x1);
subplot 212;
index=find(x==0);
plot(t,x,'b-');hold on;
plot(t,x1,'k--');grid;
hold on;
legend('Reference','Signal');title('Normalization');

num=500;  %数据数量
l=0:1:num;
x=x;
y=x1;
Ix=sum(x.^2)/num;
Iy=sum(y.^2)/num;
Ixy=sum(x.*y)/num;
c=180*acos(Ixy/(Ix*Iy)^0.5)/pi;  %计算相位差
%disp(['相位差为',num2str(c), '度']);
output = sprintf('相位差为%f度', c);
output