?? sdir2cas.m
字號:
function [C,B,A] = sdir2cas(b,a);
% s平面中直接型到級聯型的型式轉換(復數對型)
% -------------------------------------------------
% [C,B,A] = sdir2cas(b,a)
% C = 增益系數
% B = 包含各bk的K乘3維實系數矩陣
% A = 包含各ak的K乘3維實系數矩陣
% b = 直接型的分子多項式系數
% a = 直接型的分母多項式系數
%
Na = length(a)-1; Nb = length(b)-1;
% 計算增益系數 C
b0 = b(1); b = b/b0;
a0 = a(1); a = a/a0;
C = b0/a0;
%
%:
p= cplxpair(roots(a)); K = floor(Na/2);
if K*2 == Na % 當Na為偶時計算
A = zeros(K,3);
for n=1:2:Na
Arow = p(n:1:n+1,:);
Arow = poly(Arow);
A(fix((n+1)/2),:) = real(Arow);
end
elseif Na == 1 % 當 Na = 1 時計算
A = [0 real(poly(p))];
else % 當 Na 為奇及 > 1時計算
A = zeros(K+1,3);
for n=1:2:2*K
Arow = p(n:1:n+1,:);
Arow = poly(Arow);
A(fix((n+1)/2),:) = real(Arow);
end
A(K+1,:) = [0 real(poly(p(Na)))];
end
% 分子的二階因子部分:
z = cplxpair(roots(b)); K = floor(Nb/2);
if Nb == 0 % 當 Nb = 0 時計算
B = [0 0 poly(z)];
elseif K*2 == Nb % 當 Nb 為偶時計算
B = zeros(K,3);
for n=1:2:Nb
Brow = z(n:1:n+1,:);
Brow = poly(Brow);
B(fix((n+1)/2),:) = real(Brow);
end
elseif Nb == 1 % 當 Nb = 1 時計算
B = [0 real(poly(z))];
else % 當 Nb 為奇且 > 1時計算
B = zeros(K+1,3);
for n=1:2:2*K
Brow = z(n:1:n+1,:);
Brow = poly(Brow);
B(fix((n+1)/2),:) = real(Brow);
end
B(K+1,:) = [0 real(poly(z(Nb)))];
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -