中国汽车工程师之家--聚集了汽车行业80%专业人士 

论坛口号:知无不言,言无不尽!QQ:542334618 

本站手机访问:直接在浏览器中输入本站域名即可 

  • 1723查看
  • 3回复

[代码] 写了一段代码,实在是找不出哪里错了,请大家帮忙一起看看

[复制链接]

该用户从未签到

发表于 25-7-2016 20:59:09 来自手机 | 显示全部楼层 |阅读模式

汽车零部件采购、销售通信录       填写你的培训需求,我们帮你找      招募汽车专业培训老师


%1.数据准备
newData = importdata('Data160722123758.txt', '\t', 2);
data=newData.data;
sum_length=length(data(:,2));%产生的数据长度都是一样的
%sum_length=length(data(:,3));
%sum_length=length(data(:,5));
%sum_length=length(data(:,6));

%对每列产生的数据进行求和,为求中心点做铺垫
sum_ax=sum(data(:,2));
sum_ay=sum(data(:,3));
sum_anglex=sum(data(:,5));
sum_angley=sum(data(:,6));

%命名新的数组,方便以后画图
ax=data(:,2);
ay=data(:,3);
anglex=data(:,5);
angley=data(:,6);
%加速度的中心点
pointax=sum_ax/sum_length;
pointay=sum_ay/sum_length;
%角度的中心点
pointanglex = sum_anglex/sum_length;
pointangley = sum_angley/sum_length;

%角加速度的极径
p_a=sqrt(ax.^2+ay.^2);

%角速度的极径
p_angle=sqrt(anglex.^2+angley.^2);

%角加速度的极角
%a_angle=array(sum_length);
for i=1:1:sum_length
    a_angle(sum_length);
    if         ax(i)>0 && ay(i)>0
               a_angle(i)=atan(ay(i)/ax(i));
    elseif     ax(i)<0 && ay(i)<0
               a_angle(i)=atan(ay(i)/ax(i))+pi;
    else     
               a_angle(i)=atan(ay(i)/ax(i));
    end
end

%角度的极角
angle_angle(sum_length);
for i=
  • TA的每日心情

    21-3-2016 20:46
  • 签到天数: 1 天

    [LV.1]初来乍到

    发表于 25-8-2016 09:07:52 | 显示全部楼层
    &#30475;&#30475;
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    奋斗
    15-11-2017 08:40
  • 签到天数: 9 天

    [LV.3]偶尔看看II

    发表于 8-11-2017 13:46:00 | 显示全部楼层
    楼上看了就没然后了。。。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    15-3-2019 11:07
  • 签到天数: 143 天

    [LV.7]常住居民III

    发表于 6-3-2018 16:38:59 | 显示全部楼层
    后面的文件呢?怎么没有了。
    回复 支持 反对

    使用道具 举报

    快速发帖

    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    QQ|手机版|小黑屋|Archiver|汽车工程师之家 ( 渝ICP备18012993号-1 )

    GMT+8, 29-3-2024 01:22 , Processed in 0.294832 second(s), 44 queries .

    Powered by Discuz! X3.5

    © 2001-2013 Comsenz Inc.