画频谱.docVIP

插值---先插零值，再低通滤波（dsp.p.129） 时域上对某种波形的采样频率是一样的，所以频域一样，数字域变密。 b=zeros(1,1024); b(25:50)=1;%---矩形 figure(1); plot(abs(fftshift(fft(b))));%---中间一个sinc函数 c=zeros(1,4096) c(1:4:4096)=b;%----四倍插值 figure(2); plot(abs(fftshift(fft(c))));%--4000点fft----四个sinc函数 从时频域来看，采样率变为原来4倍，w=W*T，新图中1500—2500中间部分相当于原图，新旧无差别：频域无区别，数字频域由于采样率变化而变化。。。 上图经过低通滤波器，将1500—2500外的频谱滤掉，实现实际插值过程。 figure(3); plot(abs(fftshift(fft(c,512))));%--512点fft 512点和4096点都是画出整个-pi到+pi之间的频谱。 figure(4) spectrum(c); % --对数形式画频谱函数，只画一半。。。 过滤波器： fd=1;%信号符号速率 T=1/fd;%信号符号间隔 fs=8;%滤波器采样速率=信号符号速率*插值率 R=0.5;%滚降系数 ---大了边沿冲击小 del=3;%边带长度，即延时时间 h_lpf3_temp=rcosflt(1,fd,fs,normal/Fs,R,del);%输入冲击信号，产生冲击响应 h_lpf3=h_lpf3_temp;%转置 stem(h_lpf3); figure; plot(abs(fftshift(fft(h_lpf3)))); N_lpf=conv(h_lpf3, c); figure %卷积 plot(abs(fftshift(fft(N_lpf)))); plot(N_lpf); 抽取---先将几个值变零，然后去掉零值(dsp.p.129) 数字频域和频域都会变化 时域上对信号有效采样率变小，频域重复周期变小。 d=zeros(1,1024); d(25:50)=1;%---矩形 figure(1); plot(abs(fftshift(fft(d))));%---中间一个sinc函数 e=zeros(1,1024); e(25:4:50)=1;%---矩形 figure plot(abs(fftshift(fft(e))));% f=e(1:4:1024); figure plot(abs(fftshift(fft(f))));% 数字频域也是模拟频域的变化过程，采样率不变，但对信号的有效采样率降低。 频域重复频率增加。