测控13《工测试信号分析与处理》实验指导书新.docVIP

《测试信号分析与处理》实验指导书 实验一 离散系统的时域分析离散系统的时域分析Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于 信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对离散系统的时域分析/link?url=G9lkzhNeDxCOjJ6in61zJFm6fVPRpMiJUXA0WjiYbQrxRDFymy227UtT_4QXQtI6qalbqMPHuWQDKcVHb7xFx9lUnQqDNXokBW19plJTpDe） 相关MATLAB函数 sin?and?cos?函数 x = 3*sin(pi/8*n) x = -cos(5*pi/7*n-pi/3) 2.求数字卷积 y = conv(x,h)%x为数字信号，h为滤波器的脉冲响应 或 y = conv(h,x) 3.滤波：filter()函数　　y=filter(b,a,x) 由向量b和a组成的系统对输入x进行滤波　　. 求解单位响应：impz()函数　　（1）h=impz(b,a)　　 计算单位响应的序列值，取样点个数由MATLAB自动选取　　（2）h=impz(b,a,n)　　计算指定范围内（0: n-1）的单位响应的序列值　　（3）impz(b,a)　　绘制单位响应的时域波形　　. 求解频率响应：freqz()函数　　（1）[h,w]=freqz(b,a,n)　　可得到n点频率响应，这n个点均匀地分布在上半单位圆（即 ），并将这n点频率记录在w中，相应的频率响应记录在h中。n最好能取2的幂次方，如果缺省，则n=512。　　（2）[h,w]=freqz(b,a,n,whole)　　均匀选取n个点计算频率响应。　　（3）[h,w]=freqz(b,a,n,Fs)　　Fs为采样频率（以Hz为单位），在0Fs/2频率范围内选取n个频率点，计算相应的频率响应。　　（4）[h,w]=freqz(b,a,n,whole,Fs)　　 在0Fs之间均匀选取n个点计算频率响应。　　（5）freqz(b,a)　　可以直接得到系统的幅频和相频特性曲线。其中幅频特性以分贝的形式给出，频率特性曲线的横轴采用的是归一化频率，即Fs/2=1。　　. 系统函数的零极点图的绘制：zplane(b,a) Conv.m% 计算两个序列的线性卷积； clear;%清屏 N=5; M=6; L=N+M-1; x=[1,2,3,4,5];%数字信号x的采样值 h=[6,2,3,6,4,2];%数字滤波器的脉冲响应h y=conv(x,h);%求x和h的数字卷积,输出为y nx=0:N-1;%设置输入信号x的波形图横坐标为0~N-1,间隔为1 nh=0:M-1;%设置h的波形图横坐标为0~M-1,间隔为1 ny=0:L-1;%设置输出信号y的波形图横坐标为0~L-1,间隔为1 subplot(131);%绘图,图示格式为1行3列共3张图,这里显示第一幅图,其位置在最上方, stem(nx,x,.k);xlabel(n);ylabel(x(n));grid on;%绘x的图.横坐标n,纵坐标x(n),线条颜色黑色’k’指代’black’,grid on指显示网格线 subplot(132);%绘图,图示格式为1行3列共3张图,这里显示第二幅图,其位置在中间 stem(nh,h,.k);xlabel(n);ylabel(h(n));grid on;%绘h的图.横坐标n,纵坐标h(n),线条颜色黑色’k’指代’black’,grid on指显示网格线 subplot(133);%绘图,图示格式为1行3列共3张图,这里显示第三幅图,其位置在最下方 stem(ny,y,.k);xlabel(n);ylabel(y(n));grid on;%绘y的图.横坐标n,纵坐标y(n),线条颜色黑色’k’指代’black’,grid on指显示网格线 （2）filter.m;%求利用数字滤波器对信号进行滤波处理后的输出； clear; x=ones(100);%设置一维数组x,共100个1 t=1:100;%从1~100,间隔为1 b=[.001836,.007344,.011016,.007374,.001836];%设置数字滤波器的bk系数 a=[1,-3.0544,3.8291,-2.2925,.55075];%设置数字滤波器的ak系数 y=filter(b,a,x);%