23栗国齐 - 数字信号处理实验报告.docVIP

 PAGE \* MERGEFORMAT 26 指导老师：杨永双,何燕,张彦,李娜 院 系：计算机与通信工程学院 班 级：通信工程09-1班 姓 名：栗国齐 学 号：540907040123 数字信号处理实验报告 实验一 信号、系统及系统响应 实验目的： 熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对时域采样定理的理解。 熟悉离散信号和系统的时域特性。 熟悉线性卷积的计算编程方法：利用卷积的方法，观察，分析系统响应的时域特性。 掌握序列傅立叶变换的计算机实验方法，利用序列的傅里叶变换对离散信号、系统及系统响应进行频域分析。 实验代码： 信号产生子程序： 代码： 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(相位); 运行结果： 图表  SEQ 图表 \* ARABIC 1 图表  SEQ 图表 \* ARABIC 2 离散信号、系统和系统响应的分析： 代码： 1) clc; close all; clear all; n=-50:1:50; Xb=(n==0); hb=(n=