dtft code does not work...please help me

clc;
clear;
clear all;
m=3;
omegac=(pi/4);
omega=(-pi:0.01:pi);
Hm_dtft=(sin(omegac*m)./(pi*m)).*exp(-1j*omega*m);

magHm_dtft=abs(Hm_dtft); 

plot(omega, abs(Hm_dtft),'k');grid;
xlabel('\Omega '); ylabel ('|Hm(e^{j\Omega})|')
title('Magnitude Part')

[Matt J]

what does it do?

Answer 1

[Walter Roberson]

m=3;
omegac=(sym(pi)/4);
syms omega real
Hm_dtft=(sin(omegac*m)./(pi*m)).*exp(-1j*omega*m);
magHm_dtft=abs(Hm_dtft)

magHm_dtft = 

Notice the result is independent of the value of omega, which explains why your plot was a constant.

[LcySkywalker]

I just started learning Matlab. Thank you so much :)