数字信号处理的matlab实现的常用算法集合

% 用于一书 % % 离散信号和系统 % conv_m - 改进的线性卷积子程序 (第22页) % conv_tp - 用Toeplitz矩阵计算的线性卷积(第34页) % evenodd - 将实信号分解为偶和奇两部分(第15页) % impseq - 产生脉冲序列 (第6页) % sigadd - 信号相加运算(第8页) % sigfold - 信号折叠运算(第10页) % sigmult - 信号乘法运算(第9页) %

function [CBA] = dir2par（ba）;
% DIRECT-form to PARALLEL-form conversion
% --------------------------------------
% [CBA] = dir2par（ba）
%  C = Polynomial part when length（b） >= length（a）
%  B = K by 2 matrix of real coefficients containing bk‘s
%  A = K by 3 matrix of real coefficients containing ak‘s
%  b = numerator polynomial coefficients of DIRECT form
%  a = denominator polynomial coefficients of DIRECT form
%

M = length（b）; N = length（a）;

[r1p1C] = residuez（ba）;
p = cplxpair（p110000000*eps）;
I = cplxcomp（p1p）;
r = r1（I）;

K = floor（N/2）; B = zeros（K2）; A = zeros（K3）;
if K*2 == N; %N even order of A（z） odd one factor is first order
        for i=1:2:N-2
                Brow = r（i:1:i+1:）;
                Arow = p（i:1:i+1:）;
                [BrowArow] = residuez（BrowArow[]）;
                B（fix（（i+1）/2）:） = real（Brow）;
                A（fix（（i+1）/2）:） = real（Arow）;
        end
        [BrowArow] = residuez（r（N-1）p（N-1）[]）;
        B（K:） = [real（Brow） 0]; A（K:） = [real（Arow） 0];
else
        for i=1:2:N-1
                Brow = r（i:1:i+1:）;
                Arow = p（i:1:i+1:）;
                [BrowArow] = residuez（BrowArow[]）;
                B（fix（（i+1）/2）:） = real（Brow）;
                A（fix（（i+1）/2）:） = real（Arow）;
        end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        765  1997-08-30 11:49  PWS_DSPC\AFD_BUTT.M

     文件        822  1997-08-30 11:48  PWS_DSPC\AFD_CHB1.M

     文件        818  1997-08-30 11:53  PWS_DSPC\AFD_CHB2.M

     文件        893  1997-09-01 15:28  PWS_DSPC\AFD_ELIP.M

     文件        915  1997-11-24 21:40  PWS_DSPC\AMPL_RES.M

     文件        214  1997-08-07 16:54  PWS_DSPC\BLACKMAN.M

     文件        386  1997-08-07 17:21  PWS_DSPC\CAS2DIR.M

     文件        436  1997-08-07 17:22  PWS_DSPC\CASFILTR.M

     文件        830  1997-08-30 11:59  PWS_DSPC\CHEB1HPF.M

     文件        319  1997-08-30 12:00  PWS_DSPC\CIRCEVOD.M

     文件        697  1997-08-30 12:03  PWS_DSPC\CIRCONVT.M

     文件        464  1997-08-30 12:04  PWS_DSPC\CIRSHFTT.M

     文件       4887  1995-08-13 16:26  PWS_DSPC\CONTENTS.M

     文件        376  1997-08-21 22:07  PWS_DSPC\CONV_M.M

     文件        377  1997-08-07 18:25  PWS_DSPC\CONV_TP.M

     文件        355  1997-08-07 18:31  PWS_DSPC\CPLXCOMP.M

     文件        299  1997-08-07 18:39  PWS_DSPC\DB2DELTA.M

     文件        357  1997-09-17 15:04  PWS_DSPC\DELTA2DB.M

     文件        569  1997-08-30 14:56  PWS_DSPC\DFS.M

     文件        520  1997-08-30 14:58  PWS_DSPC\DFT.M

     文件        893  1997-08-07 17:20  PWS_DSPC\DIR2CAS.M

     文件        807  1997-08-30 15:02  PWS_DSPC\DIR2FS.M

     文件        708  1997-08-30 15:04  PWS_DSPC\DIR2LADR.M

     文件        423  1997-08-30 15:15  PWS_DSPC\DIR2LATC.M

     文件       1334  1998-03-14 13:16  PWS_DSPC\Dir2par.m

     文件       1031  1997-08-30 15:17  PWS_DSPC\Dir2paro.m

     文件        409  1997-08-21 16:26  PWS_DSPC\EVENODD.M

    ..A..H.      6328  1996-04-23 11:43  PWS_DSPC\FINDER.DAT

     文件        564  1997-08-07 20:05  PWS_DSPC\FREQS_M.M

     文件        781  1997-08-07 20:15  PWS_DSPC\FREQZ_M.M

............此处省略43个文件信息