《数字信号处理课程设计实验报告》.docxVIP

指导老师：***院 系：***班 级：***姓 名：***学 号：***88数字信号处理实验报告实验一 信号、系统及系统响应实验目的：熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对时域采样定理的理解。熟悉离散信号和系统的时域特性。熟悉线性卷积的计算编程方法：利用卷积的方法，观察，分析系统响应的时域特性。掌握序列傅立叶变换的计算机实验方法，利用序列的傅里叶变换对离散信号、系统及系统响应进行频域分析。实验代码：信号产生子程序：代码：clc;close all;clear all;A=444.128;a=50*pi*sqrt(2.0);W0=50*pi*sqrt(2.0);n=-50:1:50;T=1/1000;Xa=A*(exp(a*n*T)).*(sin(W0*n*T));%产生理想采样信号序列subplot(3,1,1);stem(n,Xa,r);title(理想采样信号序列);xlabel(n);ylabel(Xa);Xb=(n==0);%产生单位脉冲序列subplot(3,1,2);stem(n,Xb,r);title(单位脉冲序列);xlabel(n);ylabel(Xb);Xc=(n=0)-(n=9);%产生矩形序列subplot(3,1,3);stem(n,Xc,r); title(矩形序列);xlabel(n);ylabel(Xc);运行结果：分析理想采样信号序列的特性：代码：clc;close all;clear all;A=444.128; a=50*pi*sqrt(2.0);W0=50*pi*sqrt(2.0);n=-50:50;T1=1/1000;Xa=A*(exp(a*n*T1)).*(sin(W0*n*T1)); % T1=1/1000时得Xa序列stem(n,Xa,r);title(Xa序列);xlabel(n);ylabel(Xa);k=-25:25;X1=Xa*(exp(-j*pi/12.5)).^(n*k); %X1=FT[Xa]pause %暂停，观察波形subplot(3,2,1); stem(k,abs(X1),r);% Xa的幅度谱及相位谱title(Xa的幅度谱);ylabel(幅度);subplot(3,2,2);stem(k,angle(X1),r); title(Xa的相位谱);ylabel(相位);T2=1/300;%T2=1/300时得Xb序列Xb=A*(exp(a*n*T2)).*(sin(W0*n*T2));X2=Xb*(exp(-j*pi/12.5)).^(n*k);%X2=FT[Xb]subplot(3,2,3);stem(k,abs(X2),r); % Xb的幅度谱及相位谱title(Xb的幅度谱);ylabel(幅度);subplot(3,2,4);stem(k,angle(X2),r); title(Xb的相位谱);ylabel(相位);T3=1/200;% T3=1/200时得Xc序列Xc=A*(exp(a*n*T3)).*(sin(W0*n*T3));X3=Xc*(exp(-j*pi/12.5)).^(n*k);%X3=FT[Xc] subplot(3,2,5); stem(k,abs(X3),r);% Xc的幅度谱及相位谱title(Xc的幅度谱);ylabel(幅度);subplot(3,2,6);stem(k,angle(X3),r); title(Xc的相位谱);ylabel(相位);运行结果：图表 1 图表 2离散信号、系统和系统响应的分析：代码：1) clc;close all;clear all;n=-50:1:50;Xb=(n==0); hb=(n==0)+2.5*(n==1)+2.5*(n==2)+(n==3); Yb=conv(Xb,hb);subplot(3,1,1); stem(n,Xb,r);%产生单位脉冲序列title(单位脉冲序列);xlabel(n);ylabel(Xb);k=-25:25;X=Xb*(exp(-j*pi/12.5)).^(n*k);%X=FT[Xb]subplot(3,1,2); stem(k,abs(X),r);%产生Xb的幅度谱及相位谱title(Xb的幅度谱); ylabel(幅度);subplot(3,1,3); stem(k,angle(X),r); title(Xb的相位谱); ylabel(相位);pause%暂停观察波形subplot(3,1,1);stem(n,hb,r); title(hb序列);xlabel(n);ylabel(hb);H=hb*(exp(-j*pi/12.5)).^(n*k);%H=FT[hb]subplot(3,1,2);stem(k,abs(H),r)