Why am I getting Index in position 1 is invalid. Array indices must be positive integers or logical values.

Question

Aijalon Marsh 2022-2-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1648135-why-am-i-getting-index-in-position-1-is-invalid-array-indices-must-be-positive-integers-or-logical

回答： Image Analyst 2022-2-12

clear

syms t;

Fs = 1.5; % Sampling frequency

T = 1/Fs; % Sampling period

L = 1500; % Length of signal

t = (0:L-1)*T;

Yo=4+4.005*(sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7);

n = 2^nextpow2(L);

y = fft(Yo);

f = 2*pi*Fs*(0:(L/2))/L; % <- multiplying

p2 = abs(y/L);

p1=p2(1:L/2+1);

p1(2:end-1) = 2*p1(2:end-1);

plot(0:(Fs/n):(Fs/2-Fs/n),p1(i,1:n/2))

title ('Amplitude Spectrum')

xlabel('time')

ylabel('y(t)')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

the cyclist 2022-2-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1648135-why-am-i-getting-index-in-position-1-is-invalid-array-indices-must-be-positive-integers-or-logical#answer_894315

在 MATLAB Online 中打开

clear 
syms t; 
Fs = 1.5;            % Sampling frequency                    
T = 1/Fs;             % Sampling period       
L = 1500;             % Length of signal
t = (0:L-1)*T;
Yo=4+4.005*(sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7); 
n = 2^nextpow2(L);
y = fft(Yo); 
f = 2*pi*Fs*(0:(L/2))/L; % <- multiplying
p2 = abs(y/L); 
p1=p2(1:L/2+1); 
p1(2:end-1) = 2*p1(2:end-1); 
 
plot(0:(Fs/n):(Fs/2-Fs/n),p1(i,1:n/2))
Index in position 1 is invalid. Array indices must be positive integers or logical values.
title ('Amplitude Spectrum') 
xlabel('time') 
ylabel('y(t)')

In the expression

p1(i,1:n/2)

you have not defined i, so MATLAB is trying to use the complex value sqrt(-1) as an index.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

David Hill 2022-2-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1648135-why-am-i-getting-index-in-position-1-is-invalid-array-indices-must-be-positive-integers-or-logical#answer_894320

编辑：David Hill 2022-2-12

在 MATLAB Online 中打开

syms t; 
Fs = 1.5;            % Sampling frequency                    
T = 1/Fs;             % Sampling period       
L = 1500;             % Length of signal
t = (0:L-1)*T;
Yo=4+4.005*(sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7); 
n = 2^nextpow2(L);
y = fft(Yo); 
f = 2*pi*Fs*(0:(L/2))/L; % <- multiplying
p2 = abs(y/L); 
p1=p2(1:L/2+1); 
p1(2:end-1) = 2*p1(2:end-1); 
 
plot(t(1:length(p1)),p1);%what is the i? The length of your t array must be equal to the length of p1 
title ('Amplitude Spectrum') 
xlabel('time') 
ylabel('y(t)')