实验三 离散傅里叶变换及性质

一，实验名称： DFT 的频谱分析 二，实验目的：1. 加深对 DFT 原理的理解，熟悉DFT 的性质。

2. 掌握离散傅里叶变换的有关性质，利用Matlab 实现DFT 变换3. 深刻理解利用 DFT 分析信号频谱的原理，分析实现过程中出现的现象及解决方法三，实验原理：所谓信号的频谱分析就是计算信号的傅里叶变换。

连续信号与系统的傅里叶分析显然不便于直接用计算机进行计算，使其应用受到限制，而DFT 是一种时域和频域均离散化的变换，适合数值运算，成为分析离散信号和系统的有力工具。

工程实际中，经常遇到的连续信号Xa(t),其频谱函数Xa(jW)也是连续函数。

数字计算机难于处理，因而我们采用DFT 来对连续时间信号的傅里叶变换进行逼近，进而分析连续时间信号的频谱。

离散傅里叶变换是有限长序列的傅里叶变换，它相当于把信号的傅里叶变换进行等频率间隔采样，并且有限长序列的离散傅里叶变换和周期序列的离散傅里叶级数本质是一样的。

快速傅里叶变换（FFT ）并不是一种新的变换，它是离散傅里叶变换的一种快速算法，并且主要是基于这样的思路而发展起来的：（1）把长度为N 的序列的DFT 逐次分解成长度较短的序列的DFT 来计算。

（2）利用WN(nk)的周期性和对称性，在DFT 运算中适当的分类，以提高运算速度。

（对称性nkNnk NW W N-=+2，12-=NN W ;周期性nk N nk N nrN N k rN n N W W W W ---==)(，r 为任意整数,1=nrNNW ） 离散傅里叶变换的推导：离散傅里叶级数定义为 nk j N k p p ek x Nn x N21)(1)(π∑-== （1-1） 将上式两端乘以nm j Ne π2-并对n 在0~N-1求和可得 ⎥⎦⎤⎢⎣⎡==∑∑∑∑∑-=---=-=-=---=-10)(110101)(10N2N2N2)()(1)(N n m k n j N N k p N n N k m k n j pN n nm j p e k X ek XNen x πππ 因为{m k 1mk 0)(N )(1)(N 2N2N2-1-1N 11=≠---=-==∑m k j m k j N n m k n je eeNπππ所以∑∑-=-=--=11)()()(N2N k p N n nm j p m k k X en x δπ 这样∑-=-=10N2)()(N n nm j p p en x m X π用k代替m 得 ∑-=-=10N2)()(N n nk j p P e n x k X π （1-2）令N2πj N eW -=，则（1-2）成为DFS []∑-===10)()()(N n nk N p p p W n x k X n x （1-3）（1-1）成为 IDFS []∑-=-==1)(1)()(N n nkNpp p W k XNn x k X （1-4） 式（1-3）、（1-4）式构成周期序列傅里叶级数变换关系。

第三章离散傅里叶变换DFT: Discrete Fourier Transform第三章学习目标z理解傅里叶变换的几种形式z掌握离散傅里叶变换(DFT)及性质，圆周移位、共轭对称性，掌握圆周卷积、线性卷积及两者之间的关系z掌握频域抽样理论z掌握DFT的应用引言DFT要解决两个问题：一是频谱的离散化；二是算法的快速计算(FFT)。

这两个问题都是为了使计算机能够实时处理信号。

Fourier变换的几种可能形式时间函数频率函数连续时间、连续频率—傅里叶变换连续时间、离散频率—傅里叶级数离散时间、连续频率—序列的傅里叶变换可以得出一般的规律：一个域的离散对应另一个域的周期延拓；一个域的连续必定对应另一个域的非周期。

−jwndw e jwn 时域离散、非周期频域连续、周期z 时域周期化→频域离散化z 时域离散化→频域周期化离散连续周期性非周期性引言Fourier变换的几种可能形式时间函数频率函数连续时间、连续频率—傅里叶变换连续时间、离散频率—傅里叶级数离散时间、连续频率—序列的傅里叶变换离散时间、离散频率—周期序列的傅里叶级数由DTFT到DFS离散时间、离散频率的傅立叶级数（DFS）由上述分析可知，对DTFT，要想在频域上离散化，那么在时域上必须作周期延拓。

对长度为M的有限长序列x(n)，以N为周期延拓（N≥M）。

注意：周期序列的离散傅里叶级数（DFS）只对有限长序列作周期延拓或周期序列成立。

……四种傅里叶变换形式的归纳时间函数频率函数连续和非周期非周期和连续连续和周期(T0)非周期和离散(Ω=2π/T)离散(T)和非周期周期(Ωs=2π/T)和连续离散(T)和周期(T0)周期(Ωs=2π/T)和离散(Ω=2π/T)在进行DFS 分析时，时域、频域序列都是无限长的周期序列周期序列实际上只有有限个序列值有意义长度为N 的有限长序列可以看成周期为N 的周期序列的一个周期（主值序列）借助DFS 变换对，取时域、频域的主值序列可以得到一个新的变换—DFT ，即有限长序列的离散傅里叶变换3.1 离散傅里叶变换（DFT ）的定义及物理意义——有限长序列的离散频域表示x(n)的N 点DFT 是¾x(n)的z 变换在单位圆上的N 点等间隔抽样；¾x(n)的DTFT 在区间[0，2π）上的N 点等间隔抽样。

实验三用FFT对信号进行频谱分析和MATLAB程序实验三中使用FFT对信号进行频谱分析的目的是通过将时域信号转换为频域信号，来获取信号的频谱信息。

MATLAB提供了方便易用的函数来实现FFT。

首先，我们需要了解FFT的原理。

FFT（快速傅里叶变换）是一种快速计算离散傅里叶变换（DFT）的算法，用于将离散的时间域信号转换为连续的频域信号。

FFT算法的主要思想是将问题划分为多个规模较小的子问题，并利用DFT的对称性质进行递归计算。

FFT算法能够帮助我们高效地进行频谱分析。

下面是一个使用MATLAB进行频谱分析的示例程序：```matlab%生成一个10秒钟的正弦波信号，频率为1Hz，采样率为100Hzfs = 100; % 采样率t = 0:1/fs:10-1/fs; % 时间范围f=1;%正弦波频率x = sin(2*pi*f*t);%进行FFT计算N = length(x); % 信号长度X = fft(x); % FFT计算magX = abs(X)/N; % 幅值谱frequencies = (0:N-1)*(fs/N); % 频率范围%绘制频谱图figure;plot(frequencies, magX);xlabel('频率(Hz)');ylabel('振幅');title('信号频谱');```上述代码生成了一个10秒钟的正弦波信号，频率为1 Hz，采样率为100 Hz。

通过调用MATLAB的fft函数计算信号的FFT，然后计算每个频率分量的幅值谱，并绘制出信号频谱图。

在频谱图中，横轴表示频率，纵轴表示振幅。

该实验需要注意以下几点：1.信号的采样率要与信号中最高频率成一定比例，以避免采样率不足导致的伪频谱。

2.FFT计算结果是一个复数数组，我们一般只关注其幅值谱。

3.频率范围是0到采样率之间的频率。

实验三的报告可以包含以下内容：1.实验目的和背景介绍。

实验三 z 变换及分析、DTFT 实验一、 实验目的（1） 学会运用MATLAB 求离散时间信号的z 变换和z 反变换； （2） 学会运用MATLAB 分析离散时间系统的系统函数的零极点；（3） 学会运用MATLAB 分析系统函数的零极点分布与其时域特性的关系； （4） 学会运用MATLAB 进行离散时间系统的频率特性分析。

二、实验原理及实例分析2.1 z 正反变换序列()n x 的z 变换定义为()()[]()∑∞-∞=-==n nzn x n x z X Z （1）其中，符号Z 表示取z 变换，z 是复变量。

相应地，单边z 变换定义为()()[]()∑∞=-==0n n z n x n x z X Z （2）MATLAB 符号数学工具箱提供了计算离散时间信号单边z 变换的函数ztrans 和z 反变换函数iztrans ，其语句格式分别为Z=ztrans(x) x=iztrans(z)上式中的x 和Z 分别为时域表达式和z 域表达式的符号表示，可通过sym ：函数来定义。

【实例1】 试用ztrans 函数求下列函数的z 变换。

（1）)()cos()(n u n a n x nπ=； （2）)(])2(2[)(11n u n x n n ----=。

解：（1）z 变换MATLAB 源程序为 >>x=sym('a^n*cos(pi*n)'); >>Z=ztrans(x);>>simplify(Z) ans=z/(z+a)（2）z 变换MATLAB 源程序为 >>x=sym('2^(n-1)-(-2)^(n-1)'); >>Z=ztrans(x); >>simplify(Z) ans=z^2/(z-2)/(z+2)【实例2】 试用iztrans 函数求下列函数的z 反变换。

（1）65198)(2+--=z z z z X （2）32)2)(1()12112()(--+-=z z z z z z X解：（1）z 反变换MATLAB 源程序为 >>Z=sym('(8*z-19)/(z^2-5*z+6)'); >>x=iztrans(Z); >>simplify(x) ans=-19/6*charfcn[0](n)+5*3^(n-1)+3*2^(n-1)其中，charfcn[0](n)是)(n δ函数在MATLAB 符号工具箱中的表示，反变换后的函数形式为)()2335()(619)(11n u n n x n n --⨯+⨯+-=δ。

FFT算法分析实验实验报告一、实验目的快速傅里叶变换（Fast Fourier Transform，FFT）是数字信号处理中一种非常重要的算法。

本次实验的目的在于深入理解 FFT 算法的基本原理、性能特点，并通过实际编程实现和实验数据分析，掌握 FFT 算法在频谱分析中的应用。

二、实验原理FFT 算法是离散傅里叶变换（Discrete Fourier Transform，DFT）的快速计算方法。

DFT 的定义为：对于长度为 N 的序列 x(n)，其 DFT 为X(k) ＝∑n=0 到 N-1 x(n) e^（j 2π k n ／ N) ，其中 j 为虚数单位。

FFT 算法基于分治法的思想，将 N 点 DFT 分解为多个较小规模的DFT，从而大大减少了计算量。

常见的 FFT 算法有基 2 算法、基 4 算法等。

三、实验环境本次实验使用的编程语言为 Python，主要依赖 numpy 库来实现 FFT 计算和相关的数据处理。

四、实验步骤1、生成测试信号首先，生成一个包含不同频率成分的正弦波叠加信号，例如100Hz、200Hz 和 300Hz 的正弦波。

设定采样频率为 1000Hz，采样时间为 1 秒，以获取足够的采样点进行分析。

2、进行 FFT 计算使用 numpy 库中的 fft 函数对生成的测试信号进行 FFT 变换。

3、频谱分析计算 FFT 结果的幅度谱和相位谱。

通过幅度谱确定信号中各个频率成分的强度。

4、误差分析与理论上的频率成分进行对比，计算误差。

五、实验结果与分析1、幅度谱分析观察到在 100Hz、200Hz 和 300Hz 附近出现明显的峰值，对应于生成信号中的频率成分。

峰值的大小反映了相应频率成分的强度。

2、相位谱分析相位谱显示了各个频率成分的相位信息。

3、误差分析计算得到的频率与理论值相比，存在一定的误差，但在可接受范围内。

误差主要来源于采样过程中的量化误差以及 FFT 算法本身的近似处理。

信号实验一离散傅里叶变换及其快速算法一、实验目的1、掌握计算序列的离散傅里叶变换（FFT）的方法；2、掌握实现时间抽取快速傅里叶变换（FFT）编程方法；3、加深对DFT与序列的傅里叶变换和Z变换之间的关系的理解；4、复习复数序列的运算方法。

二、程序设计框图1.码位倒置程序框图2.蝶形图运算程序框图三、实验程序实验程序的源代码如下:#include"math.h"#include"stdio.h"/*------------------------------------------------------------------------------------------子函数部分------------------------------------------------------------------------------------------*/ void swap(float *a,float *b)//交换变量子函数{float T;T=*a;*a=*b;*b=T;}void fft (float A [],float B [],unsigned M)//数组A为序列的实部, 数组B为序列的虚部{unsigned long N,I,J,K,L,LE,LE1,P,Q,R;float Wr,Wi,W1r,W1i,WTr,WTi,theta,Tr,Ti;N=1<<M;J=0;for(I=0;I<N-1;I++){if(J>I){swap(&A [I],&A [J]);swap(&B [I],&B [J]);}K=N>>1;while(K>=2&&J>=K){J-=K;K>>=1;}J+=K;}for(L=1;L<=M;L++){LE=1<<L;LE1=LE/2;Wr=1.0;Wi=0.0;theta=(-1)*3.1415926536/LE1;W1r=cos (theta);W1i=sin (theta);for(R=0;R<LE1;R++){for(P=R;P<N-1;P+=LE){Q=P+LE1;//基本蝶形图的复数运算Tr=Wr*A[Q]-Wi*B[Q];Ti=Wr*B[Q]+Wi*A[Q];A[Q]=A[P]-Tr;B[Q]=B[P]-Ti;A[P]+=Tr;B[P]+=Ti;}WTr=Wr;WTi=Wi;Wr=WTr*W1r-WTi*W1i;Wi=WTr*W1i+WTi*W1r;}}return;}/*------------------------------------------------------------------------------------------主函数部分------------------------------------------------------------------------------------------*/ void main(){float A[20],B[20];char t1,t2,file_name[20];int M,N,i,iiff;FILE *fp;/*************************************数据读取部分************************************/ printf("请输入文件名:");//输入数据文件名scanf("%s",file_name);printf("FFT变换还是IFFT变换?(FFT:1,IFFT:-1):");//输入变换方式, 1为FFT, -1为IFFTscanf("%d",&iiff);while(iiff!=1&&iiff!=-1)//检错: 检验上一步的输入是否有错, 有错则重新输入{printf("输入错误, 请重新输入！ ");printf("FFT or IFFT?(FFT:1,IFFT:-1):");scanf("%d",&iiff);}fp=fopen(file_name,"r");//打开文件并读入数据fscanf(fp,"%d",&M);N=pow(2,M);//计算序列总数for(i=0;i<N;i++)//读取文件中的数据{fscanf(fp,"%f%c%c%f",&A[i],&t1,&t2,&B[i]);if(iiff==-1)//根据FFT或IFFT修正BB[i]=B[i]*-1;if(t2!='j')//检错: 检验读取格式是否有错{printf("输入格式错误\n");break;}if(t1=='+')//判断虚部的正负号B[i]=B[i];else if(t1=='-')B[i]=-B[i];}/****************************************变换部分****************************************/ fft(A,B,M);//FFT变换/**************************************数据输出部分**************************************/ fp=fopen("fft_result.txt","w"); //输出结果if(iiff==-1)fprintf(fp,"IFFT变换的输出结果是: \n");elsefprintf(fp,"FFT变换的输出结果是: \n");for(i=0;i<N;i++){if(iiff==-1) //根据FFT或IFFT修正B{B[i]=B[i]*-1/N;A[i]=A[i]/N;}if(B[i]>=0)//修正虚部的输出格式fprintf(fp,"%f+j%f\n",A[i],B[i]);else if(B[i]<0)fprintf(fp,"%f-j%f\n",A[i],-B[i]);else if(B[i]==0)fprintf(fp,"%f\n",A[i]);}fclose(fp);}四、程序运行结果检验（1） 1.对序列进行FFT变换输入文件fft_input.txt:21+j02+j0-1+j04+j0控制台输入:请输入文件名: fft_input.txtFFT变换还是IFFT变换?(FFT:1,IFFT:-1): 1输出文件fft_result.txt:FFT变换的输出结果是:6.00000+j0.000002.00000+j2.00000-6.00000+j0.000002.00000+j-2.00000运行结果分析:程序运行输出结果与计算结果相同, 表示傅里叶正变换（FFT）成功。

实验二
连续非周期信号的傅里叶变换（FT）及其性质一、实验目的
在理论学习的基础上，通过本实验熟悉常见信号的傅里叶变换及掌握连续时间傅里叶变换的性质。

二、相关知识
常见信号的傅里叶变换和连续时间傅里叶变换(CTFT)的性质
1、常见连续时间非周期信号及其傅里叶变换列表如下：
在本实验中可以可以对以上信号采取以下常见运算，运算结果表达式列表如下：
三、思考问题
1、X(w)和C k在量纲上分别有什么区别？
2、C k和X(w)是否分别代表周期信号和非周期信号各频率分量的振幅？
3、如果对X(w)在频域进行抽样，即令X(w)用X(KW0)代替，那么在时域对信号会产生什么影响？。

第1篇实验目的1. 理解并掌握信号变换的基本原理和方法，包括傅里叶变换、拉普拉斯变换和Z 变换。

2. 通过实验，加深对信号变换在实际应用中的理解，如信号处理、系统分析等。

3. 学习使用MATLAB等工具进行信号变换的计算和分析。

实验原理信号变换是信号处理领域的重要工具，它可以将信号从时域转换为频域或其他域，便于分析、处理和设计。

本实验主要涉及以下几种信号变换：1. 傅里叶变换（FT）：将时域信号转换为频域信号，揭示信号的频率成分。

2. 拉普拉斯变换（LT）：将时域信号转换为复频域信号，适用于分析线性时不变系统。

3. Z变换（ZT）：将离散时间信号转换为Z域信号，适用于分析离散时间系统。

实验器材1. MATLAB软件2. 示波器3. 信号发生器4. 连接线实验步骤1. 傅里叶变换实验- 使用MATLAB生成一个正弦波信号，采样频率为1000Hz，采样点数为1024。

- 对该信号进行傅里叶变换，观察频谱图。

- 改变采样频率和采样点数，观察频谱图的变化。

2. 拉普拉斯变换实验- 使用MATLAB生成一个指数衰减信号。

- 对该信号进行拉普拉斯变换，观察复频域图。

- 分析拉普拉斯变换的结果，解释信号的特性。

3. Z变换实验- 使用MATLAB生成一个离散时间信号。

- 对该信号进行Z变换，观察Z域图。

- 分析Z变换的结果，解释信号的特性。

4. 信号变换在系统分析中的应用- 使用MATLAB设计一个简单的模拟滤波器。

- 对滤波器进行傅里叶变换，观察滤波器的频率响应。

- 分析滤波器的性能，如通带、阻带、截止频率等。

实验结果与分析1. 傅里叶变换实验- 频谱图显示正弦波的频率成分与采样频率有关。

- 改变采样频率和采样点数，频谱图发生相应的变化。

2. 拉普拉斯变换实验- 复频域图显示指数衰减信号的极点和零点。

- 分析结果，可以确定信号的衰减速度和稳态值。

3. Z变换实验- Z域图显示离散时间信号的极点和零点。

- 分析结果，可以确定信号的稳定性、收敛速度等特性。

实验一时域离散信号的产生及时域处理实验目的：了解Matlab软件数字信号处理工具箱的初步使用方法。

掌握其简单的Matlab语言进行简单的时域信号分析。

实验内容：[1.1]已知两序列x1=[0,1,2,3,4,3,2,1,0];n1=[-2:6];x2=[2,2,0,0,0,-2,-2],n2=[2:8].求他们的和ya及乘积yp. 程序如下：x1=[0,1,2,3,4,3,2,1,0];ns1=-2;x2=[2,2,0,0,0,-2,-2];ns2=2;nf1=ns1+length(x1)-1;nf2=ns2+length(x2)-1;ny=min(ns1,ns2):max(nf1,nf2);xa1=zeros(1,length(ny));xa2=xa1;xa1(find((ny>=ns1)&(ny<=nf1)==1))=x1;xa2(find((ny>=ns2)&(ny<=nf2)==1))=x2;ya=xa1+xa2yp=xa1.*xa2subplot(4,4,1),stem(ny,xa1,'.')subplot(4,1,2),stem(ny,xa2,'.')line([ny(1),ny(end)],[0,0])subplot(4,1,3),stem(ny,ya,'.')line([ny(1),ny(end)],[0,0])subplot(4,1,4),stem(ny,yp,'.')line([ny(1),ny(end)],[0,0])[1.2]编写产生矩形序列的程序。

并用它截取一个复正弦序列，最后画出波形。

程序如下：clear;close alln0=input('输入序列起点:n0=');N=input('输入序列长度:N=');n1=input('输入位移:n1=');n=n0:n1+N+5;u=[(n-n1)>=0];x1=[(n-n1)>=0]-[(n-n1-N)>=0];x2=[(n>=n1)&(n<(N+n1))];x3=exp(j*n*pi/8).*x2;subplot(2,2,1);stem(n,x1,'.');xlabel('n');ylabel('x1(n)');axis([n0,max(n),0,1]);subplot(2,2,3);stem(n,x2,'.');xlabel('n');ylabel('x2(n)');axis([n0,max(n),0,1]);subplot(2,2,2);stem(n,real(x3),'.'); xlabel('n');ylabel('x3(n)的实部');line([n0,max(n)],[0,0]);axis([n0,max(n),-1,1]);subplot(2,2,4);stem(n,imag(x3),'.'); xlabel('n');ylabel('x3(n)的虚部');line([n0,max(n)],[0,0]);axis([n0,max(n),-1,1]);[1.3]利用已知条件，利用MATLAB生成图形。

傅里叶光学实验傅里叶光学原理的发明最早可以追溯到1893年阿贝（Abbe ）为了提高显微镜的分辨本领所做的努力。

他提出一种新的相干成象的原理，以波动光学衍射和干涉的原理来解释显微镜的成像的过程，解决了提高成像质量的理论问题。

1906年波特（Porter ）用实验验证了阿贝的理论。

1948年全息术提出，1955年光学传递函数作为像质评价兴起，1960年由于激光器的出现使相干光学的实验得到重新装备，因此从上世纪四十年代起古老的光学进入了“现代光学”的阶段，而现代光学的蓬勃发展阶段是从上世纪六十年代起开始。

由于阿贝理论的启发，人们开始考虑到光学成像系统与电子通讯系统都是用来收集、传递或者处理信息的，因此上世纪三十年代后期起电子信息理论的结果被大量应用于光学系统分析中。

两者一个为时间信号，一个是空间信号，但都具有线性性和不变性，所以数学上都可以用傅立叶变换的方法。

将光学衍射现象和傅立叶变换频谱分析对应起来，进而应用于光学成像系统的分析中，不仅是以新的概念来理解熟知的物理光学现象，而且使近代光学技术得到了许多重大的发展，例如泽尼克相衬显微镜，光学匹配滤波器等等，因此形成了现代光学中一门技术性很强的分支学科—傅里叶光学。

实验原理：我们知道一个复变函数f(x,y)的傅立叶变换为⎰⎰+-=ℑ=dxdy vy ux 2i y x f y x f v u F )](exp[),()},({),(π ( 1 )F (u,v)叫作f(x,y)的变换函数或频谱函数。

它一般也为复变函数，f(x,y)叫做原函数，也可以通过求F(u,v)逆傅立叶变换得到原函数f(x,y)， ⎰⎰+=ℑ=-dudv vy ux 2i v u F v u F y x f 1)](exp[),()},({),(π （2）在光学系统中处理的是平面图形，当光波照明图形时从图形反射或透射出来的光波可用空间两维复变函数（简称空间函数）来表示。

在这些情况下一般都可以进行傅里叶变换或广义的傅里叶变换。

实验三用FFT对信号作频谱分析_实验报告一、实验目的1.学习使用FFT（快速傅里叶变换）对信号进行频谱分析；2.掌握频谱分析的基本原理和方法；3.熟悉使用MATLAB进行频谱分析的操作。

二、实验原理FFT是一种基于傅里叶变换的算法，可以将时域信号转换为频域信号，并将信号的频谱特征展示出来。

在频谱分析中，我们通过分析信号的频谱可以获得信号的频率、幅值等信息，从而对信号的性质和特征进行研究。

对于一个连续信号，我们可以通过采样的方式将其转换为离散信号，再利用FFT算法对离散信号进行频谱分析。

FFT算法可以将信号从时域转换到频域，得到离散的频谱，其中包含了信号的频率分量以及对应的幅值。

MATLAB中提供了fft函数，可以方便地对信号进行FFT分析。

通过对信号进行FFT操作，可以得到信号的频谱图，并从中提取出感兴趣的频率信息。

三、实验步骤1.准备工作：（2）建立新的MATLAB脚本文件。

2.生成信号：在脚本中，我们可以通过定义一个信号的频率、幅值和时间长度来生成一个信号的波形。

例如，我们可以生成一个频率为1000Hz，幅值为1的正弦波信号，并设置信号的时间长度为1秒。

3.对信号进行FFT分析：调用MATLAB中的fft函数，对信号进行FFT分析。

通过设置采样频率和FFT长度，可以得到信号的频谱。

其中，采样频率是指在单位时间内连续采样的次数，FFT长度是指离散信号的样本点数。

4.绘制频谱图：调用MATLAB中的plot函数，并设置x轴为频率，y轴为幅值，可以绘制出信号的频谱图。

频谱图上横坐标表示信号的频率，纵坐标表示信号的幅值，通过观察可以得到信号的频率分布情况。

四、实验结果在实验过程中，我们生成了一个频率为1000Hz，幅值为1的正弦波信号，并对其进行FFT分析。

通过绘制频谱图，我们发现信号在1000Hz处有最大幅值，说明信号主要由这一频率成分组成。

五、实验总结本实验通过使用FFT对信号进行频谱分析，我们可以方便地从信号的波形中提取出频率分量的信息，并绘制出频谱图进行观察。

实验二离散信号的频谱分析一、[实验目的]（1）加深对采样定理的理解和掌握，以及对信号恢复的必要性；（2）掌握对连续信号在时域的采样与重构的方法（3）理解和加深傅里叶变换的概念及其性质。

（4）离散时间傅里叶变换(DTFT)的计算和基本性质。

（5）离散傅里叶变换(DFT)的计算和基本性质。

二、[实验内容]1.实验原理验证（一）．采样定理及采样后信号的频谱对Sa(t)的采样后信号的频谱（二）．信号重建对cos(t)的采样与重建信号cos(t) cos(t)重建信号与原信号的比较及误差（三）．离散时间信号的傅立叶变换及频谱分析（1））离散时间傅里叶变换的概念及其性质。

有限长序列x(n)={1,2,3,4,5}（2）离散傅里叶变换的概念及其性质x(n)=sin(n*pi/8)+sin(n*pi/4)，N=16的序列傅里叶变换。

2. 选取信号f(t)= cos(t)作为被采样信号（最高频率为f=8Hz），取理想低通的截止频率wc=1/2*ws。

实现对信号f(t)= cos(t)的采样及由该采样信号的恢复重建，按要求完成以下内容:(1) 分别令采样角频率ws=1.5*wm 及ws=3*wm，给出在欠采样及过采样条件下冲激取样后信号的频谱，从而观察频谱的混叠现象。

答：实验程序如下clc,cleardt=0.01;t=0:dt:1;cos(t)的3倍采样信号频谱ωF (j w )f=8; %信号频率wm=2*pi*f; %信号角频率 ft=cos(wm*t); %时域信号%bs=1.5; %采样角频率，欠采样 bs=3; %采样角频率，大于两倍采样ws=bs*wm;Ts=2*pi/ws; %采样时间间隔 wc=1/2*ws; %理想低通截止频率nTs=0:Ts:1;Tf=0.01;nTf=-10:Tf:10; f_nTs=cos(wm*nTs); %时域采样信号Fs=funexer4_1(f_nTs,nTs,Ts,nTf); figure(1); plot(nTf,Fs);title('cos(t)的3倍采样信号频谱'); xlabel('ω'); ylabel('F(jw)'); grid on%//////////////////1.5倍采样 figure(2)bs=1.5; %采样角频率，大于两倍采样ws=bs*wm;Ts=2*pi/ws; %采样时间间隔wc=1/2*ws; %理想低通截止频率nTs=0:Ts:1; Tf=0.01; nTf=-10:Tf:10;Fs=funexer4_1(f_nTs,nTs,Ts,nTf); plot(nTf,Fs); title('cos(t)的1.5倍采样信号频谱');xlabel('ω');ylabel('F(jw)'); grid on(2) 若采样角频率取为ws=3*wm ，欲使输出信号与输入信号一致为cos(t)，试根据采样信号恢复信号的误差，确定理想低通滤波器H ( jw)的截止角频率Wc 的取值范围应为多大？cos(t)的1.5倍采样信号频谱ωF (j w )Sa(t)采样后的奈奎斯特采样频谱图(4倍）ωF (j ω)答：截止频率wc 应满足： wm<wc ≤ws/2。

一、实验目的1. 理解离散信号的概念及其特点。

2. 掌握离散信号的表示方法。

3. 掌握离散信号的基本运算方法。

4. 熟悉离散系统响应的求解方法。

5. 利用MATLAB进行离散信号分析。

二、实验原理离散信号是指时间上不连续的信号，与连续信号相比，具有以下特点：1. 采样性：离散信号是在时间上等间隔取样的信号。

2. 有限性：离散信号在时间上有限，即在有限的时间内存在。

3. 线性时不变性：离散系统具有线性时不变性，即系统对信号的时延和幅度变换保持不变。

离散信号的表示方法主要有以下几种：1. 序列表示法：用括号括起来的序列表示，如x[n]。

2. 图形表示法：用坐标轴表示，横轴为时间，纵轴为信号幅度。

3. Z变换表示法：用Z变换表示，如X(z)。

离散信号的基本运算方法包括：1. 加法运算：两个离散信号相加，结果为它们的序列对应元素相加。

2. 乘法运算：两个离散信号相乘，结果为它们的序列对应元素相乘。

3. 移位运算：将离散信号沿时间轴左移或右移。

4. 展平运算：将离散信号沿时间轴展平，即将信号序列展开成矩阵形式。

离散系统响应的求解方法主要有以下几种：1. 离散卷积法：用离散卷积运算求解离散系统响应。

2. Z变换法：用Z变换求解离散系统响应。

3. 快速傅里叶变换（FFT）法：用FFT求解离散系统响应。

三、实验内容及步骤1. 实验一：离散信号的表示方法（1）在MATLAB中，创建一个离散信号序列x[n]，并绘制其图形表示。

（2）利用Z变换，将离散信号序列转换为Z变换表示。

2. 实验二：离散信号的基本运算（1）在MATLAB中，创建两个离散信号序列x[n]和y[n]，并进行加法运算、乘法运算、移位运算和展平运算。

（2）绘制运算结果，并分析运算结果的特点。

3. 实验三：离散系统响应的求解（1）在MATLAB中，创建一个离散信号序列x[n]，并设计一个离散系统。

（2）利用离散卷积法、Z变换法和FFT法求解离散系统响应。

fft实验报告傅里叶变换（Fast Fourier Transform, FFT）是一种重要的数学工具，广泛应用于信号处理、图像处理、通信系统等领域。

本文将从理论和实验两个方面，介绍FFT的原理、应用以及实验结果。

一、FFT的原理FFT是一种将时域信号转换为频域信号的算法，它基于傅里叶级数展开的思想。

傅里叶级数展开可以将一个周期信号表示为一系列正弦和余弦函数的叠加，而FFT则能够将非周期信号分解成一系列频率成分。

FFT的核心思想是将一个N点的离散信号变换为N/2个频率分量，其中前一半为正频率分量，后一半为负频率分量。

通过分别计算正频率和负频率的离散傅里叶变换（DFT），再利用对称性质进行合并，最终得到频域信号。

二、FFT的应用1. 信号处理：FFT在信号处理中有广泛应用，例如音频信号的频谱分析、滤波、降噪等。

通过将信号转换到频域，可以方便地分析信号的频率成分，从而实现各种信号处理算法。

2. 图像处理：FFT在图像处理中也有重要应用。

通过对图像进行二维FFT变换，可以将图像转换为频域表示，从而实现图像增强、去噪、压缩等操作。

例如，图像的频域滤波可以有效地去除图像中的噪声，提高图像的质量。

3. 通信系统：FFT在通信系统中也扮演着重要角色。

例如，在OFDM（正交频分复用）系统中，FFT用于将多个子载波的频域信号转换为时域信号进行传输。

这种技术能够提高信号的传输效率和抗干扰能力。

三、FFT实验结果为了验证FFT算法的正确性和效果，我们进行了一系列实验。

首先，我们使用MATLAB编程实现了FFT算法，并将其应用于音频信号处理。

通过对一段音频信号进行FFT变换，我们成功地获得了该信号的频谱图，并观察到不同频率成分的存在。

接下来，我们将FFT算法应用于图像处理。

我们选择了一张包含噪声的图像，并对其进行FFT变换。

通过对频域图像进行滤波操作，我们成功去除了图像中的噪声，并获得了清晰的图像。

最后，我们将FFT算法应用于通信系统中的OFDM技术。

第1篇一、实验背景随着信息技术的飞速发展，数字信号处理（DSP）技术已成为通信、图像处理、语音识别等领域的重要工具。

本实验旨在通过一系列实验，加深对数字信号处理基本原理和方法的理解，提高实际应用能力。

二、实验目的1. 理解数字信号处理的基本概念和原理。

2. 掌握常用信号处理算法的MATLAB实现。

3. 培养分析和解决实际问题的能力。

三、实验内容本实验共分为五个部分，具体如下：1. 离散时间信号的基本操作（1）实验目的：熟悉离散时间信号的基本操作，如加法、减法、乘法、除法、延时、翻转等。

（2）实验步骤：- 使用MATLAB生成两个离散时间信号。

- 对信号进行基本操作，如加法、减法、乘法、除法、延时、翻转等。

- 观察并分析操作结果。

2. 离散时间系统的时域分析（1）实验目的：掌握离散时间系统的时域分析方法，如单位脉冲响应、零状态响应、零输入响应等。

（2）实验步骤：- 使用MATLAB设计一个离散时间系统。

- 计算系统的单位脉冲响应、零状态响应和零输入响应。

- 分析系统特性。

（1）实验目的：掌握离散时间信号的频域分析方法，如快速傅里叶变换（FFT）、离散傅里叶变换（DFT）等。

（2）实验步骤：- 使用MATLAB生成一个离散时间信号。

- 对信号进行FFT和DFT变换。

- 分析信号频谱。

4. 数字滤波器的设计与实现（1）实验目的：掌握数字滤波器的设计与实现方法，如巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器等。

（2）实验步骤：- 使用MATLAB设计一个低通滤波器。

- 使用窗函数法实现滤波器。

- 对滤波器进行性能分析。

5. 信号处理在实际应用中的案例分析（1）实验目的：了解信号处理在实际应用中的案例分析，如语音信号处理、图像处理等。

（2）实验步骤：- 选择一个信号处理应用案例。

- 分析案例中使用的信号处理方法。

- 总结案例中的经验和教训。

四、实验结果与分析1. 离散时间信号的基本操作实验结果表明，离散时间信号的基本操作简单易懂，通过MATLAB可以实现各种操作，方便快捷。

快速傅里叶变换实验报告一、实验目的（一）加深对几个特殊概念的理解：“采样”——“混叠”；“窗函数”（截断）——“泄漏”；非整周期截取”——“栅栏”。

（二）加深理解如何才能避免“混叠”，减少“泄漏”，防止“栅栏”的方法和措施以及估计这些因素对频谱的影响。

（三）对利用通用微型计算机及相应的FFT 软件，实现频谱分析有一个初步的了解。

二、实验原理为了实现信号的数字化处理，利用计算机进行频谱分析——计算信号的频谱。

由于计算机只能进行有限的离散计算（即DFT ），因此就要对连续的模拟信号进行采样和截断。

而这两个处理过程可能引起信号频谱的畸变，从而使DFT 的计算机过于信号的实际频谱有误差。

有时由于采样和截断的处理不当，使计算出来的频谱完全失真。

因此在时域处理信号时要格外小心。

在信号数字化处理中应十分注意以下几点：（一）为了避免“混叠”，要求在采样时必须满足采样定理。

（二）为了减少“泄漏”，应适当增加截断长度和选择合适的窗。

（三）对信号进行整周期截取，则能消除“栅栏效应”。

（四）增加截断长度，则可提高频率分辨率。

三、实验内容及步骤（一）基本信号的FFT 变换1、()000sin(sin 2cos36x t w t w t w t π=+++第1组：采样频率,截断长度N=1608s f f =程序清单：n=16;%截取长度multi=8;%采样频率倍数x=0;%初始化向量，维数待定w0=2*pi;%设定基准频率for var=1:1:nx(var)=sin(w0/multi*(var-1)+pi/6)+sin(2*w0/multi*(var-1))+cos(3*w0/multi*(var-1));endy=fft(x);y=fftshift(y);ang=angle(y)/pi*180;altitude=abs(y)/n;var=1:1:n;subplot(1,2,1);bar(var,altitude,0.3);title('幅频图');xlabel('w');ylabel('幅值');subplot(1,2,2);bar(var,ang,0.3);colormap ([0 1 1]);title('相频图');xlabel('w');ylabel('相位');ang =Columns 1 through 80 -24.1998 0.0000 -155.8595 90.0000 81.2485 60.0000 157.9674 Columns 9 through 16180.0000 -157.9674 -60.0000 -81.2485 -90.0000 155.8595 -0.0000 24.1998altitude=Columns 1 through 80.0000 0.0000 0.5000 0.0000 0.5000 0.0000 0.5000 0.0000 Columns 9 through 160.0000 0.0000 0.5000 0.0000 0.5000 0.0000 0.5000 0.0000分析：1、频率分辨率是00.5sf f Nδ==2、x(t)的信号频率成分中的最高频率，满足采样定理，所以DFT 结果没有频率032sf f <混叠现象。