matlab数字信号处理求DTFT(离散时间傅里叶变换)图像和arg图像

问题
matlab运行报错,希望能得出正确图像

img

题目

img


代码


```c
clc;
clear;
n=0:29;
x1=ones(1,n);
x2=[zeros(1,15),1,ones(1,29)];

x=(0.6).^n*(x1-x2);
k=0:28000;w=(pi/28000)*k;
X=x*(exp(-j*pi/28000)).^(n'*k);
magX=abs(X);angX=angle(X);
subplot(1,2,1);stem(n,x,'.');
title('序列图');
ylabel('x(n)');axis([0,5,0,6]);grid;
subplot(1,2,2);plot(w/pi,angX);grid;
title('arg');ylabel('弧度');