[A,I]=sort(y1,'descend');%频谱高到低排列，显示各次谐波。I为返回索引 序列，表示A中元素与y1中元素的对应。m=0:799; fn=(m/800).*fs; %频率C={0,0,0,0,0,0};%表示谐波数组for i=2:7 C{i-1}=A(i)/A(1)*100;%谐波畸变率 if C{i-1}>0.05*100 disp('F') else disp('T') end end

这是一段 MATLAB 代码，用于计算一个信号的频域信息，并检查其谐波畸变率是否超过了阈值。具体来说，代码首先对信号的频谱进行了高到低的排列，并返回了排列后的索引序列。然后，代码计算了信号的各次谐波的幅值，并将其存储在一个数组中。接下来，代码计算了各次谐波的畸变率，并根据阈值判断是否存在谐波畸变。如果存在谐波畸变，则输出 'F'，否则输出 'T'。