asymptotic magnitude bode plot
显示 更早的评论

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)')
采纳的回答
更多回答(0 个)
类别
在 帮助中心 和 File Exchange 中查找有关 Frequency-Domain Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!