数字信号处理实验指导书(

《数字信号处理实验指导书(》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、“数字信号处理”实验指导书（二）一、实验课程编码：105003二、实验课程名称：数字信号处理三、实验项目名称：快速傅里叶变换的MATLAB实现四、实验目的掌握应用MATLAB实现快速傅里叶的方法，即熟悉应用MATLAB实现快速傅里叶的函数。五、主要设备安装有MATLAB软件的电脑六、实验内容编写MATLAB程序，实现下面题目：已知序列x[k]=(0.8)k,(0=

2、的MATLAB函数（参考附录1）；2、通过运行附录2中提供的例题，熟悉用MATLAB实现快速傅里叶变换的基本方法；3、根据“六、实验内容”中各个题目的要求，编写MATLAB程序代码，调试程序，分析并保存结果。八、实验结果对实验练习题编写MATLAB程序并运行，在计算机上输出仿真结果。附录1主要的相关MATLAB函数y=czt[x,M,W,A]附录2例题例1、设x(n)是由两个正弦信号及白噪声的叠加，试用fft函数对其做频谱分析。[MATLAB程序]：%用fft做谱分析%产生两个正弦加白噪声N=25

3、6;f1=0.1;f2=0.2;3fs=1;a1=5;a2=3;w=2*pi/fs;x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N);%用FFT求频谱subplot(2,2,1);plot(x(1:N/4));title('原始信号');f=-0.5:1/N:0.5-1/N;X=fft(x);y=ifft(X);subplot(2,2,2);plot(f,fftshift(abs(X)));title('频域信号');subplot(2,

4、2,3);plot(real(x(1:N/4)));title('时域信号');[运行结果]：例2、利用ChirpZ变换计算滤波器h在100Hz～200Hz的频率特性。[MATLAB程序]：%利用ChirpZ变换计算滤波器h在100Hz～200Hz的频率特性h=fir1(30,125/500,boxcar(31));Fs=1000;f1=100;f2=200;m=1024;w=exp(-j*2*pi*(f2-1)/(m*Fs));a=exp(j*2*pi*f1/Fs);y=fft(h,m);z=c

5、zt(h,m,w,a);3fy=(0:length(y)-1)'*Fs/length(y);fz=(0:length(z)-1)'*(f2-f1)/length(z)+f1;subplot(2,1,1);plot(fy(1:500),abs(y(1:500)));title('FFT');subplot(2,1,2);plot(fz,abs(z));title('CZT');[运行结果]：3