南通大学DSP课程设计实验及习题标准答案完整版精Word格式.docx
- 文档编号:14532253
- 上传时间:2022-10-23
- 格式:DOCX
- 页数:28
- 大小:523.45KB
南通大学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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 南通 大学 DSP 课程设计 实验 习题 标准答案 完整版