Материал: 893

68

if n1~=0; t(n1,n1)=t(n1,n1)+1/zn; end;

if n2~=0; t(n2,n2)=t(n2,n2)+1/zn; end;

if (n1~=0)&(n2~=0); t(n1,n2)=t(n1,n2)-1/zn;t(n2,n1)=t(n2,n1)-1/zn; end;

end;

%Индуктивность-L if pr=='L';

n1=ceil(dd{k+2});n2=ceil(dd{k+3});zn=abs(dd{k+1}); if n1~=0;

t(n1,n1)=t(n1,n1)+1/(j*zn*w);

end;

if n2~=0; t(n2,n2)=t(n2,n2)+1/(j*zn*w); end;

if (n1~=0)&(n2~=0); t(n1,n2)=t(n1,n2)-1/(j*zn*w);t(n2,n1)=t(n2,n1)-1/(j*zn*w); end;

end;

%Источник тока-I if pr=='I';

n1=ceil(dd{k+2});n2=ceil(dd{k+3});zn=abs(dd{k+1}); if n1~=0;

ww(n1,1)=ww(n1,1)+zn;

end;

if n2~=0; ww(n2,1)=ww(n2,1)-zn; end;

end;

%Переход к следующему элементу k=k+kp;

end;

disp('Проводимость-G : в схеме номинал задавать в омах!!!')

disp('Емкость-С') disp('Индуктивность-L') disp('Источник тока-I')

%

%Инт. цепь с С : fn=1e5 fk=1e7 ni=1 no=2 dd={'I' 1 1 0 'G' 1e3 1 2 'C' 25e-12 2 0};

%

ni=input('Номер входного узла= '); no=input('Номер выходного узла= '); fn=input('Нач. частота= ');

69

fk=input('Кон. частота= '); kf=input('Количество точек = '); df=(fk-fn)/kf;

[n,vv]=ytop(dd)

f=fn;k1=1; while f<=fk;

w=2*pi*f;

[t,ww]=ymatr(dd,n,w);

uv=t\ww;

uv1=uv(ni);uv2=uv(no);

ku=uv2/uv1;mf(k1)=f;mku(k1)=abs(ku);

ph(k1)=angle(ku)*180/pi;

k1=k1+1;

f=f+df;

end; %Вывод ЧХ k2=1:k1-1;

subplot(2,1,1);

semilogx(mf(k2),mku(k2)); title('АЧХ');

xlabel(' f,Гц');

ylabel('Ku,отн.ед.'); grid

subplot(2,1,2);

semilogx(mf(k2),ph(k2)); title('ФЧХ');

xlabel(' f,Гц');

ylabel('Phase,grad'); grid

disp(' End Mein ')