《DSP using MATLAB》示例Example5.10
代码:
n = 0:10; x = 10*(0.8) .^ n;
[xec, xoc] = circevod(x);
%% -----------------------------------------------------------------
%% START a
%% -----------------------------------------------------------------
figure('NumberTitle', 'off', 'Name', 'Exameple5.10a x sequence')
set(gcf,'Color','white');
subplot(3,1,1); stem(n,x); title('Original sequence x(n)'); %axis([0,10,-2.5,2.5]);
xlabel('n'); ylabel('x(n)'); grid on;
subplot(3,1,2); stem(n,xec); title('Circular-even component'); axis([-0.5,10.5,-1,11]);
xlabel('n'); ylabel('xec(n)'); grid on;
subplot(3,1,3); stem(n,xoc); title('Circular-odd component'); axis([-0.5,10.5,-4,4]);
xlabel('n'); ylabel('xoc(n)'); grid on; %% -----------------------------------------------------------------
%% END a
%% ----------------------------------------------------------------- %% -----------------------------------------------------------------
%% START b
%% -----------------------------------------------------------------
X_DFT = dft(x,11); % DFT of x
Xec_DFT = dft(xec,11); Xoc_DFT = dft(xoc,11); % DFT of xec and xoc realX_DFT = real(X_DFT); imagX_DFT = imag(X_DFT); magX_DFT = abs(X_DFT(1:1:6)); phaX_DFT = angle(X_DFT(1:1:6))*180/pi;
realXec_DFT = real(Xec_DFT); imagXec_DFT = imag(Xec_DFT); magXec_DFT = abs(Xec_DFT(1:1:6)); phaXec_DFT = angle(Xec_DFT)*180/pi;
realXoc_DFT = real(Xoc_DFT); imagXoc_DFT = imag(Xoc_DFT); magXoc_DFT = abs(Xoc_DFT(1:1:6)); phaXoc_DFT = angle(Xoc_DFT)*180/pi; figure('NumberTitle', 'off', 'Name', 'Exameple5.10b ')
set(gcf,'Color','white');
subplot(2,2,1); stem(n,realX_DFT); title('Real {DFT[x(n)]}'); axis([-0.5,10.5,-5,50]);
xlabel('k'); grid on;
subplot(2,2,2); stem(n,imagX_DFT); title('Imag {DFT[x(n)]}'); axis([-0.5,10.5,-20,20]);
xlabel('k'); grid on;
subplot(2,2,3); stem(n,realXec_DFT); title('Real {DFT[xec(n)]}'); axis([-0.5,10.5,-5,50]);
xlabel('k'); grid on;
subplot(2,2,4); stem(n,imagXoc_DFT); title('Imag {DFT[xoc(n)]}'); axis([-0.5,10.5,-20,20]);
xlabel('k'); grid on; figure('NumberTitle', 'off', 'Name', 'Exameple5.10b X DFT info')
set(gcf,'Color','white');
subplot(2,2,1); stem(n,realX_DFT); title('Real {DFT[x(n)]}'); axis([-0.5,10.5,-5,50]);
xlabel('k'); grid on;
subplot(2,2,2); stem(n,imagX_DFT); title('Imag {DFT[x(n)]}'); axis([-0.5,10.5,-20,20]);
xlabel('k'); grid on; k = 0:1:5; w = 2*pi/10*k;
subplot(2,2,3); stem(w/pi,magX_DFT); title('Magnitude DTFT'); %axis([-0.5,10.5,-5,50]);
xlabel('frequency in \pi units'); grid on;
subplot(2,2,4); stem(w/pi,phaX_DFT); title('Angle {DFT[x(n)]}'); %axis([-0.5,10.5,-20,20]);
xlabel('frequency in \pi units'); grid on; %% -----------------------------------------------------------------
%% END b
%% -----------------------------------------------------------------
运行结果:
由上图看出,性质(5.34)成立:
《DSP using MATLAB》示例Example5.10的更多相关文章
- DSP using MATLAB 示例 Example3.10
用到的性质 上代码: n = -5:10; x = rand(1,length(n)) + j * rand(1,length(n)); k = -100:100; w = (pi/100)*k; % ...
- DSP using MATlAB 示例Example2.10
上代码 % noise sequence 1 x = [3, 11, 7, 0, -1, 4, 2]; nx = [-3:3]; % given signal x(n) [y,ny] = sigshi ...
- 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.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 , ...
- DSP using MATLAB 示例 Example3.11
用到的性质 上代码: n = -5:10; x = rand(1,length(n)); k = -100:100; w = (pi/100)*k; % freqency between -pi an ...
- DSP using MATLAB 示例Example3.8
代码: x = rand(1,11); n = 0:10; k = 0:500; w = (pi/500)*k; % [0,pi] axis divided into 501 points. X = ...
随机推荐
- Base64原理
一.Base64编码由来 为什么会有Base64编码呢?因为有些网络传送渠道并不支持所有的字节,例如传统的邮件只支持可见字符的传送,像ASCII码的控制字符就 不能通过邮件传送.这样用途就受到了很大的 ...
- [Java] JSP笔记 - Java Bean
一. Java Bean 是什么? Java Bean 其实就是一个符合特定规则的 Java Class.这些规则包括: 使用 public 声明的公用类 属性使用public 的get,set方法访 ...
- APP的消息推送(极光推送)
APP的消息推送,使用的第三方平台是极光推送 简单案例(以Thinkphp为例): 1.下载下载PHPSDK 2.把PHPSDK目录下的jpush-api-php-client-3.5.1\src\J ...
- Windows平台下Qt中glut库的使用
用Qt中的QGLWidget窗体类中是不包括glut工具库的,难怪在myGLWidget(在我的程序中是QGLWidget的派生类)中绘制实心球体是说“glutSolidSphere”: 找不到标识符 ...
- C# 发送电子邮件
网上找到的发送邮件的类,改了一点点,在此谢谢原作者的奉献. 1.源码: public class CSendMail { private MailMessage mailMessage; privat ...
- jq的$()里面 一定要是字符串类型!!!!!!!!!!!!!!!!!!!!!!!!
var s = "[value="+uid+"]"; $(s).attr("checked",'true');
- PHP数组函数: array_map()
定义和用法 array_map() 函数返回用户自定义函数作用后的数组.回调函数接受的参数数目应该和传递给 array_map() 函数的数组数目一致. 语法 array_map(function,a ...
- multipart/form-data和application/x-www-form-urlencoded的区别
在Form元素的语法中,EncType表明提交数据的格式 用 Enctype 属性指定将数据回发到服务器时浏览器使用的编码类型. 下边是说明: application/x-www-form-urle ...
- CMake命令/函数汇总(翻译自官方手册)
查看官方文档 cmake命令 选项 CMake变量 CMake命令汇总 / add_custom_command add_custom_target/add_definitions/add_depen ...
- mongDB-- 3. 查询操作
1. 准备工作 (1)启动mongo 进入mongo安装目录的bin/ 目录 , ./mongod (2)启动mongo客户端 ./mongo (3) 查看所有库 show dbs; (4) 切换到l ...