实验一 Matlab软件的使用Word文档格式.docx
《实验一 Matlab软件的使用Word文档格式.docx》由会员分享，可在线阅读，更多相关《实验一 Matlab软件的使用Word文档格式.docx（15页珍藏版）》请在冰豆网上搜索。
x=[zeros（1,5）,ones（1,8）,zeros（1,23）];
gridon；
程序3：
n=-15:
20;
x=[zeros（1,15）,1,2,3,4,5,zeros（1,16）];
gridon;
程序4：
N=10;
M=10;
L=N+M-1;
x=[1,2,3,4,5,zeros（1,5）];
h=[1,2,1,2,zeros（1,6）];
y=conv（x,h）;
n=0:
L-1;
stem（n,y,'
*k'
gridon;
实验结果：
1、
2、
3、
4、
实验二DTFT和Z变换的Matlab实现
掌握DTFT和Z变换的Matlab实现的方法；
进一步了解离散时域系统的时频特性。
DTFT和z变换的计算和基本性质；
离散时域系统的仿真；
利用函数实现离散时域系统的频率响应，分析零、极点的分布，进一步了解离散时域系统的时域特性。
3、实验原理
利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统频率响应等的图像并于笔算结果进行比较，找出异同。
编译合适程序能计算取值范围不同的离散卷积。
x=zeros（1,length（n））;
h=zeros（1,length（n））;
x（[find（（n>
=0）&
（n<
=4））]）=1;
h（[find（（n>
=8））]）=0.5;
figure
（1）
subplot（3,1,1）;
stem（n,x）;
axis（[0,30,0,2]）;
title（'
输入序列'
xlabel（'
n'
ylabel（'
x（n）'
subplot（3,1,2）;
stem（n,h）;
冲激响应序列'
h（n）'
subplot（3,1,3）;
length（y）-1;
stem（n,y）;
输出响应'
y（n）'
b=[0.0181,0.0543,
0.0543,0.0181];
a=[1.000,-1.76,
1.1829,-0.2781];
w=pi*freqspace（500）;
H=freqz（b,a,w）;
MH=abs（H）;
AH=angle（H）;
subplot（2,1,1）;
plot（w/pi,MH）;
grid;
axis（[0,1,0,1]）;
w（pi）'
|H|'
幅度、相位响应'
subplot（2,1,2）;
plot（w/pi,AH）;
angle（H）'
实验三DFS和DFT的Matlab实现
加深DFS和DFT的算法原理和基本性质的理解；
熟悉利用FFT对典型数字信号进行频谱分析。
使用Matlab实现程序周期序列的DFS；
利用FFT和IFFT编程实现两个有限长序列的循环卷积；
对有限长序列使用离散Fouier变换（DFT）可以很好的反映序列的频谱特性，而且易于用快速算法在计算机上实现。
functionshiyan22（）
F=50;
N=64;
T=0.000625;
n=1:
N;
x=cos（2*pi*F*n*T）;
plot（n,x）;
holdon
X=fft（x）;
plot（n,X）;
grid
FFT|X|'
f（pi）'
%输入x（n）和冲激响应h（n）
=4））]）
=1;
=8））]）
=0.5;
stem（x）;
stem（h）;
H=fft（h）;
Y=X.*H;
y=ifft（Y）;
实验四IIR数字滤波器的设计和实现
掌握冲击响应不变法和双线形变换法设计IIR数字滤波器的原理和方法；
通过观察对实际信号的滤波，获得对数字滤波的感性认识。
巴特沃思滤波器设计实现模拟低通滤波器；
利用冲击响应不变法和双线形变换法对模拟滤波器进行A/D转换；
用数字滤波器的单位脉冲响应序列
模仿模拟滤波器的冲激响应
让
正好等于
的采样值，即
。
设计一个巴特沃思数字低通滤波器，设计指标如下：
通带内
幅度衰减不大于1dB；
阻带
幅度衰减不小于15dB。
编制实现该数字滤波器程序并且实现数字滤波。
functionshiyan3（）
fp=200*pi;
fs=300*pi;
wp=0.2*pi;
ws=0.3*pi;
Rp=1;
As=15;
T=1;
Fs=1000;
[N,fc]=buttord（fp,fs,Rp,As,'
s'
[b,a]=butter（N,fc,'
w=[0:
1000*2*pi];
[hf,w]=freqs（b,a,1000）;
[d,c]=impinvar（b,a,Fs）;
wd=[0:
512]*pi/512;
hw1=freqz（d,c,wd）;
[f,e]=bilinear（b,a,Fs）;
hw2=freqz（f,e,wd）;
OmegaP=（2/T）*tan（wp/2）;
OmegaS=（2/T）*tan（ws/2）;
ep=sqrt（10^（Rp/10）-1）;
Ripple=sqrt（1/（1+ep.^2））;
Attn=1/10^（As/20）;
plot（wd/pi,abs（hw1）/abs（hw1
（1）））;
幅度响应（冲击响应不变法）'
H'
axis（[0,1,0,1.1]）;
set（gca,'
XTickmode'
'
manual'
XTick'
[0,0.2,0.35,1.1]）;
YTickmode'
YTick'
[0,Attn,Ripple,1]）;
plot（wd/pi,20*log10（（abs（hw2）+eps）/abs（hw2
（1））））;
幅度响应（双线性变换法（dB））'
axis（[0,1,-40,5]）;
[-50,-15,-1,0]）;
plot（wd/pi,-angle（hw2））;
相位响应'
piunit'
%axis（[0,1,0,1.1]）;
[-1,0,1]）;
实验五FIR数字滤波器的设计和实现
掌握用窗函数法设计FIR数字滤波器的原理和方法；
熟悉线性相位FIR滤波器的幅频特性和相频特性；
了解不同窗函数对滤波器性能的影响。
用Hanning窗设计线性相位带通滤波器；
分别改用矩形窗和Blackman窗设计线性带通滤波器；
确定数字滤波器的性能要求，根据性能要求，合理选择单位脉冲响应
的奇偶对称性，从而确定理想频率响应
的幅频特性和相频特性；
求理想单位脉冲响应
，选择适当的窗函数
，根据
求所需设计的FIR滤波器单位脉冲响应；
设计一个FIR数字滤波器，设计指标如下：
幅度衰减不小于15dB；
编制计算设计的数字滤波器幅度特性和相位特性的程序。
tr_width=ws-wp;
M=ceil（6.6*pi/tr_width）+1;
M-1;
wc=（ws+wp）/2;
alpha=（M-1）/2;
m=n-alpha+eps;
hd=sin（wc*m）./（pi*m）;
w_ham=（hamming（M））'
;
[mag,db,pha,w]=freqz_m（h,[1]）;
delta_w=2*pi/1000;
Rp=-（min（db（1:
wp/delta_w+1）））;
As=-round（max（db（ws/delta_w+1:
501）））;
subplot（2,2,1）;
stem（n,hd）;
理想冲激响应'
axis（[0,M-1,-0.1,0.3]）;
hd（n）'
subplot（2,2,2）;
实际冲激响应'
subplot（2,2,
升级会员 