找回密码
 注册
关于网站域名变更的通知
查看: 786|回复: 3
打印 上一主题 下一主题

模拟信号采样及离散时间傅里叶变换(DTFT)案例分析

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-6-16 13:48 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
. _3 G  ~# U/ [8 R
这篇中使用的模拟信号为上篇:使用 MATLAB 实现模拟信号的近似及其连续傅里叶变换' ]# b; q. q7 t/ ?
0 E) L- W0 ?# _* ?+ X/ q0 s3 d* O
中使用的模拟信号:0 A% }* g+ F3 `0 }. ^( u/ F

0 A3 l' g2 p  L) @  p ! p! F! e2 B* K5 Y
7 y1 N8 M8 R0 u$ T) y+ g
为了研究在频域数量上的采样效果,对该信号使用两种不同的采样频率采样。
& @7 a- q4 s% }1 o' }* |* p
! F$ B! e* J# Qa. 在 fs = 5000 对信号进行采样,求出并画出其离散时间傅里叶变换;
) C* F5 z" ^0 R4 s
3 Y2 ?5 h( P, Fb. 在 fs = 1000 对信号采样,求出并画出其离散时间傅里叶变换。0 g$ D5 u6 e3 w2 w
$ D) J; [, Y& }* R/ }  m0 E
题解:
; j, Q, I8 ]( R0 n! Z
: _4 B3 t. v) z* Y( i  A1 [) f( z上篇也分析了,信号的带宽为2kHz,奈奎斯特频率就为 4000 样本/s,它小于第一问给出的采样频率,所以频谱混叠几乎不存在。3 `2 z% d( @  M1 h
; A( G& T, J( Y% X+ A
我们通过MATLAB验证下:- n! {) W8 M: M) C! f, T1 z

9 e- |: e. t8 p2 S. ^
  • clc
  • clear
  • close all
  • % Analog signal
  • Dt = 0.00005;
  • t = - 0.005:Dt:0.005;
  • xa = exp(-1000 * abs(t));
  • % Discrete-time signal
  • Ts = 0.0002;
  • n = -25:25;
  • x = exp(-1000*abs(n*Ts));
  • %Discrete-time Fourier Transform
  • K = 500;
  • k = 0:1:K;
  • w = pi*k/K;
  • X = x*exp(-j*n'*w);
  • X = real(X);
  • w = [-fliplr(w),w(2:K+1)];
  • X = [fliplr(X),X(2:K+1)];
  • subplot(2,1,1);
  • plot(t*1000,xa);
  • xlabel('t in msec');ylabel('x(n)');
  • title('Discrete signal');
  • hold on
  • stem(n*Ts*1000,x);
  • gtext('Ts = 0.2msec');
  • hold off
  • subplot(2,1,2);
  • plot(w/pi,X);
  • xlabel('Frequency in pi units');ylabel('X(w)');
  • title('Discrete - time Fourier Transform');) I2 ^' m1 i3 l  s  {; S* s2 W6 `
   
( k2 R/ Y, ^6 U6 A9 H9 Q
- X3 _+ f+ T' Q9 W9 ~; p) Q0 n' O
* k7 D2 i! I- ~* o继续作答第二问,将采样间隔改为0.001s,10ms时间区间内能采样11个点。! W( Y/ @0 G* s
$ ]3 m7 r2 v7 I3 w, E5 h
  • clc
  • clear
  • close all
  • % Analog signal
  • Dt = 0.00005;
  • t = - 0.005:Dt:0.005;
  • xa = exp(-1000 * abs(t));
  • % Discrete-time signal
  • Ts = 0.001;
  • n = -5:5;
  • x = exp(-1000*abs(n*Ts));
  • %Discrete-time Fourier Transform
  • K = 500;
  • k = 0:1:K;
  • w = pi*k/K;
  • X = x*exp(-j*n'*w);
  • X = real(X);
  • w = [-fliplr(w),w(2:K+1)];
  • X = [fliplr(X),X(2:K+1)];
  • subplot(2,1,1);
  • plot(t*1000,xa);
  • xlabel('t in msec');ylabel('x(n)');
  • title('Discrete signal');
  • hold on
  • stem(n*Ts*1000,x);
  • gtext('Ts = 1msec');
  • hold off
  • subplot(2,1,2);
  • plot(w/pi,X);
  • xlabel('Frequency in pi units');ylabel('X(w)');
  • title('Discrete - time Fourier Transform');
    2 H, }# h* u3 a4 n  `5 {* \, m
   4 m/ K' o  B6 M! B1 x

% |, H: y) B; ^- E  G0 X很明显,这个DTFT7已经变了,也就是发生了频谱混叠。7 _3 {" H5 a( p/ n
0 E! L! J4 ~& I/ z
* j& q+ o% g3 K9 [8 R
  • TA的每日心情

    2019-11-29 15:37
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2021-6-16 14:54 | 只看该作者
    模拟信号采样及离散时间傅里叶变换(DTFT)案例分析

    该用户从未签到

    3#
    发表于 2021-6-16 17:10 | 只看该作者
    模拟信号采样及离散时间傅里叶变换(DTFT)案例分析
  • TA的每日心情

    2019-11-19 15:34
  • 签到天数: 1 天

    [LV.1]初来乍到

    4#
    发表于 2021-6-16 17:11 | 只看该作者
    模拟信号采样及离散时间傅里叶变换(DTFT)案例分析
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-30 16:05 , Processed in 0.156250 second(s), 27 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表