Data Converter里的DFT理解

具体描述:
1.使用DFT时候,要确保输入样本是N周期的,否则要用加窗调制的方法来缓解最后一个数据点和第一个点的不连续性。
下面进行Matlab代码验证:

clear 
clc
f=50;                 %feature frequency
Fs = 1000;            % Sampling frequency               
T = 1/Fs;             % Sampling period       
L = 20*Fs/f;           % Length of signal x(保证样本序列是N周期)
t = (0:L-1)*T;        % Time vector
X= 0.7*sin(2*pi*f*t);
Y = myDFT(X);
P1 = abs(Y/L);
f = Fs*(0:L-1)/L;
figure(1)
plot(f,P1) 
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')

在这里插入图片描述
只改动L = 20.1*Fs/f;
在这里插入图片描述
2.加入噪声时候,增加采样点数目可以增大处理增益,使得信号和噪声更容易分离。

clear 
clc
f=50;                 %feature frequency
Fs = 1000;            % Sampling frequency               
T = 1/Fs;             % Sampling period       
L = 20*Fs/f;           % Length of signal x(保证样本序列是N周期)
t = (0:L-1)*T;        % Time vector
X= 0.7*sin(2*pi*f*t)+1*randn(size(t));
Y = myDFT(X);
P1 = abs(Y/L);
f = Fs*(0:L-1)/L;
figure(1)
plot(f,P1) 
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')

在这里插入图片描述
只将20改成200;
在这里插入图片描述

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值