找到相应的文件

findfile

%1 打开文件夹
%2 拷贝第一个文件
%3 关闭当前文件，再次拷贝新的文件，直到文件末尾结束
clc
clear
DST_PATH_t = 'C:\Users\Manual\Desktop\航次11\1234567890';%目的文件目录
directory = uigetdir('','选择指定文件夹：');
day_char = '0730'; %时间 天
hour_char = '00'; %时 00 - 23
minute_char = '00'; % 取值范围 00 - 60
num_char = '00'; %文件的序号 取值范围 00 - 60
k = 0;
hour_k = 0;
minute_k = 0;
num_k = 0;
sum_file = 0;
for hour_k = 0:9 %时hour_char = num2str(hour_k);for minute_k = 0 : 60 %分minute_char = num2str(minute_k);for num_k = 0 :60 %秒num_char = num2str(num_k);if minute_k < 10 %分if num_k < 10  %秒filepath = fullfile(directory,[day_char '0' hour_char '0' minute_char '0' num_char])endif num_k >= 10filepath = fullfile(directory,[day_char '0' hour_char '0' minute_char num_char])endend  %end if minute_k < 10 if minute_k >= 10if num_k < 10  %秒filepath = fullfile(directory,[day_char '0' hour_char  minute_char '0' num_char]) % 08 01 00 11 09endif  num_k >= 10filepath = fullfile(directory,[day_char '0' hour_char  minute_char num_char])endend %end if if minute_k >= 10if exist(filepath,'dir') == 7
%                  fprintf('%d',exist('H:\2307-8.31-10.2\RadarFile\0924002707','dir'))filepathreadfiles = dir(filepath);filename = readfiles(3).name;
%                     while 1
%                     endcd(filepath);copyfile(filename,DST_PATH_t);sum_file = sum_file + 1;break;endendend
endfor hour_k = 10:23 %时hour_char = num2str(hour_k);for minute_k = 0 : 60 % 分minute_char = num2str(minute_k);for num_k = 0 : 60 %秒num_char = num2str(num_k);
%                 filepath = fullfile(directory,[day_char hour_char minute_char num_char]);if minute_k < 10 %分if num_k < 10  %秒filepath = fullfile(directory,[day_char hour_char '0' minute_char '0' num_char])endif num_k >= 10filepath = fullfile(directory,[day_char hour_char '0' minute_char num_char])endend  %end if minute_k < 10 if minute_k >= 10if num_k < 10  %秒filepath = fullfile(directory,[day_char hour_char  minute_char '0' num_char]) % 08 01 00 11 09endif  num_k >= 10filepath = fullfile(directory,[day_char hour_char  minute_char num_char])endend %end if if minute_k >= 10if exist(filepath,'dir') == 7 %路径存在
%                     fprintf('%d',exist(filepath,'dir'))filepathreadfiles = dir(filepath);filename = readfiles(3).name;cd(filepath);copyfile(filename,DST_PATH_t);sum_file = sum_file + 1;break;endendend
end
fprintf('当前未查找到的文件个数 %d 个\n',24 - sum_file)

筛选出响应时间的谱

%1 根据大概的文件时间，找到近似的文件
%2 将文件保存到指定文件夹下
% addpath('H:\2307_20180822\D\RadarFile');
clc;
clear;
k = 9;
char = '53'
DST_PATH_t = 'C:\Users\Manual\Desktop\航次11\pu\2018_07_30';%目的文件目录
directory = uigetdir('','选择指定文件夹：');
for i = 0:2 %共24个文件夹for j = 0 : 9filepath = fullfile(directory,['CalcReslut' num2str(i,'%d') num2str(j,'%d') '_' char '.dat'])%判断文件是否存在
%         while k < 10
%             if exist(filepath,'file') == 2 %如果文件存在
%                 break;
%             end
%             filepath = fullfile(directory,['CalcReslut' num2str(i,'%d') num2str(j,'%d') '_' '0' num2str(k,'%d') '.dat'])
%             k = k + 1;
%         end
%         k = 9;while k < 60if exist(filepath,'file') == 2break;endfilepath = fullfile(directory,['CalcReslut' num2str(i,'%d') num2str(j,'%d') '_' char '.dat'])k = str2num(char);k = k + 1;char = num2str(k);endchar = '53';k = str2num(char);if ((i == 2) && (j == 4))break;end if exist(filepath,'file') == 2 copyfile(filepath,DST_PATH_t);endend
end

画雷达图像，类似600点程序

%我有一个目标：根据数据绘制雷达图像
%1°要搞清楚数据协议
%2°作图的原理
%当前程序的功能：
%1°了解绘图的机制
%2°了解我不知道的关于数据的东西clear
clc
close all
cd('C:\Users\Manual\Desktop\航次11')
%month=3;
day1=1;
day2=31;%取5月2号到5月4号的时间
mode=0;%只绘中脉冲
% DiffNum=250;%间隔100文件采样
DiffNum=1;%200个文件采样，大约10分钟一个图
mkdir('Pic3');
path_save=[cd '\Pic3\'];[filename_ogr] = uigetdir;%选择数据文件位置
FolderInformation=dir(filename_ogr);
% FolderInformation=dir('RadarFile');%当前文件夹信息
FolderNames={FolderInformation.name}';%当前文件夹下所有名字
FolderIndex=[];
for i=1:length(FolderNames)if(length(FolderNames{i,1})==10)FolderIndex=[FolderIndex i];%记录雷达文件序号；每次运行列的维度加1，等价于有逗号，如果用分号则是行的维度加1end
end
FolderNum=length(FolderIndex);%文件夹数目
for j=1:FolderNumfileFolder=[filename_ogr '\' FolderNames{FolderIndex(j),1}];%一种拼接路径的方式dirOutput=dir(fullfile(fileFolder,'*.dat'));%dir函数可以返回一个类似于结构体的东西，其中包含name项fileNames={dirOutput.name}';fileNum=length(fileNames);%文件夹下文件数目filePath=cell(fileNum,1);%%    i=0;for kk=1:DiffNum:fileNumi=i+1;%i+=1filePath=[fileFolder '\' fileNames{kk,1}];%当前文件路径文件夹路径+.dat文件名Name=fileNames{kk,1};%mm=str2num(Name(9:10));%月dd=str2num(Name(11:12));%日%if(mm==month)if(dd>=day1&&dd<=day2)SavePic(filePath,Name,mode,path_save);end%end                     end  disp(['完成' num2str(j/FolderNum*100) '%']);
end
disp(['Finished！']);

保存图像

function SavePic( pathname,filename,mode,path_save)
%UNTITLED2 此处显示有关此函数的摘要
%   此处显示详细说明R_num=600;%600R_num1=400;dr=7.5;%drfid_0=fopen(pathname,'r','l');%小端模式fseek(fid_0,10312,'bof');linenum=fread(fid_0,1,'ubit16');fclose(fid_0);if(linenum>7000)Mode=0;elseif(linenum>5000&&linenum<=7000)Mode=1;elseMode=2;end
%     if(mode~=Mode)
%         return
%     endfid_0=fopen(pathname,'r','ieee-be');%上文是小端模式，这里却又是大端，是否有问题？fseek(fid_0,18,'bof');%%_pg暂时不明原因shipDir=fread(fid_0,1,'uint32')*180/2^30;%%_pg船艏角shipZong=fread(fid_0,1,'uint32')*90/2^30;%纵摇角shipHeng=fread(fid_0,1,'uint32')*90/2^30;%横摇角fseek(fid_0,3020,'bof');%%lat=fread(fid_0,1,'int32')*90/2^30;%4个字节，32位 %%_pg经度long=fread(fid_0,1,'int32')*90/2^30;             %%_pg纬度seaDeep=fread(fid_0,1,'uint32')*0.01;%%%%     水深为0时，赋值50%vel2flow=fread(fid_0,1,'int16')*0.01*1852/3600;%%_相对航速vel2flow=fread(fid_0,1,'int16')*0.01;%VC程序上是如此算，航速velReal=fread(fid_0,1,'int16')*0.01;% wamos的数据这一位有为-9的时候，所以由 uint16 改成 int16 2011-11-14%%绝对航速velEast=fread(fid_0,1,'int16')*0.01;%%东向速度velNorth=fread(fid_0,1,'int16')*0.01;%%北向速度velWinDir=fread(fid_0,1,'ubit16')*180/2^14;%%_pg平均真风向_和协议上有些不一样velWind=fread(fid_0,1,'ubit16')*0.01;%%平均真风速fclose(fid_0);fid_0=fopen(pathname,'r');status = fseek(fid_0,3114,'bof');temp_data=fread(fid_0,[1,3601],'ubit16');%%角度信息_pgn_line=temp_data(1);%%总线数data_all=fread(fid_0,'ubit16');%%图像数据_pgdata_all=bitand(data_all,16383);  %%%128*128 = 16384  %%_pg应该是为了使大于2000H的数变为2000Hfclose(fid_0);%% 数据处理Gray=reshape(data_all,[R_num,n_line]);Gray_temp=Gray(1:400,:);if seaDeep==0seaDeep=50;endif shipHeng>180shipHeng=shipHeng-360;elseif shipHeng<-180shipHeng=shipHeng+360;endif shipZong>180shipZong=shipZong-360;elseif shipZong<-180shipZong=shipZong+360;endshipHeng=roundn(shipHeng,-2); %保留小数点后两位shipZong=roundn(shipZong,-2);str_sH=num2str(shipHeng);str_sZ=num2str(shipZong);wind_angle=360-shipDir+velWinDir;%风向角point1_x=R_num1*dr*sind(wind_angle);%风向标识起点point1_y=R_num1*dr*cosd(wind_angle);point2_x=R_num1*dr*1.1*sind(wind_angle);%风向标识终点point2_y=R_num1*dr*1.1*cosd(wind_angle);wind_x=[point1_x point2_x];wind_y=[point1_y point2_y];wind_length=((point1_x-point2_x)^2+(point1_y-point2_y)^2)^0.5;jiantou_angle=40;%风向标识箭头的角度point3_xt=wind_length*0.5*sind(wind_angle+jiantou_angle);point3_yt=wind_length*0.5*cosd(wind_angle+jiantou_angle);point4_xt=wind_length*0.5*sind(wind_angle-jiantou_angle);point4_yt=wind_length*0.5*cosd(wind_angle-jiantou_angle);point3_x=point1_x+point3_xt;%风向标识的箭头坐标point3_y=point1_y+point3_yt;point4_x=point1_x+point4_xt;point4_y=point1_y+point4_yt;wind_x1=[point3_x point1_x];%为了使用plot连线wind_y1=[point3_y point1_y];wind_x2=[point4_x point1_x];wind_y2=[point4_y point1_y];%% 画图配置alpha=linspace(0,2*pi,n_line);%极坐标，角r=linspace(0,400*dr,400);%极坐标，半径 %.600[alpha,r]=meshgrid(alpha,r);x=r.*cos(alpha);y=r.*sin(alpha);%% 为了使输出图像好看set(gcf, 'PaperPositionMode', 'manual');set(gcf, 'PaperUnits', 'inches');%     set(gcf, 'PaperPosition', [5 5 4.5 4]);set(gcf, 'PaperPosition', [5 5 5.5 4]);%% 画图和输出mesh(y(20:400,:),x(20:400,:),Gray_temp(20:400,:));%.600view(0,90);colorbarhold on
%     plot(wind_x,wind_y,'r','LineWidth',2);    %风向标识
%     plot(wind_x1,wind_y1,'r','LineWidth',2);
%     plot(wind_x2,wind_y2,'r','LineWidth',2);text(-(R_num1*(dr-1900/600+1900/600)),R_num1*(dr+500/600),['经度:',num2str(lat),'°E'],'Color','k');text(-(R_num1*(dr-1900/600+1900/600)),R_num1*(dr+100/600),['纬度:',num2str(long),'°N'],'Color','k');text(R_num1*(dr-1400/600),R_num1*(dr+500/600),['航速:',num2str(vel2flow),'节'],'Color','k');text(R_num1*(dr-1400/600),R_num1*(dr+100/600),['航向:',num2str(roundn(shipDir,-2)),'°'],'Color','k');text(R_num1*(dr-1400/600),R_num1*(dr-300/600),['风速:',num2str(velWind),'m/s'],'Color','k');text(R_num1*(dr-1400/600),R_num1*(dr-800/600),['风向:',num2str(roundn(velWinDir,-2)),'°'],'Color','k');%   %dr-2500/600 dr-2200/600 dr-1900/600text(R_num1*(dr-1400/600),-(R_num1*(dr-900/600)),['东速:',num2str(velEast),'节'],'Color','k');text(R_num1*(dr-1400/600),-(R_num1*(dr-500/600)),['北速:',num2str(velNorth),'节'],'Color','k');text(R_num1*(dr-1400/600),-(R_num1*(dr-100/600)),['地速:',num2str(velReal),'节'],'Color','k');%1900 2200 2500text(-(R_num1*(dr-1900/600+1900/600)),-(R_num1*(dr-900/600)),['纵摇:',str_sZ,'°'],'Color','k');text(-(R_num1*(dr-1900/600+1900/600)),-(R_num1*(dr-500/600)),['横摇:',str_sH,'°'],'Color','k');text(-(R_num1*(dr-1900/600+1900/600)),-(R_num1*(dr-100/600)),['海深:',num2str(seaDeep),'m'],'Color','k');%1900/600+2200/600  2200/600+2200/600  2500/600+2200/600name=filename(1:end-4);%纯粹的文件名，不带任何后缀的name_txt=name;name_txt(find(name(1,:)=='_'))='-';%text中‘_’有别的含义，在此替换为‘-’text(-(R_num1*(dr-5400/600+2100/600)),R_num1*(dr+100/600),name_txt,'Color','r');hold offprint(1, '-dpng',[path_save name]);end

绘制一维谱和二维谱

%%
%把设备输出谱文件放入pu文件夹，运行即可
%Date：2017/12/5
%Author：pg
%%
clear
clc
cd('C:\Users\Manual\Desktop\航次11')
FolderInformation=[cd '\pu'];
FolderNames=dir(FolderInformation);%当前文件夹下所有名字
mkdir('OneDimension');
mkdir('TwoDimension');
mkdir('result');useless1=0;
for DayNum=1:length(FolderNames)if(isempty(regexpi(FolderNames(DayNum).name,'^[0-9]{4}_[0-9]{2}_[0-9]{2}$')))continue;endToday=FolderNames(DayNum).name;%当前日期    fileFolder=[FolderInformation '\' Today];dirOutput=dir(fullfile(fileFolder,'*.dat'));fileNames={dirOutput.name}';fileNum=length(fileNames);%文件夹下文件数目fid=fopen([cd '\result\' num2str(str2num(Today(6:7))) num2str(Today(9:10)) '.txt'],'w');%写结果文件for kk=1:fileNumfilename_data=[fileFolder '\' fileNames{kk,1}];%     filename_data='Result2017_03_18_829.dat';%信息处理计算机生成的历史数据信息%% 数据读取TodayTime=fileNames{kk,1}(end-8:end-4);%end表示的是最后一个，等价于end-0；ShowTime=[Today  '_' TodayTime];ShowTime(find(ShowTime=='_'))='-';fid_0=fopen(filename_data,'r','ieee-be');%读取方式应当使用这一种status = fseek(fid_0,4,'bof');temp=fread(fid_0,6,'ubit16');WaveH=temp(2)*0.1;%波高WaveD=temp(3)*2^(-14)*180;%波向WaveT=temp(4)*0.1;%波周期FlowV=temp(5)*0.1;%流速，m/sFlowD=temp(6)*2^(-14)*180;%流向temp=fread(fid_0,8,'ubit8');%.可能已经没用了
%         time=[num2str(temp(1)) num2str(temp(2)) '.' num2str(temp(3)) '.' num2str(temp(4)) '-' num2str(temp(6)) ':' num2str(temp(7)) ':' num2str(temp(8))];%时间long=fread(fid_0,1,'int32')*90/2^30;%经度fat=fread(fid_0,1,'int32')*90/2^30;%纬度fclose(fid_0);%% 一维频谱fid_0=fopen(filename_data,'r');%这里的读取方式应当使用这一种而非"ieee-be"status = fseek(fid_0,32,'bof');Num=fread(fid_0,1,'ubit16');%频率个数OneDimensional_Frequency=zeros(1,Num);OneDimensional_Spec=zeros(1,Num);for i=1:NumOneDimensional_Frequency(i)=fread(fid_0,1,'int16')*0.001;OneDimensional_Spec(i)=fread(fid_0,1,'int16')/10000;endfigure(1)title(filename_data(6:end-4));plot(OneDimensional_Frequency,OneDimensional_Spec,'-*');title([ShowTime]);%% 二维频率谱Num2=fread(fid_0,1,'ubit16');%频率个数\% Num2=121;TwoDimensional_FrequencyX=zeros(1,Num2^2);TwoDimensional_FrequencyY=zeros(1,Num2^2);TwoDimensional_Spec=zeros(1,Num2^2);for i=1:Num2^2TwoDimensional_FrequencyX(i)=fread(fid_0,1,'int16')*0.001;TwoDimensional_FrequencyY(i)=fread(fid_0,1,'int16')*0.001;TwoDimensional_Spec(i)=fread(fid_0,1,'int16')/40000;end%% 信噪比、模式读取mode=fread(fid_0,1,'ubit8');Sn=fread(fid_0,1,'ubit16')*0.001;%%figure(2)fclose(fid_0);%% 二维频谱绘制处理fx_neg=reshape(TwoDimensional_FrequencyX,Num2,Num2);fy_neg=reshape(TwoDimensional_FrequencyY,Num2,Num2);sp30_neg=reshape(TwoDimensional_Spec,Num2,Num2);v0=max(max(sp30_neg));v=v0*.1:v0*0.1:v0;C=contour(fx_neg,fy_neg,sp30_neg,v);%_pg 增加一个输出，为了输出谱文件 2016.1.20shading interp;hold onangle = linspace(0,2*pi,100);r = [12,10,8,6,4];       %% 10-05-13 peach 画图显示频率R = 1./r;for ti=1:length(R)px = R(ti)*cos(angle);py = R(ti)*sin(angle);plot(px,py,'-');endR0=0.25;plot([-R0,R0],[0,0],':');plot([0,0],[-R0,R0],':');axis([-R0 R0 -R0 R0]);set(gca,'XTick',-R0:0.05:R0);set(gca,'YTick',-R0:0.05:R0);axis square%%%%%%%%%%  %% 10-05-13 peach 画图显示频率text( -R(5)/1.4,R(5)/1.4,'4 s','color','k','Fontsize',10,'fontweight','bold');text( -R(4)/1.4,R(4)/1.4,'6 s','color','k','Fontsize',10,'fontweight','bold');text( -R(3)/1.4,R(3)/1.4,'8 s','color','k','Fontsize',10,'fontweight','bold');text( -R(2)/1.4,R(2)/1.4,'10 s','color','k','Fontsize',10,'fontweight','bold');text( -R(1)/1.4,R(1)/1.4,'12 s','color','k','Fontsize',10,'fontweight','bold');newangle=0:(2*pi/12):2*pi;myx1=-R0:0.01:R0;myy1=-R0:0.01:R0;myx=sqrt(myx1.^2+ myy1.^2);for i=1:length(newangle)plot(myx.*cos(newangle(i)),myx.*sin(newangle(i)),':')hold onendtext(-0.0077, R0*0.97,'North ','color','k','Fontsize',12,'fontweight','bold');xlabel('f_x (1/s)','fontweight','bold');ylabel('f_y (1/s)','fontweight','bold');hold offcolorbartitle(ShowTime);print(1, '-dpng',['OneDimension\' 'RadarOneDimension' ShowTime ]);%'A'print(2, '-dpng',['TwoDimension\' 'RadarTwoDimension' ShowTime ]);%'B'timenow=str2num(TodayTime(1:2))*100+str2num(TodayTime(4:5));fprintf(fid,'%d   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f   %4.2f\n ',timenow,Sn,WaveH,WaveD,WaveT,mode,useless1,useless1,useless1,useless1,useless1,useless1,useless1);       disp([num2str(kk*100*DayNum/fileNum/length(FolderNames)) '%']);clfendfclose(fid);
end
fclose all

删除文件

%删除文件
clc;
clear;
cd('C:\Users\Manual\Desktop\航次11\TwoDimension');
DST_PATH_t = 'C:\Users\Manual\Desktop\航次11\NewTwoDimension';%目的文件目录
filepath = 'C:\Users\Manual\Desktop\航次11\Pic3'; %放置Pic3图像的文件夹
readfiles = dir(filepath);intmintime = 0;
inthourtime = 0;
intdaytime = 0;pic3filenumber = 0;for i = 3: length(readfiles)monthtime = readfiles(i).name(9:10); %月daytime  = readfiles(i).name(11:12); %天hourtime = readfiles(i).name(13:14); %时mintime = readfiles(i).name(15:16); %分intdaytime = str2num(daytime);inthourtime = str2num(hourtime);intmintime = str2num(mintime);if intmintime < 10FileRadarTwoDimension = fullfile(['RadarTwoDimension2018-' monthtime '-' daytime '-' hourtime '-' '0' mintime '.png']);end if intmintime >= 10FileRadarTwoDimension = fullfile(['RadarTwoDimension2018-' monthtime '-' daytime '-' hourtime '-' mintime '.png']);endwhile intmintime > 0if exist(FileRadarTwoDimension,'file') == 2%执行拷贝等功能FileRadarTwoDimensioncopyfile(FileRadarTwoDimension,DST_PATH_t);pic3filenumber = pic3filenumber + 1;break;endintmintime = str2num(mintime);intmintime = intmintime - 1;mintime = num2str(intmintime);if intmintime < 10FileRadarTwoDimension = fullfile(['RadarTwoDimension2018-' monthtime '-' daytime '-' hourtime '-' '0' mintime '.png']);endendendfprintf('Pic3图像数为 %d\n',pic3filenumber);

参考资料
1 https://ww2.mathworks.cn/help/matlab/matlab_prog/debugging-process-and-features.html 调试 MATLAB 程序
2 https://blog.csdn.net/s151506879/article/details/52074559 matlab-fullfile使用
3 https://blog.csdn.net/qq_18995069/article/details/80324588 【Matlab编程实现常见小问题之一】Matlab如何读取指定文件夹下的指定类型文件
4 https://blog.csdn.net/yeyang911/article/details/78119779 matlab 移动（复制）文件到另一个文件夹
5 https://blog.csdn.net/qikaihuting/article/details/70244416 matlab中dir函数以及sort函数的用法
6 https://www.ilovematlab.cn/thread-536023-1-1.html 如何依次循环读取多个文件夹下的所有txt文件