有关matlab命令.pptVIP

离散时间系统分析（续） 8. roots 用来求多项式的根。 调用格式为：r=roots(c) c:按降幂排列的多项式的系数； r:以c为系数的多项式的根组成的列向量； * * 《数字信号处理》中用到的相关MATLAB 命令简介 信号的产生、卷积和相关 1. rand 用来产生均值为0.5，幅度在0－1之间均匀分布的伪随机数， 在信号处理中用它近似均匀分布的白噪声信号。 举例：u=rand(N) u=rand(M,N) 信号的产生、卷积和相关(续) 2. randn 用来产生均值为零，方差为1，服从高斯（正态）分布的 白噪声信号。调用格式与rand相同。 例：产生零均值，功率为0.1且服从高斯分布的白噪声信号。 % to test randn p=0.1;N=500000; u=randn(1,N);a=sqrt(p); u=u*a;power_u=var(u); subplot(2,1,1);plot(u(1:100)); Subplot(2,1,2);plot(hist(u,50)); 信号的产生、卷积和相关(续) 3. sinc 用来产生在信号处理理论中经常用到的sinc函数 4.chirp 用来产生一中叫chirp的信号。该信号的表达式为： 调用格式为：x=chirp(T,F0,T1,F1) T横轴的时间范围向量; F0为起始频率; F1为T1时刻所具有的频率，缺省时， F0＝0, T1＝1, F1＝100 信号的产生、卷积和相关(续) 例：产生一chirp信号。 t=0:0.001:1; y=chirp(t,0,1,125);plot(t,y) 其它与信号产生相关的文件： diric:产生周期的sinc信号； gauspuls产生高斯信号； pulstran产生脉冲串信号； tripuls产生三角波脉冲信号； 信号的产生、卷积和相关(续) 5.Conv 用来实现两个离散序列的线性卷积。 调用格式为：y=conv(x,h) X,h：分别为进行线性卷积的两个序列。 6.xcorr: 用来求两个信号的互相关或一个信号的自相关。 调用格式为（1）rxy=xcorr(x,y) （2）rx=xcorr(x,Mlag,’flag’) 格式（1）为求两个序列x,y的互相关 格式（2）为求序列x的自相关； Mlag：表示rx的单边长度，总长度为2 Mlag＋1； flag：定标标志，若flag＝biased,表示是有偏估计；若 flag＝unbiased,表示是有偏估计. 信号的产生、卷积和相关(续) 离散时间系统分析 1 filter 在已知B(z),A(z)，但不知道h(n)的情况下，求一个离散系统 的输出。 调用格式为：y=filter(b,a,x) b:离散系统系统函数分子B(z)多项式系数向量； a:离散系统系统函数分母A(z)多项式系数向量； x:离散系统输入序列。 2.impz 在已知B(z),A(z)，但不知道h(n)的情况下，求一个离散系统的单位抽样响应h(n)。 调用格式为：（1）h=impz(b,a,N) （2）[h,t]=impz(b,a,N) b为离散系统系统函数分子B(z)多项式系数向量； a为离散系统系统函数分母A(z)多项式系数向量； N为所需的h(n)的长度。 格式（1）绘图时n从1开始，格式（2）绘图时n从0开始。 离散时间系统分析（续） 离散时间系统分析（续） 3.freqz 在已知B(z),A(z)的情况下，求系统的频率响应 调用格式为：[H,w]=freqz(b,a,N,’whole’,Fs) N:频率轴的分点数，建议N为2的整次幂； w:返回频率轴坐标向量； Fs:取样频率，若Fs＝1，频率轴为归一化频率 whole:指定计算的频率范围是0－Fs,缺省时为0－ Fs/2； 4.grpdelay 用来求一个系统的群延迟。 调用格式为：[gd w]=grpdelay(B,A,N)或 [gd w]=grpdelay(B,A,N，FS) B,A:H(z)的分子、分母多项式系数向量； gd:求出的群延迟； w:角频率，F为频率，两者的长度均为N； Fs:取样频率。 离散时间系统分析（续） 离散时间系统分析（续） 5. tf2zp 用来求离散时间系统系统函数H(z) =B(z)/A(z) 的零极点和增益； 调用格式为： [z,p,k] = tf2zp(b,a) b,a:分别为系统函数分子，分母多