清楚的
[文件,文件夹]=uigetfile('*.xls;*。xlsx;*。“csv”);
filename=fullfile(文件夹,文件);
Nbook=readtable(文件名);
d=方向('*.wcp');
Num_files=高度(d);
名称=单元格(0);
对于nf=1:Num\u文件
Fname=fullfile(d(nf).name);
名称{nf}=Fname;
Names=名称';
结束
对于nn=1:高度(名称)
数据文件=字符(名称(nn));
DataTable=BatchProcess(数据文件,Nbook)
结束
数据=导入wcp(数据文件);
EPSPdata=data.S{1,1};
num_samp=大小(EPSPdata);
间隔=data.t\u间隔;
t=0:间隔:(num\u samp-1)*间隔;
t=t'*1000;
[~,startI]=min(abs(t-13));
[~,endI]=min(abs(t-35));
窗口=开始时间:结束时间;
平滑因子=2;
smData=EPSPdata;
numTrace=大小(EPSPdata,2);
对于n=1:numTrace
smData(:,n,1)=平滑(EPSPdata(:,n,1),平滑因子);
结束
[峰值,~]=min(smData(窗口:);
峰值=峰值';
PeakData=峰值;
CONTROLeq=(Nbook.CADO(1)-1);
CADOeq=(Nbook.CADO\u DPCPX(1)-1);
DPCPXeq=(Nbook.CADO\u DPCPX\u NBQX(1)-1);
NBQXeq=(长度(数据记录索引)-1);
CONTROLData=smData(:,CONTROLeq);
CADOData=smData(:,CADOeq);
DPCPXData=smData(:,DPCPXeq);
NBQXData=smData(:,NBQXeq);
[~,startG]=min(abs(t-7));
[~,endG]=min(abs(t-35));
windowG=startG:endG;
子地块(2,4,1)
控制记录=绘图(t(windowG),控制数据(windowG,:);
ylim([-0.8 0.2]);
xlabel公司'时间(毫秒)'
伊莱贝尔'电压(mv)'
职务(“我”)
子地块(2,4,2)
CADOrecording=绘图(t(windowG),CADOData(windowG,:);
ylim([-0.8 0.2]);
职务(“我”)
子地块(2、4、3)
DPCPXrecording=绘图(t(windowG),DPCPXData(windowG,:);
ylim([-0.8 0.2]);
职务(“三”)
子地块(2,4,4)
NBQXrecording=绘图(t(windowG),NBQXData(windowG,:);
ylim([-0.8 0.2]);
职务(“四”)
PlotT=0:10:10*((长度(数据记录索引)-1));
PlotT=PlotT'/60;
子批次(2,4,[5,6,7,8])
TimeCoursePlot=散点图(plot,PeakData);
xlabel公司'时间(分钟)'
伊莱贝尔'第一个fEPSP振幅(mv)'
xline([PlotT(CONTROLeq)PlotT(CADOeq)PlotT(DPCPXeq)PlotT(NBQXeq)],'-',{“我”,“我”,“三”,“四”},“标签定向”,'水平');
DataTable=table(convertCharsToStrings(DataFile)、PeakData(CONTROLeq)、PeakData(CADOeq)、PeakData(DPCPXeq)、PeakData(NBQXeq),'可变名称', {'文件名','控制峰值','卡多峰','DPCPX峰值','NBQX峰值'});
结束