代码：

%% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

%%            Output Info about this m-file

fprintf('\n***********************************************************\n');

fprintf('        <DSP using MATLAB> Problem 5.38 \n\n');

banner();

%% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

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

%       1  x(n) = [2sin(0.04pi*n)+5cos(0.08pi*n)]

%                    N=40    N=50   N=60

%

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

N = 40;

nn1 = [0:N-1];

xx1 = 2*sin(0.04*pi*nn1)+5*cos(0.08*pi*nn1);

NN1 = length(xx1);                                % length 

%m = mod_1(nn1, NN1);

%x = [xx1 zeros(1, 0)];                        % padding zeros

%n = [nn1 max(nn1)+1:max(nn1)+6];

x = xx1;

n = nn1;

str = fprintf('P5.38.1 x(n) N=%d', N);

figure('NumberTitle', 'off', 'Name', str);

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

subplot(2,1,1); stem(nn1, xx1);

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

title('x(n) ori sequence');  grid on;

subplot(2,1,2); stem(n, x);

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

title('x(n) padding 0 zeros');  grid on;

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

%%                 DTFT X(w) of xn sequence, w=[0:2pi],

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

MM = 500;

[Xw_DTFT, w] = dtft1(x, n, MM);

 magXw_DTFT = abs(Xw_DTFT);   angXw_DTFT = angle(Xw_DTFT)/pi;

realXw_DTFT = real(Xw_DTFT); imagXw_DTFT = imag(Xw_DTFT);

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

%%        START X_DTFT's  mag ang real imag

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

figure('NumberTitle', 'off', 'Name', 'P5.38.1 X(w) DTFT of x(n) N=40');

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

subplot(2,2,1); plot(w/pi,magXw_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, angXw_DTFT); grid on;  % axis([-2,2,-1,1]);

title('Angle Part');

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

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

title('Real Part');

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

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

title('Imaginary Part');

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

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

%%        END X_DTFT's  mag ang real imag

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

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

%%                 DFT(k) of xn sequence, k=[0:N-1]

%%                      w=2pi*k/N       k=Nw/(2pi)

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

N1 = length(x);

k1 = [0 : N1-1];

%k2 = [-N : N-1];

%k3 = [-N/2 : N/2];

Xk_DFT = dft(x, N1);                                         % DFT

    magXk_DFT = abs( [ Xk_DFT ] );                          % DFT magnitude

    angXk_DFT = angle( [Xk_DFT] )/pi;                       % DFT angle

   realXk_DFT = real(Xk_DFT); imagXk_DFT = imag(Xk_DFT);

figure('NumberTitle', 'off', 'Name', 'P5.38.1 DFT(k) of x(n) N=40')

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

subplot(2,1,1); stem(k1, magXk_DFT); hold on; plot(N1*w/(2*pi), magXw_DTFT,'r--'); hold off;

%axis([-N/2, N/2, -0.5, 50.5]);

xlabel('k'); ylabel('magnitude(k)');

title('DFT magnitude of x(n)');  grid on;

subplot(2,1,2); stem(k1, angXk_DFT);  hold on; plot(N1*w/(2*pi), angXw_DTFT,'r--'); hold off;

%axis([-N/2, N/2, -0.5, 50.5]);

xlabel('k'); ylabel('angle(k)');

title('DFT angle of x(n)'); grid on;

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

%%                 DTFT X(w) of xn sequence, w=[-2pi:2pi],

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

MM = 500;

[Xw_DTFT, w] = dtft2(x, n, MM);

 magXw_DTFT = abs(Xw_DTFT);   angXw_DTFT = angle(Xw_DTFT)/pi;

realXw_DTFT = real(Xw_DTFT); imagXw_DTFT = imag(Xw_DTFT);

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

%%        START X_DTFT's  mag ang real imag

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

figure('NumberTitle', 'off', 'Name', 'P5.38.1 X(w) DTFT of x(n) N=40');

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

subplot(2,2,1); plot(w/pi,magXw_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, angXw_DTFT); grid on;  % axis([-2,2,-1,1]);

title('Angle Part');

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

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

title('Real Part');

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

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

title('Imaginary Part');

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

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

%%        END X_DTFT's  mag ang real imag

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

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

%%                 DFT(k) of xn sequence, k=[-N:N-1]

%%                    w=2pi*k/N       k=Nw/(2pi)

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

N2 = length(x);

%k1 = [0 : length(x)-1];

k2 = [-N2 : N2-1];

%k3 = [-N/2 : N/2];

Xk_DFT = dft(x, N2);                                                % DFT

    magXk_DFT = abs( [ Xk_DFT Xk_DFT] );                           % DFT magnitude

    angXk_DFT = angle( [Xk_DFT Xk_DFT] )/pi;                       % DFT angle

   realXk_DFT = real([Xk_DFT Xk_DFT]);

   imagXk_DFT = imag([Xk_DFT Xk_DFT]);

figure('NumberTitle', 'off', 'Name', 'P5.38.1 DFT(k) of x(n) N=40')

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

subplot(2,1,1); stem(k2, magXk_DFT); hold on; plot(N2*w/(2*pi), magXw_DTFT,'r--'); hold off;

%axis([-N/2, N/2, -0.5, 50.5]);

xlabel('k'); ylabel('magnitude(k)');

title('DFT magnitude of x(n)');  grid on;

subplot(2,1,2); stem(k2, angXk_DFT); hold on; plot(N2*w/(2*pi), angXw_DTFT,'r--'); hold off;

%axis([-N/2, N/2, -0.5, 50.5]);

xlabel('k'); ylabel('angle(k)');

title('DFT angle of x(n)');  grid on;

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

%%                 DTFT X(w) of xn sequence, w=[-pi:pi],

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

MM = 500;

[Xw_DTFT, w] = dtft3(x, n, MM);

 magXw_DTFT = abs(Xw_DTFT);   angXw_DTFT = angle(Xw_DTFT)/pi;

realXw_DTFT = real(Xw_DTFT); imagXw_DTFT = imag(Xw_DTFT);

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

%%        START X_DTFT's  mag ang real imag

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

figure('NumberTitle', 'off', 'Name', 'P5.38.1 X(w) DTFT of x(n) N=40');

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

subplot(2,2,1); plot(w/pi,magXw_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, angXw_DTFT); grid on;  % axis([-2,2,-1,1]);

title('Angle Part');

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

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

title('Real Part');

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

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

title('Imaginary Part');

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

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

%%        END X_DTFT's  mag ang real imag

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

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

%%                 DFT(k) of xn sequence, k=[-N/2:N/2]

%%                    w=2pi*k/N       k=Nw/(2pi)

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

N3 = length(x);

%k1 = [0 : length(x)-1];

%k2 = [-N : N-1];

k3 = [-N3/2 : N3/2];

Xk_DFT = dft(x, N3);                                                   % DFT

    magXk_DFT = abs( [ Xk_DFT(N3/2+1 : N3) Xk_DFT(1 : N3/2+1)] );        % DFT magnitude

    angXk_DFT = angle( [Xk_DFT(N3/2+1 : N3) Xk_DFT(1 : N3/2+1)] )/pi;    % DFT angle

%   realXk_DFT = real([Xk_DFT(N/2+1 : N) Xk_DFT(1 : N/2+1)]);

%   imagXk_DFT = imag([Xk_DFT(N/2+1 : N) Xk_DFT(1 : N/2+1)]);

figure('NumberTitle', 'off', 'Name', 'P5.38.1 DFT(k) of x(n) N=40')

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

subplot(2,1,1); stem(k3, magXk_DFT); hold on; plot(N3*w/(2*pi), magXw_DTFT,'r--'); hold off;

%axis([-N/2, N/2, -0.5, 50.5]);

xlabel('k'); ylabel('magnitude(k)');

title('DFT magnitude of x(n)');  grid on;

subplot(2,1,2); stem(k3, angXk_DFT); hold on; plot(N3*w/(2*pi), angXw_DTFT,'r--'); hold off;

%axis([-N/2, N/2, -0.5, 50.5]);

xlabel('k'); ylabel('angle(k)');

title('DFT angle of x(n)');  grid on;

　　运行结果：

下面是N=50的情况，附图

《DSP using MATLAB》Problem 5.38的更多相关文章

  1. 《DSP using MATLAB》Problem 8.38

    代码: function [wpLP, wsLP, alpha] = bp2lpfre(wpbp, wsbp) % Band-edge frequency conversion from bandpa ...

  2. 《DSP using MATLAB》Problem 7.38

    代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...

  3. 《DSP using MATLAB》Problem 7.27

    代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...

  4. 《DSP using MATLAB》Problem 7.26

    注意:高通的线性相位FIR滤波器,不能是第2类,所以其长度必须为奇数.这里取M=31,过渡带里采样值抄书上的. 代码: %% +++++++++++++++++++++++++++++++++++++ ...

  5. 《DSP using MATLAB》Problem 7.25

    代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...

  6. 《DSP using MATLAB》Problem 7.24

    又到清明时节,…… 注意:带阻滤波器不能用第2类线性相位滤波器实现,我们采用第1类,长度为基数,选M=61 代码: %% +++++++++++++++++++++++++++++++++++++++ ...

  7. 《DSP using MATLAB》Problem 7.23

    %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output Info a ...

  8. 《DSP using MATLAB》Problem 7.16

    使用一种固定窗函数法设计带通滤波器. 代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ...

  9. 《DSP using MATLAB》Problem 7.15

    用Kaiser窗方法设计一个台阶状滤波器. 代码: %% +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ...

随机推荐

  1. Win10系列：VC++媒体播放控制3

    (5)添加视频进度条 视频进度条可以用来显示当前视频的播放进度,并可以通过拖动视频进度条来改变视频的播放进度.接下来介绍如何实现视频进度条,首先打开MainPage.xaml文件,并在Grid元素中添 ...

  2. POJ 3436 ACM Computer Factory 最大流,拆点 难度:1

    题目 http://poj.org/problem?id=3436 题意 有一条生产线,生产的产品共有p个(p<=10)零件,生产线上共有n台(n<=50)机器,每台机器可以每小时加工Qi ...

  3. day15-python常用内置模块的使用

    在日常的开发工作中,我们要写很多的python代码,如果都写在一个文件中,会导致代码特别难维护,为了拓展代码的可维护性,我们把函写在不同的文件里,这样每个文件包含的文件就比较少,逻辑更加清楚.在pyt ...

  4. tomcat 服务器线程问题

    http://blog.csdn.net/wtopps/article/details/71339295 http://blog.csdn.net/wtopps/article/details/713 ...

  5. vue-router-5-命名路由

    创建 Router 实例的时候,在 routes 配置中给某个路由设置名称 const router = new VueRouter({ routes: [ { path: '/user/:userI ...

  6. Linux系统命令行中vim编辑器取消高亮显示

    由于在使用vim编辑代码的时候不小心忘记首先输入i(insert)模式,导致写的代码出现了棕黄色的阴影显示 摸索了很久终于找到了解决方法: 1.退出vim编译器 2.在在命令行下输入:nohl,回车 ...

  7. Linux学习 ：多线程编程

    1.Linux进程与线程() 进程:通过fork创建子进程与创建线程之间是有区别的:fork创建出该进程的一份拷贝,创建时额外申请了新的内存空间以及存储代码段.数据段.BSS段.堆.栈空间,     ...

  8. ssh框架中，工具类调用service层方法（参考https://www.cnblogs.com/l412382979/p/8526945.html）

    代码如下: package common.dataService; import javax.annotation.PostConstruct; import org.springframework. ...

  9. DevExpress WPF v18.2新版亮点（六）

    买 DevExpress Universal Subscription  免费赠 万元汉化资源包1套! 限量15套!先到先得,送完即止!立即抢购>> 行业领先的.NET界面控件2018年第 ...

  10. 解决NPM无法安装任何包的解决方案（npm ERR! code MODULE_NOT_FOUND）

    前言 今天突然发现npm无法使用了,执行任何命令都报如下错误: npm ERR! code MODULE_NOT_FOUND npm ERR! Cannot find module 'internal ...