神经网络探第二章.ppt

第三章 Z变换 . 在MATLAB语言中有专门对信号进行正反Z变换的函数ztrans( ) 和itrans( )。其调用格式分别如下： F=ztrans( f ) ?????对f(n)进行Z变换，其结果为F(z) F=ztrans(f,v)? ?对f(n)进行Z变换，其结果为F(v) F=ztrans(f,u,v) ???对f(u)进行Z变换,其结果为F(v) f=itrans ( F )? ??对F(z)进行Z反变换,其结果为f(n) f=itrans(F,u)? ?对F(z)进行Z反变换，其结果为f(u) f=itrans(F,v,u ) ?对F(v)进行Z反变换,其结果为f(u) 注意： 在调用函数ztran( )及iztran( )之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量 Z变换 例①.求数列 fn=e-n的Z变换及其逆变换。命令如下： syms n z fn=exp(-n); Fz=ztrans(fn,n,z) %求fn的Z变换 f=iztrans(Fz,z,n) %求Fz的逆Z变换 Z变换 例③ ．已知一离散信号的Z变换式为 , 求出它所对应的离散信号f(k). MATLAB程序如下： syms k z Fz=2*z/(2*z-1);?????? %定义Z变换表达式 fk=iztrans(Fz,k)??????? %求反Z变换 运行结果如下：fk = (1/2)^k 例6.7.1(例6.5.1) clear all; fp=100;fs=300;Fs=1000;rp=3;rs=20; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; Fs=Fs/Fs; % let Fs=1 % Firstly to finish frequency prewarping ; wap=tan(wp/2);was=tan(ws/2); % [n,wn]=buttord(wap,was,rp,rs,s) % Note: s! [z,p,k]=buttap(n); % [bp,ap]=zp2tf(z,p,k) % [bs,as]=lp2lp(bp,ap,wap) % % Note: s=(2/Ts)(z-1)/(z+1);Ts=1,that is 2Fs=1,Fs=0.5; [bz,az]=bilinear(bs,as,Fs/2) % [h,w]=freqz(bz,az,256,Fs*1000); plot(w,abs(h));grid on; 例6.7.1(例6.5.1) clear all; wp=.2*pi;ws=.6*pi;Fs=1000; rp=3;rs=20;% % Firstly to finish frequency prewarping; wap=2*Fs*tan(wp/2);was=2*Fs*tan(ws/2); [n,wn]=buttord(wap,was,rp,rs,s);% Note: s! [z,p,k]=buttap(n); [bp,ap]=zp2tf(z,p,k); [bs,as]=lp2lp(bp,ap,wap) w1=[0:499]*2*pi; h1=freqs(bs,as,w1); [bz,az]=bilinear(bs,as,Fs) % Note: z=(2/ts)(z-1)/(z+1); [h2,w2]=freqz(bz,az,500,Fs); plot(w1/2/pi,abs(h1),w2,abs(h2),k);grid on; 例6.7.1(例6.5.1) clear all; wp=.2*pi;ws=.6*pi;Fs=1000; rp=3;rs=20; [n,wn]=buttord(wp/pi,ws/pi,rp,rs); [bz,az]=butter(n,wp/pi) [bz1,az1]=butter(n,wn) [h,w]=freqz(bz,az,128,Fs); [h1,w1]=freqz(bz1,az1,128,Fs); plot(w,abs(h),w1,abs(h1),g.);grid on; 例4.7.2 设x(n)由三个实正弦组成，频率分别是8HZ， 8.22HZ 和9HZ， 抽样频率是40HZ ，时域取128点。 clear; % 生成滤波器系数h和混有高斯白噪的正弦信号x h=fir1(10,0.3,hanning(11)); N=500;p=0.05;f=1/16; u=randn(1,N)*sqrt(p);% s=sin(2*pi*f*[0:N-1]); x=u(1:N)+s; % 将x分为长度为L的小段 L=2