采样率转换的MATLAB程序实现

随着数字信号处理技术的发展，信号处理系统中的处理、编码、传输和存储等工作量越来越大。为了节省计算工作量和存储空间，一个信号处理系统通常需要不同的抽样率并可在相互之间转换。抽取与内插是多抽样信号处理的基本操作，可以实现采样率的整数倍转换。两者结合还可以实现采样率的非整数倍转换，即先进行L倍的内插然后再进行M倍抽取，就可实现采样率的L/M倍转换。

%% 抽取，d倍
close all;
clear all;
clc;
N=100; 
n=0:0.5:N-1; 
x=cos（0.125*pi*n）;
subplot（211） 
stem（nx）;
axis（[-133-1.11.1]）;
xlabel（‘t‘）;title（‘x（t）‘）; 
d=5;
% y=zeros（1length（x））; 
% for i=1:d:length（x）; 
% y（i）=x（i）; %采样值 
% end 
y=x（1:d:end）
subplot（212） 
stem（y）;
axis（[-133-1.11.1]）;
xlabel（‘t‘）;title（‘y（t）‘）;
z1=fft（x）;%内插前频谱图
figure;
plot（abs（fftshift（z1）））;
z2=fft（y）;%内插后频谱图
figure;
 w=0:length（z2）-1;
plot（w/piabs（fftshift（z2）））;
% plot（abs（fftshift（z2）））

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         525  2020-06-07 11:44  chouqu.m
     文件         647  2020-06-07 11:43  neicha.m
     文件         517  2020-06-07 11:45  neicha_chouqu.m