作者:范二小姐儿 | 来源:互联网 | 2023-10-11 16:09
本文由编程笔记#小编为大家整理,主要介绍了《DSP using MATLAB》Problem 6.21相关的知识,希望对你有一定的参考价值。
![技术分享图片](https://img6.php1.cn/3cdc5/c64b/1c8/d630577ad95d3cdb.jpeg)
代码:
%% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%% Output Info about this m-file
fprintf(‘
***********************************************************
‘);
fprintf(‘ Problem 6.21
‘);
banner();
%% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%format long;
format short;
fprintf(‘
FIR filter All-Zero LATTICE-form:
‘);
Klc = [2, 0.6, 0.3, 0.5, 0.9]
% ----------------------------------------------------------
% NOTE: linear-phase can not use LATTICE-form
% ----------------------------------------------------------
fprintf(‘
Convert All-Zero LATTICE-form to DIRECT-form :
‘);
[b] = latc2dir(Klc)
a = 1
fprintf(‘
Convert DIRECT-form to PARALLEL-form :
‘);
[C, Bp, Ap] = dir2par(b, a)
if size(C)==0
C = 0;
end
fprintf(‘
Convert DIRECT-form to CASCADE-form :
‘);
[b0, Bc, Ac] = dir2cas(b, a)
fprintf(‘
Convert TF-form to SOS-form :
‘);
[sos, g] = tf2sos(b, a)
fprintf(‘
Convert DIRECT-form to FREQUENCY-SAMPLE-form 1 :
‘);
[Cfs, Bfs, Afs] = dir2fs(b)
fprintf(‘
Convert DIRECT-form to FREQUENCY-SAMPLE-form 2 :
‘);
r = 0.99;
[Cfs_r, Bfs_r, Afs_r, rM] = dir2fs_r(b, r)
% -----------------------------------------
% START check
% -----------------------------------------
n = [0:7];
delta = impseq(0, 0, 7)
%format long
format short
hcas = casfiltr(b0, Bc, Ac, delta)
hltc = latcfilt(Klc, delta)
%hladr = ladrfilt(Klr, Clr, delta)
hdir = filter(b, a, delta)
% -------------------------------------------
% END check
% -------------------------------------------
运行结果:
得到的直接形式系数
![技术分享图片](https://img6.php1.cn/3cdc5/c64b/1c8/0c5886dea91f94c5.jpeg)
串联形式系数
![技术分享图片](https://img6.php1.cn/3cdc5/c64b/1c8/0c952da7d92d2011.jpeg)
频率采样形式1系数
![技术分享图片](https://img6.php1.cn/3cdc5/c64b/1c8/4f47a106557e6cff.jpeg)
频率采样形式2系数
![技术分享图片](https://img6.php1.cn/3cdc5/c64b/1c8/82e22268b1064461.jpeg)
Lattice和直接形式流图
![技术分享图片](https://img6.php1.cn/3cdc5/c64b/1c8/73c9319328642005.jpeg)