离散时间系统时域.docVIP

实验一、离散时间系统时域、变换域分析 班级：通信C113 姓名：冉启锋 学号：115237 一、实验目的 1．熟悉信号处理软件MATLAB的使用。 2．离散信号的基本运算实现。 3．了解基本序列及复杂序列的产生方法。 4．运用卷积方法观察系统的时域特性。 5. 熟悉信号的采样与恢复。 6．掌握线性时不变系统的变换域表示方法。 7. 掌握序列傅里叶变换性质。 二、实验内容 （一）验证性实验内容 1．熟悉扩展函数 2．运行例题程序 （1）序列的基本运算实验结果 序列加的结果 y1 = 1 5 9 13 17 10 n1 = 1 2 3 4 5 6 序列乘的结果 y2 = 0 6 20 42 72 0 n2 = 1 2 3 4 5 6 序列移位结果 y3 = 1 3 5 7 9 n3 = 4 5 6 7 8 序列翻褶结果 y4 = 10 8 6 4 2 n4 = -6 -5 -4 -3 -2 序列和结果 y5 = 25 序列积结果 y6 = 945 序列能量结果 ex = 165 （2）产生序列并画图 （3）产生复信号 -10=n=10，并画出复序列的实部、虚部、幅值和相位离散图。 （4）线性时不变系统的变换域表示 （5）验证采样定理 1、请写出时的程序 Ts=0.001; n=-25:1:25; x=exp(-1000*abs(n*Ts)); K=500;k=0:1:K;w=pi*k/K; X=x*exp(-j*n*w); X=real(X); w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; subplot(2,2,3); stem(n*Ts*1000,x); xlabel(time(millisecond)); gtext(Ts=0.2ms); ylabel(x1(n)); title(discrete signal); subplot(2,2,4); plot(w/pi,X); xlabel(frequency(radian)); ylabel(x1(jw));title(DTFT); 2、将实验结果画出 （二）提高综合运用能力部分内容 1．编写利用扩展函数产生下列序列并画图的程序 (a) -5=n=5 n=-5:5; x=2*impseq(-2,-5,5)-impseq(4,-5,5); stem(n,x) xlabel(n); ylabel(x(n)); (b) 和 0=n=50 w(n)为白噪声，函数w=randn(size(n)) n=0:50; x=cos(0.04*n*pi); y=x+0.2*randn(size(n)); subplot(2,1,1); stem(n,x); xlabel(n); ylabel(x(n)); subplot(2,1,2); stem(n,y); xlabel(n); ylabel(y(n)); 2.编写程序求系统输出 y(n)的零状态响应 设线性移不变系统的抽样响应为 ， 输入序列为 提示: 输出为输入和抽样响应的卷积 n=0:30; h=(0.9.^n); x=stepseq(0,0,20)-stepseq(10,0,20); [y,ny]=conv_m(x,n,h,n); h=[h zeros(1,(length(ny)-length(h)))]; x=[x zeros(1,(length(ny)-length(x)))]; subplot(3,1,1);stem(ny,h); subplot(3,1,2);stem(ny,x); subplot(3,1,3);stem(ny,y); 3. 编写程序实现序列傅里叶变换频移性质、时域卷积后的频域特性（时域卷积定理的程序设计实现）。 时域卷积定理： w=-pi:2*pi/255:pi; x1=[2 4 6 8 10 12 14 16 18]; x2=[1 -2 3 -2 1]; y=conv(x1,x2); h1=freqz(x1,1,w); h2=freqz(x2,1,w); hp=h1.*h2; h3=freqz(y,1,w); subplot(2,2,1); plot(w/pi,abs(hp));