function?[Xk?w]?=?DFT（x?M?str）
%離散傅里葉變換及x（n）FT[x（n）]的繪圖函數
????n?=?0?:?M?-?1;
????%單位圓離散譜點數
????N?=?200;
????%畫2個周期
????k?=?-N?:?N?-?1;
????%離散頻率橫坐標
????w?=?（2?*?pi?/?N）?*?k;
????%x（e^jw）=∑x（n）e^（jwn），使用矩陣乘法實現
????Xk?=?x?*?exp（?-j?*?（2?*?pi?/?N）?）?.^?（?n‘?*?k?）;
????magX?=?abs（Xk）;
????
????subplot（2?2?1）;
????t?=?max（x）;
????n?=?0?:?M?-?1;
????stem（n?x?‘.‘）;
????xlabel（‘n‘）;
????ylabel（‘x（n）‘）;
????title（‘時域信號波形‘）;
????text（?（0.3?*?M）?（0.8?*?t）?str?）;
????
????subplot（2?2?3）
????plot（w/pi?magX?‘b‘）;
????xlabel（‘ω/π‘）;
????ylabel（‘|X（e^j^w）|‘）;
????title（‘信號的幅度譜‘）;

%?????%使用FFT算法
%?????subplot（2?2?4）
%?????Xk?=?fft（x?N）;
%?????Xk?=?[Xk?Xk];%補足2個周期
%?????magX?=?abs（Xk）;
%?????stem（w/pi?magX?‘r‘）;
%?????xlabel（‘ω/π‘）;
%?????ylabel（‘|X（e^j^w）|‘）;
%?????title（‘信號的FFT幅度譜‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????872??2008-10-15?18:01??實驗1信號系統及其響應參考代碼\DFT.m

?????文件?????????84??2006-11-05?09:11??實驗1信號系統及其響應參考代碼\impseq.m

?????文件???????5617??2008-10-15?18:01??實驗1信號系統及其響應參考代碼\main.m

?????文件????????143??2008-10-15?17:14??實驗1信號系統及其響應參考代碼\mainseq.m

?????文件?????????90??2008-10-15?17:17??實驗1信號系統及其響應參考代碼\stepseq.m

?????目錄??????????0??2008-10-15?18:04??實驗1信號系統及其響應參考代碼

-----------?---------??----------?-----??----

?????????????????6806????????????????????6