实验四(华电版).docxVIP

实验四 离散系统分析一、实验目的深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。二、 实验原理MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。例题：离散系统的单位脉冲响应N=40; a=[1,0.4,-0.12];b=[1,2]; y=impz(b,a,N); stem(y) xlabel(k);title(h[k])N=100;b=[1,2]; a=[1,0.4,-0.12]; x=ones(1,N); y=filter(b,a,x)y = Columns 1 through 12 1.0000 2.6000 2.0800 2.4800 2.2576 2.3946 2.3131 2.3621 2.3327 2.3504 2.3398 2.3461 Columns 13 through 24 2.3423 2.3446 2.3432 2.3441 2.3436 2.3439 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 Columns 25 through 36 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 Columns 37 through 48 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 Columns 49 through 60 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 Columns 61 through 72 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3438 2.3438 Columns 73 through 84 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 Columns 85 through 96 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 Columns 97 through 100 2.3438 2.3438 2.3438 2.3438 N=100; b=[1,2]; a=[1,0.4,-0.12]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi)y = Columns 1 through 12 0.8400 2.7840 1.9872 2.5392 2.2228 2.4156 2.3005 2.3697 2.3282 2.3531 2.3381 2.3471 Columns 13 through 24 2.3417 2.3450 2.3430 2.3442 2.3435 2.3439 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 Columns 25 through 36