《DSP using MATLAB》Problem 5.17

1、代码

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

%%            Output Info about this m-file

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

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

banner();

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

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

%          X(k) is 10-point DFTs of real-valued sequence x(n)

%          X(k) = [10, -2+3j, 3+4j, 2-3j, 4+5j, 12, 4-5j, 2+3j, 3-4j, -2-3j]

%               N = 10       k=[0:9]

%             1   y(n) = x((2-n))10

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

    k1 = [0:9];

Xk_DFT = [10, -2+3j, 3+4j, 2-3j, 4+5j, 12, 4-5j, 2+3j, 3-4j, -2-3j];

    N1 = length(Xk_DFT);                                              % length is 10

    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.17.1 X(k), DFT of x(n)')

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

subplot(2,2,1); stem(k1, magXk_DFT);

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

title('magnitude DFT of x(n), N=10');  grid on;

subplot(2,2,3); stem(k1, angXk_DFT);

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

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

title('angle DFT of x(n), N=10');  grid on;

subplot(2,2,2); stem(k1, realXk_DFT);

xlabel('k'); ylabel('real (k)');

title('real DFT of x(n), N=10');  grid on;

subplot(2,2,4); stem(k1, imagXk_DFT);

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

xlabel('k'); ylabel('imag (k)');

title('imag DFT of x(n), N=10');  grid on;

[xn] = real(idft(Xk_DFT, N1));                     % real-valued sequence

n = [0 : N1-1];

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

%         1st way to get y(n)-----circular shifft

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

xn_cirfold = xn(mod(-n,N1)+1);

m = 2;                                         % shift

yn1 = cirshftt(xn_cirfold, m, length(xn));

Xfk_DFT = dft(xn_cirfold, N1);

    magXfk_DFT = abs( [ Xfk_DFT ] );                          % DFT magnitude

    angXfk_DFT = angle( [Xfk_DFT] )/pi;                       % DFT angle

   realXfk_DFT = real(Xfk_DFT); imagXfk_DFT = imag(Xfk_DFT);

figure('NumberTitle', 'off', 'Name', 'P5.17.1 X((-k)), DFT of x((-n))')

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

subplot(2,2,1); stem(k1, magXfk_DFT);

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

title('magnitude DFT of x((-n)), N=10');  grid on;

subplot(2,2,3); stem(k1, angXfk_DFT);

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

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

title('angle DFT of x((-n)), N=10');  grid on;

subplot(2,2,2); stem(k1, realXfk_DFT);

xlabel('k'); ylabel('real (k)');

title('real DFT of x((-n)), N=10');  grid on;

subplot(2,2,4); stem(k1, imagXfk_DFT);

axis([0, 10, -5, 5.05]);

xlabel('k'); ylabel('imag (k)');

title('imag DFT of x((-n)), N=10');  grid on;

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

%       2ed way to get y(n)-----IDFT of Y(k)

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

k1 = [0:9];

Yk_DFT = exp(-j*2*pi*(2*k1)/10) .* Xk_DFT(mod(-k1, N1)+1);

    N1 = length(Yk_DFT);                                    % length is 10

    magYk_DFT = abs( [ Yk_DFT ] );                          % DFT magnitude

    angYk_DFT = angle( [Yk_DFT] )/pi;                       % DFT angle

   realYk_DFT = real(Yk_DFT); imagYk_DFT = imag(Yk_DFT);

figure('NumberTitle', 'off', 'Name', 'P5.17.1 DFT(k) of y(n)')

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

subplot(2,2,1); stem(k1, magYk_DFT);

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

title('magnitude DFT of y(n), N=10');  grid on;

subplot(2,2,3); stem(k1, angYk_DFT);

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

title('angle DFT of y(n), N=10');  grid on;

subplot(2,2,2); stem(k1, realYk_DFT);

xlabel('k'); ylabel('real (k)');

title('real DFT of y(n), N=10');  grid on;

subplot(2,2,4); stem(k1, imagYk_DFT);

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

xlabel('k'); ylabel('imag (k)');

title('imag DFT of y(n), N=10');  grid on;

[yn2] = real(idft(Yk_DFT, N1));

n = [0 : N1-1];

figure('NumberTitle', 'off', 'Name', 'P5.17.1 x(n) & y(n)')

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

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

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

title('x(n), IDFT of X(k)');  grid on;

subplot(2,2,2); stem(n, yn1);

xlabel('n'); ylabel('y(n)');

title('y(n) by circular shift x((2-n))_N N=10');  grid on;

subplot(2,2,3); stem(n, xn_cirfold);

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

title('x((-n)) N=10');  grid on;

subplot(2,2,4); stem(n, yn2);

xlabel('n'); ylabel('y(n)');

title('y(n) by IDFT of Y(k)');  grid on;

　　运行结果：

4、代码：

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

%%            Output Info about this m-file

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

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

banner();

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

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

%          X(k) is 10-point DFTs of real-valued sequence x(n)

%          X(k) = [10, -2+3j, 3+4j, 2-3j, 4+5j, 12, 4-5j, 2+3j, 3-4j, -2-3j]

%               N = 10       k=[0:9]

%             4   y(n) = x(n) Circular-Conv x((-n))10

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

    k1 = [0:9];

Xk_DFT = [10, -2+3j, 3+4j, 2-3j, 4+5j, 12, 4-5j, 2+3j, 3-4j, -2-3j];

    N1 = length(Xk_DFT);                                              % length is 10

    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.17.4 X(k), DFT of x(n)')

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

subplot(2,2,1); stem(k1, magXk_DFT);

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

title('magnitude DFT of x(n), N=10');  grid on;

subplot(2,2,3); stem(k1, angXk_DFT);

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

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

title('angle DFT of x(n), N=10');  grid on;

subplot(2,2,2); stem(k1, realXk_DFT);

xlabel('k'); ylabel('real (k)');

title('real DFT of x(n), N=10');  grid on;

subplot(2,2,4); stem(k1, imagXk_DFT);

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

xlabel('k'); ylabel('imag (k)');

title('imag DFT of x(n), N=10');  grid on;

[xn] = real(idft(Xk_DFT, N1));                     % real-valued sequence

n = [0 : N1-1];

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

%         1st way to get y(n)-----circular

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

xn_cirfold = xn(mod(-n,N1)+1);

       yn1 = circonvt(xn, xn_cirfold, N1);

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

%       2ed way to get y(n)-----IDFT of Y(k)

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

k1 = [0:9];

Xfk_DFT = dft(xn_cirfold, N1);

Yk_DFT = Xk_DFT .* Xfk_DFT;

    N1 = length(Yk_DFT);                                    % length is 10

    magYk_DFT = abs( [ Yk_DFT ] );                          % DFT magnitude

    angYk_DFT = angle( [Yk_DFT] )/pi;                       % DFT angle

   realYk_DFT = real(Yk_DFT); imagYk_DFT = imag(Yk_DFT);

figure('NumberTitle', 'off', 'Name', 'P5.17.4 DFT(k) of y(n)')

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

subplot(2,2,1); stem(k1, magYk_DFT);

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

title('magnitude DFT of y(n), N=10');  grid on;

subplot(2,2,3); stem(k1, angYk_DFT);

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

title('angle DFT of y(n), N=10');  grid on;

subplot(2,2,2); stem(k1, realYk_DFT);

xlabel('k'); ylabel('real (k)');

title('real DFT of y(n), N=10');  grid on;

subplot(2,2,4); stem(k1, imagYk_DFT);

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

xlabel('k'); ylabel('imag (k)');

title('imag DFT of y(n), N=10');  grid on;

[yn2] = real(idft(Yk_DFT, N1));

n = [0 : N1-1];

figure('NumberTitle', 'off', 'Name', 'P5.17.4 x(n) & y(n)')

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

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

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

title('x(n), IDFT of X(k)');  grid on;

subplot(2,2,2); stem(n, yn1);

xlabel('n'); ylabel('y(n)');

title('y(n) by x(n) Circular-Conv x((-n))_N N=10');  grid on;

%subplot(2,2,3); stem(n, xn_cirfold);

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

%title('x((-n)) N=10');  grid on;

subplot(2,2,4); stem(n, yn2);

xlabel('n'); ylabel('y(n)');

title('y(n) by IDFT of Y(k)');  grid on;

　　运行结果：

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

《DSP using MATLAB》Problem 6.17
代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...
《DSP using MATLAB》Problem 3.17
用差分方程两边进行z变换,再变量带换得到频率响应函数(或转移函数,即LTI系统脉冲响应的DTFT). 代码: %% ------------------------------------------ ...
《DSP using MATLAB》Problem 2.17
1.代码: %% ------------------------------------------------------------------------ %% Output Info abo ...
《DSP using MATLAB》Problem 8.17
代码: %% ------------------------------------------------------------------------ %% Output Info about ...
《DSP using MATLAB》Problem 4.17
《DSP using MATLAB》Problem 5.22
代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% O ...
《DSP using MATLAB》Problem 5.15
代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...
《DSP using MATLAB》Problem 2.18
1.代码: function [y, H] = conv_tp(h, x) % Linear Convolution using Toeplitz Matrix % ----------------- ...
《DSP using MATLAB》Problem 7.28
又是一年五一节,朋友圈都是晒名山大川的,晒脑袋的,我这没钱的待在家里上网转转吧频率采样法设计带通滤波器,过渡带中有一个样点代码: %% ++++++++++++++++++++++++++++++ ...

随机推荐

【转载】offer经验
http://www.cnblogs.com/figure9/archive/2013/01/09/2853649.html
linux内核initcall
include/linux/init.h #define pure_initcall(fn) __define_initcall(fn, 0) #define core_initcall(fn) __ ...
.NET接入UnionPay银联支付（一）手机wap支付
最近呢,比较忙,公司在接入银联全渠道支付,博主接手的wap支付,发表一下博主在接入的时候遇到的坑和注意事项,方便大家学习接入,爬坑的路上更顺利一点~ 开发步骤 1. 以表单的方式组装要发送给银联全渠道 ...
unity3D 中的C#脚本一个类调用另一类中简单方法
案例展示 SubMenuManage类中的实例化代码如下: static SubMenuManage sub_this; public static SubMenuManage Instance() ...
day12 生成器和各种推导式
今天主要学习了 1.生成器 2.生成器函数 3.各种推导式(比较诡异,理解了很简单,不理解很难) 4.生成器表达式(重点) 一.生成器 def func(): print'我叫周润发' return ...
解决WDCP3环境gbk网站编码程序乱码问题
因为默认WDCP V3版本环境编码格式是UTF-8版本,如果我们程序采用的是GBK编码肯定都会有乱码问题. 我们到WDCP后台,"网站管理"-"PHP设置",看 ...
控制台程序读取WIKI形式的TXT文件并一表格的形式显示在Word中
'Imports System.Collections.Generic 'Imports System.Text 'Imports System.IO 'Imports office = Micros ...
C++基础知识：类的静态成员
1.普通成员变量通过对象名能够访问public成员变量每个对象都可以有只属于自己的成员变量成员变量不能在对象之间共享 2.从命名空间的角度:类的静态成员只是类这个命名空间中的全局变量和全局函数不同之处 ...
Linux系统管理常用命令用法总结（2）
1.free指令会显示内存的使用情况,包括实体内存,虚拟的交换文件内存,共享内存区段,以及系统核心使用的缓冲区等. 语法:free [-bkmotV][-s <间隔秒数>] 参数说明: - ...
用swagger生成接口文档代码
1.Swagger2类: package com.example.demo; import com.google.common.base.Predicate; import io.swagger.an ...

《DSP using MATLAB》Problem 5.17

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

随机推荐

热门专题