资源描述：
《北航测试信号与处理技术实验》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、测试信号与处理技术实验报告报告人学号及姓名：实验地点实验时间实验名称：信号分析处理综合练习一、实验目的要求学生会用MATLAB语言进行编程，绘出所求波形，并且运用FFT求对连续信号进行分析。二、实验要求1、用Matlab产生正弦波，矩形波，并显示各自的时域波形图；2、进行FFT变换，显示各自频谱图，其中采样率、频率、数据长度自选，要求注明；3、绘制三种信号的均方根图谱；4、用IFFT回复信号，并显示恢复的正弦信号时域波形图。三、实验原理用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行频谱分析的信号是模拟信号和时域离散信号。频谱分
2、辨率直接和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对信号进行谱分析。四、实验步骤1.打开matlab并设定采样频率、正弦信号频率、采样点等数据；2.输入源
3、程序绘制正弦信号图形；3.进行FFT变换，并绘制频谱图和方均根图；1.进行IFFT变换恢复信号；2.输入源程序绘制矩形信号图形；3.进行FFT变换，并绘制频谱图和方均根图；4.进行IFFT变换恢复信号。一、实验结果二、结果分析说明及实验感想通过对本次应用MATLAB对信号进行频谱分析及滤波设计，使我们更加系统的理解课本所学的相关知识。如果序列x(n)主值序列长度为M，最高频率为fc只有当频率采样点N≥M，采样频率fs≥2fc时，才有x(n)=IDFT[X(k)]=x(n)可有频域采样X（k）恢复原序列x(n),否则产生时域混叠失真现象。一、源程序
4、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);subpl
5、ot(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:lengt
6、h(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;
7、%进行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('均方根谱');title('矩形波均方根谱')
8、;grid;%用IFFT恢复原始信号xifft=ifft(y);magx=real(xifft);ti=[0:length(xifft)
