matlab中生成基四FFT然后输入信号频率为2KHz的方波信号，并且生成至少5次谐波，计算谐波失真度和谐波的归一化振幅

好的，我可以回答你的问题。首先，我们需要生成基四FFT，可以使用matlab中的fft函数，具体代码如下：

N = 1024; % 采样点数
fs = 8e3; % 采样率
Ts = 1/fs; % 采样间隔
t = (0:N-1)*Ts; % 时间向量

f0 = fs/N; % 基频
f = (0:N-1)*f0; % 频率向量

W = exp(-1j*pi/2); % 基四根号

X = zeros(1,N); % 生成基四FFT的系数
for k = 0:N/4-1
    X(1+4*k) = 1;
    X(1+4*k+1) = W^(k*1);
    X(1+4*k+2) = 1j;
    X(1+4*k+3) = W^(k*3);
end

接下来，我们可以生成频率为2kHz的方波信号，并使用fft函数进行傅里叶变换，得到频域波形。具体代码如下：

f1 = 2e3; % 方波信号频率
A = 1; % 方波信号幅值
x = A*square(2*pi*f1*t); % 生成方波信号

Xf = fft(x); % 使用fft函数进行傅里叶变换

现在，我们可以计算谐波失真度和谐波的归一化振幅了。由于方波信号的频率为2kHz，我们可以通过观察频域波形的谐波分量来计算失真度和归一化振幅。具体代码如下：

harm = 5; % 计算前5个谐波

% 计算失真度和归一化振幅
d = zeros(1,harm+1);
A_norm = zeros(1,harm+1);
for k = 1:harm+1
    d(k) = abs(Xf(k))/abs(Xf(1));
    A_norm(k) = abs(Xf(k))/(N/2);
end

% 输出结果
for k = 1:harm+1
    fprintf('谐波 %d 失真度为 %f\n',k-1,d(k));
    fprintf('谐波 %d 归一化振幅为 %f\n',k-1,A_norm(k));
end

上述代码中，我们计算了前5个谐波的失真度和归一化振幅，并输出了结果。你可以根据需要修改harm的值来计算更多的谐波。