专业课程设计Sat函数程序波形图.doc

程序：
持续信号Sa(t)：
t=-20:0.5:20；
f=sinc(t/pi)；
plot(t,f)；
xlabel('t');
ylabel('x(t)');
title('持续信号sa(t)=sinc(t/pi)波形')；
grid;

Sa(t)频谱：
t=-20:0.5:20;
f=sinc(t/pi)；
N=1000;
k=-N:N;
w1=10；
w=k*w1/N；
F=f*exp(-j*t'*w)*0.5；
plot(w,F)；
xlabel('x')；
ylabel('x(w)')；
title(' sa(t)=sinc(t/pi)信号频谱图')

Sa(t)临界抽样及重构：
wm=1；
wc=wm；
Ts=pi/wm；
ws=2*pi/Ts；
n=-100:100；
nTs=n*Ts
f=sinc(nTs/pi)；
Dt=0.005;
t=-20:Dt:20；
fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))))；
t1=-20:0.5:20;
f1=sinc(t1/pi);
subplot(211)；
stem(t1,f1);
xlabel('kTs');
ylabel('f(kTs)')；
title('sa(t)=sinc(t/pi)临界采样信号')；
subplot(212)；
plot(t,fa);
xlabel('t');
ylabel('fa(t)');
title('由sa(t)=sinc(t/pi)临界采样信号重构sa(t)');
grid;

Sa(t)过抽样及重构：
wm=1;
wc=1.1*wm；
Ts=0.7*pi/wm;
ws=2*pi/Ts；
n=-100:100;
nTs=n*Ts
f=sinc(nTs/pi);
Dt=0.005;
t=-10:Dt:10； fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));
error=abs(fa-sinc(t/pi));
t1=-10:0.5:10;
f1=sinc(t1/pi)；
m=-20:20；
N=41；
Xw=abs(fft(f,N))；

subplot(2,2,1);
plot(m,Xw)；
axis([-20 20 1.1*min(Xw) 1.1*max(Xw)])；
xlabel('w')；
ylabel('|Xw|')；
title('抽样后频谱波形图');
subplot(2,2,2)；
stem(t1,f1);
xlabel('kTs');
ylabel('f(kTs)')；
title('sa(t)=sinc(t/pi)采样信号')；
subplot(2,2,3);
plot(t,fa)
xlabel('t')；
ylabel('fa(t)');
title('由sa(t)过采样信号重构sa(t)')；
grid；
subplot(2,2,4)；
plot(t,error)；
xlabel('t');
ylabel('error(t)');
title('过采样信号与原信号误差')

Sa(t)欠抽样及重构：
wm=1;
wc=wm；
Ts=1.5*pi/wm;
ws=2*pi/Ts；
n=-100:100;
nTs=n*Ts
f=sinc(nTs/pi);
Dt=0.005;
t=-10:Dt:10；
fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));
error=abs(fa-sinc(t/pi));
t1=-10:0.5:10;
f1=sinc(t1/pi)；
m=-20:20；
N=41；
Xw=abs(fft(f,N))；
subplot(2,2,1);
plot(m,Xw)；
axis([-20 20 min(Xw) max(Xw)])；
xlabel('w')；
ylabel('|Xw|')；
title('抽样后频谱波形图');
subplot(2,2,2)；
stem(t1,f1);
xlabel('kTs');
ylabel('f(kTs)')；
title('sa(t)=sinc(t/pi)采样信号')；
subplot(2,2,3);
plot(t,fa)
xlabel('t')；
ylabel('fa(t)');
title('由sa(t)欠采样信号重构sa(t)')；
grid；
subplot(2,2,4)；
plot(t,error)；
xlabel('t');
ylabel('error(t)');
title('欠采样信号与原信号误差')