数字信号处理实验指导

最新推荐文章于 2023-12-27 17:59:23 发布

lwz36

最新推荐文章于 2023-12-27 17:59:23 发布

阅读量676

点赞数

分类专栏：数字信号处理文章标签：信号处理

本文链接：https://blog.csdn.net/m0_55944825/article/details/130148456

版权

本文是一份关于数字信号处理的实验指导，涵盖抽样定理、FFT在频谱分析中的应用、C语言实现FFT及IFFT算法、IIR与FIR滤波器设计等内容。通过一系列实验，读者将深入理解信号处理的基本概念和方法，包括信号的重建、频谱分析、滤波器设计等。

摘要由CSDN通过智能技术生成

实验一、抽样定理

（一）实验目的：
1、验证乃奎斯特采样定理，加深对时域采样信号频谱变化的认识。
2、掌握序列傅里叶变换的计算机实现方法，利用序列傅里叶变换对连续信号、离散信号及系统响应进行频域分析。
（二）实验内容：
1、设x(n)=(0.8.) nu(n)。
（1）利用conv函数求出x(n)*x(n)的前50个样本，并作图。
相关m文件：
nb=0:50;
n=0:100;
x=(0.8).^nb;
yb=conv(x,x);
stem(n,yb);
xlabel(‘n’);
ylabel(‘yb(n)’);
title(‘convlution output’);
在这里插入图片描述

(2) 利用filter函数求出x(n)*x(n)的前50个样本，并作图。
相关m文件：
nb=0:50;
x=(0.8).^nb;
yb=filter(1,[1,-0.8],x);
stem(nb,yb);
axis([-1,51,-1,3])
xlabel(‘n’);
ylabel(‘y_b(n)’);
title(‘Filter output’);
在这里插入图片描述

2、考虑模拟信号在Ts=0.01s, 0.05s和0.1s间隔采样得到x(n)。
（1）对每一个Ts，画出x(n)。
（2）采用sinc内插（用Δt=0.001）用样本x(n)重建模拟信号ya(t)，并从图中求出在ya(t)中的频率（不管末端效果）。
（3）采用三次样条内插从样本x(n)重建模拟信号，并从图中求出在ya(t)中的频率（不管末端效果）。
相关m文件：
（1）
t=0:0.001:1;xa=cos(20pit);
Ts=0.01;N1=round(1/Ts);n1=0:N1;x1=cos(20pin1Ts);
subplot(3,1,1);plot(t,xa,n1Ts,x1,‘o’);
ylabel(‘x_1(n)’);title(‘Sampling of x_a(t) using Ts=0.01’);
Ts=0.05;N2=round(1/Ts);n2=0:N2;x2=cos(20pin2Ts);
subplot(3,1,2);plot(t,xa,n2Ts,x2,‘o’);
ylabel(‘x_2(n)’);title(‘Sampling of x_a(t) using Ts=0.05’);
Ts=0.1;N3=round(1/Ts);n3=0:N3;x3=cos(20pin3Ts);
subplot(3,1,3);plot(t,xa,n3Ts,x3,‘o’);
ylabel(‘x_3(n)’);title(‘Sampling of x_a(t) using Ts=0.1’);
在这里插入图片描述

（2）
Ts=0.01;Fs=1/Ts;
xa1=x1sinc(Fs(ones(length(n1),1)t-(n1Ts)'ones(1,length(t))));
subplot(3,1,1);plot(t,xa1);axis([0,1,-1.1,1.1]);
ylabel(‘x_a(t)’);title(‘Reconstruction of x_a(t) when Ts=0.01’);
Ts=0.05;Fs=1/Ts;
xa2=x2sinc(Fs*(ones(length(n2),1)t-(n2Ts)‘ones(1,length(t))));
subplot(3,1,2);plot(t,xa2);axis([0,1,-1.1,1.1]);
ylabel(‘x_a(t)’);title(‘Reconstruction of x_a(t) when Ts=0.05’);
Ts=0.1;Fs=1/Ts;
xa3=x3sinc(Fs*(ones(length(n3),1)t-(n3Ts)’*ones(1,length(