《DSP using MATLAB》示例Example5.7

代码：

x = [1, 1, 1, 1, zeros(1,4)];

N = 8;                                                 % zero-padding operation

X_DFT = dft(x,N);                                      % DFT of x(n)

magX_DFT = abs(X_DFT)

phaX_DFT = angle(X_DFT)*180/pi                         % degrees

realX_DFT = real(X_DFT); imagX_DFT = imag(X_DFT);

angX_DFT = angle(X_DFT);                               % radias

n = 0:(N - 1);

k = 0:1000; w = (pi/500)*k;                            % [0,2pi] axis divided into 501 points.

%k = 0:500; w = (pi/500)*k;                            % [0,pi] axis divided into 501 points.

X_DTFT = x * (exp(-j*pi/500)) .^ (n'*k);               % DTFT of x(n)

magX_DTFT = abs(X_DTFT); angX_DTFT = angle(X_DTFT); realX_DTFT = real(X_DTFT); imagX_DTFT = imag(X_DTFT);

figure('NumberTitle', 'off', 'Name', 'Exameple5.7 x sequence')

set(gcf,'Color','white');

stem(n,x); title('x sequence'); axis([0,N,-0.5,1.5]);

xlabel('n'); ylabel('x(n)'); grid on;

%% --------------------------------------------------------------

%%        START X_DTFT's  mag ang real imag

%% --------------------------------------------------------------

figure('NumberTitle', 'off', 'Name', 'X_DTFT its Magnitude and Angle, Real and Imaginary Part');

set(gcf,'Color','white');

subplot(2,2,1); plot(w/pi,magX_DTFT); grid on;  % axis([-2,2,0,15]);

title('Magnitude Part');

xlabel('frequency in \pi units'); ylabel('Magnitude  |X\_DTFT|');

subplot(2,2,3); plot(w/pi, angX_DTFT*180/pi); grid on;  % axis([-2,2,-1,1]);

title('Angle Part');

xlabel('frequency in \pi units'); ylabel('Degrees'); %axis([-200,200,0,2]);

subplot('2,2,2'); plot(w/pi, realX_DTFT); grid on;

title('Real Part');

xlabel('frequency in \pi units'); ylabel('Real');

subplot('2,2,4'); plot(w/pi, imagX_DTFT); grid on;

title('Imaginary Part');

xlabel('frequency in \pi units'); ylabel('Imaginary');

%% --------------------------------------------------------------

%%        END X_DTFT's  mag ang real imag

%% --------------------------------------------------------------

%% --------------------------------------------------------------

%%        START X_DFT's  mag ang real imag

%% --------------------------------------------------------------

figure('NumberTitle', 'off', 'Name', 'Example5.7 X_DFT its Magnitude and Angle');

set(gcf,'Color','white');

subplot(2,1,1); stem(n,magX_DFT); hold on; plot(4*w/pi,magX_DTFT,'--'); hold off;

grid on;   axis([-0.5,8.2,-1,5]);

title('Magnitude Part of the DFT: N = 8');

xlabel('k'); ylabel('Magnitude  |X\_DFT|');

subplot(2,1,2); stem(n, phaX_DFT); hold on; plot(4*w/pi,angX_DTFT*180/pi,'--'); hold off;

grid on;   axis([-0.5,8.2,-200,200]);

title('Angle Part of the DFT: N = 8');

xlabel('k'); ylabel('Degrees'); %axis([-200,200,0,2]);

%subplot('2,2,2'); stem(n, realX_DFT); grid on;

%title('Real Part');

%xlabel('frequency in \pi units'); ylabel('Real');

%subplot('2,2,4'); stem(n, imagX_DFT); grid on;

%title('Imaginary Part');

%xlabel('frequency in \pi units'); ylabel('Imaginary');

%% --------------------------------------------------------------

%%        END X_DFT's  mag ang real imag

%% --------------------------------------------------------------

　　结果：

将序列末尾补12个零，长度达到16位

代码：

x = [1, 1, 1, 1, zeros(1,12)];                         % append 12 zeros to the end of x(n)

N = 16;                                                % zero-padding operation

X_DFT = dft(x,N);                                      % DFT of x(n)

magX_DFT = abs(X_DFT)

phaX_DFT = angle(X_DFT)*180/pi                         % degrees

realX_DFT = real(X_DFT); imagX_DFT = imag(X_DFT);

angX_DFT = angle(X_DFT);                               % radias

n = 0:(N - 1);

k = 0:1000; w = (pi/500)*k;                            % [0,2pi] axis divided into 501 points.

%k = 0:500; w = (pi/500)*k;                            % [0,pi] axis divided into 501 points.

X_DTFT = x * (exp(-j*pi/500)) .^ (n'*k);               % DTFT of x(n)

magX_DTFT = abs(X_DTFT); angX_DTFT = angle(X_DTFT); realX_DTFT = real(X_DTFT); imagX_DTFT = imag(X_DTFT);

figure('NumberTitle', 'off', 'Name', 'Exameple5.7 x sequence')

set(gcf,'Color','white');

stem(n,x); title('x sequence, N = 16'); axis([0,N,-0.5,1.5]);

xlabel('n'); ylabel('x(n)'); grid on;

%% --------------------------------------------------------------

%%        START X_DTFT's  mag ang real imag

%% --------------------------------------------------------------

figure('NumberTitle', 'off', 'Name', 'X_DTFT its Magnitude and Angle, Real and Imaginary Part');

set(gcf,'Color','white');

subplot(2,2,1); plot(w/pi,magX_DTFT); grid on;  % axis([-2,2,0,15]);

title('Magnitude Part');

xlabel('frequency in \pi units'); ylabel('Magnitude  |X\_DTFT|');

subplot(2,2,3); plot(w/pi, angX_DTFT*180/pi); grid on;  % axis([-2,2,-1,1]);

title('Angle Part');

xlabel('frequency in \pi units'); ylabel('Degrees'); %axis([-200,200,0,2]);

subplot('2,2,2'); plot(w/pi, realX_DTFT); grid on;

title('Real Part');

xlabel('frequency in \pi units'); ylabel('Real');

subplot('2,2,4'); plot(w/pi, imagX_DTFT); grid on;

title('Imaginary Part');

xlabel('frequency in \pi units'); ylabel('Imaginary');

%% --------------------------------------------------------------

%%        END X_DTFT's  mag ang real imag

%% --------------------------------------------------------------

%% --------------------------------------------------------------

%%        START X_DFT's  mag ang real imag

%% --------------------------------------------------------------

figure('NumberTitle', 'off', 'Name', 'Example5.7 X_DFT its Magnitude and Angle');

set(gcf,'Color','white');

subplot(2,1,1); stem(n,magX_DFT); hold on; plot(8*w/pi,magX_DTFT,'--'); hold off;

grid on;   axis([-0.5,16.2,-1,5]);

title('Magnitude Part of the DFT: N = 16');

xlabel('k'); ylabel('Magnitude  |X\_DFT|');

subplot(2,1,2); stem(n, phaX_DFT); hold on; plot(8*w/pi,angX_DTFT*180/pi,'--'); hold off;

grid on;   axis([-0.5,16.2,-200,200]);

title('Angle Part of the DFT: N = 16');

xlabel('k'); ylabel('Degrees'); %axis([-200,200,0,2]);

%subplot('2,2,2'); stem(n, realX_DFT); grid on;

%title('Real Part');

%xlabel('frequency in \pi units'); ylabel('Real');

%subplot('2,2,4'); stem(n, imagX_DFT); grid on;

%title('Imaginary Part');

%xlabel('frequency in \pi units'); ylabel('Imaginary');

%% --------------------------------------------------------------

%%        END X_DFT's  mag ang real imag

%% --------------------------------------------------------------

　　结果：

继续补零，长度N=128，

代码：

x = [1, 1, 1, 1, zeros(1,124)];                        % append 124 zeros to the end of x(n)

N = 128;                                               % zero-padding operation

X_DFT = dft(x,N);                                      % DFT of x(n)

magX_DFT = abs(X_DFT)

phaX_DFT = angle(X_DFT)*180/pi                         % degrees

realX_DFT = real(X_DFT); imagX_DFT = imag(X_DFT);

angX_DFT = angle(X_DFT);                               % radias

n = 0:(N - 1);

k = 0:1000; w = (pi/500)*k;                            % [0,2pi] axis divided into 501 points.

%k = 0:500; w = (pi/500)*k;                            % [0,pi] axis divided into 501 points.

X_DTFT = x * (exp(-j*pi/500)) .^ (n'*k);               % DTFT of x(n)

magX_DTFT = abs(X_DTFT); angX_DTFT = angle(X_DTFT); realX_DTFT = real(X_DTFT); imagX_DTFT = imag(X_DTFT);

figure('NumberTitle', 'off', 'Name', 'Exameple5.7 x sequence')

set(gcf,'Color','white');

stem(n,x); title('x sequence, N = 128'); axis([0,N,-0.5,1.5]);

xlabel('n'); ylabel('x(n)'); grid on;

%% --------------------------------------------------------------

%%        START X_DTFT's  mag ang real imag

%% --------------------------------------------------------------

figure('NumberTitle', 'off', 'Name', 'X_DTFT its Magnitude and Angle, Real and Imaginary Part');

set(gcf,'Color','white');

subplot(2,2,1); plot(w/pi,magX_DTFT); grid on;  % axis([-2,2,0,15]);

title('Magnitude Part');

xlabel('frequency in \pi units'); ylabel('Magnitude  |X\_DTFT|');

subplot(2,2,3); plot(w/pi, angX_DTFT*180/pi); grid on;  % axis([-2,2,-1,1]);

title('Angle Part');

xlabel('frequency in \pi units'); ylabel('Degrees'); %axis([-200,200,0,2]);

subplot('2,2,2'); plot(w/pi, realX_DTFT); grid on;

title('Real Part');

xlabel('frequency in \pi units'); ylabel('Real');

subplot('2,2,4'); plot(w/pi, imagX_DTFT); grid on;

title('Imaginary Part');

xlabel('frequency in \pi units'); ylabel('Imaginary');

%% --------------------------------------------------------------

%%        END X_DTFT's  mag ang real imag

%% --------------------------------------------------------------

%% --------------------------------------------------------------

%%        START X_DFT's  mag ang real imag

%% --------------------------------------------------------------

figure('NumberTitle', 'off', 'Name', 'Example5.7 X_DFT its Magnitude and Angle');

set(gcf,'Color','white');

subplot(2,1,1); stem(n,magX_DFT); hold on; plot(64*w/pi,magX_DTFT,'--'); hold off;

grid on;   axis([-0.5,128.2,-1,5]);

title('Magnitude Part of the DFT: N = 128');

xlabel('k'); ylabel('Magnitude  |X\_DFT|');

subplot(2,1,2); stem(n, phaX_DFT); hold on; plot(64*w/pi,angX_DTFT*180/pi,'--'); hold off;

grid on;   axis([-0.5,128.2,-200,200]);

title('Angle Part of the DFT: N = 128');

xlabel('k'); ylabel('Degrees'); %axis([-200,200,0,2]);

%subplot('2,2,2'); stem(n, realX_DFT); grid on;

%title('Real Part');

%xlabel('frequency in \pi units'); ylabel('Real');

%subplot('2,2,4'); stem(n, imagX_DFT); grid on;

%title('Imaginary Part');

%xlabel('frequency in \pi units'); ylabel('Imaginary');

%% --------------------------------------------------------------

%%        END X_DFT's  mag ang real imag

%% --------------------------------------------------------------

　　运行结果：

《DSP using MATLAB》示例Example5.7的更多相关文章

DSP using MATLAB 示例Example3.21
代码: % Discrete-time Signal x1(n) % Ts = 0.0002; n = -25:1:25; nTs = n*Ts; Fs = 1/Ts; x = exp(-1000*a ...
DSP using MATLAB 示例 Example3.19
代码: % Analog Signal Dt = 0.00005; t = -0.005:Dt:0.005; xa = exp(-1000*abs(t)); % Discrete-time Signa ...
DSP using MATLAB示例Example3.18
代码: % Analog Signal Dt = 0.00005; t = -0.005:Dt:0.005; xa = exp(-1000*abs(t)); % Continuous-time Fou ...
DSP using MATLAB 示例Example3.23
代码: % Discrete-time Signal x1(n) : Ts = 0.0002 Ts = 0.0002; n = -25:1:25; nTs = n*Ts; x1 = exp(-1000 ...
DSP using MATLAB 示例Example3.22
代码: % Discrete-time Signal x2(n) Ts = 0.001; n = -5:1:5; nTs = n*Ts; Fs = 1/Ts; x = exp(-1000*abs(nT ...
DSP using MATLAB 示例Example3.17
DSP using MATLAB示例Example3.16
代码: b = [0.0181, 0.0543, 0.0543, 0.0181]; % filter coefficient array b a = [1.0000, -1.7600, 1.1829, ...
DSP using MATLAB 示例 Example3.15
上代码: subplot(1,1,1); b = 1; a = [1, -0.8]; n = [0:100]; x = cos(0.05*pi*n); y = filter(b,a,x); figur ...
DSP using MATLAB 示例 Example3.13
上代码: w = [0:1:500]*pi/500; % freqency between 0 and +pi, [0,pi] axis divided into 501 points. H = ex ...
DSP using MATLAB 示例 Example3.12
用到的性质代码: n = -5:10; x = sin(pi*n/2); k = -100:100; w = (pi/100)*k; % freqency between -pi and +pi , ...

随机推荐

tensrflow python [defunct]
在ubuntu上面安装了GPU版本的tensorflow后,很容易碰到zombie thread 的问题,无法正常关闭tensorflow的线程,用ps aux|grep python可以看到 pyt ...
<<< tomcat启动报错StandardServer.await: create[8005]
启动tomcat的时候出现异常严重: StandardServer.await: create[8005]: java.net.BindException: Address already in u ...
图解JAVA对象的创建过程
前面几篇博文分别介绍了JAVA的Class文件格式.JVM的类加载机制和JVM的内存模型,这里就索性把java对象的创建过程一并说完,这样java对象的整个创建过程就基本上说明白了(当然你要有基础才能 ...
jquery实现简单瀑布流布局
jquery实现简单瀑布流布局是开头都会说的原理瀑布流布局有两种,一种是固定列,一种是非固定列.在此主要记述第一种的实现. 固定列的特征是:无论页面如何缩放,每行的总列数都一致. 一行4列的瀑布流 ...
HBase中批量修改
先随便写写..做个随笔记录使用Rest连接操作Hbase.. 是微软提供的 Microsoft.Hbase.Client 类库.. 版本是0.4.1.0 一直知道 client.StoreCe ...
Spring系列之AOP
一.什么是AOPAOP(Aspect-OrientedProgramming,面向方面编程),可以说是OOP(Object-Oriented Programing,面向对象编程)的补充和完善.OOP引 ...
Android ViewPager sharedpreferences
http://www.cnblogs.com/dwinter/archive/2012/02/27/AndroidViewPager%E5%A4%9A%E9%A1%B5%E9%9D%A2%E6%BB% ...
tyvj1172 自然数拆分Lunatic版
背景话说小小鱼看了P1171(自然数拆分)之后感觉异常不爽,于是异常邪恶地将题目加强. 描述输入自然数n,然后将其拆分成由若干数相加的形式,参与加法运算的数可以重复. 输入格式输入只有一个整数n ...
UIScrollView的代理（delegate）
很多时候,我们想在UIScrollView正在滚动或滚动到某个位置或者停止滚动时做一些特定的操作要想完成上述功能,前提条件就是能够监听到UIScrollView的整个滚动过程当UIScr ...
URL、URN、URI的区别？
URL.URN.URI区别既然Web应用程序的文件等资源是放在服务器上,而服务器是因特网(Internet)上的主机,当然必须要有个方法,告诉浏览器到哪里取得文件等资源.通常会听到有人这么说:“你要 ...

《DSP using MATLAB》示例Example5.7

《DSP using MATLAB》示例Example5.7的更多相关文章

随机推荐

热门专题