x1n=[ones（14）];?????%產生R4（n）序列向量
X1k8=fft（x1n8）;?????%計算x1n的8點DFT
X1k16=fft（x1n16）;???%計算x1n的16點DFT
%以下繪制幅頻特性曲線
N=8;
f=2/N*（0:N-1）;
figure（1）;
subplot（121）;stem（fabs（X1k8）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（1a）?8點DFT[x_1（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=16;
f=2/N*（0:N-1）;
subplot（122）;stem（fabs（X1k16）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（1a）?16點DFT[x_1（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
%x2n?和?x3n
M=8;xa=1:（M/2）;??xb=（M/2）:-1:1;?
x2n=[xaxb];????%產生長度為8的三角波序列x2（n）
x3n=[xbxa];
X2k8=fft（x2n8）;
X2k16=fft（x2n16）;
X3k8=fft（x3n8）;
X3k16=fft（x3n16）;
figure（2）;
N=8;
f=2/N*（0:N-1）;
subplot（221）;stem（fabs（X2k8）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（2a）?8點DFT[x_2（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（223）;stem（fabs（X3k8）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（3a）?8點DFT[x_3（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=16;
f=2/N*（0:N-1）;
subplot（222）;stem（fabs（X2k16）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（2a）?16點DFT[x_2（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（224）;stem（fabs（X3k16）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（3a）?16點DFT[x_3（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
%x4n?和?x5n
N=8;n=0:N-1;
x4n=cos（pi*n/4）;
x5n=cos（pi*n/4）+cos（pi*n/8）;
X4k8=fft（x4n8）;
X4k16=fft（x4n16）;
X5k8=fft（x5n8）;
X5k16=fft（x5n16）;
figure（3）;
N=8;
f=2/N*（0:N-1）;
subplot（221）;stem（fabs（X4k8）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（4a）?8點DFT[x_4（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（223）;stem（fabs（X5k8）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（5a）?8點DFT[x_5（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=16;
f=2/N*（0:N-1）;
subplot（222）;stem（fabs（X4k16）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（4a）?16點DFT[x_4（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（224）;stem（fabs（X5k16）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（5a）?16點DFT[x_5（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;

%x8n
Fs=64;?T=1/Fs;?
N=16;n=0:N-1;?%對于N=16的情況
nT?=?n*T;
x8n=cos（8*pi*nT）+cos（16*pi*nT）+cos（20*pi*nT）
X8k16=fft（x8n16）;
N=16;
f=2/N*（0:N-1）;
figure（4）;
subplot（221）;stem（fabs（X8k16）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（8a）?16點DFT[x_8（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=32;n=0:N-1;?%對于N=16的情況
nT?=?n*T;
x8n=cos（8*pi*nT）+cos（16*pi*nT）+cos（20*pi*nT）
X8k32=fft（x8n32）;
N=32;
f=2/N*（0:N-1）;
subplot（222）;stem（fabs（X8k32）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（8a）?32點DFT[x_8（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=64;n=0:N-1;?%對于N=16的情況
nT?=?n*T;
x8n=cos（8*pi*nT）+cos（16*pi*nT）+cos（20*pi*nT）
X8k64=fft（x8n64）;
N=64;
f=2/N*（0:N-1）;
subplot（223）;stem（fabs（X8k64）‘.‘）;?%繪制8點DFT的幅頻特性圖
title（‘（8a）?64點DFT[x_8（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-06-16?13:52??數字信號處理實驗\
?????目錄???????????0??2018-12-31?10:14??數字信號處理實驗\實驗一\
?????文件????????1390??2018-12-16?19:15??數字信號處理實驗\實驗一\shiyan1_1.asv
?????文件????????2718??2018-12-16?19:22??數字信號處理實驗\實驗一\shiyan1_1.m
?????文件??????184320??2018-12-16?19:40??數字信號處理實驗\實驗一\實驗一截圖.doc
?????文件??????193024??2018-12-16?19:21??數字信號處理實驗\實驗一\實驗一：用FFT對信號作頻譜分析_實驗報告.doc
?????目錄???????????0??2018-12-31?10:14??數字信號處理實驗\實驗三\
?????文件?????????269??2018-12-16?19:56??數字信號處理實驗\實驗三\myplot.m
?????文件?????????132??2018-12-27?11:14??數字信號處理實驗\實驗三\shiyan3_1.m
?????文件?????????149??2018-12-27?13:57??數字信號處理實驗\實驗三\shiyan3_2.m
?????文件?????????368??2018-12-27?14:35??數字信號處理實驗\實驗三\shiyan3_3.m
?????文件?????????224??2018-12-16?19:54??數字信號處理實驗\實驗三\tplot.m
?????文件????????1317??2018-12-17?13:22??數字信號處理實驗\實驗三\xtg.m
?????文件???????73728??2018-12-17?13:17??數字信號處理實驗\實驗三\實驗三--FIR數字濾波器設計與軟件實現.doc
?????文件??????201216??2018-12-27?14:37??數字信號處理實驗\實驗三\實驗三截圖.doc
?????目錄???????????0??2018-12-31?10:14??數字信號處理實驗\實驗二\
?????文件?????????559??2018-12-16?19:54??數字信號處理實驗\實驗二\mstg.m
?????文件?????????269??2018-12-16?19:56??數字信號處理實驗\實驗二\myplot.m
?????文件?????????873??2018-12-16?20:25??數字信號處理實驗\實驗二\shiyan2.m
?????文件?????????224??2018-12-16?19:54??數字信號處理實驗\實驗二\tplot.m
?????文件??????203776??2018-12-16?20:26??數字信號處理實驗\實驗二\實驗二截圖.doc
?????文件??????228238??2018-12-16?20:26??數字信號處理實驗\實驗二\實驗二：IIR數字濾波器設計及軟件實現.doc
?????目錄???????????0??2018-12-31?10:14??數字信號處理實驗\實驗四\
?????文件?????????287??2018-12-27?15:57??數字信號處理實驗\實驗四\shiyan4_1.m
?????文件?????????149??2018-12-27?16:22??數字信號處理實驗\實驗四\shiyan4_2.m
?????文件???????79872??2018-12-27?16:23??數字信號處理實驗\實驗四\實驗四截圖.doc