数字信号处理上机答案(含程序及图片)第三版高西全著Word格式.doc
《数字信号处理上机答案(含程序及图片)第三版高西全著Word格式.doc》由会员分享，可在线阅读，更多相关《数字信号处理上机答案(含程序及图片)第三版高西全著Word格式.doc（14页珍藏版）》请在冰点文库上搜索。
y1（n）'
;
stem（n,y1n,'
（b）系统对R8（n）的响应y1（n）'
ylabel（yn）;
y2n=filter（B,A,x2n,xi）;
length（y2n）-1;
subplot（2,2,4）;
y2（n）'
stem（n,y2n,'
（c）系统对u（n）的响应y2（n）'
内容二
x1n=[11111111];
h1n=[ones（1,10）zeros（1,10）];
h2n=[12.52.51zeros（1,10）];
y21n=conv（h1n,x1n）;
y22n=conv（h2n,x1n）;
M1=length（y21n）-1;
M2=length（y22n）-1;
n1=0:
1:
M1;
n2=0:
M2;
n11=0:
length（h1n）-1;
n22=0:
length（h2n）-1;
tstem（n11,h1n）;
（d）系统单位脉冲响应h1（n）'
ylabel（h1（n））;
stem（n1,y21n,'
fill'
（e）h1（n）与R8（n）的卷积y21（n）'
ylabel（y21（n））;
subplot（2,2,3）;
tstem（n22,h2n）;
（f）系统单位脉冲响应h2（n）'
ylabel（h2（n））;
stem（n1,y22n,'
（g）h2（n）与R8（n）的卷积y22（n）'
ylabel（y22（n））;
内容三
谐振器对u（n）的响应
xn=[1,zeros（1,250）];
B=[1/100.49,-1/100.49];
A=[1,-1.8237,0.9801];
yn=filter（B,A,xn,xi）
length（yn）-1;
subplot（1,1,1）;
谐振器对正弦信号的响应
xsin=sin（0.014*n）+sin（0.4*n）;
yn=filter（B,A,xsin,xi）
数字信号处理实验三
实验
（1）
x1n=[ones（1,4）];
X1k8=fft（x1n,8）;
X1k16=fft（x1n,16）;
N=8;
f=2/N*（0:
N-1）;
figure
（1）;
subplot（1,2,1）;
stem（f,abs（X1k8）,'
（la）8µ
ã
DFT[x_1（n）]'
\omega/\pi¡
¯
ylabel（¡
®
|（e^j^\omega）|'
N=16;
subplot（1,2,2）;
stem（f,abs（X1k16）,'
（la）16µ
\omega/\pi'
ylabel（'
实验（1-2，1-3）
M=8;
xa=1:
（M/2）;
xb=（M/2）:
-1:
1;
x2n=[xa,xb];
x3n=[xb,xa];
X2k8=fft（x2n,8）;
X2k16=fft（x2n,16）;
X3k8=fft（x3n,8）;
X3k16=fft（x3n,16）;
figure
（2）;
stem（f,abs（X2k8）,'
（2a）8µ
DFT[x_2（n）]'
stem（f,abs（X3k8）,'
（3a）8µ
DFT[x_3（n）]'
stem（f,abs（X2k16）,'
（2a）16µ
'
stem（f,abs（X3k16）,'
（3a）16µ
实验（2-1，2-2）
N-1;
x4n=cos（pi*n/4）;
x5n=cos（pi*n/4）+cos（pi*n/8）;
X4k8=fft（x4n,8）;
X4k16=fft（x4n,16）;
X5k8=fft（x5n,8）;
X5k16=fft（x5n,16）;
figure（3）;
stem（f,abs（X4k8）,'
（4a）8µ
DFT[x_4（n）]'
stem（f,abs（X5k8）,'
（5a）8µ
DFT[x_5（n）]'
stem（f,abs（X4k16）,'
（4a）16µ
stem（f,abs（X5k16）,'
（5a）16µ
实验（3）
Fs=64;
T=1/Fs;
nT=n*T;
x8n=cos（8*pi*nT）+cos（16*pi*nT）+cos（20*pi*nT）;
X8k16=fft（x8n,16）;
figure（4）;
stem（f,abs（X8k16）,'
（8a）16µ
DFT[x_8（n）]'
N=32;
X8k32=fft（x8n,32）;
stem（f,abs（X8k32）,'
（8a）32µ
N=64;
X8k64=fft（x8n,64）;
stem（f,abs（X8k64）,'
（8a）64µ
数字信号处理实验四
functionst=mstg
N=800
Fs=10000;
Tp=N*T;
t=0:
T:
（N-1）*T;
k=0:
f=k/Tp;
fc1=Fs/10;
fm1=fc1/10;
fc2=Fs/20;
fm2=fc2/10;
fc3=Fs/40;
fm3=fc3/10;
xt1=cos（2*pi*fm1*t）.*cos（2*pi*fc1*t）;
xt2=cos（2*pi*fm2*t）.*cos（2*pi*fc2*t）;
xt3=cos（2*pi*fm3*t）.*cos（2*pi*fc3*t）;
st=xt1+xt2+xt3;
fxt=fft（st,N）;
subplot（3,1,1）
plot（t,st）;
grid;
t/s'
s（t）'
axis（[0,Tp/8,min（st）,max（st）]）;
（a）s（t）的波形'
subplot（3,1,2）
stem（f,abs（fxt）/max（abs（fxt））,'
（b）s（t）的频谱'
axis（[0,Fs/5,0,1.2]）;
f/Hz'
幅度'
st=mstg;
%低通滤波器设计与实现
fp=280;
fs=450;
wp=2*fp/Fs;
ws=2*fs/Fs;
rp=0.1;
rs=60;
[N,wp]=ellipord（wp,ws,rp,rs）;
[B,A]=ellip（N,rp,rs,wp）;
y1t=filter（B,A,st）;
subplot（3,1,1）;
[H,w]=freqz（B,A,1000）;
m=abs（H）;
plot（w/pi,20*log（m/max（m）））;
grid
on;
低通滤波损耗函数曲线'
w/pi'
axis（[0,1,0,1.2*max（H）]）yt='
y1（t）'
subplot（3,1,2）;
plot（t,y1t）;
低通滤波后的波形'
ylabel（y1（t））;
%带通滤波器设计与实现
fpl=440;
fpu=560;
fsl=275;
fsu=900;
wp=[2*fpl/Fs,2*fpu/Fs];
ws=[2*fsl/Fs,2*fsu/Fs];
y2t=filter（B,A,st）;
带通滤波损耗函数曲线'
y2（t）'
plot（t,y2t）;
带通滤波后的波形'
ylabel（y2（t））;
%高通滤波器设计与实现
fp=890;
fs=600;
[B,A]=ellip（N,rp,rs,wp,'
high'
y3t=filter（B,A,st）;
高通滤波损耗函数曲线'
y3（t）'
plot（t,y3t）;
高通滤波后的波形'
ylabel（y3（t））;
低通滤波器损耗函数及其分离出的调幅信号y1（t）
带通滤波器损耗函数及其分离出的调幅信号y2（t）
高通滤波器损耗函数及其分离出的调幅信号y3（t）
数字信号处理实验五
1、functionxt=xtg（N）
Fs=1000;
fc=Fs/10;
f0=fc/10;
mt=cos（2*pi*f0*t）;
ct=cos（2*pi*fc*t）;
xt=mt.*ct;
nt=2*rand（1,N）-1;
fp=150;
fs=200;
Rp=0.1;
As=60;
fb=[fp,fs];
m=[0,1];
dev=[10^（-As/20）,（10^（Rp/20）-1）/（10^（Rp/20）+1）];
[n,fo,mo,W]=remezord（fb,m,dev,Fs）;
hn=remez（n,fo,mo,W）;
yt=filter（hn,1,10*nt）;
xt=xt+yt;
fst=fft（xt,N）;
plot（t,xt）;
x（t）'
axis（[0,Tp/5,min（xt）,max（xt）]）;
（a）信号加噪声波形'
plot（f,abs（fst）/max（abs（fst）））;
（b）信号加噪声的频谱'
axis（[0,Fs/2,0,1.2]）;
2、xt=xtg;
N=1000;
fp=120;
fs=150;
wc=（fp+fs）/Fs;
B=2*pi*（fs-fp）/Fs;
M=ceil（11*pi/B）;
hn=fir1（M-1,wc,blackman（M））;
Hw=abs（fft（hn,N））;
ywt=fftfilt（hn,xt,N）;
figure;
subplot（2,1,1）;
plot（f,20*log10（Hw）/max（Hw））;
gridon
幅度（dB）'
（a）低通滤波器的幅频特性'
axis（[0,500,-160,5]）;
subplot（2,1,2）;
plot（t,ywt）;
y_1（t）'
（b）滤除噪声后的信号波形'
升级会员 