با توجه به اینکه X/D به ازای هر مقدار از g^2 و w یه عدد مختلط نتیجه میده پس باید واسه محور عمودی "قدر مطلق X/D" یا "قسمت حقیقی X/D" و یا "قسمت موهومی X/D" که یه عدد صحیحه در نظر بگیری.
در این صورت میتونی با نوشتن یه حلقه For تو در تو به راحتی نمودار این فرمول رو رسم کنی!
مثلا اگه منظور از X/D قدر مطلق X/D باشه کدش این میشه:
در این صورت میتونی با نوشتن یه حلقه For تو در تو به راحتی نمودار این فرمول رو رسم کنی!
مثلا اگه منظور از X/D قدر مطلق X/D باشه کدش این میشه:
کد پیاچپی:
clear all
clc
m=1;
XD=zeros(3,201);
absXD=zeros(3,201);
for g2=[.5 .1 .05]
n=1;
for w=0:.01:2
XD(m,n)=.62*w*1i/((1-g2)*w^4-.62*w^3*1i-2.012*w^2+.62*w*1i+1);
absXD(m,n)=abs(XD(m,n));
n=n+1;
end
m=m+1;
end
plot(0:.01:2,absXD(1,: ),'k')
hold on
plot(0:.01:2,absXD(2,: ),'r')
hold on
plot(0:.01:2,absXD(3,: ))
xlabel('\omega')
ylabel('ABS(X/D)')
legend('g^2 = 0.5','g^2 = 0.1','g^2 = 0.05')