Maltab作业31.doc

一、实验任务及要求： 实验任务 设某模拟信号由3种频率成分组成，f1=19KHz,f2=20KHz,f3=21KHz.x(t)=2cos(2pif1t)+cos(2pif2t)+3cos(2pif3t) 对其进行频谱分析，要求分辨出信号所包含的3种频率成分。 实验要求： 1）分析并选择合理的采样频率 2）分析并选择合适的fft点数，要求点数为 2的整数次幂。 3）分析并仿真点数的选择对频率分辨率的 影响。 4）将实验过程和结果写成实验报告，打印 上交。 实验分析： 选择合理的采样频率并选择合适的fft点数,要求点数为2的整数次幂编写程序时选择了fs=40000*m;N=2^(7+m); 实验中取四组数据观察图像m=1,2,3,4; 实验结果 实验结论 采样频率越大，选择的点数越多，分辨率越高。为了实验简洁，只取了几组数据，但在粗略的分析下我们依然可以得到上述正确结论 程序代码： clc;clear all;close all; for m=1:4 fs=40000*m; N=2^(7+m); %采样频率和数据点数 n=0:N-1;t=n/fs; %时间序列 x=2*cos(2*pi*19000*t)+cos(2*pi*20000*t)+3*cos(2*pi*21000*t); %信号 y1=fft(x,N); %对信号进行快速Fourier变换 y2=fftshift(y1); mag1=abs(y1); %求得Fourier变换后的振幅 mag2=abs(y2); f1=n*fs/N; %频率序列 f2=n*fs/N-fs/2;%频率序列 figure(m); subplot(2,1,1),plot(f1,mag1,r); %绘出随频率变化的振幅 xlabel(频率/Hz); ylabel(振幅);title(usual FFT,color,r);grid on; subplot(2,1,2),plot(f2,mag2,k); %绘出随频率变化的振幅 xlabel(频率/Hz); ylabel(振幅);title([FFT after fftshift,color,k, N=,num2str(N), fs=,num2str(fs)]);grid on; end MATLAB信号处理上机作业 2012年11月6日 学号：XXXX 姓名：XXXX 班级：XXXX - 1 -