Matlab中的信号处理函数讲课.ppt

* 与本章内容有关的MATLAB文件主要是fft, ifft和 czt.m。顾名思义，fft实现快速傅立叶变换，ifft实现快速傅立叶反变换，czt.m 用来实现线性调频Z变换。 fft的调用格式是： X=fft(x), 或 X=fft(x, N)。 以高频pi为频谱中心 fftshift（X）：移位使以零频为频谱中心 与第八、九章有关的 MATLAB 文件 * X=fft(x) fftshift(X) 信号x * fftfilt.m 用叠接相加法实现长序列卷积。格式是: y=fftfilt(h,x) 或 y=fftfilt(h, x，N) 记 的长度为 ， 的长度为 。 若采用第一个调用方式，程序自动地确定对 分段的长度 及做FFT的长度 ， 显然， 是最接近 的2的整次幂。分的段数为 。 采用第二个调用方式，使用者可自己指定做FFT的长度。建议使用第一个调用方式。 * 例1. 设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(4,1,1); plot(x(1:N/4)); axis tight f1=0:2*pi/N:2*pi-pi/N; f=-pi:2*pi/N:pi-pi/N; X=fft(x); y=ifft(X); subplot(4,1,2); plot(f1, abs(X));axis tight subplot(4,1,3); plot(f, fftshift(abs(X)));axis tight subplot(4,1,4); plot(real(y(1:N/4))); axis tight * X=fft(x) fftshift(X) 信号x x’=ifft(X) * clear; % 用叠接相加法，计算滤波器系数h和输入信号x的卷积 % 其中h为10阶hanning窗，x是带有高斯白噪的正弦信号 h=fir1(10,0.3,hanning(11)); % h: is the impulse response of a N=500;p=0.05;f=1/16; % low-pass filter. u=randn(1,N)*sqrt(p); % u:white noise s=sin(2*pi*f*[0:N-1]); % s:sine signal x=u(1:N)+s; % x: a long sequence; y=fftfilt(h,x); % Overlap-add method of FIR filtering using FFT % y=filter(h,1,x); % y=x*h subplot(211); plot(x); subplot(212); plot(y); figure; subplot(211); plot(abs(fft(x))); subplot(212); plot(abs(fft(y))); 例2 令x(n)为一正弦加白噪声信号，长度为500，h(n)是用fir1.m文件设计的一个低通FIR滤波器，长度为11.试用fftfilt实现长序列的卷积 * 信号x 滤波输出信号y 信号x频谱 信号y频谱 * clear; % 按照overlap-add算法编程实现长信号滤波输出 % 生成滤波器系数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