数字信号处理实验四-李马元.docVIP

电子信息工程学系实验报告 课程名称：数字信号处理 实验项目名称：用DFT（FFT）对连续信号进行频谱分析 实验时间：2012.5.5 实 班级：电信092 姓名: 李马元 学号：910706203 一、实 验 目 的 ①掌握用DFT（FFT）对模拟信号进行谱分析的方法，理解可能出现的分析误差及其原因， 以便在实际中正确应用FFT。 熟悉应用FFT实现两个序列的线性卷积的方法。 计算机混叠现象截断效应栅栏效应 2、①编写程序如下，结果如图一： clc;clear;fs=20003; %**************************** tp1=1.5/1000;t10=0:1/fs:tp1;x1=ones(1,length(t10)); t11=(tp1+1/fs):1/fs:2*tp1;x2=zeros(1,length(t11)); t1=[t10,t11];x1t=[x1,x2]; subplot(3,1,1);plot(t1,x1t);axis([0,3/1000,0,1.5]); xlabel(t1-李马元);ylabel(x1(t));title(limayuan); %****************************************** f=50;tp2=1/10;t2=0:1/fs:tp2;x2t=sin(2*pi*f*t2+pi/8); subplot(3,1,2);plot(t2,x2t); xlabel(t2-李马元);ylabel(x2(t));title(limayuan); %****************************************** tp3=1;t3=0:1/fs:tp3;x3t=cos(8*pi*t3)+cos(16*pi*t3)+cos(20*pi*t3); 图一 个实验信号图 subplot(3,1,3);plot(t3,x3t); xlabel(t3-李马元);ylabel(x3(t));title(limayuan); 3、(1)编写程序如下，结果如图二所示： clc;clear;fs=64000;tp1=3/1000;t1=0:1/fs:tp1;t11=1.5/1000;n1=ones(1,t11*fs);n11=zeros(1,fs*(tp1-t11)); x1=[n1,n11];fs1=4000;x1t=x1(1:fs/fs1:tp1*fs);subplot(3,1,1);stem(0:1/fs1:(tp1*fs1-1)/fs1,x1t,.); xlabel(t1-李马元);ylabel(x1(t));title(limayuan); %*********************** tp1=3/1000;t1=0:1/fs:tp1;t11=1.5/1000; n1=ones(1,t11*fs);n11=zeros(1,fs*(tp1-t11)); x1=[n1,n11];fs1=8000;x1t=x1(1:fs/fs1:tp1*fs); subplot(3,1,2);stem(0:1/fs1:(tp1*fs1-1)/fs1,x1t,.); xlabel(t1-李马元);ylabel(x1(t));title(limayuan); %************************** tp1=3/1000;t1=0:1/fs:tp1;t11=1.5/1000; n1=ones(1,t11*fs);n11=zeros(1,fs*(tp1-t11)); x1=[n1,n11];fs1=16000;x1t=x1(1:fs/fs1:tp1*fs); subplot(3,1,3);stem(0:1/fs1:(tp1*fs1-1)/fs1,x1t,.); xlabel(t1-李马元);ylabel(x1(t));title(limayuan); 图2 不同采样频率下的频谱 (2)编写程序如下，结果如图三： clear;clc;fs=64000;f=203;T=1/f;fs1=4*f;tp2