课件：神经网络第二章.ppt

* * 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=20;M=length(h); y=zeros(1,N+M-1); tempy=zeros(1,M+L-1); tempX=zeros(1,L); for k=0:N/L-1 tempx(1:L)=x(k*L+1:(k+1)*L); tempy=conv(tempx,h); y=y+[zeros(1,k*L),tempy,zeros(1,N-(k+1)*L)]; end subplot(211);plot(x) subplot(212);plot(y(1:N)) * * hilbert.m 文件用来计算信号Hilbert变换。调用 的格式是: y=hilbert(x),y的实部就 是 ，虚部是的Hilbert变换 。 所以，y 实际上是 x 的解析信号。 * czt.m 调用格式是： X＝czt(x, M, W, A) 。x是待变换的时域信号，其长度设为N，M是变换的长度，W确定变换的步长，A确定变换的起点。若M=N, A=1, 则CZT变成DFT。 A=exp(j*2*pi*f0/fs); W=exp(-j*2*pi*DELf/fs); * 例4.7.2 设x(n)是两个正弦信号及白噪声的叠加，进行频谱分析。 clear all; % 产生两个正弦加白噪声； N=256; f1=.1;f2=.2;fs=1; a1=5;a2=3; w=2*pi/fs; x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N); % 应用FFT 求频谱； subplot(3,1,1); plot(x(1:N/4)); f=-0.5:1/N:0.5-1/N; X=fft(x); y=ifft(X); subplot(3,1,2); plot(f,fftshift(abs(X))); subplot(3,1,3); plot(real(y(1:N/4))); * * 例4.7.2 设x(n)由三个实正弦组成，频率分别是8HZ， 8.22HZ 和9HZ， 抽样频率是40HZ ，时域取128点。 用CZT计算的DFT 用FFT计算的DFT 7~(7+M×0.05) HZ 的CZT变换 * 程序 clear all； % 构造三个不同频率的正弦信号的叠加作为试验信号 N=128; f1=8;f2=8.22;f3=9;fs=40; stepf=fs/N; n=0:N-1; t=2*pi*n/fs; n1=0:stepf:fs/2-stepf; x=sin(f1*t)+sin(f2*t)+sin(f3*t); M=N; W=exp(-j*2*pi/M); % A=1时的czt变换 A=1; Y1=czt(x,M,W,A); subplot(311) plot(n1,abs(Y1(1:N/2)));grid on; * % DTFT Y2=abs(fft(x)); subplot(312) plot(n1,abs(Y2(1:N/2)));grid on; % 详细构造A后的czt M=60; f0=7.2; DELf=0.05; A=exp(j*2*pi*f0/fs); W=exp(-j*2*pi*DELf/fs); Y3=czt(x,M,W,A); n2=f0:DELf:f0+(M-1)*DELf; subplot(313);plot(n2,abs(Y3));grid on; * 1．filtfilt.m 本文件实现零相位滤波。其调用格式是：y=filtfilt(B, A, x) 。式中B是 的分子多项式，A是分母多项式，x是待滤波信号，y是滤波后的信号。 clear; N=32; n=-N/2:N+N/2; w=0.1*pi; x=cos(w*n)+cos(2*w*n); subplot(311);stem(n,x,.);grid on; xlabel(n); b=[0.06745 0.1349 0.06745]; a=[1 -1.143 0.4128]; y=filtfilt(b,a,x); % 用给定系统(b,a)对信号 x 作零相位滤波； y1=filter(b,a,x); % 用给定系统(b,a