北航测试信工号与处理技术实验.docVIP

测试信号与处理技术 实验报告 报告人学号及姓名： 实验地点 实验时间 实验名称：信号分析处理综合练习 实验目的 要求学生会用MATLAB语言进行编程，绘出所求波形，并且运用FFT求对连续信号进行分析。 实验要求 1、用Matlab产生正弦波，矩形波，并显示各自的时域波形图； 2、进行FFT变换，显示各自频谱图，其中采样率、频率、数据长度自选，要求注明； 3、绘制三种信号的均方根图谱； 4、用IFFT回复信号，并显示恢复的正弦信号时域波形图。 实验原理 用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行频谱分析的信号是模拟信号和时域离散信号。频谱分辨率直接和FFT的变换区间N有关，因为FFT能够实现频率分辨率是2π/N。 x(n)是一个长度为M的有限长序列，则x(n)的N点离散傅立叶变换为：X(k)=DFT[x(n)]=,k=0,1,...,N-1 逆变换：x(n) =IDFT[X(k)]= ,k=0,1,...,N-1 但FFT是一种比DFT更加快速的一种算法，提高了DFT的运算速率,为数字信号处理技术应用于各种信号处理创造了条件，大大提高了数字信号处理技术的发展。本实验就是采用FFT，IFFT对信号进行谱分析。 实验步骤 打开matlab并设定采样频率、正弦信号频率、采样点等数据； 输入源程序绘制正弦信号图形； 进行FFT变换，并绘制频谱图和方均根图； 进行IFFT变换恢复信号； 输入源程序绘制矩形信号图形； 进行FFT变换，并绘制频谱图和方均根图； 进行IFFT变换恢复信号。 实验结果 结果分析说明及实验感想 通过对本次应用MATLAB对信号进行频谱分析及滤波设计，使我们更加系统的理解课本所学的相关知识。如果序列x(n)主值序列长度为M，最高频率为fc只有当频率采样点N≥M，采样频率fs≥2fc时，才有x(n)=IDFT[X(k)]=x(n)可有频域采样X（k）恢复原序列x(n),否则产生时域混叠失真现象。 源程序 fs=100;%设定采样频率 N=64; n=0:N-1; t=n/fs; f0=5;%设定正弦信号频率 %生成正弦信号 x=sin(2*pi*f0*t); figure(1); subplot(221); plot(t,x);%作正弦信号的时域波形 xlabel(t); ylabel(y); title(正弦信号y=2*pi*6t时域波形); grid; %进行FFT变换并做频谱图 y=fft(x,N);%进行fft变换 mag=abs(y);%求幅值 f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换 figure(1); subplot(222); plot(f,mag);%做频谱图 axis([0,100,0,30]); xlabel(频率(Hz)); ylabel(幅值); title(正弦信号y=2*pi*6t幅频谱图N=64); grid; %均方根谱 sq=abs(y); figure(1); subplot(223); plot(f,sq); xlabel(频率(Hz)); ylabel(均方根谱); title(正弦信号y=2*pi*6t均方根谱); grid; %IFFT恢复原始信号 xifft=ifft(y); magx=real(xifft); ti=[0:length(xifft)-1]/fs; figure(1); subplot(224); plot(ti,magx); xlabel(t); ylabel(y); title(IFFT转换的正弦信号波形); grid; %矩形 fs=10;%设定采样频率 t=-10:0.1:10; x=square(t); x=x(1:190); figure(2); subplot(221); plot(t(1:190),x);%作矩形波的时域波形 axis([-10,10,-1.5,1.5]); xlabel(t); ylabel(y); title(矩形波时域波形); grid; %进行FFT变换并做频谱图 y=fft(x);%进行fft变换 mag=abs(y);%求幅值 f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换 figure(2); subplot(222); plot(f,mag);%做频谱图 xlabel(频率(Hz)); ylabel(幅值); title(矩形波幅频谱图); grid; %求均方根谱 sq=abs(y); figure(2); subplot(223); plot(f,sq); xlabel(频率(Hz)); ylabel(均方根谱);