《DSP using MATLAB》Problem 5.7
代码:
%% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 5.7 \n\n');banner(); %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++% -------------------------------------------------- % 1 x(n) = [ 2,4,6,1,3,5 ] N--period % | % -------------------------------------------------- x = [2, 4, 6, 1, 3, 5]; N = 6; n = [0 : N-1]; % periodic sequence with period N y = x( mod(-n, N)+1 ); % circular foldfigure('NumberTitle', 'off', 'Name', 'P5.7 x(n) and x(-n)') set(gcf,'Color','white'); subplot(2,1,1); stem(n, x); xlabel('n'); ylabel('x(n)'); title('$\tilde x(n) sequence, period N=6$', 'interpreter', 'latex', 'fontsize', 14); grid on; subplot(2,1,2); stem(n, y); xlabel('n'); ylabel('y(n)'); title('x(-n) sequence'); grid on;%% ------------------------------------------------------------------ %% DFS(k) of xn sequence, k=[0:N-1], N=6 %% ------------------------------------------------------------------ k1 = [0 : length(x)-1]; %k2 = [-N : N-1]; %k3 = [-N/2 : N/2]; Xk_1 = dfs(x, N); % DFS magXk = abs( [ Xk_1 ] ); % DFS magnitude angXk = angle( [Xk_1] )/pi; % DFS anglefigure('NumberTitle', 'off', 'Name', 'P5.7 DFS(k) of x(n)') set(gcf,'Color','white'); subplot(2,1,1); stem(k1, magXk); %axis([-N/2, N/2, -0.5, 50.5]); xlabel('k'); ylabel('magnitude(k)'); title('DFS magnitude of x(n), N=6'); grid on; subplot(2,1,2); stem(k1, angXk); %axis([-N/2, N/2, -0.5, 50.5]); xlabel('k'); ylabel('angle(k)'); title('DFS angle of x(n), N=6'); grid on;%% ------------------------------------------------------------------ %% DFS(k) of yn sequence, k=[0:N-1], N=6 %% ------------------------------------------------------------------ k1 = [0 : length(y)-1]; %k2 = [-N : N-1]; %k3 = [-N/2 : N/2]; Yk_1 = dfs(y, N); % DFS magYk = abs( [ Yk_1 ] ); % DFS magnitude angYk = angle( [Yk_1] )/pi; % DFS anglefigure('NumberTitle', 'off', 'Name', 'P5.7 DFS(k) of y(n)') set(gcf,'Color','white'); subplot(2,1,1); stem(k1, magYk); %axis([-N/2, N/2, -0.5, 50.5]); xlabel('k'); ylabel('magnitude(k)'); title('DFS magnitude of y(n), N=6'); grid on; subplot(2,1,2); stem(k1, angYk); %axis([-N/2, N/2, -0.5, 50.5]); xlabel('k'); ylabel('angle(k)'); title('DFS angle of y(n), N=6'); grid on;
运行结果:
周期序列及其圆周折叠序列:
原始序列的DFS系数
圆周折叠序列的DFS系数
转载于:https://www.cnblogs.com/ky027wh-sx/p/9314201.html
《DSP using MATLAB》Problem 5.7相关推荐
- 《DSP using MATLAB》Problem 6.24
代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...
- 《DSP using MATLAB》Problem 6.12
代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...
- 《DSP using MATLAB》Problem 6.20
先放子函数: function [C, B, A, rM] = dir2fs_r(h, r);% DIRECT-form to Frequency Sampling form conversion % ...
- 《DSP using MATLAB》Problem 7.26
注意:高通的线性相位FIR滤波器,不能是第2类,所以其长度必须为奇数.这里取M=31,过渡带里采样值抄书上的. 代码: %% +++++++++++++++++++++++++++++++++++++ ...
- 《DSP using MATLAB》Problem 6.6
代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...
- 《DSP using MATLAB》Problem 7.36
代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output In ...
- 《DSP using MATLAB》Problem 7.16
使用一种固定窗函数法设计带通滤波器. 代码: %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ...
- 《DSP using MATLAB》Problem 8.22
时光飞逝,亲朋会一个一个离我们远去,孤独漂泊一阵子后,我们自己也要离开, 代码: %% -------------------------------------------------------- ...
- 《DSP using MATLAB》Problem 8.1
代码: %% ------------------------------------------------------------------------ %% Output Info about ...
最新文章
- How to enable javascript in windows server 2008 R2 enterprise
- Android Service使用方法--简单音乐播放实例
- Linux Graphic DRI Wayland 显示子系统
- nginx php默认首页,iis、apache、nginx设置默认首页index.html方法图解 电脑维修技术网...
- 常用IOC注解按照作用分类
- Hadoop之资源调度器与任务推测执行
- 前端传递多个数组以及其他参数,springboot后台如何接收
- Visual Stdio 2017显示UTF-8字符串(中文的情况下),字符串中的字符无效
- java.util.regex.PatternSyntaxException: Dangling meta character '?' near index 0特殊字符表达
- 帆软报表属性设计基础
- Openssl s_server命令
- 射手网字幕打包下载(73.16G)
- Good, then we can start
- 数据结构之红黑树简介
- linux多线程编程 实验,linux操作系统-实验五-linux 多线程编程.docx
- Mac 远程 屏幕共享 screen sharing
- 计数器代码php,php 的计数器程序代码示例
- Sql Server Report Builder 计算标准偏差
- C++入门:让计算机“开口说话”
- python 高级部分