数字信号处理——基于计算机的方法第四版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-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediaapplications
     文件     1310943  2010-07-15 13:29  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediaapplicationsApplications2.pdf
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademos
     文件       13584  2010-07-20 09:20  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosdemos.htm
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_demo_1
     文件      228842  2004-09-18 06:16  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_demo_1lackwhite.png
     文件      463972  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_demo_1 ulips.PNG
     文件      225814  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_demo_1 ulips_BLUE.PNG
     文件      221998  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_demo_1 ulips_GREEN.PNG
     文件      210874  2004-09-18 03:18  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_demo_1 ulips_RED.PNG
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_enhance
     文件      138349  2004-09-19 08:57  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_enhancefiltered_image.png
     文件      155232  2004-09-19 08:58  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_enhance
oisy_image.png
     文件       17408  2004-09-20 09:04  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosimage_enhanceThumbs.db
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesis
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesisother
     文件     3102964  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesisother
onLinear_bells.wav
     文件     2431792  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesisotherphyMod_string_major.wav
     文件     1232452  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesisotherspecMod_major.wav
     文件     1904676  2004-09-18 06:21  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesisotherwavetable_piano_oldsong.wav
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesis winkle
     文件     2123088  2004-09-17 07:40  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesis winkle
onLinear_horn_twinkle.wav
     文件     2044296  2004-09-17 07:41  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesis winklephyMod_guitar_twinkle.wav
     文件     1980708  2004-09-17 07:41  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesis winklespecMod_twinkle.wav
     文件     1351716  2004-09-17 07:39  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademosmusic_synthesis winklewavetable_piano_twinkle.wav
     目录           0  2011-11-10 10:28  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademossound_aliasing
     文件        2066  2004-09-20 08:12  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademossound_aliasingaliasing.m
     文件        2063  2004-09-17 08:17  14597-数字信号处理——基于计算机的方法（第四版）（英文改编版）mediademossound_aliasingaliasing.m.bak
............此处省略284个文件信息