3．实验过程
3.1整数倍抽取实验
已知一连续时间信号x(t)的表达式为x(t)=Acos(2πf1t)+Bcos(2πf2t) 其中f1=50Hz,f2=100Hz,A=1.5,B=1。
按2整数倍对信号进行抽取，程序如下： Fs=1000; A=1.5; B=1; f1=50; f2=100; t=0:1/Fs:1
x=A*cos(2*pi*f1*t)+cos(2*pi*f2*t); y=decimate(x,2); subplot(212); stem(x(1:25),'.'); xlabel('时间，nT'); ylabel('输入信号'); grid on; subplot(222); stem(y(1:25),'.'); xlabel('时间,nT'); ylabel('输出抽取信号'); grid on;
3.2整数倍插值实验
已知一连续时间信号x(t)的表达式为x(t)=Acos(2πf1t)+Bcos(2πf2t) 其中f1=50Hz,f2=100Hz,A=1.5,B=1。
按4整数倍对信号进行插值，程序如下： Fs=1000;
15
A=1.5; B=1; f1=50; f2=100; t=0:1/Fs:1
x=A*cos(2*pi*f1*t)+cos(2*pi*f2*t); y=interp(x,4); subplot(221); stem(x(1:25),'.'); xlabel('时间，nT'); ylabel('输入信号'); grid on; subplot(222); stem(y(1:100),'.'); xlabel('时间,4nT'); ylabel('输出插值信号'); grid on;
3.3用有理因子I/D的采样率转换进行的实验
已知x(n)=1.5cos(0.3πn),对x(n)按有理因子3/5进行采样速率转换，程序如下： n=0:24;
x=1.5*cos(0.3*pi*n); [y,h]=resample(x,3,5); figure(1); stem(n,x); xlabel('时间，n'); ylabel('输入信号'); ny=0:length(y)-1; figure(2);
16
stem(ny,y); xlabel('时间，n'); ylabel('输出变换信号'); w=(0:511)*2/512;
H=20*log10(abs(fft(h,512))); figure(3); plot(w,H); grid on xlabel('频率'); ylabel('幅值');
17
4．实验结果
4.1信号的整数倍抽取
4.2信号的整数倍插值
18
4.3用有理因子I/D的采样速率转换
19