clear;clc; format compact; close all; % Minden törlése

     D=4;pt=50;s=0;Y=0;ff=0;
     R=27;L=52e-3;C=100e-6;

     for m=1:pt*D;k=m/pt;
     f=10^(s+k);w=2*pi*f;
     
     Y(m)=j*w*L/(j*3*w*L+R)-R/(R+1/(j*w*C));
     
     ff(m)=f;
     end;

figure;
     subplot(211);plot(Y);title('Nyquist diagram');
     %axis([-.5 0 -.15 .15]);
     grid
     
     subplot(212);Y;semilogx(ff,abs(Y));grid;
     axis([1 10000 0 .8]);
     ylabel('Amplitúdó abs(Y)'),xlabel('frekvencia [1/s]')
     title('Kimeneti jel');

figure;
     subplot(211);Ydb=20*log10(abs(Y));semilogx(ff,Ydb)
     axis([1 10000 -50 0]);ylabel('Amplitúdó [dB]');grid
     title('Amplitúdó Bode diagram')
     
     subplot(212);fi=180/pi*angle(Y);semilogx(ff,fi)
     axis([1 10000 -190 -70]);grid
     ylabel('fáziszsög  [fok]'),xlabel('frekvencia [1/s]')
     title('Fázis Bode diagram')