新·Matlab 在信号处理中的应用.docVIP

Matlab 在信号处理中的应用 产生锯齿波或三角波信号函数sawtooth 格式一：x=sawtooth（t） 功能：产生周期为2pi，振幅从-1到1的锯齿波。在2pi的整数倍处值为-1到1，这一段波形斜率为1/pi。 格式二：sawtooth（t，width） 功能：产生三角波，width在0到1 之间。 例：产生周期为0.02的三角波。 Fs=10000；t=0:1/Fs：1 X1=sawtooth（2*pi*50*t，0）； X2=sawtooth（2*pi*50*t，1）； Subplot（2,1,1），plot（t，x1），axis（[0,0.2，-1,1]）； Subplot（2,1,2），plot（t，x2），axis（[0,0.2，-1,1]）； 产生方波信号函数square 格式：x=square（t） 功能：产生周期为2pi，振幅从-1到1的方波。 格式：x=square（t，duty） 功能：产生指定周期的方波，duty为正半周期的比例。 例：产生周期为0.02的方波。 Fs=10000；t=0:1/Fs：1 X1=square（2*pi*50*t，20）； X2=square（2*pi*50*t，80）； Subplot（2,1,1），plot（t，x1），axis（[0,0.2，-1.5,1.5]）； Subplot（2,1,2），plot（t，x2），axis（[0,0.2，-1.5,1.5]）； 产生sinc函数波形函数sinc 格式：y=sinc（x） 功能：sinc（x）用于计算sinc函数，即： 1 t=0 Sinc（t）= sin（pi*t）/（pi*t） Sinc函数十分重要，它的傅里叶变换正好是振幅为1的矩阵脉冲 例：产生sinc函数波形。 X=linspace（-4,4） Y=sinc（x） Plot（x，y） 产生非周期方波信号函数rectpuls 格式：y=rectpuls（t） 功能：产生非周期方波信号，方波的宽度为时间轴的一半 格式：y=rectpuls（t，w） 功能：产生指定宽度为w的非周期方波。 例：t=0:0.01:1 Y1=rectpuls（t）；y2=rectpuls（t，0.6） Subplot（2,1,1），plot（t，y1），grid Subplot（2,1,2），plot（t，y2），grid 5.产生非周期三角波信号函数tripuls 格式：y=tripuls（t） 功能：产生非周期三角波信号，三角波的宽度为时间轴的一半。 格式：y=tripuls（t，w，s） 功能：产生周期为w的非周期方波，斜率为s（-1s1）。 例：t=0:0.01:1; y1=tripuls(t)； y2=tripuls(t,0.6)； Subplot(2,1,2),plot(t,y1),grid； subplot(2,1,2),plot(t,y2),grid; 6.产生线性调频扫频信号函数chirp 格式：y=chirp（t，f0，t1，f1） 功能：产生一个线性（频率随时间线性变化）信号，其时间轴设置由数组t定义。时刻0的瞬间频率为f0，时刻t1的瞬间频率为f1.默认情况，f0=0Hz，t1=1，f1=100Hz。 格式：y=chirp（t，f0，t1，f1，‘method’） 功能：指定改变扫频的方法。可用的方法有‘linear’（线性扫频）、‘quadratic’（二次扫频）和‘logarithmic’（对数扫频）；默认时为‘linear’。注意：对于对数扫频，必须有f1f0。 格式：y=chirp（t，f0，t1，f1，‘method’，phi） 功能：指定信号的初始相位为phi（单位为度），默认时phi=0 例：绘制一线性调频信号。 T=0:0.01:2; Y=chirp(t,0,1,150)； Figure(1)； Plot(t,y)?; Axis([0,0.5,0,1])?; 7.产生冲击串信号函数pulstran 格式：y=pulstran（t，d，‘func’） 功能：在指定的时间范围t，对连续函数func，按向量d提供的平移量进行平移后抽样产生的冲击信号y=func（t-d（1））+ func（t-d（2））+...。其中函数func必须是t的函数，且可以用函数句柄形式调用，即y=pulstran（t，d，@func）。 例：产生一不对称的锯齿冲击串信号，要求锯齿宽度为0.1s，波形间隔为1/3s t=0:1/1E3:1;抽样频率1kHz，连续时间为1s d=0?:1/3?:1?;3Hz重复频率 ?y=pul