Index out of bounds error very unusual

1 次查看（过去 30 天）
Nicola
Nicola 2014-2-7
评论： Walter Roberson 2015-8-12
Good evening my friends, i have an unusual error; if my code loads these values from file it gives me error Attempted to access head(3); index out of bounds because numel(head)=1.
If it loads other values like these it works
The txt file is like that
1 1994 1 8 3 59 27 0 15 54 11
8.9 1
9.2 1
10.2 1
12.9 1
11.4 1
10.7 1
12.4 1
8.3 1
7.7 1
2.2 6
2.0 1
etc etc
I hope you have some ideas. Thanks so much for your help!
  3 个评论
显示 1更早的评论
Nicola
Nicola 2014-2-9
The problem is i don't touch anything, with some line of txt it works. Any idea? Because
1 1994 1 8 3 59 27 0 15 54 11
doen't work
but
5 1994 4 8 9 14 40 0 31 16 5
works
Michael
Michael 2015-8-12
I had the same problem. Turns out somehow sum2str got declared as a variable. Subsequently matlab tries to access the 1st entry of the variable (which didn't exist). Hence the error.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（3 个）

Walter Roberson
Walter Roberson 2014-2-12
The error could be further down, such as near the end of file.
Try switching to
thisline = fgetl(head);
if ~ischar(thisline); break; end %end of file
head = sscanf(thisline, '%d', 11);
if length(head) < 11
fprintf('warning: line has too few elements: |%s|', thisline);
continue;
end
  1 个评论
Nicola
Nicola 2014-2-12
编辑：Nicola 2014-2-12
Thanks for your help How can an error at the end affect on the beginning of the code if i use debugger?
The complete code is (freq=32)
function ssttest(freq)
%Whitesands
%elev is the elevation angle (deg)
%site latitude (deg)
lat=32.54%Withesands
%site altitude above mean sea level (m)
Hamsl=1463%Whitesands
%polarization angle (deg)
pola=45;
arr=2;
z=2;
temp=2;
camp=6;
temp2=1;
partenza=1;
n_elev=1;
part=0;
ultimo=0;
fine=0;
cond=0;
evento_pioggia=0;
ffne=0;
d=1;
scrivi=0;
colori=cool(10);
cas=randi(100);
delete('Eventi_anno/*.txt');
year=['1994';'1995';'1996';'1997';'1998'];
res='60';
[ny,ans]=size(year); %ny=5 ans=4
[nr,ans]=size(res); %nr=1 ans=2
pathin='BancaDati/';
pathout='RisSSTtest/';
TtotRain=0;
Ttot=0;
for ir=1:nr,
deltat=str2num(res(ir,:));
for iy=1:ny,
nev=0;
name=[pathin,'r',year(iy,3:4),res(ir,:),'.txt'];
fid=fopen(name,'r');
while ~feof(fid)
head=fscanf(fid,'%d',11);
if ~isempty(head)
giorno=head(3);
mese=head(4); %prende il mese
ora=head(5);
minuti=head(6);
durata_h=head(8);
durata_m=head(9);
m=head(4);
tsR=[];
%vedere altezza sole in base ad ora
nname=['Risultatigiorno/',num2str(giorno),'_',num2str(mese),'.txt'];
nfid=fopen(nname,'r');
while ~feof(nfid)
hhead=fscanf(nfid,'%f',2);
if ~isempty(hhead)
tempo=hhead(1);
oora(arr)=fix(tempo);
mmminuti(arr)=abs(tempo-oora(arr));
mmminuti(arr)=mmminuti(arr)*100;
mminuti(arr)=int32(mmminuti(arr));
if mminuti(arr)==60 %risolve problema minuti .06
mminuti(arr)=mminuti(arr)/10;
end
gggradi(arr)=hhead(2);
ggradi(arr)=fix(gggradi(arr));
% if cond==0 && durata_h>=0 && durata_h<24; %guardo quante ore e minuti dura
% n_inc=ceil((durata_m/camp)+(10*durata_h)); %approssimo per eccesso
% end
if cond==1 && temp<=fine %gradi utili successivi al primo
g_elev(n_elev)=ggradi(arr-1); %forse non mi prende l'ultimo valore
on_elev=n_elev;
n_elev=n_elev+1;
if temp==fine
cond=0; %nel caso mi accorgessi che salta ultimo valore, metto qui un g_elev(n_elev)=ggradi(arr-1); n_elev=n_elev+1;
ultimo=1;
scrivi=1;
end
temp=temp+1;
temp2=temp2+1;
end
% if temp<=fine && n_elev==2 %evita il problema se ho solo 1 campione 1? evento
% cond=0;
% end
%
if ora==oora(arr) && minuti >= mminuti(arr-1) && minuti < mminuti(arr) %inizio evento
g_elev(n_elev)=ggradi(arr-1);
on_elev=n_elev;
n_elev=n_elev+1;
cond=1;
scrivi=0;
one=1;
temp2=1;
evento_pioggia=evento_pioggia+1;
fnamea=['Eventi_anno/',year(iy,3:4),'.txt'];
fidone=fopen(fnamea,'a');
fprintf(fidone,num2str(evento_pioggia));
fprintf(fidone,'\n');
fclose(fidone);
%temp=temp+1; %toglie il problema se ho solo 1 campione
end
if cond==1 && one==1
for ii=1:head(11), %fa scorrere per tutti gli eventi di pioggia?d
aux1=fscanf(fid,'%f',1); %prende la colonna dell'intensit? della pioggia
aux2=fscanf(fid,'%d',1); %prende la colonna del numero di eventi
tsR=[tsR zeros(1,aux2)+aux1]; %crea un array di 1 riga, lungo aux2 e ci mette dentro tante volte il valore di aux1
end
[A,t,deltaTa]=SyntStormEM(tsR,str2num(res),lat,Hamsl,freq,pola,g_elev(on_elev)); %mi serve lanciarlo per vedere n_inc
if n_elev==2
li=length(t);
n_inc=ceil(t(li)/(camp*60));
end
fine=partenza+n_inc;
if partenza==1
fine=fine-1;
end
one=0;
end
if cond==1 || ultimo==1
ultimo=0;
disp([res(ir,:),' - y: ',year(iy,:),' - m: ',num2str(m),' / ',num2str(head(1))]);
[A,t,deltaTa]=SyntStormEM(tsR,str2num(res),lat,Hamsl,freq,pola,g_elev(on_elev));%usa la nuova Hr di ITU-R
% %Grafico attenuazione x controllo
% figure(n_elev)
% leg=g_elev(on_elev);
% legend(num2str(leg))
% hold on;
% plot(t,A,'Color', colori(temp,:)) %gi? in secondi
% xlabel('Tempo (s)')
% ylabel('Attenuazione (dB)')
% pause(1)
TtotRain=TtotRain+length(tsR);
Ttot=Ttot+length(t);
fname=[pathout,'SST_WHS19_',year(iy,3:4),'XX',num2str(evento_pioggia,'%04.0f'),'.dat'];
fido=fopen(fname,'a');
fprintf(fido,'%%Channel model: SST, Sampling rate: %4.2f sample/s, Station name: White Sands, Link elevation: 37.8?.\n',1/(t(2)-t(1)));
fprintf(fido,'\n');
fne=(camp*(temp-1)*60);
% if fne>length(t)
% fne=length(t);
% end
for c=1:1:length(A)
if t(1,c)>=part && t(1,c)<=fne
fprintf(fido,'%6.3f\n',A(c));
y(d)=A(c);
x(d)=(t(c));
ffne=ffne+1;
d=d+1; %devo vedere dove azzerarlo a 1 /fatto
end
end
% %stampa grafico pezzi
% figure(999)
% plot(x,y)
% axis square
% xlabel('Tempo (s)')
% ylabel('Attenuazione (dB)')
if scrivi==1 %scrivilo su file
f = figure('Visible', 'off');
plot(x,y,'-','linewidth',1);
xlabel('Tempo (s)')
ylabel('Attenuazione (dB)')
saveas(f,['Grafici/Totale/Grafico_totale_evento',num2str(head(1)), '_del_' num2str(head(3)),'-',num2str(head(4)),'-',num2str(head(2)),'.jpeg']);
scrivi=0;
clear x
clear y
d=1;
end
part=fne;
fclose(fido);
end
arr=arr+1;
end
end
fclose(nfid);
if temp2<fine %tramonta ma continua a piovere %verifica se spostare
manca=fine-temp2;
fine=fine-manca;
end
if fine~=0 %controlla se serve e va qui
partenza=1;
end
% %funzione che salva in automatico i grafici
% f = figure('Visible', 'off');
% plot(t,A,'r');
% saveas(f,['Grafici/',num2str(elev),'_Evento_',num2str(head(1)), '_del_' num2str(head(3)),'-',num2str(head(4)),'-',num2str(head(2)),'.jpeg']);
end
end
fclose(fid);
end
end
deltaTA=t(2)-t(1);
save SSTTimeInfo9400 TtotRain Ttot deltaTA
if you want to try the code create a txt called r9460.txt and paste this text
1 1994 1 8 3 59 27 0 15 54 11
8.9 1
9.2 1
10.2 1
12.9 1
11.4 1
10.7 1
12.4 1
8.3 1
7.7 1
2.2 6
2.0 1
2 1994 2 8 0 32 30 0 22 22 15
17.2 1
3.3 4
2.3 1
2.3 5
15.3 1
56.9 1
17.9 1
18.7 1
24.3 1
12.5 1
9.1 1
6.3 2
7.1 1
7.5 1
2.7 1
3 1994 2 8 22 56 18 0 4 5 5
15.8 1
12.2 1
33.6 1
28.0 1
1.8 1
4 1994 3 8 0 41 7 0 45 59 34
31.3 1
90.6 1
119.4 1
126.4 1
93.9 1
110.6 1
119.3 1
116.2 1
86.4 1
104.1 1
114.3 1
144.8 1
94.9 1
72.7 1
80.3 1
84.3 1
76.7 1
63.0 1
59.2 1
46.1 1
33.5 1
20.6 1
15.9 1
7.9 1
7.6 1
5.4 1
5.3 1
5.1 1
3.0 4
3.0 1
2.7 5
2.5 1
2.1 6
2.0 1
5 1994 4 8 9 14 40 0 31 16 5
0.6 1
0.6 25
0.9 1
3.5 4
0.9 1
6 1994 14 8 23 37 16 0 12 28 12
15.9 1
12.9 1
18.9 1
24.9 1
22.6 1
15.9 1
10.0 1
4.8 1
4.3 2
5.9 1
8.2 1
3.8 1
7 1994 15 8 13 35 24 1 4 25 5
0.4 1
0.4 37
0.4 1
0.6 25
0.2 1
8 1994 16 8 1 48 41 0 2 3 3
7.4 1
7.4 1
0.4 1
9 1994 23 8 19 26 30 0 20 22 10
19.5 1
22.7 1
17.2 1
3.9 1
2.5 5
2.8 1
2.8 4
2.7 1
2.6 5
1.0 1
10 1994 30 8 0 36 14 3 25 50 56
23.6 1
24.3 1
25.7 1
33.4 1
21.2 1
19.2 1
21.8 1
27.6 1
17.0 1
16.2 1
9.3 1
5.3 1
1.7 8
8.8 1
14.1 1
18.4 1
24.0 1
12.3 1
2.3 6
5.3 1
8.1 1
9.0 1
9.4 1
4.9 3
0.6 1
0.3 56
9.3 1
19.5 1
30.8 1
28.1 1
30.9 1
33.7 1
26.8 1
24.2 1
34.5 1
36.5 1
35.4 1
45.2 1
51.3 1
27.7 1
27.6 1
56.7 1
44.5 1
31.8 1
3.5 4
1.9 1
1.0 14
0.8 1
0.4 35
0.5 1
1.1 14
1.5 1
1.5 9
1.3 1
1.2 11
1.0 1
11 1994 31 8 4 30 38 0 33 29 11
40.9 1
42.9 1
39.4 1
20.0 1
12.1 1
10.6 1
5.0 1
4.0 2
3.9 1
0.6 23
0.3 1
12 1994 3 9 23 44 38 1 8 35 15
1.2 1
1.2 11
1.0 1
0.9 17
3.6 1
3.8 3
2.9 1
2.9 4
1.8 1
1.2 12
1.4 1
1.7 8
1.9 1
2.1 6
1.2 1
13 1994 4 9 22 18 53 0 42 49 9
6.5 1
6.5 1
4.0 1
2.6 5
0.9 1
0.7 20
0.7 1
1.2 12
1.0 1
14 1994 8 9 21 30 41 0 18 23 19
11.6 1
13.6 1
17.4 1
26.3 1
26.6 1
26.5 1
23.5 1
22.7 1
34.1 1
31.6 1
32.2 1
33.5 1
37.7 1
35.3 1
34.7 1
34.2 1
38.1 1
18.4 1
5.1 1
15 1994 14 10 11 37 35 0 36 23 7
0.6 1
0.6 25
2.1 1
9.8 1
4.9 1
1.9 7
0.7 1
16 1994 15 10 2 48 1 2 10 30 41
3.7 1
3.7 3
2.6 1
2.5 5
1.6 1
1.5 9
6.4 1
6.5 1
3.8 1
2.2 6
0.9 1
0.3 47
23.4 1
11.6 1
6.1 1
5.1 1
1.8 8
17.7 1
18.2 1
23.0 1
29.0 1
26.5 1
23.7 1
15.7 1
10.5 1
31.9 1
35.4 1
28.6 1
33.3 1
16.9 1
18.5 1
25.5 1
28.7 1
20.2 1
9.1 1
1.3 11
1.5 1
1.5 9
29.6 1
39.6 1
7.6 1
17 1994 3 11 1 26 1 0 29 59 15
22.4 1
24.9 1
15.8 1
18.1 1
20.3 1
23.1 1
16.4 1
10.3 1
17.2 1
5.5 2
4.0 1
1.5 9
1.6 1
1.9 7
1.8 1
18 1994 8 11 14 30 15 0 33 11 3
0.5 1
0.5 32
0.1 1
19 1994 11 11 17 36 51 1 17 54 5
0.5 1
0.5 29
0.4 1
0.3 46
0.3 1
20 1994 11 11 21 35 55 0 19 10 3
0.8 1
0.8 18
0.1 1
21 1994 11 11 23 31 42 0 58 51 33
1.4 1
1.4 9
2.1 1
3.1 4
2.9 1
2.8 5
5.0 1
5.2 2
6.6 1
6.6 1
7.1 1
7.3 1
7.1 1
7.0 1
7.5 1
8.1 1
7.5 1
7.0 1
5.8 1
3.8 3
4.1 1
4.8 2
4.8 1
5.1 2
5.0 1
5.0 2
5.0 1
4.9 3
4.0 1
4.0 2
4.0 1
3.8 3
3.2 1
22 1994 12 11 3 54 31 1 7 20 9
2.1 1
2.1 6
1.3 1
0.8 17
0.7 1
0.5 32
0.6 1
1.8 8
0.6 1
23 1994 12 11 7 31 40 0 11 46 3
1.3 1
1.3 10
1.0 1
24 1994 12 11 9 39 17 0 26 28 3
0.6 1
0.6 25
0.3 1
25 1994 12 11 12 14 52 2 49 41 41
1.7 1
1.7 8
1.0 1
1.0 14
0.6 1
0.4 37
0.5 1
0.6 27
2.4 1
2.5 5
3.1 1
3.2 3
3.2 1
3.6 4
4.1 1
4.2 2
4.3 1
4.7 3
4.9 1
4.9 2
4.8 1
4.8 2
6.6 1
7.3 1
5.5 1
4.4 2
4.5 1
5.2 2
4.8 1
3.4 4
3.4 1
3.4 3
3.3 1
2.8 5
2.6 1
2.5 5
2.2 1
2.1 6
1.7 1
1.0 14
0.7 1
26 1994 12 11 16 42 24 0 59 39 15
2.0 1
2.0 6
2.2 1
2.3 6
3.3 1
3.6 3
2.9 1
2.3 6
2.4 1
2.4 5
1.9 1
1.6 8
1.6 1
0.8 18
0.5 1
27 1994 12 11 19 42 19 1 21 25 5
0.6 1
0.6 25
0.3 1
0.3 54
0.1 1
28 1994 5 12 16 56 50 12 26 21 133
0.6 1
0.6 26
0.4 1
0.4 35
0.9 1
1.0 14
1.0 1
0.9 16
0.9 1
1.3 11
1.3 1
1.4 10
1.5 1
1.8 8
2.2 1
2.4 5
2.3 1
1.9 7
1.9 1
1.8 8
1.8 1
1.8 7
1.7 1
1.7 8
1.8 1
1.9 7
1.5 1
1.0 14
1.4 1
2.7 5
3.1 1
3.4 3
3.4 1
3.6 4
2.6 1
2.4 5
2.3 1
2.3 6
1.6 1
1.6 8
1.8 1
2.3 6
2.3 1
2.3 6
1.8 1
1.7 8
1.5 1
1.5 9
1.4 1
1.4 9
1.4 1
1.3 11
1.5 1
1.9 7
1.7 1
1.2 12
1.1 1
1.1 12
1.1 1
1.2 12
1.6 1
2.4 5
2.4 1
1.4 10
1.5 1
1.6 9
1.9 1
1.9 7
2.0 1
2.0 6
2.2 1
3.1 4
3.1 1
3.3 4
3.7 1
3.9 3
2.8 1
2.4 5
2.0 1
1.3 11
3.0 1
3.3 3
3.4 1
3.6 4
1.0 1
0.9 15
2.6 1
4.2 3
7.3 1
7.7 1
4.6 1
4.3 2
3.5 1
1.8 8
3.0 1
3.0 4
2.5 1
2.4 5
2.1 1
1.9 7
1.4 1
1.0 14
1.2 1
1.2 12
1.5 1
1.6 8
1.5 1
1.2 12
1.2 1
1.2 12
1.2 1
1.3 11
1.2 1
1.0 14
0.7 1
0.4 37
0.6 1
1.0 15
0.5 1
0.4 33
1.2 1
1.5 9
2.4 1
2.9 4
2.4 1
1.9 7
1.6 1
1.0 15
0.8 1
0.6 23
0.5 1
0.4 39
0.1 1
29 1994 6 12 8 4 2 3 40 45 70
1.4 1
1.4 9
2.5 1
4.8 2
4.6 1
3.8 3
3.5 1
1.2 12
1.6 1
1.8 7
1.9 1
4.0 3
4.3 1
5.2 2
5.8 1
7.0 1
7.0 1
7.0 2
7.5 1
7.6 1
5.7 1
5.6 1
6.0 1
7.2 1
7.0 1
6.0 2
8.6 1
10.1 1
8.3 1
7.4 1
5.0 2
4.7 1
3.9 3
3.9 1
3.9 3
3.1 1
1.9 7
1.7 1
0.8 18
1.8 1
2.3 6
1.1 1
1.0 14
1.6 1
4.2 3
5.1 1
5.9 2
4.1 1
4.0 2
4.0 1
4.2 3
4.4 1
4.5 2
4.1 1
1.8 8
1.7 1
1.7 8
2.0 1
2.2 6
1.8 1
1.5 9
1.3 1
0.9 16
1.3 1
1.9 7
1.5 1
1.1 13
1.8 1
2.0 6
1.5 1
30 1994 12 12 23 18 13 6 47 6 21
0.3 1
0.3 44
0.3 1
0.3 49
0.4 1
0.4 39
0.4 1
0.3 46
0.4 1
0.4 36
0.4 1
0.4 41
0.4 1
0.4 34
0.4 1
0.3 46
0.4 1
0.4 35
0.5 1
0.6 27
0.1 1
31 1994 26 12 10 50 5 0 8 53 7
4.5 1
4.5 2
10.3 1
13.1 1
10.2 1
4.8 2
4.2 1
32 1994 27 12 4 10 20 0 42 48 17
2.4 1
2.4 5
3.2 1
3.9 3
3.9 1
3.9 3
3.4 1
3.2 4
2.7 1
2.6 4
2.7 1
2.8 5
3.6 1
3.8 3
2.3 1
1.8 7
1.4 1
33 1994 27 12 16 41 20 0 17 42 7
1.8 1
1.8 7
2.3 1
3.0 4
3.2 1
3.8 3
2.7 1
Thanks so much for the help

请先登录，再进行评论。

Nicola
Nicola 2014-2-13
is anyone able to help me? :(
  2 个评论
Rizwana
Rizwana 2014-2-14
编辑：Rizwana 2014-2-14
Iam beginner in matlab. Best way to know why there is errors is run each command by yourself specially when the code is written by someone else. Moreover comments in the code are not in english. Iam not able to understand what it is doing.
Nicola
Nicola 2014-2-19
You are right, so i simplified the code

请先登录，再进行评论。

Nicola
Nicola 2014-2-19
Guys, i simplified the problem The code is:
name=['r9460.txt']
fid=fopen(name,'r');
while ~(feof(fid)),
head=fscanf(fid,'%d',11)
if ~(isempty(head)),
m=head(4)
end
end
fclose(fid)
The txt is
5 1994 4 8 9 14 40 0 31 16 5
0.6 1
0.6 25
0.9 1
3.5 4
0.9 1
Matlab says Attempted to access head(4); index out of bounds because numel(head)=1.
Please, help me, i'm in your hands!
  1 个评论
Walter Roberson
Walter Roberson 2015-8-12
head=fscanf(fid,'%d',11) tells MATLAB to read 11 elements even if it requires reading multiple lines to do so. The first line of the file has 11 elements and so works, but then multiple lines get read for the next fscanf(). Reading will leave off even in the middle of a text line if that is what is needed to read 11 items, and the next fscanf() would pick up in the middle of a line if need be. At some point you are ending up at the end of file with a line which is not empty but which does not have 4 items to index.
If you have a mix of line formats that imply different data structures then you need a mix of fscanf() calls. If you do not know ahead of time which data structure a line is going to be for until you read the line and count the number of items, then you need to read the line as a string, then use sscanf() to parse the string to find out how many items are on the line and react accordingly.
Any time the line format has to be examined to determine something about how the data should be interpreted, you should be reading a line at a time use fgetl() or fgets() and examining the line.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Data Type Identification 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by