搞fft时出现了索引超出数组元素的数目(6)该咋办啊qwq 提醒的是在这一步FTy_a=FTy_a(1:K/2);

Question

问题遇到的现象和发生背景

问题相关代码，请勿粘贴截图

clc;
close all;
clear all;
[y,Fs]=audioread('雷声.wav');
sound(y,Fs);
K=length(y);
t=0.1/Fs:(K-1)/Fs;
fo=Fs/K;
f=0:fo:Fs/2-fo;
FTy=fft(t);
FTy_a=abs(FTy);
FTy_a=FTy_a(1:K/2);
figure(1)
subplot(2,1,1);
plot(y);xlabel('t');
title('语音信号时域波形图');
subplot(2,1,2);
plot(f,FTy_a);xlabel('t');
title('语音信号频域波形图')

运行结果及报错内容

索引超出数组元素的数目(6)。

出错 xinxichuli (第 12 行)
FTy_a=FTy_a(1:K/2);

我的解答思路和尝试过的方法

我想要达到的结果

Answer 1

你好，代码里面有两处错误、一处警告，一一帮你修改 并加以说明了，有帮助还望给个采纳支持一下（点击下方“采纳该答案”，非常感谢）

clc;
close all;
clear;
[y,Fs]=audioread('handel.wav');
sound(y,Fs);
K=length(y);
t=0:1/Fs:(K-1)/Fs; % 这里给你修改成0:1/Fs:(K-1)/Fs
fo=Fs/K;
f=0:fo:Fs/2-fo;
FTy=fft(y); % fft是对y进行变换的
FTy_a=abs(FTy);
FTy_a=FTy_a(1:floor(K/2)); % 这里用floor保证整数下标
figure(1)
subplot(2,1,1);
plot(y);xlabel('t');
title('语音信号时域波形图');
subplot(2,1,2);
plot(f,FTy_a);xlabel('t');
title('语音信号频域波形图')