第一章 数字信号处理.ppt

图1.2.2 例题1.2.2图 其MATLAB程序演示如下： 卷积和演示 h=[1 2 1]; x=[h,2*ones(1,17)]; %x(n)只取前20个点 y=conv(x,h); stem(y);xlabel(n);ylabel(y(n)); 卷积结果如图1.2.4所示。 图1.2.4 离散序列卷积求和结果 a=[1 -0.7 0.1];b=[7 -2]; N=50; n=0:N-1; h=impz(b,a,n); %脉冲响应 x=ones(1,N); %构造输入u(n) y=filter(b,a,x); %零状态响应 subplot(3,1,1);stem(n,h);xlabel(n);ylabel(h(n)); subplot(3,1,2);stem(n,x);xlabel(n);ylabel(x(n)); subplot(3,1,3);stem(n,y);xlabel(n);ylabel(y_{zs}(n)); 其结果如下图所示： 图1.2.5 例题1.2.5脉冲响应和零状态响应 MATLAB应用举例 %采样定理演示程序 Dt=0.00005;t=-0.005:Dt:0.005;xa=exp(-1000*abs(t)); %产生模拟信号xa(t) Wmax=2*pi*2000;K=500;k=0:1:K;w=k*Wmax/K; Xa=xa*exp(-j*t*w)*Dt; Xa=real(Xa); %连续时间傅立叶变换 W=[-fliplr(w),w(2:501)]; %将频率限制在-Wmax到Wmax之间 Xa=[fliplr(Xa),Xa(2:501)]; %将Xa频率限制在-Wmax到Wmax之间 subplot(2,2,1);plot(t*1000,xa); xlabel(时间(ms)); ylabel(x_{a}(t));title(模拟信号); subplot(2,2,2);plot(W/(2*pi*1000),Xa*1000);xlabel(频率(KHz)); ylabel(|X_{a}(j\Omega)|);title(模拟信号的CFT); Ts=0.0002;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); %采样信号的CFT w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; w1=w/pi; w1=[w1,w1+5,w1+2*5]; T=0:3002; D=(0:3)*1001; Y=pulstran(T,D,X); %将X进行周期延拓 subplot(2,2,3);stem(n*Ts*1000,x); xlabel(时间(ms)); ylabel(x_{1}(n));title(采样信号); subplot(2,2,4);plot(w1,Y); xlabel(频率(KHz)); ylabel(|X_{a1}(j \Omega)|);title(采样信号的CFT(f_{s}=5KHz)) Ts=0.0005;n=-10:1:10; 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)]; h=w(end)-w(1); w=[w,w+h,w+2*h]; T=0:3002; D=(0:3)*1001; Y=pulstran(T,D,X); figure(2);plot(w/pi,Y); xlabel(频率(KHz)); ylabel(|X_{a2}(j \Omega)|);title(采样信号的CFT(f_{s}=2KHz)); x=[x zeros(1,3)]; z=x; for n0=1:3 X=[x(end-n0+1:end),x(1:end-n0)]; % x为输入 z=z+X; end y=z/4; 通过选择“深圳能源”2009年8月20日当日各时刻股票成交价格作为输入信号，观测其平均走势。 Matlab程序如下： x=[12.30 12.54 12.52 12.54 12.47 12.41 12.37 12.40 12.46 12.4