实验二20110599.docVIP

数字信号处理实验 姓名：刘书梅 学号班级：11信工2班 实验名称：Z变换及快速傅里叶变换 实验目的 熟悉对离散系统的频率响应分析方法 加深对零点和极点分布的概念理解 熟悉快速傅利叶变换的应用及画图 利用求系统的零极点，单位脉冲响应及其波形，系统的频率响应及其波形 二、 习题1：代码及图形如下 clf;clc;clear; fs=1000; %采样频率和数据点数 N=1024; n=0:N-1; t=n/fs; %时间序列 f1=50;f2=120; %信号的频率成分 x1=sin(2*pi*f1*t)+2*sin(2*pi*f2*t); %无噪声信号 figure(1) subplot(2,1,1) plot(t(1:100),x1(1:100)) title(无噪声); x2=x1+2*randn(1,length(t)); %含有噪声的信号 subplot(2,1,2) plot(t(1:100),x2(1:100)) title(有噪声); %------------------------------------------------------- y1=fft(x1,N); %对无噪声信号进行快速Fourier变换 mag1=abs(y1); %求取振幅 f=n*fs/N; %真实频率序列 figure(2) subplot(2,1,1) plot(f(1:N/2),mag1(1:N/2)*2/N); %绘出Nyquist频率之前的振幅 xlabel(频率/Hz); ylabel(振幅); title(N=1024（无噪声）); grid on; %------------------------------------------------------- y2=fft(x2,N); %对含有噪声信号进行Fourier变换 mag2=abs(y2); %得到振幅值 subplot(2,1,2) plot(f(1:N/2),mag2(1:N/2)*2/N); %绘出Nyquist频率之前的振幅 xlabel(频率/Hz); ylabel(振幅); title(N=1024（有噪声）); grid on; 说明：图1和图2中均用subplot把图形分隔成两部分，各自绘图。图1中第一部分为无噪声信号的图形，信号的频率成分有低频和高频，分别为f1=50，f2=120，由图可见，无噪声信号的图形曲线很平滑，第二部分则在第一部分的基础上用randn命令加入了随机信号，使信号产生了噪声，图中信号的曲线有较尖锐的变化。图2中为对1图中的信号进行快速傅里叶变换并分别绘出振幅随频率变化的曲线。 三、 习题2：代码与图形如下 clear;clc;clf; fs=100;N=128; %采样频率和数据个数 n=0:N-1;t=n/fs; %时间序列 x=sin(2*pi*40*t)+sin(2*pi*15*t); %时间域信号 subplot(2,2,1),plot(t,x);xlabel(时间/s); ylabel(x);title(原始信号); grid on; y=fft(x,N); %Fourier变换 mag=abs(y); %得到振幅谱 f=n*fs/N; %频率序列 subplot(2,2,2), plot(f(1:N/2),mag(1:N/2)*2/N); %绘制Nyquist频率前的振幅 xlabel(频率/Hz);ylabel(振幅); title(原始信号的快速Fourier变换);grid on; xifft=ifft(y); %进行Fourier逆变换 realx=real(xifft); %求取Fourier逆变换的实部 ti=[0:length(xifft)-1]/fs; %Fourier逆变换的时间序列 subplot(2,2,3),plot(ti,realx); xlabel(时间/s);ylabel(x); title(运用Fourier逆变换得到的信号);grid on; yif=fft(xifft,N); %将Fourier逆变换得到的时间域信号进行Fourier变换 mag=abs(yif); f=[0:length(y)-1]*fs/length(y); %频率序列 subplot(2,2,4),plot(f(1:N/2),mag(1:N/2)*2/N);%绘制Nyquist频率前的振幅 xlabel(频率/Hz);ylabel(振幅); title(运用IFFT得到信号的快速Fourier变换);grid on; %误差 e1=x-xifft; e2=y-