数字信号处理高西全实验课matlab程序

四川大学通信工程的专业课《数字信号处理》实验代码，matlab代码。配套书籍为《数字信号处理》高西全版，代码文件按“p123.m”格式命名，标明该段代码是书本第几页的实现代码。可供参考。

%实验二：时域采样与频域采样
%时域采样

%原信号
Tp=64/1000;%观测时间64ms 
Fs=1000;  T=1/Fs; 
M=Tp*Fs;%取样点数64点
n=0:M-1;t=n*T;
A=444.128;a=pi*50*2^0.5;w=pi*50*2^0.5;
xnt=A*exp（-a*t）.*sin（w*t）;
Xk=T*fft（xntM）;   %M=64点FFT[xnt]
figure（1）
subplot（421）;  plot（nxnt）;title（‘原信号波形‘）;
k=0:M-1;  wk=k/Tp; 
subplot（422）;plot（wkabs（Xk））;title（‘Fs=1KHz幅频特性‘）;
xlabel（‘f/Hz‘）;ylabel（‘幅度（H1（jf））‘）;

%采样频率1KHz
Fs=1000;T=1/Fs; %采样频率1KHz
M=Tp*Fs;%取样点数64点
n=0:M-1;t=n*T;
A=444.128;a=pi*50*2^0.5;w=pi*50*2^0.5;
xnt=A*exp（-a*t）.*sin（w*t）;
Xk=T*fft（xntM）;   %M=64点FFT[xnt]
subplot（423）;  stem（nxnt‘.‘）;title（‘Fs=1KHz采样序列‘）;
k=0:M-1; wk=k/Tp;
subplot（424）;plot（wkabs（Xk））;title（‘Fs=1KHz幅频特性‘）;
xlabel（‘f/Hz‘）;ylabel（‘幅度（H1（jf））‘）;

%采样频率300Hz
Fs=300;T=1/Fs;  M=Tp*Fs;n=0:M-1;t=n*T;
A=444.128;a=pi*50*2^0.5;w=pi*50*2^0.5;
xnt=A*exp（-a*t）.*sin（w*t）;
Xk=T*fft（xntM）; 
subplot（425）; stem（nxnt‘.‘）;title（‘Fs=300Hz采样序列‘）;
k=0:M-1; wk=k/Tp;
subplot（426）;plot（wkabs（Xk））;title（‘Fs=300Hz幅频特性‘）;
xlabel（‘f/Hz‘）;ylabel（‘（H2（jf））‘）;

%采样频率200Hz
Fs=200;T=1/Fs;  M=Tp*Fs;n=0:M-1;t=n*T;
A=444.128;a=pi*50*2^0.5;w=pi*50*2^0.5;
xnt=A*exp（-a*t）.*sin（w*t）;
Xk=T*fft（xntM）; 
subplot（427）;  stem（nxnt‘.‘）;title（‘Fs=200Hz采样序列‘）;
k=0:M-1;  wk=k/Tp;
subplot（428）;plot（wkabs（Xk））;title（‘Fs=200Hz幅频特性‘）;
xlabel（‘f/Hz‘）;ylabel（‘（H3（jf））‘）;

%=========================================================================%

%频域采样
M=27;N=32;n=0:M;
xn=（n>=0&n<=13）.*（n+1）+（n>=14&n<=26）.*（27-n）+（n<0|n>26）.*0;     %产生x（n）
Xk=fft（xn1024）;                                  %1024点FFT[x（n）] 
k=0:1023;
wk=2*k/1024;         %连续频谱图的横坐标取值
figure（2）
subplot（321）; plot（wkabs（Xk））; title（‘FT[x（n）]‘）;
xlabel（‘w/\pi‘）;ylabel（‘|X（e^j^\w）|‘）;axis（[010200]）;
subplot（322）; stem（nxn‘.‘）; title（‘x（n）‘）;
xlabel（‘n‘）;ylabel（‘x（n）‘）;axis（[032020]）

k=0:N-1;                %离散频谱图的横坐标取值
X32k=fft（xn32）;                                  %32点FFT[x（n）]
x32n=ifft（X32k）;                                  %32点IFFT[X32（k）]得到x32（n）
subplot（325）;stem（kabs（X32k）‘.‘）;title（‘32点频域采样‘）;
xlabel（‘k‘）;ylabel（‘|X_3_2（k）|‘）;axis（[0160200]）
n1=0:N-1;
subplot（326）;stem（n1x32n‘.‘）;title（‘32IDFT[X_3_2（k）]‘）;
k=0:N/2-1;             %离散频谱图的横坐标取值 
X16k=X32k（1:2:N）;       %隔点抽取X32（k）得到X16（k）
x16n=ifft（X16kN/2）;  %16点IFFT[X16（k）]得到x16（n）
subplot（323）;stem（kabs（X16k）‘.‘）;title（‘16点频域采样‘）;
xlabel（‘k‘）;ylabel（‘|X_1_6（k）|‘）;axis（[080200]）
n1=0:N/2-1;           
subplot（324）;stem（n1x16n‘.‘）;title（‘16IDFT[X_1_6（k）]‘）;
xlabel（‘n‘）;ylabel（‘x_1_6（n）‘）;axis（[032020]）
xlabel（‘n‘）;ylabel（‘x_3_2（n）‘）;axis（[032020]）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-01-19 19:01  shuzixinhaochulishiyan\
     文件        2793  2017-12-27 15:11  shuzixinhaochulishiyan\exp2.m
     文件        2728  2017-12-27 21:04  shuzixinhaochulishiyan\exp3.m
     文件        4062  2017-12-27 20:48  shuzixinhaochulishiyan\exp3new.m
     文件        1844  2017-12-27 22:01  shuzixinhaochulishiyan\exp4.m
     文件        1825  2017-12-04 10:32  shuzixinhaochulishiyan\exp42.m
     文件        1127  2017-12-28 01:38  shuzixinhaochulishiyan\exp6.m
     目录           0  2018-01-19 19:01  shuzixinhaochulishiyan\matlab章节仿真\
     文件       11167  2017-12-10 23:01  shuzixinhaochulishiyan\matlab章节仿真\AH.ps
     文件         264  2017-12-10 22:54  shuzixinhaochulishiyan\matlab章节仿真\p165.m
     文件         270  2017-12-10 22:54  shuzixinhaochulishiyan\matlab章节仿真\p171.m
     文件         263  2017-11-01 10:03  shuzixinhaochulishiyan\matlab章节仿真\p173.m
     文件         258  2017-11-01 10:26  shuzixinhaochulishiyan\matlab章节仿真\p177.m
     文件         279  2017-12-10 22:59  shuzixinhaochulishiyan\matlab章节仿真\P177a.m
     文件         134  2017-12-10 23:08  shuzixinhaochulishiyan\matlab章节仿真\p177b.m
     文件         231  2017-12-10 23:07  shuzixinhaochulishiyan\matlab章节仿真\p179.m
     文件         207  2017-11-07 21:05  shuzixinhaochulishiyan\matlab章节仿真\p180.m
     文件         482  2017-11-19 11:06  shuzixinhaochulishiyan\matlab章节仿真\p187.m
     文件         552  2017-11-07 22:46  shuzixinhaochulishiyan\matlab章节仿真\p193.m
     文件         245  2017-11-07 22:49  shuzixinhaochulishiyan\matlab章节仿真\p195.m
     文件         175  2017-11-07 22:54  shuzixinhaochulishiyan\matlab章节仿真\p196.m
     文件         332  2017-11-19 14:12  shuzixinhaochulishiyan\matlab章节仿真\p218.m
     文件         452  2017-11-19 14:59  shuzixinhaochulishiyan\matlab章节仿真\p219.m
     文件         344  2017-11-19 15:07  shuzixinhaochulishiyan\matlab章节仿真\p220.m
     文件         839  2017-11-19 15:29  shuzixinhaochulishiyan\matlab章节仿真\p225.m
     文件         724  2017-11-19 15:37  shuzixinhaochulishiyan\matlab章节仿真\p226.m
     文件         427  2017-11-19 16:25  shuzixinhaochulishiyan\matlab章节仿真\p232.m
     文件         467  2017-11-19 16:31  shuzixinhaochulishiyan\matlab章节仿真\p232b.m
     文件         301  2017-10-16 07:56  shuzixinhaochulishiyan\matlab章节仿真\p69.m
     文件         276  2017-10-16 08:00  shuzixinhaochulishiyan\matlab章节仿真\p71.m
     文件         667  2017-10-16 09:52  shuzixinhaochulishiyan\matlab章节仿真\p82.m
............此处省略27个文件信息