pam.m

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

clear all;clc;
t0=1;ts=0.001;ts1=0.01;
ts2=0.0005;
fs=1/ts;fs1=1/ts1;fs2=1/ts2;
df=0.5;
t=[-t0/2:ts:t0/2];

t2=[-t0/2:ts2:t0/2];
x2=sin(200*t2);
m2=x2./(200*t2);
w2=t0/(2*ts2)+1;
m2(w2)=1;
%m2=m2.*m2;
[M2,mn2,df2]=fftseq(x2,ts2,df);
M2=M2/fs2;
N2=[M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2];
f2=[-7*df2*length(mn2):df2:6*df2*length(mn2)-df2]-fs2/2;
subplot(211);stem(t2,m2);
xlabel('time');
subplot(212);plot(f2,abs(fftshift(N2)));