【运动学】基于matlab计步【含Matlab源码 524期】

举报
海神之光 发表于 2022/05/29 01:47:05 2022/05/29
【摘要】 一、获取代码方式 获取代码方式1: 完整代码已上传我的资源:【运动学】基于matlab计步【含Matlab源码 524期】 获取代码方式2: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得...

一、获取代码方式

获取代码方式1:
完整代码已上传我的资源:【运动学】基于matlab计步【含Matlab源码 524期】

获取代码方式2:
通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码。

备注:
订阅紫极神光博客付费专栏,可免费获得1份代码(有效期为订阅日起,三天内有效);

二、部分源代码

%对原始观测得到的数据进行获取
clc, clear;
%加速度数组
accspe_data = [];
time_record = [];
fid=fopen('sensor_normal.txt', 'r', 'n', 'utf-8');
%初始赋值使循环开始
tline = 1;  
%对文本进行循环读入
%line = 'owieo';
while tline
    %fgetl获取某一行
    tline=fgetl(fid);
    %转换编码形式
    tline = native2unicode(tline);
    %enc = detect_encoding(tline)
    %转换为char类型
    line = char(tline); 
    %对值进行输出
    %disp(line);
    %提取的数据
   % data = regexp(line, ':', 'split');
    %disp(line);
    %disp(length(line));
%     res = '';
% %     for len = 1:length(line)
% %         res = strcat(res,line(len));         %字符串相加
% %         %disp('oiwei');
%     end
    
    data = strsplit(line, ':');
    %从字符串中提取所需要的数据
    if strcmp(data(1), '加速度')
        volecity = data(2);
        vole = strsplit(char(volecity));
        %disp(vole);
        vole_x = str2double(vole(2));
        vole_y = str2double(vole(3));
        vole_z = str2double(vole(4));
        vole_res = sqrt(power(vole_x, 2) + power(vole_y, 2) + power(vole_z, 2));
        accspe_data = [accspe_data, vole_res];
    end
    
    if strcmp(data(1), '时间')
        time = strcat(data(2), data(3), data(4));
        time_data = strsplit(char(time));
        time_record = [time_record, (time_data(3))];
    end
   % end
    %print('oiweoi')
    %print(class(tline));
    %提取其中的数据 
end
%判断是否缺少数据
len = length(time_record);
time_check = ones(1, len);
for num = 1:len
    time_num = str2double(char(time_record(num)));
    
    hour = floor(time_num / 10000);
    minute = floor((time_num - hour * 10000) / 100);
    second = floor((time_num - hour * 10000 - minute * 100));
    
    %将时间转换为秒
    time_check(num) = hour * 3600 + minute * 60 + second;
    if num >= 2 && (time_check(num) - time_check(num - 1)) > 1
        disp('数据采集有误,部分数据缺失!!!');
        disp(time_check(num) - time_check(num - 1)) ;
    end
end
fclose(fid);

  
 
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71

三、运行结果

在这里插入图片描述

在这里插入图片描述

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社,2013.

文章来源: qq912100926.blog.csdn.net,作者:海神之光,版权归原作者所有,如需转载,请联系作者。

原文链接:qq912100926.blog.csdn.net/article/details/114853973

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。