Matlab 在信号处理中的应用.doc

Matlab 在信号处理中的应用1. 产生锯齿波或三角波信号函数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：1X1=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）；2. 产生方波信号函数square格式：x=square（t）功能：产生周期为2pi，振幅从-1到1的方波。格式：x=square（t，duty）功能：产生指定周期的方波，duty为正半周期的比例。例：产生周期为0.02的方波。Fs=10000；t=0:1/Fs：1X1=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）；3. 产生sinc函数波形函数sinc格式：y=sinc（x）功能：sinc（x）用于计算sinc函数，即： 1 t=0Sinc（t）= sin（pi*t）/（pi*t） Sinc函数十分重要，它的傅里叶变换正好是振幅为1的矩阵脉冲 例：产生sinc函数波形。X=linspace（-4,4） Y=sinc（x）Plot（x，y）4. 产生非周期方波信号函数rectpuls格式：y=rectpuls（t）功能：产生非周期方波信号，方波的宽度为时间轴的一半格式：y=rectpuls（t，w）功能：产生指定宽度为w的非周期方波。例：t=0:0.01:1Y1=rectpuls（t）；y2=rectpuls（t，0.6）Subplot（2,1,1），plot（t，y1），gridSubplot（2,1,2），plot（t，y2），grid5.产生非周期三角波信号函数tripuls格式：y=tripuls（t）功能：产生非周期三角波信号，三角波的宽度为时间轴的一半。格式：y=tripuls（t，w，s）功能：产生周期为w的非周期方波，斜率为s（-1sf0。格式：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/3st=0:1/1E3:1;抽样频率1kHz，连续时间为1sd=0:1/3:1;3Hz重复频率y=pulstran(t,d,tripuls,0.1,1);调用tripuls函数实现冲击串plot(t，y) ；8.产生Dirichlet信号的函数diric格式：y=diric（x，n）功能：用于产生x的dirichlet函数，即：Dirichlet（x）= （-1）（k*（n-1）x=2*pi*k k=整数 Sin（nx/2）/（n*sin（x/2）其他9.产生高斯正弦脉冲信号函数guaspuls格式一：yi=guaspuls（T,FC,BW,BWR）功能：返回持续时间为T，中心频率为FC(Hz)，宽带为BW的幅度为1的高斯正弦脉冲（RF）信号的抽样。脉冲宽度是指信号幅度下降（相对于信号包络峰值）到BWR(dB)时所对应宽度的100*BW%。注意BW0，BER0，默认时，FC=1000Hz，BW=0.5，BER=-6dB。格式二：TC= guaspuls(cutoff，FC,BE,BER,TPB)功能：返回按参数TEP(dB)计算所对应的截断时间TC。参数TPE(TPE0)是指脉冲拖尾幅度相对包络最大幅度的下降程度，默认时FPE=-60dB。例：产生频率为50kHz，宽度为60%的高斯RF脉冲信号。要求脉冲包络幅度下降到40dB处截断，采样频率为1MHz。tc=gauspuls(cutoff,50E3,0.6,-40);计算截断时间，参数BER用空代替t=-tc:1E-6:tc;按采样频率，对时间t赋值yi=gauspuls(t,50E3,0.6);按要求产生高斯脉冲信号plot(t,yi)10.产生高斯单脉冲信号函数gmonopuls格式：y=gmonopuls（t，fc）功能：返回最大振幅值为1的高斯单脉冲信号，时间数轴由t给定，fc为中心频率。默认情况下，fc=1000Hz。格式：tc=gmonopuls（cutoff，fc）功能：返回信号的最大值和最小值之间的持续时间。例：产生一个2GHz的高斯单脉冲信号，采样频率为100GHz。fc=2e9;fs=100e9; tc=gmonopuls( cutoff,fc);t=-2*tc:1/fs:2*tc;y=gmonopuls(t,fc);plot(t,y)11.电压控制振荡器函数vco格式：y=vco（x，fc，fs）功能：产生一个采样频率为fs的振荡信号。其震荡频率由输入向量或数组x指定。fc为载波或参考频率，如果x=0，则y的采样频率为fs（Hz）、振幅为1、频率为fc（Hz）的余弦信号。X的取值范围为-11，如果x=-1，输出y的频率为0：如果x=0，输出y的频率为fc（Hz）：如果x=1，输出y的频率为2*fc。输出y与x的维数一样。默认情况下，fs=1，fc=fs/4。如果x是矩阵，vco函数按列产生一个振荡信号矩阵，x与它列对应。格式：y=vco（x，fmin，fmax，fs）功能：可调整频率调制的范围，使得x=-1时产生频率为fmin（Hz）的振荡信号；x=1时产生频率为fmax（Hz）的振荡信号。为了得到最好的结果，fmin和fmax的取值范围应该在0fs/