数字实验报告2离散傅立叶变换
时间:2020-11-11 16:47:46 来源:勤学考试网 本文已影响 人
附 页
第 页
陕西科技大学实验报告
第 页
班级 学号 姓名 实验组别
实验日期 室温 报告日期 成绩
实验名称:离散时间傅里叶变换
一、实验目的
1.复习离散时间傅立叶正反变换
2.复习DTFT的两个重要特性
3.复习DTFT的其它特性
4.离散LTI系统的频率响应
5.采样及重构信号
二、实验原理
1、信号的离散时间傅立叶变换(DTFT)
2、DTFT的两个重要特性
周期性:离散时间傅立叶变换是w的周期函数,其周期为2π。
对称性:对于实值的X(n),是共扼对称的。即实部为偶对称,虚部为奇对称。
3、DTFT的其他特性
线性;
时移:
共扼:
折叠:
卷积:
乘法:
能量:
4、LTI系统的频率响应
5、模拟信号的采样与重构
采样定理
重构: 步骤如下
(a)先把样本集转换成一个加权脉冲串列
(b)再将此脉冲串列通过一个带宽为F的低通滤波器进行滤波。
以上 两个步骤可用插值公式来描述:
实验内容
1.求信号的离散时间傅立叶变换并分析其周期性和对称性;
给定正弦信号x(t)=2*cos(2*pi*10*t),fs=100HZ,求其DTFT。
close all
clc
f=10;T=1/f;w=-10:0.2:10;
t1=0:0.0001:1;t2=0:0.01:1;
n1=-2;n2=8;n0=0;n=n1:0.01:n2;
x5=[n>=0.01];
x1=2*cos(2*f*pi*t1);
x2=2*cos(2*f*pi*t2);
x3=(exp(-j).^(t2'*w));
x4=x2*x3;
subplot(2,2,1);plot(t1,x1);
axis([0 1 1.1*min(x2) 1.1*max(x2)]);
xlabel('x(n)'); ylabel('x(n)');
title('原信号x1');
xlabel('t');ylabel('x1');
subplot(2,2,3);stem(t2,x2);
axis([0 1 1.*min(x2) 1.1*max(x2)]);
title('原信号采样结果x2');
xlabel('t');ylabel('x2');
subplot(2,2,2);stem(n,x5);
axis([0 1 1.*min(x5) 1.1*max(x5)]);
xlabel('n');ylabel('x2');
title('采样函数x2');
subplot(2,2,4);stem(t2,x4);
axis([0 1 -0.2+1.1*min(x4) 1.1*max(x4)]);
xlabel('t');ylabel('x4');
2.用以下两个有限长序列来验证DTFT的线性、卷积和共轭特性;
x1(n)=[1 2 3 4 5 6 7 8 9 10 11 12];x2(n)=R10(n)
1)线性:
w=linspace(-8,8,10000);
nx1=[0:11]; nx2=[0:9];
x1=[1 2 3 4 5 6 7 8 9 10 11 12];
x2=[1 1 1 1 1 1 1 1 1 1 ];
x3=[x2,zeros(1,(length(x1)-length(x2)))];
x4=2*x1+3*x3;
X1=x1*exp(-1i*nx1'*w);
X3=x3*exp(-1i*nx1'*w);
X4=x4*exp(-1i*nx1'*w);
subplot(5,3,4),stem(nx1,x1),axis([-1,13,0,20]);
subplot(5,3,2),stem(nx2,x2),axis([-1,13,0,5]);
subplot(5,3,3),stem(nx1,x4),axis([-1,13,0,26]);
subplot(5,3,1),plot(w,abs(X1));
subplot(5,3,7),plot(w,abs(X1));
subplot(5,3,10),plot(w,abs(X1));
subplot(5,3,13),plot(w,abs(X1));
subplot(5,3,5),plot(w,abs(X3));
subplot(5,3,8),plot(w,abs(X3));
subplot(5,3,11),plot(w,abs(X3));
subplot(5,3,14),plot(w,abs(X3));
subplot(5,3,6),plot(w,abs(X4));
subplot(5,3,9),plot(w,angle(X4));
subplot(5,3,12),plot(w,angle(X4));
subplot(5,3,15),plot(w,angle(X4));
2)卷积性:
close all
clc
nx1=0:11; nx2=0:9; nx3=0:20;
w=linspace(-8,8,40); %w=[-8,8]分10000份
x1=[1 2 3 4 5 6 7 8 9 10 11 12];
x2=[1 1 1 1 1 1 1 1 1 1];
x3=conv(x1,x2);% x1卷积x2
x4=x1*exp(-j*nx1'*w);% x1频率特性
x5=x2*exp(-j*nx2'*w);% x2频率特性
x6=x3*exp(-j*nx3'*w);% x1卷积x2频率特性
x7=x4.*x5;
subplot(2,2,1),stem(nx1,x1),axis([-1,15,0,15]),title('x1');
subplot(2,2,2),stem(nx2,x2),axis([-1,15,0,5]),title('x2');
subplot(2,1,2),stem(nx3,x3),axis([-1,25,0,80]);title('x1卷积x2结果x3');
figure,subplot(2,2,1),stem(x4,'filled'),title('x1的DTFT结果x4');
subplot(2,2,2),stem(x5,'filled'),title('x2的DTFT结果x5');
subplot(2,2,3),stem(x6,'filled'),title('x3的DTFT结果x6');
subplot(2,2,4),stem(x7,'filled'),title('x4的DTFT结果x7');
figure,subplot(3,2,1),stem(w,abs(x6)), ylabel('幅度'),title('x1卷积x2的DTFT');
subplot(4,2,3),stem(w,angle(x6)),ylabel('相位')
subplot(4,2,5),stem(w,real(x6)),ylabel('实部')
subplot(4,2,7),stem(w,imag(x6)),ylabel('虚部')
subplot(4,2,2),stem(w,abs(x7)), title('x1与x2的DTFT的乘积');
subplot(4,2,4),stem(w,angle(x7));
subplot(4,2,6),stem(w,real(x7));
subplot(4,2,8),stem(w,imag(x7))
3)共轭性
x1n=[1 2 3 4 5 6 7 8 9 10 11 12];
w=-10:10;
N1=length(x1n);
n1=0:N1-1;
x1=real(x1n);
x2=imag(x1n);
x2n=x1-j*x2;
X1=x2n*(exp(-j).^(n1'*w));
X2=x1n*(exp(j).^(n1'*w));
x3=real(X2); x4=imag(X2);
X2=x3-j*x4;
figure,subplot(211);stem(w,X1,'.');title('x1n共轭的DTFT');
subplot(212);stem(w,X2,'.');title('x1n的DTFT取共轭且反折')
3.求LTI系统的频率响应
给定系统H(Z)=B(Z)/A(Z),A=[0.98777 -0.31183 0.0256];B=[0.98997 0.989 0.98997],求系统的幅频响应和相频响应。(要求使用filter(B,A,δ(n))求解。
程序代码如下:
A=[0.98777 -0.31183 0.0256];
B=[0.98997 0.989 0.98997];
C=[1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
y=filter(B,A,C);
subplot(2,2,1);stem(y,'.');title('原始序列');
mag=abs(y);
ph=angle(y);
ph=ph*180/pi;
subplot(2,2,2);stem(mag,'.');title('幅频特性');
xlabel('时间信号n');ylabel('信号幅度');
subplot(2,2,3);stem(ph,'.');title('相频特性');
xlabel('时间信号n');ylabel('信号相位');
4. 采样和频谱混叠
给定信号x(t)=100*exp(-100*t)*cos(2*pi*500*t),求该信号的频谱;当采样频率分别为fs1=2000HZ,fs2=1000HZ;fs3=500HZ;fs4=200HZ,时输出序列的DTFT。
x=100*exp(-100*t).*cos(2*pi*500*t);
t=-2:0.1:-1;w=-5:0.1:5;
x=100*exp(-100*t).*cos(2*pi*500*t);
y=x*(exp(-j).^(t'*w));
subplot(2,1,1),plot(t,x);
subplot(2,1,2),plot(w,y);title('原始信号的频谱');
figure,fs1=2000;Ts1=1/fs1;n1=-2;Ts1:2;
fs2=1000;Ts2=1/fs2;n2=-2;Ts2:2;
fs3=500;Ts3=1/fs3;n3=-2;Ts3:2;
fs4=200;Ts4=1/fs4;n4=-2;Ts4:2;
x1=100.*exp(-100*n1).*cos(2*pi*500*n1);y1=x1*(exp(-j).^(n1'*w));
subplot(221);plot(w,y1);title('经2000Hz采样后信号的DTFT');
x2=100.*exp(-100*n2).*cos(2*pi*500*n2);y2=x2*(exp(-j).^(n2'*w));
subplot(222);plot(w,y2);title('经1000Hz采样后信号的DTFT');
x3=100.*exp(-100*n3).*cos(2*pi*500*n3);
y3=x3*(exp(-j).^(n3'*w));
subplot(223);plot(w,y3);title('经500Hz采样后信号的DTFT');
x4=100.*exp(-100*n4).*cos(2*pi*500*n4);y4=x4*(exp(-j).^(n4'*w));
subplot(224);plot(w,y4);title('经200Hz采样后信号的DTFT');
实验结论
通过这次实验,让我复习了离散时间傅立叶正反变换的各种基本性质及DTFT的两个重要特性和其它特性,离散LTI系统的频率响应和采样及重构信号,离散傅里叶变换实现了频域离散化,使数字信号处理可以在频域采用数值运算的方法进行,这样就增加了数字信号处理的灵活性。