资源描述：
《同样具有同样复件.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、1.FT与DFT关系的MATLAB仿真任务要求：（1）任设一个长度为N的实序列；（2）编写MATLAB代码，实现该序列的傅立叶变换和离散傅立叶变换；（3）调试运行代码，得到仿真图形；（4）通过仿真图形观察二者的关系，并得出结论。M文件的代码如下：%x为任意的序列,N变换的点数%X为x的傅里叶变换的结果,XK为x离散傅里叶变换的结果function[X,Xk]=wjl_ft_dft(x,N>w=2*pi*(0:127>/128。X=x*exp(-j*[1:length(x>]'*w>%x的傅立叶变换FTb5E2RGbCAPsubplot(211>plot(w
2、,abs(X>>title('x的傅里叶变换FT'>。xlabel('omega/pi'>。ylabel('
3、X
4、'>p1EanqFDPwk=0:N-1。XK=x*exp(-j*[1:length(x>]'*(2*pi*k>/N>%x的离散傅立叶变换DFTDXDiTa9E3dsubplot(212>stem(k,abs(XK>>title('x的离散傅里叶变换DFT'>。xlabel('k'>。ylabel('
5、XK
6、'>end在命令窗口中输入如下指令：>>x=[0804451233]。>>N=32。>>wjl_ft_dft(x,N>得到如下运行结果：
7、2.DFT的快速算法FFT的仿真任务要求：10/10（1）在熟悉FFT算法原理的基础上，对所给的时域信号进行分析；（2）对信号做时域分析时要注意混叠现象的发生；（3）对混合信号做快速傅立叶变换FFT运算，其中、、自设，变换点数根据所设的情况自定；（4）对变换结果做分析。M文件代码如下：%x是原始输入信号,xfft是x的FFT结果%f1,f2,f3是信号x的频率含量,fs是采样频率%n是变换的点数function[x,xfft]=wjl_fft(f1,f2,f3,fs,n>t=0:1/fs:1。x=sin(2*pi*f1*t>+sin(2*pi*f2*t>+
8、sin(2*pi*f3*t>。subplot(211>plot(t,x>%信号x的时域信号图形title('x的时域波形'>。xlabel('t'>。ylabel('x'>subplot(212>xfft=fft(x,n>f=fs*(0:n/2-1>/n。plot(f,abs(xfft(1:n/2>>>%信号x的n点FFT图形title('x的FFT波形'>。xlabel('omega/pi'>。ylabel('
9、X
10、'>RTCrpUDGiTend在命令窗口输入如下代码：>>f1=10。f2=15。f3=25。fs=100。n=1024。>>wjl_f
11、ft(f1,f2,f3,fs,n>得到如下图形：3.利用MATLAB信号处理工具箱设计FIR滤波器要求如下：(1>熟悉FIR滤波器的工作原理；(2>设计一个N阶高通数字滤波器，阶数和截止频率自定；10/10(3>设计一个N阶带通数字滤波器，阶数和截止频率自定；(4)设计一个N阶低通数字滤波器，阶数和截止频率自定；(5)有完整的MATLAB代码。1、FIR高通数字滤波器设计M文件代码如下：%hisimpulseresponseofFIRfilter,sisoriginalsignal,sfisfilteredsignal5PCzVD7HxA%f1andf2a
12、refrequencyofsignals,fsissamplefrequencyjLBHrnAILgfunction[h,s,sf]=wjl_FIR_hpf(f1,f2,fs>m=(0.3*f1>/(fs/2>。%definetransitionbandwidthxHAQX74J0XM=round(8/m>。%definelengthofhammingwindowLDAYtRyKfEN=M-1。%defineorderofFIRfilterZzz6ZB2Ltkb=fir1(N,0.9*f2/(fs/2>,'high'>。[h,f]=freqz(b,1,51
13、2>。%frequencyresponseoffilterdvzfvkwMI1plot(f*fs/(2*pi>,20*log10(abs(h>>>%drawmagnitude_frequencyoffilterrqyn14ZNXIxlabel('frequency/Hz','fontsize',8,'fontweight','bold'>EmxvxOtOcoset(gca,'fontsize',8,'fontweight','bold'>figuresubplot(211>t=0:1/fs:1。s=sin(2*pi*f1*t>+sin(2*pi*f2*t>
14、。plot(t,s>。%draworiginalsignalsSi
