常用matlab脚本

%% 采集数据出现的次数sequencedistance1=sort(Data1); %对距离从小到大排列起来,排下序，方便后面获取每个数的位置, B=sort(A) 对一维或二维数组进行升序排序c = unique(sequencedistance1); %% 得到sequencedistance1的不重复数列，并获取每个不重复数的位置count = hist(sequencedistance1,unique(sequencedistance1)); %得到每个元素重复的次数；figure(1);stem(c,count,'MarkerFaceColor',[ 1 0 0 ]);%% 找出一行中全是某值的数据，仅有这些行的数值加起来为某个值for i=1:1500 if data2(i,1)+data2(i,2)+data2(i,3)+data2(i,4)+data2(i,5)+data2(i,6)+data2(i,7)+data2(i,8)==8*50; endend%% 将16进制N*8的矩阵转化为10进制N*8的矩阵data=textread('E:\OriginalSignal.txt','%s')';%以字符形式打开文件,打开之后就为一行或者一列data1=hex2dec(data)'; %16进制转化为10进制数，存入data1矩阵data2=reshape(data1',8,1500)';%换成1500*8的矩阵%% 将10进制1500*8矩阵转化为16进制以字符串格式呈现的1500*8矩阵，导入到txt格式后可以空格去掉引号，得到我们想要的txt格式%方法：先将10进制的数通过dec2hex转化为元胞矩阵，在通过元胞矩阵转化为字符串）data3=reshape(data2',8,1500);data4=dec2hex(data3);data5=cellstr(data4);data6=reshape(data5,8,1500)';%% 查找矩阵中满足条件的某些值对应的横纵坐标yt3_100index=find(yt3>51);%选取大于51的数据的索引号yt3_100=yt3(yt3_100index);%选取大于51的数据值中%% 查找一系列数的峰值以及峰值对应的横纵坐标[peakycoordiante,peakxcoordiante]=findpeaks(data);%% 行为（134,0x000005e)文件格式文件的读取fp_Read = fopen('E:\data\ljfxy\21.549.csv','r'); %读取数据文件，修改改文件即可 index = 1; while ~feof(fp_Read) %判断是不是文件末尾 tline=fgetl(fp_Read); %读取一行，注意，读文一行后，光标就会自动到下一行 Length = length(tline); s = tline(1,Length-7:Length); %去除前面的0X部分 Data(1,index)=hex2dec(s); %把十六进制换成乘10进制 index = index + 1; end%% 把一个矩阵直接变为1行或1列reshapesamplingdata=a(:)'; %% 去除某值所在的行[row,column]=find(data2==306); data2(row,:)=[];%% 查找峰值（可查找峰值大于60的峰值） [peaksy,peaksx]=findpeaks(data4,'minpeakheight',60);%% 输入一行或一列16进制的数（AF 4D 36 5A..........................）data=textread('E:\data\gao_data_case_5.28_12\data_d21.086m_t12.txt','%s')';%以字符形式打开文件 data1=hex2dec(data)'; %16进制转化为10进制数，存入data1矩阵,并转化为1行%% 寻找函数的波峰与波谷值 IndMin=find(diff(sign(diff(data)))>0)+1;IndMax=find(diff(sign(diff(data)))<0)+1;%其中的，IndMin, data(IndMin)对应的是波谷点的数据,IndMax,data(IndMax)对应的是波峰点的数据%% 通过增加最小间隔条件可以过滤掉单周期内的干扰极值 [maxv,maxl] = findpeaks(data,'minpeakdistance',mpd); %mpd 设定两峰值间的最小间隔数%% 同时读入多个文件的办法A = dir(fullfile('C:\Users\admin\Desktop','*.txt'));% 读取后A的格式为 % name -- filename% date -- modification date% bytes -- number of bytes allocated to the file% isdir -- 1 if name is a directory and 0 if notA = struct2cell(A);num = size(A);for k =0:num(2)-1 x(k+1) = A(5*k+1);% 找出name序列endfor k = 1:num(2) newpath = strcat('C:\Users\admin\Desktop','\',x(k)); data{k} = load(char(newpath));end%% 求矩阵元素减去某个值后的最小值所在的位置[minv num]=min(abs(seljitt(:,4)-ug)+abs(seljitt(:,5)-dg));%% 2个byte的存储空间作为存储两个变量的值，其中前一个byte的高四位存储一个变量，前一个byte的低四位和后一个byte所有位存储另一个变量var1=bitshift(bitand(high4,240),-4);var2index=bitand(low4,15);var2=low4*256+all;%% 2个byte的存储空间作为存储两个变量的值，前一个byte所有位的和后一个byte高四位存储一个变量，后一个byte的低四位存储一个变量angleIndex(cnt)=bitshift(bitand(b(cnt),240),-4); %取后一个byte高四位存储的值angle(cnt)=a(cnt)*16+angleIndex(cnt); %前一个byte的所有位+后一个byte的高四位为其中一个变量值channel(cnt)=bitand(b(cnt),15); %后一个byte的后四位为其中另一个变量值；