Добавил:
debilX2
t.me
Фулл всегда есть! А если нет, то Вы плохо его ищите! ИиКГ, СКДИКТ, ОКИТПЭС и тд https://t.me/whitedevil752rnТак же веду разработку КД (конструкторской документации) согласно ГОСТ. Имеется опыт работы при производстве на одном из ведущих в области радиэлектроники предприятии. Пишите)
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:pr3
.m clc
clear all
close all
%%%%%%%%%%%%%%%%%%%%
%Расчет согласования
%%%%%%%%%%%%%%%%%%%%
l=0.075;
a=0.003;
f=0.6*10^9:5*10^6:1.2*10^9;
c=3*10^8;
lambda=c./f;
k=2*pi./lambda;
Zist=50+1j*0;
for i=1:length(k)
fun= @(teta)((cos(k(i)*l*cos(teta))-cos(k(i)*l)).^2)/sin(teta);
R(i)=60*integral(fun,0,pi, 'arrayvalued',true);
alfa=R(i)/(120*l*(log(l/a)-1)*(1-sin(2*k(i)*l)/(2*k(i)*l)));
beta=k(i);
gamma(i)=alfa+j*beta;
Zv=120*(log(l/a)-1)*(1-j*alfa/beta);
Zvh(i)=Zv.*coth(gamma(i)*l);
Kotr=(real(Zvh)-Zist)/(real(Zvh)+Zist);
KSV(i)=(1+abs(Kotr))/(1-abs(Kotr));
end
figure(100)
plot(f,real(Zvh),f,imag(Zvh))
grid
xlabel('f'), ylabel('Zвх')
title('График зависимости вх. сопр. от f')
figure(200)
plot(f,KSV)
grid
xlabel('f'), ylabel('КСВ')
title('График зависимости КСВ от f')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Расчет ДН симметричного вибратора
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
f2=1*10^9;
lambda2=(3*10^8)/f2;
l2=[0.25*lambda2 0.5*lambda2 0.625*lambda2];
thetaindegrees=0:1:360;
theta=thetaindegrees*(pi/180);
k2=2*pi/lambda2;
figure(300)
for ii=l2
F2=abs((cos(k2*ii.*cos(theta))-cos(k2*ii))/(1-cos(k2*ii)).*sin(theta));
polar(theta,F2)
hold on
end
legend('l=0.25*lambda','l=0.5*lambda','l=0.625*lambda')
title('График ДН симметричного вибратора')
clear all
close all
%%%%%%%%%%%%%%%%%%%%
%Расчет согласования
%%%%%%%%%%%%%%%%%%%%
l=0.075;
a=0.003;
f=0.6*10^9:5*10^6:1.2*10^9;
c=3*10^8;
lambda=c./f;
k=2*pi./lambda;
Zist=50+1j*0;
for i=1:length(k)
fun= @(teta)((cos(k(i)*l*cos(teta))-cos(k(i)*l)).^2)/sin(teta);
R(i)=60*integral(fun,0,pi, 'arrayvalued',true);
alfa=R(i)/(120*l*(log(l/a)-1)*(1-sin(2*k(i)*l)/(2*k(i)*l)));
beta=k(i);
gamma(i)=alfa+j*beta;
Zv=120*(log(l/a)-1)*(1-j*alfa/beta);
Zvh(i)=Zv.*coth(gamma(i)*l);
Kotr=(real(Zvh)-Zist)/(real(Zvh)+Zist);
KSV(i)=(1+abs(Kotr))/(1-abs(Kotr));
end
figure(100)
plot(f,real(Zvh),f,imag(Zvh))
grid
xlabel('f'), ylabel('Zвх')
title('График зависимости вх. сопр. от f')
figure(200)
plot(f,KSV)
grid
xlabel('f'), ylabel('КСВ')
title('График зависимости КСВ от f')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Расчет ДН симметричного вибратора
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
f2=1*10^9;
lambda2=(3*10^8)/f2;
l2=[0.25*lambda2 0.5*lambda2 0.625*lambda2];
thetaindegrees=0:1:360;
theta=thetaindegrees*(pi/180);
k2=2*pi/lambda2;
figure(300)
for ii=l2
F2=abs((cos(k2*ii.*cos(theta))-cos(k2*ii))/(1-cos(k2*ii)).*sin(theta));
polar(theta,F2)
hold on
end
legend('l=0.25*lambda','l=0.5*lambda','l=0.625*lambda')
title('График ДН симметричного вибратора')