matlba中扩展函数dir2fsfunction[C,B,A]=dir2fs(h)%直接型到频率采样型的转换%---------------------------------%[C,B,A]=dir2fs(h)%C=包含各并行部分增益的行向量%B=包含按行排列的分子系数矩阵%A=包

Question

matlba中扩展函数dir2fs
function [C,B,A] = dir2fs(h)
% 直接型到频率采样型的转换
% ---------------------------------
% [C,B,A] = dir2fs(h)
% C = 包含各并行部分增益的行向量
% B = 包含按行排列的分子系数矩阵
% A = 包含按行排列的分母系数矩阵
% h = FIR滤波器的脉冲响应向量
%
M = length(h);
H = fft(h,M);
magH = abs(H); phaH = angle(H)';
% check even or odd M
if (M == 2*floor(M/2))
L = M/2-1; % M为偶数
A1 = [1,-1,0;1,1,0];
C1 = [real(H(1)),real(H(L+2))];
else
L = (M-1)/2; % M is odd
A1 = [1,-1,0];
C1 = [real(H(1))];
end
k = [1:L]';
% 初始化 B 和 A 数组
B = zeros(L,2); A = ones(L,3);
% 计算分母系数
A(1:L,2) = -2*cos(2*pi*k/M); A = [A;A1];
% 计算分子系数
B(1:L,1) = cos(phaH(2:L+1));
B(1:L,2) = -cos(phaH(2:L+1)-(2*pi*k/M));
% 计算增益系数
C = [2*magH(2:L+1),C1]';

Answer 1

没有问题啊!结果是：
Optimization terminated:magnitude of directional derivative in search
direction less than 2*options.TolFun and maximum constraint violation
is less than options.TolCon.
Active inequalities (to within options.TolCon = 1e-006):
lower upper ineqlin ineqnonlin
3 2 9
8 6
9 7
10 12
13
14
x =
Columns 1 through 12
0.1164 1.0000 -0.0000 0.6726 0.8836 1.0000 1.0000 0.0000 -0.0000 0.0000 0.2641 1.0000
Columns 13 through 16
0.0000 -0.0000 0.7359 0.3274
fval =
104.8100
105.0249
exitflag =
105.0249
V =
0.0325
0.0324
0.0338
0.0463
0.0326
0.0342
0.0697