MATLAB求余弦函数的dtft,【 MATLAB 】用 MATLAB 实现离散时间傅里叶变换(DTFT)的两个案例分析...

先给出离散时间傅里叶变换的简单介绍:

如果 x(n) 是绝对可加的,即

c3b85c5bebf8c6156d7a638409cf4c78.gif

那么它的离散时间傅里叶变换给出为:

3ed5118970745ca929ed91d0e10d56c2.gif

w 称为数字频率,单位是每样本 rad(弧度)或 (弧度/样本)(rad/sample)

案例1:

9e1c2225f7ca274df0d87dbec59a9da6.gif的离散时间傅里叶变换

933497a5ec96aa0e448ac32aeb6c1aef.gif,并用MATLAB将

d687e3106796329510da66d314250ade.gif

edda1c04931777707f823363edb7a81c.gif之间的501个等分点上求值,并画出它的幅度、相位、实部和虚部。

题解:

由于x(n)是无限长的序列,所以不能直接用MATLAB直接从x(n)得到

66ac6aeaa0a802696c4f63c904cc2e01.gif。然而,我们可以用它对表达式

0bd50d93fa3c0d65df10c307a69ed19c.gif

5579a1f53a03d42f817a0b50e99be631.gif频率点上求值,然后画出它的幅度和相位(实部或虚部)。

e5b985be1349253c19a0d669a1876017.gif

0d1a8d5627f0764d02bb710858ef53cf.gif

脚本:

clc

clear

close all

w = [0:500]*pi/500; %[0,pi] axis divided into 501 points

X = exp(j*w)./( exp(j*w) - 0.5*ones(1,501) );

magX = abs(X);

angX = angle(X);

realX = real(X);

imagX = imag(X);

subplot(2,2,1)

plot(w/pi,magX);

grid;

title('Magnitude Part');

xlabel('frequency in pi units');ylabel('Magnitude');

subplot(2,2,2)

plot(w/pi,angX);

grid;

title('Angle Part')

xlabel('frequency in pi units');ylabel('Radians');

subplot(2,2,3)

plot(w/pi,realX);

grid

title('Real Part');

xlabel('frequency in pi units');ylabel('Real');

subplot(2,2,4);

plot(w/pi,imagX);

grid;

title('Imaginary Part');

xlabel('frequency in pi units');ylabel('Imaginary');

aa3d788870db9d866eb6c766f54a91f1.png

案例2:

求下面有限长序列的离散时间傅里叶变换:

e6d426083301599c63640d03ca94e2bb.gif

在[0,pi]之间的501个等分频率上进行数值求值。

题解:

cad3dbc91b2141937b951251ccc16d25.gif

我们可以直接对上式进行MATLAB编程,但是这种方法在有限长序列的DTFT中不是太方便,我们可以直接由 x(n) 来求它的DTFT。

我们使用向量化的编程方法,最后得到一个通用的公式。推导如下:

12ce5bbb00214b141604d57933a0ce76.png

用MATLAB实现如下:

k = [0:M];

n = [n1:n2];

X = x * (exp(-j * pi/M)).^(n'*k);

给出MATLAB脚本语言如下:

clc

clear

close all

n = -1:3;

x = 1:5;

k = 0:500;

w = (pi/500)*k;

X = x * (exp(-j * pi/500)).^(n' * k);

magX = abs(X);

angX = angle(X);

realX = real(X);

imagX = imag(X);

subplot(2,2,1);

plot(w/pi,magX);

title('Magnitude Part');

xlabel('w/pi');ylabel('Magnitude');

subplot(2,2,2);

plot(w/pi,angX);

title('Angle Part');

xlabel('w/pi');ylabel('Radians');

subplot(2,2,3);

plot(w/pi,realX);

title('Real part');

xlabel('w/pi');ylabel('Real');

subplot(2,2,4);

plot(w/pi,imagX);

title('Imaginary Part');

xlabel('w/pi');ylabel('Imaginary');

87da79826a838cd8f7d45bd49664c400.png

本文同步分享在 博客“李锐博恩”(CSDN)。

如有侵权,请联系 support@oschina.cn 删除。

本文参与“OSC源创计划”,欢迎正在阅读的你也加入,一起分享。

  • 0
    点赞
  • 28
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值