无法在 MATLAB 中生成系数矩阵

Question

我有以下一段代码，我正在尝试计算系数矩阵a_k 来求解线性系统以获得h[n]，即脉冲响应。我正在使用快速傅里叶逆变换。

N = 9; % period is chosen to be 16
n = 0:N-1; %vector for x must start at n = 0
for k = 1:9
    y3 = zeros(1,9);
    y3(k+1) = N/2;
    y3(N - k + 1) = N/2;  
end
x3 = ifft(y3);

figure;
subplot(2,2,1);stem(n,real(x3));xlabel('n'); //line 52
ylabel('real(x3)');axis([0 N-1 -1 1]);
subplot(2,2,2);stem(n,imag(x3));xlabel('n');
ylabel('imag(x3)');axis([0 N-1 -1 1]);
subplot(2,2,3);stem(n,real(y3)/N);xlabel('k');
ylabel('real(a_k)');axis([0 N-1 -1 1]);
subplot(2,2,4);stem(n,imag(y3)/N);xlabel('k');
ylabel('imag(a_k)');axis([0 N-1 -1 1]);

但是，当我运行此代码时，我收到以下错误：

Error using stem (line 43)
X must be same length as Y.

Error in fft_examples (line 52)
subplot(2,2,1);stem(n,real(x3));xlabel('n');

我不确定我在哪里出错。我知道k 的矩阵是从1 到9。因此，我做了一个for 循环。 y 值变得不匹配。

Answer 1

size(real(x3))  % --> 1 10
size(n)         % --> 1 9

所以它们的大小不一样。您正在增加 y3(k+1) = N/2; 中 y3 的大小

另外，为什么要在每次迭代中创建矩阵 y3：y3 = zeros(1,9);