第7章2012matlab在信号处理中的应用-学生.pptVIP

MATLAB在信号处理中的应用 信号的MATLAB表示及运算 本节首先讨论用MATLAB表示和分析连续信号的问题 严格说来，只有用符号推理的方法才能分析连续系统。用数值方法是不能表示连续信号的，因为它给出的是各个样本点的数据。只有当样本点取得很密时才可看成连续信号。所谓密，是相对于信号变化的快慢而言。以下均假定相对于采样点密度而言，信号变化足够慢。 7.1 常见连续信号的表示 正弦型信号 内部函数cos(w*t ) 和sin(w*t ) 周期方波 square(w*t,占空比) 单位冲激函数 单位阶跃 由于内存有限，MATLAB无法表示无限序列。 在MATLAB中可用一个行向量来表示一个有限长序列。然而，这样一个向量并没有任何有关样本位置k的信息，因此，x(k)的准确表示要求有两个向量；一个对x，另一个对k。 二.常用的典型序列 单位脉冲序列 单位阶跃序列 复指数序列 常用的典型序列及生成 1.单位冲激序列 7.3序列的运算和变换 序列的基本运算为加法、乘法、移位、反折。 1.序列的加法、乘法 要求x1(k)和x2(k)的长度必须相同。如果序列长度不相等，或者即使长度相等的序列而样本位置不同，也不能直接用算符 为此首先必须对x1(k)和x2(k)扩大或延长以使它们具有相同位置的向量k（从而也有相同的长度）。 例：给出两个序列x1和x2。 x1=[1,2,3,4]; 起始位置ns1=-2 x2=[2,2,0,0]， 起始位置ns2＝0； ws=8; %采样角频率 Ts=2*pi/ws; %采样间隔 n=-15:15;nTs =n*Ts; %t的时域采样 Sa_nTs=sinc(nTs/pi); %Sa(nTs)采样信号 stem(n,Sa_nTs); xlabel(nTs); ylabel(Sa(nTs)); title([\omegas=,num2str(ws), 的Sa(t)的采样信号]); grid; ，对该模拟信号采样得到 ， 设采样角频率为ws=8，试确定采样间隔Ts，绘制 已知信号x(t)= Sa(t) nTs =(-15:15)*Ts的Sa信号的抽样图。 随机序列 弦波序列 方波序列 stem函数 迟延的单位冲激序列 单位脉冲序列的MATLAB表示 % unit puls sequence k=-5:5; delta=[zeros(1,5),1,… zeros(1,5)]; stem(k,delta) 方法一： 方法二： k=-5:5; delta=[k==0]; stem(k,delta) 2.单位阶跃序列 迟延的单位阶跃序列 % unit step sequence k=-5:5; uk=[zeros(1,5), ones(1,6)]; stem(k,uk) 单位阶跃序列的MATLAB表示 k=-5:5; f=[k=0]; stem(k,f) 方法一： 方法二： % x(n) = {...,5,4,3,2,1,5,4,3,2,1,...};-10=n=9 n=-10:9; %20个元素 x=[5,4,3,2,1]; xn=repmat(x,1,4); subplot(2,2,4);stem(n,xn); title(例 2.1d 的序列图) ylabel(x(n));axis([-10,9,-1,6]); text(10,-1,n) 若给定超始位置，如何编程 y1 = 1 2 3 4 0 0 y2 = 0 0 2 2 0 0 * * 由内部函数产生 一、周期信号 画出角频率为w的信号 周期方波（square） (1)square(w*t) %产生角频率为w的周期性方波 (2)square(w*t,占空比) %产生角频率为w，占空比＝τ/T*100的周期性方波 A=1;f=4; w=2*pi*f;T=1/f t=0:0.001:10*T; y=A*square(w*t,50); plot(t,y); xlabel(Time(s)); title(square wave) axis([0,2.6,-1.5,1.5]) 例如，产生一个幅度为1，基频为4Hz，占空比为20%的周期方波 矩形波(方波) y = rectpuls(t,width) 三角波 y = tripuls(t, width,skew) 指数信号 y = A*exp(a*t) 取样信号Sa(x) sinc(t) 由内部函数 产生 单位冲激信号 单位阶跃信号 自编程序产生 二.非