chenpc%2f文件下载%2f信号处理%2f离散时间信号与系统(第一章).pptVIP

解本题的MATLAB程序如下： x=[3 4 5];nx=[-1:1]; h=[1 2 3];nh=[-1:1]; [y,ny]=myconv(x,nx,h,nh) subplot(221) stem(nx,x) xlabel(n) ylabel(x(n)) title(第一个序列) subplot(223) stem(nh,h) xlabel(n) ylabel(h(n)) title(第二个序列) subplot(222) stem(ny,y) xlabel(n) ylabel(y(n)) title(卷积和序列) 上机作业：给出两个序列： 求线性卷积： 。 答案： 三、因果性：某时刻的输出只取决于该时刻及其以前的输入。 四、稳定性：有界输入产生有界输出。 例题：一个线性时不变系统的单位脉冲响应为： 求当输入信号为矩形序列 时的 输出 。 解： 作业：P.42 2 8 例题： 求连续信号及其对应的离散信号的频谱，观察采用不同的抽样频率对函数 的频谱的影响。 （1）以 对 抽样得到 ，求 ； （2）以 对 抽样得到 ，求 。 解：连续信号的频谱： 离散信号的频谱： 用MATLAB语言编程绘出连续信号和离散信号的频谱并验证抽样定理的正确性。 方法一： %本程序绘出连续信号及其对应的离散信号的频谱， %以说明Shannon抽样定理的正确性。 % chouyang.m fs=5000;T=1/fs;n=-25:25; xa=exp(-1000*abs(n*T)); %抽样信号 W=-7000:0.01:7000; Xa=2000./(1000^2+W.^2); subplot(321) plot(1000*n*T,xa) %绘连续信号 xlabel(t in msec) %时间以毫秒为单位 ylabel(x_a(t)) subplot(322) plot(W,Xa) xlabel(\Omega) ylabel(X_a(\Omega)) title(连续信号频谱) axis([-1e4 1e4 0 2.2e-3]) x1=xa; subplot(323) stem(n*T*1000,x1) xlabel(t in msec) ylabel(x_1(n)) w=W*T; X1=(1-exp(-2000*T))./(1-2*cos(w)*exp(-1000*T)+exp(-2000*T)); subplot(324) plot(W,T*X1) xlabel(\Omega) ylabel(X_1(\Omega)) title(离散信号频谱(f_s=5000Hz)) axis([-1e4 1e4 0 2.2e-3]) fs=1000;T=1/fs;n=-5:5; x2=exp(-1000*abs(n*T)); w=W*T; X2=(1-exp(-2000*T))./(1-2*cos(w)*exp(-1000*T)+exp(-2000*T)); subplot(325) stem(n*T*1000,x2) xlabel(t in msec) ylabel(x_2(n)) subplot(326) plot(W,T*X2) xlabel(\Omega) ylabel(X_2(\Omega)) title(离散信号频谱(f_s=1000Hz)) 程序运行结果为： 方法二： % jiaocai2d1.m % 本程序完成门爱东著《数字信号处理》P.13页 % 例2-1 fs=5000;T=1/fs;n=-25:25; xa=exp(-1000*abs(n*T)); %抽样信号 W=-7000:0.1:7000; Xa=2000./(1000^2+W.^2); %连续信号的频谱 subplot(321) plot(1000*n*T,xa) %绘连续信号 xlabel(‘t in msec’) %时间以毫秒为单位 ylabel(x_a(t)) subplot(322) plot(W,Xa) %绘连