使用matlab軟體畫出頻譜圖
x = [zeros(1,3*m),ones(1,n),zeros(1,3*m)]; % 給出輸入序列
w = [-15:0.1:15] + 1e-10;
x = sin(0.5*n*w*t)./sin(0.5*w*t); % 給出頻譜序列
subplot(1,3,1); stem(n,x,'.'); % 畫出輸入序列
axis([-20,20,-0.1,1.1]), grid on
xlabel('n'); title('(a) 序列幅度');
subplot(1,3,2); plot(w,x),grid on % 畫出頻譜序列
xlabel('\omega'); title('(b) 幅頻特性');
subplot(1,3,3); plot(w,x),grid on % 改變橫軸比例，畫出頻譜序列
v = axis;axis([-pi/t,pi/t,v(3),v(4)]);
xlabel('\omega'); title('(c) 橫軸放大後幅頻特性');
set(gcf,'color','w'); % 置圖形背景色為白
結果如圖所示：
經過驗證，頻譜圖符合要求。
離散時間傅利葉變換（DTFT）簡記
本文適合在看完 離散訊號與模擬訊號之間的頻率關係 由模擬訊號取樣得到的離散訊號 這篇博文之後閱讀，這樣才能理解後面的前因後果，也就是我的邏輯！下面先簡單粗暴地介紹一下dtft 性質只列了一條，因為我想討論的內容和這個相關，這不是一篇科普性的博文，所以沒有正正經經地把性質列出來，背景說出來，那樣不如去...
使用python實現離散時間傅利葉變換
我們經常使用傅利葉變換來計算數碼訊號的頻譜，進而分析數碼訊號，離散時間傅利葉變換的公式為 x ejw n n x n e jwn x e jw n n x n e j wn可是自己動手實現一遍才是最好的學習。在數字分析裡面，傅利葉變換預設等時間間隔取樣，不需要時間序列，只需要訊號陣列即可分析。分析過...
離散傅利葉變換 快速傅利葉變換C 實現
傅利葉變換是將時域訊號變換為頻域訊號的一種方式，我主要用它來做兩件事情 1 求一段資料的週期性。2 通過傅利葉變換及其逆變換，進行低通濾波 去躁 首先需要做幾點說明 1.快速傅利葉變換是離散傅利葉變換的快速演算法，當資料來源較大時 大於1000 快速傅利葉變換有明顯優勢。2.快速傅利葉變換的訊號源長...