hw_pam.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

%抽样定理
%PAM编译码实验
%产生频率为1Kz的正弦波
clc;
fs=8000;ts=1/fs;
dt=ts/8;
t=0:dt:0.005;
xt=2*sin(2*1000*pi*t)
subplot(311);plot(t,xt);
title('频率为1Kz的正弦波');xlabel('t');ylabel('xt');
%抽样信号，频率为8Kz
t1=0:ts:0.005;
st=2*sin(2*1000*pi*t1)
subplot(312);stairs(t1,st)
title('频率为8Kz的抽样信号');xlabel('抽样点t1');ylabel('抽样值st')
%通过巴特沃兹滤波器后，抽样信号还原
t2=dt:dt:8*length(st)*dt
x1=zeros(1,8*length(st));
for i=1:length(st)%插值，使得波形平缓抽样信号
    x1(8*i-7)=st(i)
end
[b a]=butter(9,2*1000*dt)
y=filter(b,a,x1);
subplot(313);plot(t2,y)
axis([0 0.005 -0.2 0.2]);grid
title('抽样信号还原波形');xlabel('t2');ylabel('y');

%figure(2)
%频率混叠
%产生频率为5Kz的正弦波
clc;
figure(2);
fs=8000;ts=1/fs;
dt=ts/8;
t=0:dt:0.003;
xt=2*sin(2*5000*pi*t)
subplot(311);plot(t,xt);
title('频率为5Kz的正弦波');xlabel('t');ylabel('xt');
%抽样信号，频率为8Kz
t1=0:ts:0.003;
st=2*sin(2*5000*pi*t1)
subplot(312);stem(t1,st)
title('频率为8Kz的抽样信号');xlabel('抽样点t1');ylabel('抽样值st')
%通过巴特沃兹滤波器后，抽样信号还原
t2=dt:dt:8*length(st)*dt
x1=zeros(1,8*length(st));
for i=1:length(st)%插值，使得波形平缓抽样信号
    x1(8*i-7)=st(i)
end
[b a]=butter(9,2*1000*dt)
y=filter(b,a,x1);
subplot(313);plot(t2,y)
axis([0 0.003 -0.02 0.02]);grid
title('抽样信号还原波形');xlabel('t2');ylabel('y');