每行代码进行解释 找出报错原因

Question

clear all
path1='d:\temp\kerncountyX\';

filename=dir([path1,'*']);
filename=struct2cell(filename);
filename=char(filename(1,3:length(filename)));%5,length(filename)
filename=char(filename);
N=int16(size(filename,1));

ss=1;
count=1;

for k=1:N

    
    file_current=[path1,filename(k,:)];  %currently direct file 
    fid1=fopen(file_current);
    
        for j=1:4
        fline=fgets(fid1);
        switch j
            
            case 2
                 Dt=str2double(fline(41:44));%震源深度
           
     
        end
    end
        data=fscanf(fid1,'%f',inf);
      
        fclose(fid1);
        eval(['Rec.data' num2str(k) '=data'])    %把data数据记录到一个结构体Rec的dataK变量中,调用这个记录就直接用Rec.dataK
       
      
end


fileds = fieldnames(Rec);
k1=[1 2 3 4];
    for i=1:length(k1)
    i1=k1(i);
     kk1 = fileds(i1);
      key = kk1{1};
R1=Rec.(key);
r11=length(R1)/28;
R1=reshape(R1,28,r11);
R1=R1';
for nn=2:7
    RL=abs((R1(:,nn+1)-R1(:,nn)));
   Res(nn)=max(RL); 
end
h=[0 5.5 4.8 4.7 4.8 3.7 5.5];%height
%计算层间位移角
Cj(i,1:7)=Res./h;
    end
%求塔顶最大位移
    for i=1:N
    i1=k1(i);
     kk1 = fileds(i1);
      key = kk1{1};
R1=Rec.(key);
r11=length(R1)/8;
R1=reshape(R1,8,r11);
R1=R1';
t=R1(:,1);
Dis(i,1:r11)=R1(:,2)-R1(:,8);
    end
 Dmax=max(Dis);   
  

Answer 1

你好，我是有问必答小助手，本次您提出的有问必答问题，已经由其他小伙伴为您做出解答，并解决了您的问题。

本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。