|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
" V! D6 g' O& [% a F; _ g
这篇中使用的模拟信号为上篇:使用 MATLAB 实现模拟信号的近似及其连续傅里叶变换
3 r5 o' |" F: z1 w! Z$ N
- N7 W8 v+ }8 g+ e: ^+ i中使用的模拟信号:. Q ^6 e6 J% T% r) O
/ H- @# j# C, v2 s5 ]% T
. K& y9 N- [1 I. g' }
' E% `5 f( E1 j2 B$ ^
为了研究在频域数量上的采样效果,对该信号使用两种不同的采样频率采样。: S; o3 `& ~: V: ]& h
' Y3 V7 k+ ?) d% ?7 s/ A K! Aa. 在 fs = 5000 对信号进行采样,求出并画出其离散时间傅里叶变换;
+ f, W1 I' B& T1 o; B: b
% I" p- ?8 f2 H* _0 q1 bb. 在 fs = 1000 对信号采样,求出并画出其离散时间傅里叶变换。
; k6 J3 |1 ^/ i" g |: x. `# R% ]" ~: T/ Y5 A: S
题解:3 V8 Q* R1 \- o5 l; H
8 m9 h! C5 X6 { v5 C上篇也分析了,信号的带宽为2kHz,奈奎斯特频率就为 4000 样本/s,它小于第一问给出的采样频率,所以频谱混叠几乎不存在。
. [: w5 W, A2 x: [; P; \$ d4 O6 S* l6 [4 R& I7 T' [
我们通过MATLAB验证下:! Y6 s* W3 L7 H& A: X* K
, k& r; h( E1 c9 C3 p. F7 t8 a- 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');
$ c9 z! j: |1 G E- Q% G; X' R$ } 8 P$ Z" }' t; _2 G& ~" C
# `& x& H. y( K( |
7 O3 t, N$ G4 g- s2 ]( q+ ~3 o
继续作答第二问,将采样间隔改为0.001s,10ms时间区间内能采样11个点。
; Z5 O4 m( P! G$ y1 u& n" B$ A# o" I' F3 \% ]
- 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');) ?' h) {# R, f& T
9 c6 r: m% N9 ^, V
- z* X1 @" o' w, _$ {- F2 Q很明显,这个DTFT7已经变了,也就是发生了频谱混叠。' f$ k2 S3 A. S/ I% N6 P' c* o; U
1 G& j/ y6 p! n: w$ q+ O7 O" ?3 q
1 }) a6 A5 t+ E" {3 ~ |
|