exe0308.m

来自「Matlab程序为清华大学胡广书版信号处理的上机习题答案」· M 代码 · 共 31 行

clear;

N=50;

PP=0:1:(N-1);
k=0:1:(N-1);
for kj=0:1:(N-1)
    
    PPtemp=0;
    for ni=0:1:(N-1)
        pptemp=0;
        for m=0:1:(N/2-1)
            if (ni==(2*m))
                delt=1;
            else
                delt=0;
            end
            pptemp=pptemp+delt*exp(-j*2*pi/N*kj*ni);
        end
        PPtemp=PPtemp+pptemp;
    end
    PP(kj+1)=PPtemp;
end
PP
the=angle(PP)*180/pi
figure
PP=abs(PP);
stem(k,PP)