帮我用MATLAB代码实现DTFT

以下是MATLAB代码示例：

%生成一个长度为N的信号x
N = 1000;
n = 0:N-1;
x = sin(2*pi*0.05*n) + 0.5*sin(2*pi*0.2*n) + 0.3*sin(2*pi*0.3*n);

%计算DTFT
w = -pi:pi/100:pi; %频率序列
X = zeros(size(w));
for k = 1:length(w)
X(k) = sum(x.*exp(-1j*w(k)*n));
end

%绘制幅度和相位谱
subplot(2,1,1);
plot(w, abs(X));
xlabel('频率');
ylabel('幅度');
title('DTFT幅度谱');
subplot(2,1,2);
plot(w, angle(X));
xlabel('频率');
ylabel('相位');
title('DTFT相位谱');

请注意，此代码仅仅是一个简单的示例，您需要根据自己的需求和信号进行适当的修改。