2016-12-30 27 views
-1

我真的不明白如何绘制任意公式的传递函数。 Matlab文档说要得到波特图,我需要创建一个传递函数。他们举了一个例子:tf([]),但这只给了三级。如何在Matlab中绘制该公式的传递函数?

如何绘制特定公式的波特图?

它是:。。

Y = 1./((Z0./Zl).^2 *Ĵ罪(乙 d)/ 2 + EXP(j * B * d )*(1 + Z0。/ Z1))

其中Z1,d和B是每个频率值的向量。这是我完整的Matlab脚本。

numElem = 200; 
w = linspace(2*pi*10^8,2*pi*10^10,numElem); 
wRes = 2*pi*10^9; 

%%Insert wRes (resonant frequency) 
for n = 1:numElem 
    if w(n) > wRes 
     w = [w(1:n - 1), wRes, w(n:numElem)]; 
     resLocation = n; 
     break 
    end 
end 

freq = w ./ (2*pi); 
lambda = 299792458 ./ freq; 

%Beta 
B = pi/2 .* w ./wRes; 
B(resLocation) = pi/2; 

%Distance 
d = lambda(resLocation) ./ 4; 

Z0 = 50; %Ohms 
C = 10^(-12); 
L = 1/(wRes^2*C); 
Zl = 1 ./ (j .* w .* C + 1 ./ (j .* w .* L)); 

y = 1./((Z0./Zl).^2 .* j.*sin(B .* d)/2 + exp(j .* B .* d).*(1 + Z0 ./ Zl)) 

回答

0

我以前遇到过这个问题,它很简单。 您的幅度20log(|数据|)和你的相位是反正切(实(数据)/ IMAG(数据)),所以你可以用下面这个公式:

Amplitude=20*log10(abs(y)); 
phase = angle(y)*180/pi; 
subplot(211) 
semilogx(freq,Amplitude) 
subplot(212) 
semilogx(frq,phase) 

结果是: enter image description here