南通大学DSP课程设计实验及习题标准答案完整版精Word格式.docx
《南通大学DSP课程设计实验及习题标准答案完整版精Word格式.docx》由会员分享，可在线阅读，更多相关《南通大学DSP课程设计实验及习题标准答案完整版精Word格式.docx（28页珍藏版）》请在冰豆网上搜索。
stem（n,x1,'
subplot（4,1,3;
stem（n,x2,'
subplot（4,1,4;
stem（n,x3,'
实验二
n=-100:
100;
x1=（n>
=-2&
（n<
2;
m=（-j*0.65.*pi*n;
w=-2*pi:
0.01:
2*pi;
X1=x1*exp（-j*n'
*w;
x=eps.^m;
x2=x1.*x;
X2=x2*exp（-j*n'
x3=（n>
=0&
4;
X3=x3*exp（-j*n'
figure（1
plot（w,abs（X1;
figure（2
plot（w,abs（X2;
figure（3
plot（w,abs（X3;
实验三
N1=64;
N2=64；
N3=100;
n1=0:
N1-1;
n2=0:
N2-1;
n3=0:
N3-1;
x1=（n1>
（n1<
8;
X1=fft（x1,N1;
a11=cos（0.5*pi.*n1;
a21=cos（0.25*pi.*n1;
a31=cos（0.125*pi.*n1;
a12=cos（0.5*pi.*n2;
a22=cos（0.25*pi.*n2;
a32=cos（0.125*pi.*n2;
x21=a11+a21+a31;
x22=a12+a22+a32;
X21=fft（x21,N1;
X22=fft（x22,N2;
a4=cos（0.5*pi.*n3;
a5=cos（0.25*pi.*n3;
a6=cos（0.125*pi.*n3;
x3=a4+a5+a6;
X3=fft（x3,N3
N4=250;
n4=0:
N4-1;
u1=（n4>
u2=（n4>
=N4;
a7=cos（0.5*pi.*n4;
a8=cos（0.25*pi.*n4;
a9=cos（0.125*pi.*n4;
x4=a7+a8+a9;
y1=u1-u2;
y2=hamming（N4;
x41=x4.*y1;
x42=x4.*y2'
;
X41=fft（x41,N4;
X42=fft（x42,N4;
stem（n1,abs（X1,'
subplot（2,1,1;
stem（n1,abs（X21,'
subplot（2,1,2;
stem（n2,abs（X22,'
stem（n3,abs（X3,'
figure（4
stem（n4,abs（X41,'
stem（n4,abs（X42,'
stem（n,abs（X41,'
stem（n,abs（X42,'
实验四
dt=0.0005;
t=-0.005:
dt:
0.005;
xa=exp（-1000*abs（t;
n=-500:
500;
w=-3*pi:
0.001:
3*pi;
Ts1=1/20000;
Ts2=1/5000;
Ts3=1/1000;
x1=exp（-1000*abs（n*Ts1;
x2=exp（-1000*abs（n*Ts2;
x3=exp（-1000*abs（n*Ts3;
figure（1;
plot（t,xa;
subplot（3,1,1;
stem（w,abs（X1
subplot（3,1,2;
stem（w,abs（X2
subplot（3,1,3;
stem（w,abs（X3
实验五
A=[1,-3.1836,4.6223,-3.7795,1.8136,-0.48,0.0544];
B=[0.0004,0.0022,0.0056,0.0075,0.0056,0.0022,0.0004];
[H,w]=freqz（B,A;
plot（w,abs（H,'
'
figure（2;
h=impz（B,A;
stem（h,'
figure（3;
zplane（B,A;
z=roots（B
p=roots（A
x=[-4,-2,0,-4,-6,-4,-2,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,4,0,0,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];
figure（4;
y=filter（B,A,x;
stem（y
stem（x
实验六
Fs=5000;
[b,a]=butter（3,pi/4*Fs,'
s'
[c,d]=impinvar（b,a,Fs;
[H,w]=freqz（c,d;
subplot（2,1,1
plot（w/pi,abs（H;
ylabel（'
³
å
»
÷
Ï
ì
Ó
¦
²
±
ä
·
¨
[m,n]=bilinear（b,a,Fs;
[M,w]=freqz（m,n;
subplot（2,1,2
plot（w/pi,abs（M;
Ë
«
ß
Ð
Ô
grid;
wp=0.4*pi;
ws=0.5*pi;
ap=3;
as=20;
[N,wn]=buttord（wp/pi,ws/pi,ap,as;
[b,a]=butter（N,wn,'
low'
freqz（b,a;
n=0:
47;
x=1+cos（pi*n/4+cos（2*pi*n/3;
y=filter（b,a,x;
plot（n,x;
plot（n,y;
ws=0.4*pi;
wp=0.5*pi;
high'
figure（5
wp=[0.2*pi0.48*pi];
ws=[0.1*pi0.55*pi];
[N,wn]=buttord（wp./pi,ws./pi,ap,as;
bandpass'
figure（6
figure（7
实验七
N=25;
wn=0.6;
b=fir1（N,wn,hamming（N+1;
[H,w]=freqz（b,1;
N;
stem（b,'
grid
plot（w/pi,20*log10（abs（H;
b2=fir1（N,wn,boxcar（N+1;
[H2,w]=freqz（b2,1;
stem（b2,'
plot（w/pi,20*log10（abs（H2;
b3=fir1（N,wn,blackman（N+1;
[H3,w]=freqz（b3,1;
stem（b3,'
plot（w/pi,20*log10（abs（H3;
N2=30;
wn2=[0.30.6];
b4=fir1（N2,wn2,hamming（N2+1;
[H4,w]=freqz（b4,1;
stem（b4,'
figure（8
plot（w/pi,20*log10（abs（H4;
习题一
n=-12:
12;
x=（n>
x0=（2.^n
x1=x.*x0;
x01=（n>
=-4&
x02=（n>
-4&
x2=x02-x01;
x3=cos（0.3*pi*n;
习题二
N1=16;
N2=25;
N3=32;
N4=50;
=32;
x0=x01-x02;
x1=sin（0.25*n*pi.*x0;
X1=dtft（x1,n,w;
plot（w,angle（X1;
x2=sin（0.25.*n1*pi;
X2=fft（x2,N1;
subplot
升级会员 