基于MATLAB的信号采样与离散系统的时域分析.docVIP

基于matlab的信号采样与离散系统的时域分析 原理 所谓“取样”就是利用取样脉冲序列s(t)从连续信号f(t)中“抽取”一系列离散样本值的过程。这样得到的离散信号称为取样信号fs(t) 。 画fs(t)的频谱时，设定ωS ≥2ωm ,这时其频谱不发生混叠，因此能设法(如利用低通滤波器)从Fs (j()中取出F(j()，即从fs(t)中恢复原信号f(t); 否则将发生混叠。 设有限时宽余弦信号f(t)=cos(2Πt/3)（0≤t≤40）,用它近似理想余弦信号，用MATLAB编程求该信号和其采样信号的频谱，对比观察过采样和欠采样状态。采样周期通过键盘输入。 程序： % 时域采样定理 display(奈奎斯特周期0.25秒，Ts0.25，过采样，Ts0.2.5，欠采样); display(Please input the value of sample period); Ts = input(Ts = ); %绘制有限长余弦信号y=cos(2/3*pi*t) t = 0:0.01:40; y = cos(2/3*pi*t); subplot(221); plot(t,y); axis([0 6 -1.1 1.1]); xlabel(t 单位：s,Fontsize,8); title(f(t)); line([0 6],[0 0],color,[0 0 0]); %数值求解余弦信号的频谱 N = 300; W = 2*pi*5; k = -N:N; w = k*W/N; Y = 0.01*y*exp(-j*t*w); %求f(t)的傅里叶变换F1(ω) Y = abs(Y); subplot(222); plot(w/pi,Y) axis([-2,2,0,pi*7+0.2]); title(F(j\omega)); xlabel(\omega 单位:pi); %采样后的余弦信号 subplot(223); plot(t,y,b:); %绘制包络 hold on t2=0:Ts:40; y2=cos(2/3*pi*t2); stem(t2,y2); axis([0 6 -1.1 1.1]); xlabel(t 单位：s,Fontsize,8); title(fs(t)); hold off %采样后余弦信号的频谱 Y2 = Ts*y2*exp(-j*t2*w); Y2 = abs(Y2); subplot(224); plot(w/pi,Y,b) %蓝色绘制原信号频谱 xlabel(\omega 单位:pi); title(Fs(j\omega)); hold on plot(w/pi,Y2,r); %红色绘制采样信号频谱 axis([-2,2,0,pi*10]); hold off %end 2. 离散时间系统(1)线性时不变 (LTI) 离散时间系统用常系数线性差分方程进行描述： 其中，f[k]和y[k]分别表示系统的输入和输出，N=max(n,m)是差分方程的阶数。 系统的零状态响应就是在系统初始状态为零条件下微分方程的解。在零初始状态下，MATLAB控制系统工具箱提供了一个filter函数，可以计算差分方程描述的系统的响应，其调用形式为： y=filter(b,a,f) 其中，、分别是系统差分方程左、右端的系数向量，f表示输入向量，y表示输出向量。注意，输出序列的长度与输入序列的长度相同。 (2)离散系统的冲激响应、阶跃响应分别是输入信号为和所对应的零状态响应。MATLAB控制系统工具箱专门提供了两个函数求解离散系统的冲激响应和阶跃响应。 冲激响应：h = impz(b, a, K)，其中的h表示系统的单位序列响应，、分别是系统差分方程左、右端的系数向量，K表示输出序列的时间范围。阶跃响应：g = stepz(b, a, N)，其中的g表示系统的单位阶跃序列响应，b和a的含义与上相同，N表示输出序列的长度。 (3)卷积是信号与系统中一个最基本、也是最重要的概念之一。在时域中，对于LTI 离散时间系统，其零状态响应等于输入信号与系统冲激响应的卷积；而利用卷积定理，这种关系又对应频域中的乘积。 MATLAB信号处理工具箱提供了一个计算两个离散序列卷积和的函数conv。设向量a、b代表待卷积的两个序列，则c = conv(a, b)就是a与b卷积后得到的新序列。 我们知道两个序列卷积以后，一般而言所得新序列的时间