asymptotic magnitude bode plot

2 次查看(过去 30 天)
i want to find figure 2 for 'FRACTIONAL ZERO',i took the values as per the paper but both the graph is not crossing each other as shown in the figure,the code is given below
clc;
hold off;alpha=0.9;
Wcr=2.15845;a=2;
w1=logspace(-2,log10(Wcr),1000);w2=logspace(log10(Wcr),2,1000);
Mag1= (20*log10(abs(a)))*ones(1,numel(w1));Mag2=@(w2) 20*alpha*log10(w2);
semilogx(w1,Mag1);hold on;semilogx(w2,Mag2(w2));
grid on;xlabel('Frequency');ylabel('Magnitude dB')
clc;
clear
syms w alpha
alpha=0.9;a=2;w=logspace(-2,2,1000);
Mag1=@(w) (20*log10(abs(sqrt(((1j.*w).^alpha + a).^2))));
semilogx(w,Mag1(w));hold on;grid on;
xlabel('Frequency (rad/sec)');ylabel('Magnitude (dB)')

采纳的回答

Ameer Hamza
Ameer Hamza 2020-4-9
The graph in the paper is exaggerated to clarify the concepts about Asymptotes. But in fact, it is totally wrong. Asymptote is a line that comes infinitely closer to your curve as the frequency tends to infinity. The graph shown in the paper has started to move away from the Asymptote, and therefore, it is wrong. The graph generated by your code is correct.

更多回答(0 个)

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by