基于MATLAB的图形仿真与矫正参考.doc

武汉理工大学 学生实验报告书 实验课程名称：测试信号分析与处理 开 课 学 院 ：机 电 工 程 学 院 指导教师姓名：陈 祯 学 生 姓 名 ：乾 引 学生专业班级：测 控1303 学号班级序号：18 《测试信号分析与处理》实验 实验一 离散系统的时域分析离散系统的时域分析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,) 在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];h=[6,2,3,6,4,2];y=conv(x,h); nx=0:N-1;nh=0:M-1;ny=0:L-1; subplot(131);stem(nx,x,.k);xlabel(n);ylabel(x(n));grid on; subplot(132);stem(nh,h,.k);xlabel(n);ylabel(h(n)); subplot(133); stem(ny,y,.k);xlabel(n);ylabel(y(n));grid on; （2）filter.m;%求利用数字滤波器对信号进行滤波处理后的输出； clear;x=ones(100);t=1:100;b=[.,.,.,.,.]; a=[1,-3.0544,3.8291,-2.2925,.55075];y=filter(b,a,x);stem(t,y,.);ylabel(y)); xlabel(n); （3）impz .m% 计算滤波器的脉冲响应 clear; b=[.,.,.,.,.];a=[1,-3.0544,3.8291,-2.2925,.55075]; [h,t]=impz(b,a,40); stem(t,h,.);grid on;ylabel(h(n))xlabel(n) （4）filter.m% 计算滤波器的阶跃响应 x=ones(100);t=1:100;b=[.,.,.,.,.]; a=[1,-3.0544,3.8291,-2.2925,.55075];y=filter(b,a,x);plot(t,x,g.,t,y,k-);grid on; ylabel(x(n) and y(n)); xlabel(n); （5）dtft.m%求取非周期数字信号频谱并绘图；