请给出一段频谱扩展与加窗收缩的matlab示例代码

这是一段频谱扩展与加窗收缩的matlab示例代码：

% 频谱扩展
N = 512; % FFT长度
x = randn(N,1); % 生成随机信号
X = fft(x); % 进行FFT变换
k = 2; % 扩展倍数
X_ext = [X(1:N/k); zeros(N*(k-1),1); X(N/k+1:end)]; % 将FFT结果扩展k倍
x_ext = ifft(X_ext); % 进行IFFT变换

% 加窗收缩
M = 256; % 截断FFT长度
w = hamming(N); % 生成汉明窗
X_win = fft(x.*w); % 加窗后进行FFT变换
X_shrink = [X_win(1:M/2); X_win(N-M/2+1:N)]; % 截断FFT结果
x_shrink = real(ifft(X_shrink)); % 进行IFFT变换并取实部

希望这个示例代码能够帮到你！