上篇博文通过在理论上讨论了DFT的三个性质:【 MATLAB 】DFT性质讨论(一)线性、循环反转与共轭

分别讨论:

一、线性

给出一个例子,给出x1和x2,x3 = 0.3*x1+0.8*x2;

之后我们求x3的DFT,和x1和x2的DFT的线性组合是否一致,即可验证线性性质。

clc,clear,close all;%signal 1
n1 = 0:3;
x1 = [1,1,1,1];%signal 2
n2 = 0:5;
x2 = ones(1,6);% signal 3 = signal 1 +signal 2
n3 = 0:max(length(n1),length(n2))-1;
% modify x1 and x2
x1 = [ones(1,4),zeros(1,length(n3)-4)];
x2 = [ones(1,6),zeros(1,length(n3)-6)];x3 = 0.3*x1+0.8*x2;N = length(n3);% DFT of x1
X1 = dft(x1,N);% DFT of x2
X2 = dft(x2,N);% DFT of x3
X3 = dft(x3,N);% Linear property
X3_l = 0.3*X1 + 0.8*X2;subplot(2,1,1);
stem(n3,X3);
title('DFT of x3');subplot(2,1,2);
stem(n3,X3_l);
title('DFT of  0.3*x1+0.8*x2');

二、循环反转

给出例子:

求:

a. 求出并画出 

b. 验证循环反转性质。

题解:


a

clc,clear,close all;n = 0:10;
x = 10 * (0.8).^n;
N = 11;y = x( mod(-n,N)+1 );subplot(2,1,1)
stem(n,x);
title('Original Sequence');
xlabel('n');ylabel('x(n)');subplot(2,1,2)
stem(n,y);
title('Circularly folded Sequence');
xlabel('n');ylabel('x(-n mod 10)');


b

clc,clear,close all;n = 0:10;
x = 10 * (0.8).^n;
N = 11;y = x( mod(-n,N)+1 );% subplot(2,1,1)
% stem(n,x);
% title('Original Sequence');
% xlabel('n');ylabel('x(n)');
%
% subplot(2,1,2)
% stem(n,y);
% title('Circularly folded Sequence');
% xlabel('n');ylabel('x(-n mod 10)');
k = n;
X = dft(x,N);
Y = dft(y,N);
RealX = real(X);
ImagX = imag(X);RealY = real(Y);
ImagY = imag(Y);subplot(2,2,1);
stem(k,RealX);
title('DFT Real Part of x(n)');
xlabel('k');ylabel('Real{ DFT(x(n)) }');subplot(2,2,2);
stem(k,ImagX);
title('DFT Imag Part of x(n)');
xlabel('k');ylabel('Imag{ DFT(x(n)) }');subplot(2,2,3);
stem(k,RealY);
title('DFT Real Part of x(-n mod N)');
xlabel('k');ylabel('Real{DFT(x(-n mod N))}');subplot(2,2,4);
stem(k,ImagY);
title('DFT Imag Part of x(-n mod N)');
xlabel('k');ylabel('Imag{DFT(x(-n mod N))}');


三、共轭与时序列的对称性

我们讨论时序列的对称性,采用的例子是这篇博文中的:【 MATLAB 】模拟信号采样及离散时间傅里叶变换(DTFT)案例分析

使用的模拟信号:

在 fs = 5000 对信号进行采样,可以得到离散时间序列。

之后讨论它的循环共轭对称性。

clc
clear
close all% Discrete-time signal
Ts = 0.0002;
n = 0:25;
x = exp(-1000*abs(n*Ts));
N = length(n);x1 = x( mod(-n,N) + 1 );% DFT of x
X = dft(x,N);X1 = dft(x1,N);X1_c = conj(X1);subplot(2,2,1);
stem(n,x);
title('Original sequence');
xlabel('n');ylabel('x(n)');subplot(2,2,2);
stem(n,X);
title('DFT of x(n)');
xlabel('k');ylabel('DFT{x(n)}');subplot(2,2,3);
stem(n,x1);
title('Circularly folded sequence');
xlabel('n');ylabel('x(-n mod N)');subplot(2,2,4)
stem(n,X1);
title('DFT of x((-n mod N))');
xlabel('k');ylabel('DFT{x(-n mod N)}');figure
stem(n,X);
title('DFT of x(n)');
xlabel('k');ylabel('DFT{x(n)}');figurestem(n,X1_c);
title('Conjugation{DFT of x((-n mod N))}');
xlabel('k');ylabel('Conj{DFT{x(-n mod N)}}');

【 MATLAB 】DFT性质讨论(一)线性、循环反转、共轭与实序列的对称性的MATLAB实现相关推荐

  1. 【 MATLAB 】DFT性质讨论(一)线性、循环反转、共轭与实序列的对称性

    本篇博文通过手稿版的形式对DFT的线性,循环反转以及共轭性质给予了总结. 下篇博文设想通过MATLAB来实现这几个性质.

  2. 【数字信号处理】傅里叶变换性质 ( 序列傅里叶变换共轭对称性质 | 实序列的幅频特性偶对称 | 实序列相频特性奇对称 | 示例说明 )

    文章目录 一.实序列的 幅频特性 和 相频特性 对称性质 二.性质由来 三.示例说明 一.实序列的 幅频特性 和 相频特性 对称性质 如果 x(n)x(n)x(n) 序列是 " 实序列 &q ...

  3. 【 MATLAB 】【 MATLAB 】DFT的性质讨论(三)序列的循环卷积及其 MATLAB 实现

    这篇博文要将的是循环卷积,循环卷积和线性卷积还是有很大区别的,我们都知道,两个N点序列之间的线性卷积会得出一个更长的序列,不得不再一次要将区间限制在 0 <= n <= N -1.因此代替 ...

  4. 【 MATLAB 】DFT的性质讨论(二)序列的循环移位及其 MATLAB 实现(频域方法)

    上篇博文:[ MATLAB ]DFT的性质讨论(二)序列的循环移位及其 MATLAB 实现(时域方法) 提到了对序列x(n)做循环移位后的DFT形式为: 上篇博文已经讨论过了第一种实现循环移位的方法, ...

  5. 循环卷积的时域DFT性质推导

    循环卷积的时域DFT性质推导 最近研究了OFDM系统,借着学习循环前缀的机会,重新理解了一下循环卷积的思想,循环前缀主要利用了循环卷积的时域DFT性质,可以大幅的降低接收端的检测复杂度,对于这个性质的 ...

  6. matlab中线性变压器,用MATLAB优化计算差动变压器的线性范围

    基金项目:河南省自然科学基金资助项目(994043500) 作者简介:胡明江(1974 - ) ,男 ,河南洛阳人 ,硕士生. 收稿日期:2003 - 04 - 23 文章编号 :1672 - 687 ...

  7. Matlab变量、分支语句和循环语句

    Matlab变量.分支语句和循环语句 一.Matlab的变量 1.变量的命名 Matlab的变量名由数字.字母和下划线组成,可是仅仅能由字母开头.大写和小写敏感,最大长度为63个字符.不能使用Matl ...

  8. dft计算傅里叶级数系数_一道国外的DFT性质的题目

    由于上半年实在太忙太忙,所以导致很久没更新公众号了,特意向各位长期关注的小伙伴表示歉意.今天分享的是DFT性质的应用. 背景:DFT的对称性在解题中是非常常见的,很多同学,一看到"实序列&q ...

  9. 回归预测 | MATLAB实现CNN-GRU(卷积门控循环单元)多输入单输出

    回归预测 | MATLAB实现CNN-GRU(卷积门控循环单元)多输入单输出 目录 回归预测 | MATLAB实现CNN-GRU(卷积门控循环单元)多输入单输出 基本介绍 模型结构 CNN神经网络 G ...

最新文章

  1. 应届前端面试——看这篇就够了(一)
  2. 压缩aspx页面,移除aspx多余的空格 供学习参考
  3. Java是否支持默认参数值?
  4. 真AI用钱表达:这家少年班毕业生创办的AI公司三年就盈利,增速300%
  5. python – IOError:[Errno 28] pip install 设备上没有空间
  6. java设计模式2--抽象工厂模式(Abstract Factory)
  7. javascript中原型模式创建对象特点分析
  8. 1.Python介绍
  9. UCOS操作系统——时间片轮转调度(五)
  10. 推荐5个免费的项目管理工具
  11. 职业规划-服务器运维工程师(待完善)
  12. OKR 如何转变你的绩效管理策略
  13. IDEA进行对json数据格式化
  14. python考试分几级_全国计算机等级考试2级分几类?到底怎么考啊?
  15. 汇编语言学习篇2——MASM的环境搭建(win10与Ubuntu1604下的配置)【有问题,待更正】
  16. STM8学习笔记1:CPU简介
  17. GameFrameWork框架(Unity3D)使用笔记(八) 实现场景加载进度条
  18. GitLab针对关键账户接管漏洞发布安全补丁
  19. Redis核心技术笔记——Redis数据结构
  20. 思科下一代模拟器EVE-NG安装

热门文章

  1. python3爬取数据存入mysql_Python如何爬取51cto数据并存入MySQL
  2. 如何删除oracle用户数据库用户,oracle删除指定用户的原数据库,建立该用户的新数据库...
  3. java与与短路与_Java中短路运算符与逻辑运算符示例详解
  4. ctrl选中多个_PS快捷键/Ctrl/Alt/Shift
  5. 用计算机画函数图像教案,信息技术应用 用计算机画函数图象教案1
  6. Getting Started with OpenCV
  7. 粘贴铜箔简易实验电路制作
  8. 图像修复中的方法--AI智能.
  9. ST-3806系列单圈编码器 说明书
  10. 实验电路的放大器伴侣:LMV358