clc;
clear;
data=load('D:\桌面\中期\paper\整理\000数据库\zjy.mat');
Au=data.zjy(:,6);
Au=Au(5384:7384,1);
Au=Au-mean(Au);
Fs = 200;%采样频率
N = length(Au);%求取抽样点数
t = (0:N-1)/Fs;%显示实际时间
% t = 0:1/Fs:(N-1)/Fs; %生成信号时间
plot(t,Au);
axis([0 10,-1.5,1.5]);
xlabel('时间t/s');
set(gca,'Fontsize',14)%设置字体大小为14
ylabel('Amplitude');
关于如何把一个数据导入进来画出来相应的时域图,对于初入门的学者来说还是比较简单。最后给大家展示一下这个代码最后的呈现效果。
下面我简单讲解一下这个代码,也是一个比较简单的画时域代码
data=load('D:\桌面\中期\paper\整理\000数据库\zjy.mat');
Au=data.zjy(:,6);
Au=Au(5384:7384,1);
这个里面是把自己的数据导入进来,但是由于导入的.mat文件里面还有一层,所以需要第二行,也就是data.zjy(:,6),这个意思是data这个盒子打开,发现里面还有个盒子叫zjy,要打开zjy,必须通过一个格式,就是data.zjy,注意,后面的括号里面第一个是行,第二个是列。
(:,6)的意思是所有行,第6列。
我们继续往后看
Au=Au-mean(Au);
这是啥意思?其实和我放入的数据有关,我放一下原来的数据:
发现没有?这里的数据居然不是以0线为基本波动,我们平时处理还是建议在0上下浮动。
这里就是减去它的平均值,做了一个转换。转换之后的效果:
很明显看出来,调整好了,在0上下波动。这也是比较简单的处理方法。
Fs = 200;%采样频率
N = length(Au);%求取抽样点数
t = (0:N-1)/Fs;%显示实际时间
% t = 0:1/Fs:(N-1)/Fs; %生成信号时间
这又是啥情况?
因为我们导入数据,是没有时间的,要做时域图的x轴坐标是时间,我们只能通过设置采样频率来求时间
Fs是采样频率,也就是每秒钟采样的点数,N是数据的中共点数,t是不是一下子好求了?
这里我们的t还有一个表示
t = (1:N)/Fs;%显示实际时间
t点和Au点要对应,Au有2001个点,每个点出现对应一个时间,那么这个时间点也有2001个
最后的就不过多解释了,就是把图像呈现出来