matlab

[x1freqz, x1freqzw]=freqz(x,1,3000,"whole"); 

[X, W]=me_dtft(x1",pi,3000); figure title("my") % plot(W/pi,20*log10(abs(X))); plot(W/pi,abs(X)) ax = gca; % ax.YLim = [-40 70]; xlabel("Normalized Frequency (\times\pi rad/sample)") ylabel("Magnitude (dB)") function [X, w]=me_dtft(x,whalfrange, nsample) w= linspace(-whalfrange,whalfrange,nsample); t=0:1:size(x,2)-1; X=zeros(1,size(w,2)); for i=1:1:size(w,2) X(i)=x*exp(-t*1i*w(i))"; end end 

何か間違ったことを教えてください！これらの変換は、純粋な正弦波（図のコード）から取得されます。右から、fft、freqzを見ることができます。マナーは上と（左）前に説明したようにDTFTを説明しました。

「ジェイソンR」のコメントの後に編集：わかりました、混乱するので対数目盛も削除しました。その後、直感的にあなたは同様に、次の画像でわかるように、なぜそれらが完全に同じではないのですか（対数目盛で最後の画像を参照しますか？）？

freqz：

[x1freqz, x1freqzw]=freqz(fliplr(XX"),1,3000,"whole"); figure title("freqz") % plot((x1freqzw/pi)-1,20*log10(abs(fftshift(x1freqz)))) plot((x1freqzw/pi)-1,abs(fftshift(x1freqz))) ax = gca; % ax.YLim = [-40 70]; ax.XTick = -1:.5:2; xlabel("Normalized Frequency (\times\pi rad/sample)") ylabel("Magnitude") 

サインサンプル：

Fs=1000; Ts=1/Fs; time=0:Ts:1; Freqs=500; Xs=zeros(length(Freqs),length(time)); for i=1:length(Freqs) Xs(i,:)= cos(2*pi*Freqs(i)*time); end XX=Xs; XX=XX./ max(abs(XX)); figure;plot(time, XX); axis(([0 time(end) -1 1])); xlabel("Time (sec)"); ylabel("Singal Amp."); title("A sample audio signal"); sound(XX,Fs) 

コメント

• プロットを生成する完全なスクリプトを投稿する必要があります
• 1つの違いは、コードでDTFTを-fs / 2から+ fs / 2にプロットしているのに対し、FFTは0からfsになっていることです（ただし、 " fft shitt "）すべてのコードが表示されない場合、2つのプロットの違いは丸め誤差であり、より多くのポイントを含めるか、倍精度浮動小数点数にすることで減少すると思われます。 DTFTは、サンプリングレートの整数の倍数である周波数の正弦波を作成しない限り、FFTがサンプルとなるSinc関数になることに注意してください。また、FFTをゼロパディングすることでDTFTのサンプルを取得できることも知っておいてください。
• @ DanBoschen、fftの図では、1000はfs / 2を意味しますか？
• 1000はfs、波形を意味しませんはfs / 2にあります
• @DanBoschenなぜそれら（私のDTFTとfreqz、つまり）が完全に同じではないのか説明できますか？ fftをパディングすると、sincが表示されますか？