Добавил:
debilX2
t.me
Фулл всегда есть! А если нет, то Вы плохо его ищите! ИиКГ, СКДИКТ, ОКИТПЭС и тд https://t.me/whitedevil752rnТак же веду разработку КД (конструкторской документации) согласно ГОСТ. Имеется опыт работы при производстве на одном из ведущих в области радиэлектроники предприятии. Пишите)
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:pr5
.m clc
clear all
close all
f=(2:2:30)*10^6;
lambda=(3*10^8)./f;
SNRmin=10;
k=1.38*10^(-23);
df=3*10^3;
N=15;
T0=293;
r=30*10^3;
P=10;
e=4.4;
sigma=0.001;
E = load('fieldstrength.txt')';
E = (10.^(E/20))*1e-6;
nu = load('nu.dat')';
G = load('gain.dat')';
T_ash=[1.5e11 2.9e9 1.5e8 1e8 2.9e7 1.5e7 1.4e7 1.3e7 1.2e7 3.9e6 2.9e6 2.8e6 2.6e6 2.2e6 1.5e6];
T_pr=(10^(N/10)-1)*T0;
Psh=(k.*T_pr.*df)+(k.*T_ash.*df.*nu);
Pshdb=10*log10(Psh);
Ps=(E.^2.*lambda.^2.*G)./(480*pi^2);
Psdb=10*log10(Ps);
SNR=Ps./Psh;
SNRdb=10*log10(SNR);
figure(100)
plot(f,Psdb,f,Pshdb)
grid
title('Мощность шума и мощность сигнала от частоты')
xlabel('f, МГц')
ylabel('P, дБ/мВт (dBm)')
figure(200)
plot(f,SNRdb)
grid
title('ОСШ от частоты')
xlabel('f, МГц')
ylabel('SNR, дБ/мВт (dBm)')
clear all
close all
f=(2:2:30)*10^6;
lambda=(3*10^8)./f;
SNRmin=10;
k=1.38*10^(-23);
df=3*10^3;
N=15;
T0=293;
r=30*10^3;
P=10;
e=4.4;
sigma=0.001;
E = load('fieldstrength.txt')';
E = (10.^(E/20))*1e-6;
nu = load('nu.dat')';
G = load('gain.dat')';
T_ash=[1.5e11 2.9e9 1.5e8 1e8 2.9e7 1.5e7 1.4e7 1.3e7 1.2e7 3.9e6 2.9e6 2.8e6 2.6e6 2.2e6 1.5e6];
T_pr=(10^(N/10)-1)*T0;
Psh=(k.*T_pr.*df)+(k.*T_ash.*df.*nu);
Pshdb=10*log10(Psh);
Ps=(E.^2.*lambda.^2.*G)./(480*pi^2);
Psdb=10*log10(Ps);
SNR=Ps./Psh;
SNRdb=10*log10(SNR);
figure(100)
plot(f,Psdb,f,Pshdb)
grid
title('Мощность шума и мощность сигнала от частоты')
xlabel('f, МГц')
ylabel('P, дБ/мВт (dBm)')
figure(200)
plot(f,SNRdb)
grid
title('ОСШ от частоты')
xlabel('f, МГц')
ylabel('SNR, дБ/мВт (dBm)')