.mat文件在Matlab中的读取使用操作

已于 2023-07-10 17:48:43 修改
阅读量9.4w 收藏 553
点赞数 64
分类专栏： Matlab 文章标签： matlab
于 2022-06-30 13:28:34 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_49220054/article/details/125538440
版权

目录

前言 

1. 数据中的采样频率为100Hz, 根据所分析的信号的特点确定合适的采样频率，截断长度（记录点数）

2. 通过DFT频域分析方法，编程画出该信号的频谱与时域图形

3. 使用滤波器设计方法将这些频率成分分别获得，并以图形显示

4. 将滤得的频率成分叠加生产原信号，并与给的标准原信号进行对比

前言 

mat数据格式是Matlab的数据存储的标准格式，.mat文件的方便之处在于其可以连同数据的变量名一同保存下来，并且不需要控制数据的存储格式，matlab会自动保存并区分我们所存储的内容，可以利用load指令可以实现对.mat文件的读取与使用。

利用load指令读取.mat文件：

load('E:\文件名'); //不要加.mat后缀

以下图为例，读取并画出.mat文件中前160个点形成的波形。 左侧双击.mat文件后，下方出现了mat文件内部变量的名称与类型，处理数据或读写的时候要以变量的名称来处理。

clear all;
load('E:\MATLAB\zibiancx\yy59');

for i=1:160   %数组索引必须为正整数或逻辑值,所以i从1开始
    x(i)=yy9(i);
    t(i)=i/100;
end
t = t-0.01;  %i从1开始，这里减去1
figure(1);
plot(t,x);title('前160个点'),xlabel('时间/s');

1. 数据中的采样频率为100Hz, 根据所分析的信号的特点确定合适的采样频率，截断长度（记录点数）

         绘制出所有点的图，发现这是一个周期信号。原始信号采样频率100Hz，原始采样周期0.01s，结合前160点图进行分析发现信号周期约为0.5s，画出第50-100个点（0.5-1s）的波形如图所示：

%1.确定合适的采样频率，截断长度（记录点数）
clear all;
load('E:\MATLAB\zibiancx\yy59');
for i=1:4096
    x(i)=yy9(i);
    t(i)=i;
end
figure(2);
plot(t,x);title('yy59所有点'),xlabel('点数');

clear all;
load('E:\MATLAB\zibiancx\yy59');
Fs = 100;  % 原始信号采样频率,  100Hz
Ts = 1/Fs; % 原始信号采样周期,  0.01s
L = 50;    % 序列长度，每20点近似为一个周期
N = 50;    % 频率域采样个数，取N=L
for i=1:L+1    %取第50-100点进行处理
    x(i)=yy9(i+50);
    t(i)=i/Fs+0.5;
end
figure(3);
plot(t,x);title('yy59 0.5-1s内的信号波形'),xlabel('时间/s');

        根据香农采样定理:如果一个系统以超过信号最高频率至少两倍的速率对模拟信号进行均匀采样，那么原始模拟信号就能从采样产生的离散值中完全恢复。所以，先对信号进行DFT频域分析（2. 通过DFT频域分析方法，编程画出该信号的频谱与时域图形）获得频谱，该信号包含12Hz、16Hz、30Hz的信号，信号的最高频率分量为30Hz，所以采样频率应该大于等于60Hz，而给出的数据采样频率为100Hz，符合采样定理的要求。故采样频率沿用100Hz，截断长度（记录点数）取第51-100点，该部分点信号比较完整。

2. 通过DFT频域分析方法，编程画出该信号的频谱与时域图形

         MATLAB的fft函数，其用法为Y=fft（x，N），其中x为原信号序列，Y为DFT（也就是FFT）变换后的，FFT与DFT等价，FFT为DFT的算法优化，fft函数中有三个参数：采样频率Fs，原序列长度L，N-DFT变换的N。令Fs=100(Hz)，L=50，N=50，对50—100点内的信号进行编程处理。

%2.DFT频域分析画出信号频谱
%DFT
clear all;
load('E:\MATLAB\zibiancx\yy59');
Fs = 100;  % 原始信号采样频率,  100Hz
Ts = 1/Fs; % 原始信号采样周期,  0.01s
L = 50;    % 原信号序列长度，每50点近似为一个周期
N = 50;    % 频率域采样个数，取N=L
t=(0:L-1).*Ts;%时域自变量
for i=1:L
    x(i)=yy9(i+50);
end%原信号
figure(4)
subplot(311);
plot(t,x);
title("原信号"),xlabel("t/s");
grid on;
Y=fft(x,N);%fft变换
Y=abs(Y)./N;%实际幅值变换
f=(0:N-1)*Fs./N;%实际频率变换
subplot(312);
stem(f(1:N/2),Y(1:N./2));
title("N-DFT变换幅频响应单边"),xlabel("f/Hz");
grid on;
f=f-Fs./2;%移位
subplot(313);
stem(f,fftshift(Y));%移位
title("N-DFT变换幅频响应双边"),xlabel("f/Hz");
grid on;
%时域内各分量与合成
clear all;
load('E:\MATLAB\zibiancx\yy59');
Fs = 100;  % 原始信号采样频率,  100Hz
Ts = 1/Fs; % 原始信号采样周期,  0.01s
L = 50;   % 使用原始序列个数，接近一个周期
for i=1:L
    x(i)=yy9(i+50);
    t(i)=i/Fs;
end
figure(5)
subplot(511);
t = t-Ts;
plot(t,x);title('原信号'),xlabel('时间/s');grid on
subplot(512);
t=(0:L-1).*Ts;
x=sin(12*2*pi*t);
plot(t,x,'r');title('12Hz分量'),xlabel('时间/s');grid on
subplot(513);
t=(0:L-1).*Ts;
x=sin(16*2*pi*t);
plot(t,x, 'g');title('16Hz分量'),xlabel('时间/s');grid on
subplot(514);
t=(0:L-1).*Ts;
x=sin(30*2*pi*t);
plot(t,x, 'b');title('30Hz分量'),xlabel('时间/s');grid on
subplot(515);
for i=1:L
    x(i)=yy9(i+50);
    t(i)=i/Fs;
end
t = t-Ts;
plot(t,x);title('合成'),xlabel('时间/s');hold on;
t=(0:L-1).*Ts;
x=sin(12*2*pi*t);
plot(t,x, 'r');hold on;
t=(0:L-1).*Ts;
x=sin(16*2*pi*t);
plot(t,x, 'g');hold on;
t=(0:L-1).*Ts;
x=sin(30*2*pi*t);
plot(t,x, 'c');hold on;

         通过将一个周期的50个点进行DFT频域变换，分析发现该信号是由12Hz、16Hz和30Hz的频率信号组成的。下图为三个分量信号的时域图与叠加后的图像。

3. 使用滤波器设计方法将这些频率成分分别获得，并以图形显示

         利用MATLAB Filter Designer设计低通滤波器对信号进行滤波。

%3.使用滤波器设计方法将这些频率成分分别获得，并以图形显示。
%12Hz
clear all;
load('E:\MATLAB\zibiancx\yy59');
Fs = 100;  % 原始信号采样频率,  100Hz
Ts = 1/Fs; % 原始信号采样周期,  0.01s
L = 160;   % 使用原始序列个数
for i=1:L
    y(i)=yy9(i);
    t(i)=i/100;
end

figure(6);
t = t-0.01;
subplot(311);
plot(t,y);title('时域原始信号'),xlabel('时间/s');hold on;

output = filter(filter12, y);
subplot(312);
t = t+0.01;
plot(t,output, 'b');title('时域滤波信号'),xlabel('时间/s');hold on;

y=sin(12*2*pi*t);
subplot(313);
plot(t,y, 'r');title('拟合信号(12Hz)'),xlabel('时间/s');
%16hz
clear all;
load('E:\MATLAB\zibiancx\yy59');
Fs = 100;  % 原始信号采样频率,  100Hz
Ts = 1/Fs; % 原始信号采样周期,  0.01s
L = 160;   % 使用原始序列个数
for i=1:L
    y(i)=yy9(i);
    t(i)=i/100;
end

figure(7);
t = t-0.01;
subplot(311);
plot(t,y);title('时域原始信号'),xlabel('时间/s');hold on;

output = filter(filter16, y);
subplot(312);
t = t+0.01;
plot(t,output, 'b');title('时域滤波后信号'),xlabel('时间/s');hold on;

y=sin(16*2*pi*t);
subplot(313);
plot(t,y, 'r');title('拟合信号(16Hz)'),xlabel('时间/s');
%30hz
clear all;
load('E:\MATLAB\zibiancx\yy59');
Fs = 100;  % 原始信号采样频率,  100Hz
Ts = 1/Fs; % 原始信号采样周期,  0.01s
L = 160;   % 使用原始序列个数
for i=1:L
    y(i)=yy9(i);
    t(i)=i/100;
end

figure(8);
t = t-0.01;
subplot(311);
plot(t,y);title('时域原始信号'),xlabel('时间/s');hold on;

output = filter(filter30, y);
subplot(312);
t = t+0.01;
plot(t,output, 'b');title('时域滤波后信号'),xlabel('时间/s');hold on;

y=sin(30*2*pi*t);
subplot(313);
plot(t,y, 'r');title('拟合信号(30Hz)'),xlabel('时间/s');


        （1）12Hz低通滤波器：滤除16Hz和30Hz信号，如图所示，分别为时域原始信号，滤波后信号以及12Hz正弦产生信号用来对比。

         （2）16Hz带通滤波器

         （3）30Hz高通滤波器

        通过分析发现信号滤波之后相位会发生一定的滞后 。

4. 将滤得的频率成分叠加生产原信号，并与给的标准原信号进行对比

%4.将滤得的频率成分叠加生产原信号，并与给的标准原信号进行对比。
clear all;
load('E:\MATLAB\zibiancx\yy59');
Fs = 100;  % 原始信号采样频率,  100Hz
Ts = 1/Fs; % 原始信号采样周期,  0.01s
L = 160;   % 使用原始序列个数
for i=1:L
    z(i)=yy9(i);
    t(i)=i/100;
end
figure(9);
output1 = filter(filter12, z);
output2 = filter(filter16, z);
output3 = filter(filter30, z);
subplot(511);
plot(t,output1);title('滤波信号(12Hz)'),xlabel('时间/s');
subplot(512);
plot(t,output2, 'b');title('滤波信号(16Hz)'),xlabel('时间/s');
subplot(513);
plot(t,output3, 'g');title('滤波信号(30Hz)'),xlabel('时间/s');
subplot(514);
plot(t,output1+output2+output3, 'r');title('叠加信号'),xlabel('时间/s');
subplot(515);
plot(t,z, 'r');title('原始信号'),xlabel('时间/s');

        分析可知，从0.8s开始，叠加信号开始稳定，并且呈现出周期性，周期与原始信号一致，但是图像形状与原始信号有差距。可能是因为滤波器参数设置问题，也可能是因为滤波过程中有损失等

下图则为程序运行的全数据：

李小风风风
关注
  • 64
    点赞
  • 踩
  • 553
    收藏
    觉得还不错? 一键收藏
  • 13
    评论
专栏目录
matlab如何读取.mat文件,matlab中读取mat文件
weixin_34280468的博客
03-18 6033
Matlab 训练好的 BP 神经网络如何保存和读取方法 附实例说 (明)看到论坛里很多朋友都在提问如何存储和调用已经训练好的神经网络。 本人前几天也遇到了这样的问题,......2. 2 HDF 文件数据的读取 Mat lab 提供了三种从 H DF 或者 H DF2EOS 数据集中读取数据的方法: H DF 输入工具、 Matlab 高级输 入函数以及 Mat lab 命令......= st...
matlab.mat文件读取,matlab中mat文件的生成和读取
weixin_39548733的博客
03-18 9343
1、mat文件的生成(1)直接在Matlab中创建并保存矩阵数据打开Matlab软件，点击左上角文件(File)，然后点击新建(new)，选择变量(Variable)，就新建了一个mat文件。点击你新建的文件，即可打开一个编辑器，输入数据即可。(2)保存单一矩阵数据如图所示，输入数据，生成数组。最后输入save x1 x。该语句表示将敲入的数组x生成数据x1.mat文件。(3)如果矩阵A.B都已赋...
【Matlab】Matlab导入多个.mat文件并画图的过程详解
热门推荐
weixin_44201830的博客
11-20 1万+
matlab导入多个.mat文件，并根据文件中的数据进行画图，图像绘制包括散点图、折线图、子图、图例设置等。
.mat文件打开方式方法
qq_50988206的博客
11-03 8265
通过使用MATLAB的数据处理工具箱，可以对.mat文件中的数据进行滤波、降噪、插值等操作，为后续分析提供更准确、可靠的数据。.mat文件是MATLAB软件中使用的二进制文件格式，用于存储多维数组、矩阵、结构和其他数据类型。通过使用MATLAB的信号处理和图像处理工具箱，可以对.mat文件中的信号和图像进行滤波、压缩、特征提取等操作。5、数值计算和仿真：.mat文件存储了数值计算和仿真过程中生成的数据。.mat文件作为MATLAB的数据文件格式，在科学研究、数据分析和工程领域具有广泛的应用。
Matlab批量读取mat格式文件
02-13
针对指定文件夹内存在的较多的文件，需要批量读取。本代码提供一种批量读取数据并组合成为一个数组的方式。
mat文件的读写. mat 文件 MATLAB读写操作
09-07
在 MATLAB 中，使用 `save` 函数可以将工作空间中的变量保存到 MAT 文件中。例如，要将数组 A 保存到名为 `A.mat` 的文件中，可以按照以下方式编写代码： 1. 定义数组 A： ```matlab A = [1 2 3 4 5; 6 7 8 9 0]; `...
matlab批量读取同一路径下的mat文件.rar_mat _matlab mat_strugglewru_批量 mat_读取m
07-15
从同一路径下，批量读取mat文件到matlab中，以便进行下一步计算
matio-1.3.4.zip_.mat文件读取_matio_matio-1.3.4_matio.h_读取 .mat文件
07-13
在处理与MATLAB数据交换的过程中，有时我们需要在C/C++环境中读取和操作.mat文件。这时，matio库便成为一个重要的工具。matio-1.3.4是matio库的一个版本，它提供了方便的接口，使得开发者能够轻松地将.mat文件转换为...
MAT.jl：用于读取MATLAB文件的Julia模块
02-03
MAT.jl是Julia语言中的一个模块，专为读取和写入MATLAB文件而设计。这个库使得在Julia环境中处理由MATLAB生成的数据变得简单便捷，这对于跨平台的数据交换和分析工作尤为有用。MATLAB文件格式通常用于存储数组、矩阵...
mat文件的读与写
青锋断尘
02-08 1640
.mat文件的 创建.mat文件 mat文件简介 ​ mat文件是matlab的数据存储的标准格式。mat文件是标准的二进制文件,还可以ASCII码形式保存和加载，在MATLAB中打开显示类似于单行EXCEL表格。 优点 ​ .mat文件的方便之处在于，它可以连同数据的变量名一同保存下来（这一点在程序调试的过程中很方便），并且不需要控制数据的存储格式，matlab会自动保存并区分我们所存储的内容。...
MATLAB R2017a 安装与破解 ​编辑
Hello-FPGA
07-29 412
转载自。
MATLAB基础：函数与函数控制语句
Zjs15630471263的博客
07-26 1225
今天我们继续学习Matlab中函数相关知识。help 命令是最基本的查询方法，可查询所有目录、指定目录、命令、函数。我们直接点击帮助菜单即可查询所需的API函数。lookfor 关键字用于搜索相关的命令和函数。如，我们输入lookfor input，MATLAB就会显示所有关于input的API函数type M 文件名用于显示文件代码。如果想显示以.m扩展名结尾的就要加通配符 which M 文件名用于显示指定文件的路径。返回一个路径demo 用于演示matlab功能。如图，我们可以获取各种演示代码来解
使用MATLAB 2022 中的Lidar Camera Calibrator 实现16线LiDAR 与相机的联合标定
最新发布
u014374826的博客
08-02 523
使用MATLAB 2022 中的Lidar Camera Calibrator 实现16线LiDAR 与相机的联合标定
MATLAB学习日志DAY23
T3ARWIPE的博客
07-29 480
在这种情况下，randi 的第一个输入描述整数可能值的范围，后面两个输入描述行和列的数量。列向量为 m×1 矩阵，行向量为 1×n 矩阵，标量为 1×1 矩阵。要手动定义矩阵，使用方括号 [ ] 来表示 数组的开始和结束。有关创建和处理矩阵的详细信息，参阅“创建、串联和扩展矩阵”。MATLAB 提供了许多函数，用于创建各种类型的矩阵。在标量（1×1 矩阵）的情况下，括号不是必需的。在括号内，使用分号;
MATLAB（10）分类算法
m0_73399576的博客
08-01 218
MATLAB——分类算法篇
数值分析——埃尔米特（Hermit）插值
qq_40619699的博客
08-01 969
在之前提到的插值方法——线性插值、抛物线插值、拉格朗日插值、牛顿插值中，插值函数需要满足的条件都是在插值点与原函数数值相同。在部分实际问题中要求在插值点的一阶导数甚至高阶导数也相等，满足这种条件的插值函数就是今天的主角——埃尔米特（Hermit）插值。
MATLAB画散点密度图(附代码和测试数据的压缩包)
WZZHHH的博客
07-29 410
怀俄明探空站数据解算PWV和Tm：怀俄明探空站数据解算PWV和Tm怀俄明多线程下载探空站数据（包括检查和下载遗漏数据的代码）：怀俄明多线程下载对IGRAv2进行质量控制得到PWV和Tm的 matlab 代码：IGRAv2进行质量控制得到PWV和Tm算 IGRAv2 探空站的 Tm 和 PWV：IGRAv2计算Tm和PWV提取探空站 IGRAv2 全部数据：ERA5 解算合集（温度、气压、PWV、水汽压和 Tm）代码获取：ERA5合集。
.mat文件读取和输入matlab
07-27
您可以使用以下代码将.mat文件读取到MATLAB中： ```matlab data = load('filename.mat'); ``` 其中，'filename.mat'是您要读取的.mat文件的文件名。读取后，变量"data"将包含.mat文件中的所有数据。 要将数据输入到MATLAB中，您可以使用以下方法之一： 1. 使用命令行输入：在MATLAB命令行中直接输入您的数据。例如，要输入一个向量，您可以键入： ```matlab vec = [1, 2, 3, 4, 5]; ``` 2. 使用脚本或函数：将数据保存在一个脚本或函数文件中，并在MATLAB中运行该文件。例如，创建一个名为"my_script.m"的脚本文件，并在其中定义向量： ```matlab vec = [1, 2, 3, 4, 5]; ``` 然后，在MATLAB命令行中运行该脚本文件： ```matlab run my_script.m ``` 3. 使用MATLAB GUI：如果您使用的是MATLAB的图形用户界面（GUI）版本，您可以使用界面中的编辑器来输入和编辑您的数据。 希望这些信息对您有所帮助！如果您有任何其他问题，请随时提问。
评论 13
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值