Matlab读取txt数据的实用方法[通俗易懂]

Matlab读取txt数据的实用方法[通俗易懂]需求有个朋友需要我帮忙写个matlab脚本读取100个txt文档的实验数据,这些文档的结构相同,分为四列,从第一列到第四列依次是时间、位置、速度、加速度。读取完数据之后需要对数据进行处理,具体的处理方式是:提取以0.002为采样周期的数据,分类存储起来。文件内容是这样的:技术难点技术难点在于,这些文件中的数据是从一个软件中仿真得到的,由于采用的是变步长仿真,因此采样时间不统一,很难采用对…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

需求

有个朋友需要我帮忙写个matlab脚本读取100个txt文档的实验数据,这些文档的结构相同,分为四列,从第一列到第四列依次是时间、位置、速度、加速度。读取完数据之后需要对数据进行处理,具体的处理方式是:提取以0.002为采样周期的数据,分类存储起来。
文件内容是这样的:
在这里插入图片描述

技术难点

技术难点在于,这些文件中的数据是从一个软件中仿真得到的,由于采用的是变步长仿真,因此采样时间不统一,很难采用对时间取余的操作来实现。

解决办法

首先对于给定的文档结构,采用textread函数读取四列数据分别存放在四个变量a1,a2,a3,a4中。小伙伴儿可能会问,由于文件头的存在,读取的数据前面几行并不是实际数据,怎么办?好办,把前几行去掉不就行了,使用MATLAB很容易提取。

[a1,a2,a3,a4]=textread(filename,'%s%s%s%s','headerlines',4);

由于读入到matlab中的数据属于文本变量,对应的类型是cell,而我们需要的数据类型是double,怎么办?类型转换吧;处理过程:cell->mat->double

function res=cell2double(input)
[n,m]=size(input);
res=zeros(n,1);
for i=1:n
   temp=cell2mat(input(i));
   res(i)=str2double(temp);
end
end

重点来了,由于采样时间不是固定的(变步长采样时间有些地方大,有些地方小),因此采用

mod(time(i),0.002)==0

这样的简单判断是不行的,本人试过这种方式来提取数据,结果有的文件提取的数据多,有的文件提取的数据少,无法跟时间进行对齐,确实很让人头疼。
考虑到采样时间精度在0.002秒以下,如果只是采用时间除以0.002取余数,很多数据会被错误判定为不合适,从而丢失数据。怎么办呢?
对时间四舍五入嘛,将时间四舍五入到小数点后3位不就行了。
但是呢,问题又来了,比如采样时间0.0082,0.0084时,则这两个数据都会被加入,这是行不通的,处理的结果依然跟前面的一样。
从四舍五入后的数据中寻找特征吧。对于四舍五入后的时间点,当存在相同的采样时间是只取其中的一个即可,仔细观察可以发现这样的数据特征,每一个满足要求的数据,它的后面一个数据的时间节点总是比它大。
代码就这样写吧:

  for j=1:time_rows-5-1
          if mod(time(j),0.002)==0%能被0.002整除,条件1
              if time(j+1)>time(j)%后一个采样时间比前一个大,条件2
                   tempTime(k)=time(j);
                   tempPos(k)=pos(j);
                   tempSpeed(k)=speed(j);
                   tempAcc(k)=acc(j);
                   k=k+1;
              end
          end
      end

好啦,贴一下完整的代码,仅供参考。

参考代码

clc;
clear all;
N=100;
M=5000;
dataToGetPos=zeros(M,N+1);
dataToGetSpeed=zeros(M,N+1);
dataToGetAcc=zeros(M,N+1);
for i=1:N
str1='result';
str2='.txt';
filename=sprintf('%s%d%s',str1,i,str2);
[a1,a2,a3,a4]=textread(filename,'%s%s%s%s','headerlines',4);
[pos_rows,pos_cols]=size(a2);
[time_rows,time_cols]=size(a1);
[speed_rows,speed_cols]=size(a3);
[acc_rows,acc_cols]=size(a4);
tempPos=zeros(M,1);
tempTime=zeros(M,1);
tempSpeed=zeros(M,1);
tempAcc=zeros(M,1);
time=cell2double(a1(5:time_rows));
pos=cell2double(a2(5:pos_rows));
speed=cell2double(a3(5:speed_rows));
acc=cell2double(a4(5:acc_rows));
time=round(time,3);
k=1;
if i==1 
for j=1:time_rows-5-1
if mod(time(j),0.002)==0
if time(j+1)>time(j)
tempTime(k)=time(j);
tempPos(k)=pos(j);
tempSpeed(k)=speed(j);
tempAcc(k)=acc(j);
k=k+1;
end
end
end
[tempPosRows,tempPoscols]=size(tempPos);
if tempPosRows<M
for l=tempPosRows+1:M
tempPos(l)=0;
tempTime(l)=0;
tempSpeed(l)=0;
tempAcc(l)=0;
end
end
dataToGetPos(:,1)=tempTime;
dataToGetPos(:,2)=tempPos;
dataToGetSpeed(:,1)=tempTime;
dataToGetSpeed(:,2)=tempSpeed;
dataToGetAcc(:,1)=tempTime;
dataToGetAcc(:,2)=tempAcc;
else
for j=1:time_rows-5-1
if mod(time(j),0.002)==0
if time(j+1)>time(j)
tempPos(k)=pos(j);
tempTime(k)=time(j);
tempSpeed(k)=speed(j);
tempAcc(k)=acc(j);
k=k+1;
end
end
end
[tempPosRows,tempPoscols]=size(tempPos);
if tempPosRows<M
for l=tempPosRows+1:M
tempPos(l)=0;
tempSpeed(l)=0;
tempAcc(l)=0;
end
end
dataToGetPos(:,i+1)=tempPos;
dataToGetSpeed(:,i+1)=tempSpeed;
dataToGetAcc(:,i+1)=tempAcc;
end
disp(['Processing: ',num2str(i/N*100),'%.......']);
end
save dataToGetPos.mat dataToGetPos;
save dataToGetSpeed.mat dataToGetSpeed;
save dataToGetAcc.mat dataToGetAcc;
figure(1)
plot(dataToGetPos(:,1),dataToGetPos(:,2:101));
xlabel('时间/s');
ylabel('位置/mm');
set(gca,'FontSize',13);
figure(2)
plot(dataToGetSpeed(:,1),dataToGetSpeed(:,2:101))
xlabel('时间/s');
ylabel('速度/mm.s^{-1}');
set(gca,'FontSize',13);
figure(3)
plot(dataToGetAcc(:,1),dataToGetAcc(:,2:101))
xlabel('时间/s');
ylabel('加速度/mm.s^{-2}');
set(gca,'FontSize',13);
disp('Done!');

数据转换函数:

function res=cell2double(input)
[n,m]=size(input);
res=zeros(n,1);
for i=1:n
   temp=cell2mat(input(i));
   res(i)=str2double(temp);
end
end

运行以上代码生成下面的几幅图:
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/196688.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • Android直播开发之旅(10):AndroidUSBCamera,UVCCamera开发通用库(支持开启多路摄像头)[通俗易懂]

    Android直播开发之旅(10):AndroidUSBCamera,UVCCamera开发通用库(支持开启多路摄像头)[通俗易懂]AndroidUSBCamera基于[saki4510t/UVCCamera](https://github.com/saki4510t/UVCCamera)开发,该项目对USBCamera(UVC设备)的使用和视频数据采集进行了高度封装,能够帮助开发者通过几个简单的API实现USBCamera设备的检测、连接、预览和音视频数据采集,最重要的是手机无需root,只需支持otg功能即可驱………

  • 【Tensorflow】Windows安装tensorflow错误原因查询、卸载tensorflow与重新安装

    【Tensorflow】Windows安装tensorflow错误原因查询、卸载tensorflow与重新安装Windows安装tensorflow错误原因查询、卸载tensorflow与重新安装深度学习这个大坑的苦与甜,谁踩谁知道。上文介绍了windows如何安装tensorflow点我:Windows(64)位安装tensorflow,但因为一般电脑上会有好几个版本的python,比如我电脑上装了从2.7各版到3.6等等(具体怎么查询电脑的python版本下面会介绍到),所以在安装的时候一定要确认…

  • C++键盘钩子

    C++键盘钩子C++键盘钩子Windows系统是建立在事件驱动的机制上的,整个系统都是通过消息传递来实现的。而钩子是Windows系统中非常重要的系统接口,用它可以截获并处理送给其他应用程序的消息,来完成普通应用程序难以实现的功能。钩子可以监视系统或进程中的各种事件消息,截获发往目标窗口的消息并进行处理。这样,我们就可以在系统中安装自定义的钩子,监视系统中特定事件的发生,完成特定的功能,比如截获键盘、鼠标的…

  • css opacity属性_CSS中的opacity属性[通俗易懂]

    css opacity属性_CSS中的opacity属性[通俗易懂]cssopacity属性CSS|不透明度属性(CSS|opacityProperty)Withthegrowingneedofmakingwebsites,theneedforstylingthemhasalsoincreased.Therefore,CSShasbecomeanindispensablepartofcreating…

  • robotium android,android自动化测试框架robotium配置和使用

    robotium android,android自动化测试框架robotium配置和使用背景介绍android开发过程中,每次迭代升级都需要去回归一下之前版本功能,看看最新的修改有没有影响到之前的正常功能。然而这个过程永远都是在做一些繁琐的重复的操作,大大浪费人力,所以我们决定使用自动化来做这个事情,这就引入了我们接下来要介绍的自动化测试框架——robotium。Robotium是一款国外的Android自动化测试框架,主要针对Android平台的应用进行黑盒自动化测试,它提供了模拟…

  • SLAM技术概述_SRAM工艺

    SLAM技术概述_SRAM工艺导语随着最近几年机器人、无人机、无人驾驶、VR/AR的火爆,SLAM技术也为大家熟知,被认为是这些领域的关键技术之一。本文对SLAM技术及其发展进行简要介绍,分析视觉SLAM系统的关键问题以及在实际应用中的难点,并对SLAM的未来进行展望。1.SLAM技术SLAM(SimultaneousLocalizationandMapping),同步定位与地图构建,最早在机器人领域提出,…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号