数字信号处理——基于计算机的方法第四版MATLAB程序

本书是数字信号领域的经典教材Digital Signal Processing: A Computer-BasedApproach, FourthEdition的中文翻译版，内容涵盖了信号与信号处理、时域中的离散时间信号、频域中的离散时间信号、离散时间系统、有限长离散变换、z变换、变换域中的LTI离散时间系统、数字滤波器结构、IIR数字滤波器设计、FIR数字滤波器设计、DSP算法实现等方面。本书的特点是，在讲解上述内容的同时，给出了大量简单而实用的例子，并用MATLAB程序进行了验证，同时提供了大量的高质量习题和仿真练习。

%
% Example of sound aliasing
%
clear all;

fc = 1.8e3;               % tone freq
Fs = 1e3*8;               % sample freq
Ts = 1/Fs;                % sample time
timelen = 3;              % sample file length

n = 0:ceil（timelen/Ts）-1;
x = 0.9*cos（2*pi*fc*n*Ts）;

disp（[‘tone freq = ‘num2str（fc）‘ sampling frequency = ‘num2str（Fs）]）;
wavplay（xFs）;

%
% downsampling
%
d2=3/4;
disp（‘press any key to continue‘）;
pause
disp（[‘tone freq = ‘num2str（fc）‘ sampling frequency = ‘num2str（Fs*d2）]）;
x_d2 = 0.9*cos（2*pi*fc*n*Ts/d2）;
wavplay（x_d2（1:floor（length（x_d2）*d2））Fs*d2）;

d3=1/3;
disp（‘press any key to continue‘）;
pause
disp（[‘tone freq = ‘num2str（fc）‘ sampling frequency = ‘num2str（Fs*d3）]）;
x_d3 = 0.9*cos（2*pi*fc*n*Ts/d3）;
wavplay（x_d3（1:floor（length（x_d3）*d3））Fs*d3）;

%
% write wav files
%
%wavwrite（xFs16‘output01.wav‘）;
%wavwrite（x_d2（1:floor（length（x_d2）*d2））Fs*d216‘output02.wav‘）;
%wavwrite（x_d3（1:floor（length（x_d3）*d3））Fs*d316‘output03.wav‘）;

%
% show plot
%
Xw = fft（xlength（x））;
L = 0:length（Xw）-1;
figure（1）;
subplot（211）;
plot（Fs*L/length（Xw）/1e3 20*log10（abs（Xw）））; grid;
xlabel（‘Frequency KHz‘）; ylabel（‘Gain dB‘）;
title（[‘Tone Freq = ‘num2str（fc）‘ Sampling Freq = ‘num2str（Fs）]）;
subplot（212）;
stem（1:30x（1:30））;
xlabel（‘n‘）; ylabel（‘Magnitude‘）;

Xw_d2 = fft（x_d2length（x_d2））;
L_d2 = 0:length（Xw_d2）-1;
figure（2）;
subplot（211）;
plot（（Fs*d2）*L_d2/length（Xw_d2）/1e3 20*log10（abs（Xw_d2）））; grid;
xlabel（‘Frequency KHz‘）; ylabel（‘Gain dB‘）;
title（[‘Tone Freq = ‘num2str（fc）‘ Sampling Freq = ‘num2str（Fs*d2）]）;
subplot（212）;
stem（1:30x_d2（1:30））;
xlabel（‘n‘）; ylabel（‘Magnitude‘）;

Xw_d3 = fft（x_d3length（x_d3））;
L_d3 = 0:length（Xw_d3）-1;
figure（3）;
subplot（211）;
plot（（Fs*d3）*L_d3/length（Xw_d3）/1e3 20*log10（abs（Xw_d3）））; grid;
xlabel（‘Frequency KHz‘）; ylabel（‘Gain dB‘）;
title（[‘Tone Freq = ‘num2str（fc）‘ Sampling Freq = ‘num2str（Fs*d3）]）;
subplot（212）;
stem（1:30x_d3（1:30））;
xlabel（‘n‘）; ylabel（‘Magnitude‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\
     文件          59  2010-07-16 08:30  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\autorun.inf
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\applications\
     文件     1310943  2010-07-15 13:29  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\applications\Applications2.pdf
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\
     文件       13584  2010-07-20 09:20  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\demos.htm
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_demo_1\
     文件      228842  2004-09-18 06:16  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_demo_1\blackwhite.png
     文件      463972  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_demo_1\tulips.PNG
     文件      225814  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_demo_1\tulips_BLUE.PNG
     文件      221998  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_demo_1\tulips_GREEN.PNG
     文件      210874  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_demo_1\tulips_RED.PNG
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_enhance\
     文件      138349  2004-09-19 08:57  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_enhance\filtered_image.png
     文件      155232  2004-09-19 08:58  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_enhance\noisy_image.png
     文件       17408  2004-09-20 09:04  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\image_enhance\Thumbs.db
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\other\
     文件     3102964  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\other\nonLinear_bells.wav
     文件     2431792  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\other\phyMod_string_major.wav
     文件     1232452  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\other\specMod_major.wav
     文件     1904676  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\other\wavetable_piano_oldsong.wav
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\twinkle\
     文件     2123088  2004-09-17 07:40  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\twinkle\nonLinear_horn_twinkle.wav
     文件     2044296  2004-09-17 07:41  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\twinkle\phyMod_guitar_twinkle.wav
     文件     1980708  2004-09-17 07:41  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\twinkle\specMod_twinkle.wav
     文件     1351716  2004-09-17 07:39  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\music_synthesis\twinkle\wavetable_piano_twinkle.wav
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\sound_aliasing\
     文件        2066  2004-09-20 08:12  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\sound_aliasing\aliasing.m
     文件        2063  2004-09-17 08:17  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）\media\demos\sound_aliasing\aliasing.m.bak
............此处省略284个文件信息