matlab在数字信号处理中运用.pptVIP

第七讲 Matlab在数字信号 处理中的运用 ;本章内容：;信号;数字信号;1. 离散时间信号与系统;例7-1 产生一个32点的，右移20各单位的冲击序列。 程序如下： clear all; N=32; k=20; x=zeros(N); x(k)=1; xn=0:N-1; stem(xn,x);; 单位阶跃序列μ(n);例7-3：产生一个频率为150Hz，幅度为0.45，初始相位为35度的正弦波，信号持续时间为5s。;sound(X,Fs); disp(播放结束，下面将音频数据存盘为C:\my50HzSIN.wav); wavwrite(X,Fs,C:\my50HzSIN.wav); clear; [R,Rs]=wavread(C:\my50HzSIN.wav); sound(R,Rs); ;Matlab实现：n=0:N-1; x=a.^n; 例： x(n)= 实现程序： n=[0:10]; x=(0.9).^n; stem(n,x); 随机序列 ;例7-4 产生随机序列;Matlab在数字信号处理中运用;2）波形发生器;例如： x=square(2*pi*f*t+fai,duty) % t为时间取样序列 % f为方波的基频率 % fai为方波的初始相位 % duty为占空比，为0-100之间的数。 例如 duty=20， 即占空比为20% % x 返回的幅度为1的矩形波样值序列; 三角函数sawtooth： ; 其它波形发生函数：;3）序列的操作;倍率： 每一个采样值乘以一个常数a，表示为 a{x(n)} = {ax(n) } 用算术运算符“*”实现。 折叠： x(n)的每个样本都对n=0翻转，得到一个折叠后的序列y(n)。 y(n)= {x(-n)} 由flghr(x)实现。;样本和： 将n1和n2之间的所有样本x(n)加起来： 由sum(x(n1): x(n2))函数实现。 样本积： 将n1和n2之间的所有样本x(n)乘起来： 由prod x((n1: n2))函数实现。;4）线性系统及其Matlab实现; 时域响应工具箱函数;例如：计算 的冲激响应。 sys=tf(1,[1,1,0.5]) Transfer function: 1 ------------- s^2 + s + 0.5 impulse(sys); ; 离散系统模型时域表示 ;impz函数： 直接给出系统的单位冲激相应，调用格式： impz(b,a);例7-5：当系统的输入差分方程为： y(n)-0.8y(n-1)-0.5y(n-2)=0.7x(n)+0.3(n-1), 分别利用filter函数和impz函数求系统的单位冲激相应。 clear all; pulse=[1,zeros(1,63)]; b=[0.7,0.3]; a=[1,-0.8,-0.5]; h1=filter(b,a,pulse); h2=impz(b,a,64); subplot(2,1,1),stem(h1),title(filter fuction); subplot(2,1,2),stem(h2),title(impz fuction);; 传递函数响应;freqz函数： 用于计算并??出离散系统的幅频响应和相频响应。该函数使用基于FFT算法计算系统传递函数响应模型中的系数向量a和b。常用调用格式为： [h,f]=freqz(b,a,n,fs) 其中： b为传递函数H(z)分子多项式系数， a为分母多项式系数， n为指定计算频率点数(由于采用FFT算法，n常取 2的幂次方，以提高计算速度)， fs为离散系统的采样频率， h为对应于频率点序列f的复频率响应。如无输出变量，则自动作出幅频响应和相频响应图， f为记录频率点数。;2. 数字滤波器概述;由于信号通常夹杂噪声及无用信号成分，所以必须将这些干扰成分滤除。滤波器可以对信号进行筛选，只让特定信号通过。 一般而言，噪声信号往往是高频信号。而经典滤波器正是假定有用信号与噪声具有不同的频段，所以利用经典的滤波器可以将噪声滤除。 但是如果信号和噪声频谱相互重叠，那么利用经典的滤波器就不能发挥作用。现代滤波器的作用是从含有噪声的数据记录中估计除信号的某些特征或信号本身，那么估计出来的信号和源信号相比，就具有更高的信噪比。;滤波是信号处理的基础，滤波运算是信号处理中的基本运算。因