Поможем написать учебную работу
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Мордовский государственный университет имени Н.П.Огарёва
Специальность: автоматизированные системы обработки информации
и управления
ЛАБОРАТОРНАЯ РАБОТА №2
по теоретическим основам автоматизированного управления
Анализ системы обработки данных (СОД): оценка границ работоспособности СОД.
ЛР 02069964 ТОАУ 02 13
ВЫПОЛНИЛИ ПРОВЕРИЛ
студенты 441 группы ________01.10.2013
01.10.2013 Панкратов М. В.
Пырякина К. А.
Пигачев А. А.
Саранск 2013
Цель работы: моделирование типовой системы обработки данных (СОД) и оценка ее работоспособности при вариации параметров коэффициентов потерь.
Практическая часть.
Задание 1:
Выполнение задания:
с1 =17
с2=25
L1=19
L2=27
Т_предел=3.0207
clear,clc,close all
helpdlg('c1 - НОРМАТИВНАЯ ТРУДОЕМКОСТЬ РАБОТЫ С МАССИВАМИ В МИНУТАХ 1-Й ПОДСИСТЕМЫ; c2 - НОР-МАТИВНАЯ ТРУДОЕМКОСТЬ РАБОТЫ С МАССИВАМИ В МИНУТАХ 2-Й ПОДСИСТЕМЫ; L1 - ЧИСЛО ЭВМ В 1-Й ПОД-СИСТЕМЕ; L2 - ЧИСЛО ЭВМ ВО 2-Й ПОДСИСТЕМЕ.','СПРАВКА: НАЗНАЧЕНИЕ ПАРАМЕТРОВ')
disp('ВВЕДИТЕ ПАРАМЕТРЫ c1,c2,L1,L2 В КОМАНДНОЙ СТРОКЕ')
%- АНАЛИЗ ВВОДА НОРМАТИВНОЙ ТРУДОЕМКОСТИ РАБОТЫ 1-Й ПОДСИСТЕМЫ В МИНУТАХ (ЧАСАХ) -
c1=input('ВВЕДИТЕ c1: ');
if c1==0|c1<0|c1==inf|isreal(c1)==0
errordlg('ОШИБКА ВВОДА: c1 ДОЛЖНО БЫТЬ КОНЕЧНЫМ ДЕЙСТВИТЕЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ! ПО-ВТОРИТЕ ВВОД')
c1=input('ВВЕДИТЕ c1: ');
if c1==0|c1<0|c1==inf|isreal(c1)==0
errordlg(' ПОВТОРНАЯ ОШИБКА ВВОДА! ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%-------------------------------------------------------------------------------------------------------------------
% -АНАЛИЗ ВВОДА НОРМАТИВНОЙ ТРУДОЕМКОСТИ РАБОТЫ 2-Й ПОДСИСТЕМЫ В МИНУТАХ (ЧАСАХ) -
c2=input('ВВЕДИТЕ c2: ');
if c2==0|c2<0|c2==inf|isreal(c2)==0
errordlg('ОШИБКА ВВОДА: c2 ДОЛЖНО БЫТЬ КОНЕЧНЫМ ДЕЙСТВИТЕЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ! ПО-ВТОРИТЕ ВВОД')
c2=input('ВВЕДИТЕ c2: ');
if c2==0|c2<0|c2==inf|isreal(c2)==0
errordlg(' ПОВТОРНАЯ ОШИБКА ВВОДА! ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%---------------------------------------------------------------------------------------------------------------
% ---------------------- АНАЛИЗ ВВОДА ЧИСЛА ЭВМ 1-Й ПОДСИСТЕМЫ ------------------
L1=input('ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): ');
if L1==0|L1<0|isreal(L1)==0|L1-fix(L1)~=0
errordlg('ОШИБКА ВВОДА ЧИСЛО ЭВМ ДОЛЖНО БЫТЬ НАТУРАЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
L1=input('ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): ');
if L1==0|L1<0|L1-fix(L1)~=0
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ЧИСЛА ЭВМ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%------------------------------------------------------------------------------------------------------------
%----------------------- АНАЛИЗ ВВОДА ЧИСЛА ЭВМ 2-Й ПОДСИСТЕМЫ ---------------
L2=input('ВВЕДИТЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: ');
if L2==0|L2<0|L2-fix(L2)~=0
errordlg('ОШИБКА ВВОДА ЧИСЛО ЭВМ ДОЛЖНО БЫТЬ НАТУРАЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
L1=input('ВВЕДИТЕ ДОПУСТИМОЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: ');
if L2==0|L2<0|L2-fix(L2)~=0
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ЧИСЛА ЭВМ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%----------------------------------------------------------------------------------------------------
% -------------------------- АНАЛИЗ ИНТЕРВАЛА ВРЕМЕНИ --------------------------
disp('------------------------------------------------------------------------------')
disp('Tmin - МИНИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ')
Tmin=c1/L1+c2/L2
disp('Tmax - МАКСИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ')
Tmax=c1+c2
disp('-----------------------------------------------------------------------------------------')
T=input('ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: ');
if isempty(T)==1|T<=0|isreal(T)==0|T==inf|T<Tmin|T>Tmax
errordlg('ОШИБКА ВВОДА ВРЕМЕНИ ОБРАБОТКИ. ВРЕМЯ ДОЛЖНО БЫТЬ ДЕЙСТВИТЕЛЬНЫМ КОНЕЧНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
T0=input('ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: ');
if isempty(T)==1|T<=0|isreal(T)==0|T==inf|T<Tmin|T>Tmax
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ВРЕМЕНИ ОБРАБОТКИ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
disp('----------------------------------------------------------------------------------------')
t=rand(1);
if T==Tmax
disp('ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ')
t=T;
else
disp('ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ')
disp('t - СЛУЧАЙНО НАЗНАЧАЕМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ В СОД')
t=t*(T-Tmin)+Tmin
end
X1=(t-c2/L2)*L1/c1;
if X1>L1
X1=L1;
end
X2=(t-c2)*L1/c1;
if X2<1
X2=1;
end
ALPHA1=rand(1);
disp('ALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ В 1-Й ПОДСИСТЕМЕ')
ALPHA1=ALPHA1*(X1-X2)+X2
disp('ALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ ВО 2-Й ПОДСИСТЕМЕ')
ALPHA2=(t-c1*ALPHA1/L1)*L2/c2
ALPHA1t=[0 t*L1/c1];
ALPHA2t=[t*L2/c2 0];
disp('kvadrat - КООРДИНАТЫ ДОПУСТИМОГО РАБОЧЕГО ПРЯМОУГОЛЬНИКА НА ПЛОСКОСТИ: ')
kvadrat=[1 1;L1 L2]
kvadratx=[1 L1 L1 1 1];
kvadraty=[L2 L2 1 1 L2];
Tminx=[Tmin*L1/c1 0];
Tminy=[0 Tmin*L2/c2];
Tmaxx=[Tmax*L1/c1 0];
Tmaxy=[0 Tmax*L2/c2];
Tx=[0 T*L1/c1];
Ty=[T*L2/c2 0];
figure(1),
plot(kvadratx,kvadraty,'linewidth',2), grid on, hold on,
plot(Tx,Ty,'k:','linew',2),
plot(Tminx,Tminy,'m:','linew',2),
plot(Tmaxx,Tmaxy,'g:','linew',2)
plot(ALPHA1t,ALPHA2t,'r','linew',2),
plot(ALPHA1,ALPHA2,'r*','linew',2),
xlabel('\bf\alpha_1'),ylabel('\bf\alpha_2'),
legend('Valid square','T extreme','Tmin','Tmax','T','Random point'),
set(gcf,'color','w'),hold off
%--------------------- ВЫБОР ТОЧКИ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ ------------------------
disp('---------------------------------------------------------------------------------------------------------')
disp('ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ОБРАБОТКИ ДАННЫХ ПО ЗАДАННЫМ КОЭФФИЦИЕНТАМ ПОТЕРЬ')
TALPHA1=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: ');
if (TALPHA1<1)|(TALPHA1>L1)|TALPHA1==inf|isreal(TALPHA1)==0
errordlg('ОШИБКА ВЫБОРА: КООРДИНАТА (TALPHA1) НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПОВТОРИТЕ ВВОД ДЕЙСТВИТЕЛЬНОГО ЧИСЛА')
TALPHA1=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: ');
if TALPHA1<1|TALPHA1>L1|TALPHA1==inf|isreal(TALPHA1)==0
errordlg('ПОВТОРНАЯ ОШИБКА ВЫБОРА: КООРДИНАТА TALPHA1 НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
close all
break
end
end
TALPHA2=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: ');
if (TALPHA2<1)|(TALPHA2>L2)|TALPHA2==inf|isreal(TALPHA2)==0
errordlg('ОШИБКА ВЫБОРА: КООРДИНАТА (TALPHA2) НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПОВТОРИТЕ ВВОД ДЕЙСТВИТЕЛЬНОГО ЧИСЛА')
TALPHA2=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: ');
if TALPHA2<1|TALPHA2>L2|TALPHA2==inf|isreal(TALPHA2)==0
errordlg('ПОВТОРНАЯ ОШИБКА ВЫБОРА: КООРДИНАТА TALPHA2 НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
close all
break
end
end
%--------------------------------------------------------------------------------------------
disp(' Tpoint - ВРЕМЯ ОБРАБОТКИ')
Tpoint=c1*TALPHA1/L1+c2*TALPHA2/L2
xpoint=[0 Tpoint*L1/c1];
ypoint=[Tpoint*L2/c2 0];
disp('---------------------------------------------------------------------------------------')
%----------------------------- ГРАФИЧЕСКИЕ ПОСТРОЕНИЯ --------------------------------------
figure(2),
plot(kvadratx,kvadraty,'linew',2), grid on, hold on,
plot(Tx,Ty,'k:','linew',2),
plot(Tminx,Tminy,'m:','linew',2),
plot(Tmaxx,Tmaxy,'g:','linew',2)
plot(ALPHA1t,ALPHA2t,'r','linew',2),
plot(ALPHA1,ALPHA2,'r*','linew',2),
plot(xpoint,ypoint,'k','linew',2),
plot(TALPHA1,TALPHA2,'k*','linew',2),
xlabel('\bf\alpha_1'),ylabel('\bf\alpha_2'),
legend('Volid square','T extreme','Tmin','Tmax','T','Random point','Tpoint','Input point'),
set(gcf,'color','w'),hold off
В результате получили:
ВВЕДИТЕ ПАРАМЕТРЫ c1,c2,L1,L2 В КОМАНДНОЙ СТРОКЕ
ВВЕДИТЕ c1: 17
ВВЕДИТЕ c2: 25
ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): 19
ВВЕДИТЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: 27
------------------------------------------------------------------------------
Tmin - МИНИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ
Tmin =
1.8207
Tmax - МАКСИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ
Tmax =
42
-----------------------------------------------------------------------------------------
ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: 3.0207
----------------------------------------------------------------------------------------
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ
t - СЛУЧАЙНО НАЗНАЧАЕМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ В СОД
t =
2.7124
ALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ В 1-Й ПОДСИСТЕМЕ
ALPHA1 =
1.3909
ALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ ВО 2-Й ПОДСИСТЕМЕ
ALPHA2 =
1.5854
kvadrat - КООРДИНАТЫ ДОПУСТИМОГО РАБОЧЕГО ПРЯМОУГОЛЬНИКА НА ПЛОСКОСТИ:
kvadrat =
1 1
19 27
---------------------------------------------------------------------------------------------------------
ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ОБРАБОТКИ ДАННЫХ ПО ЗАДАННЫМ КОЭФФИЦИЕНТАМ ПОТЕРЬ
ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: 1.3909
ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: 1.5854
Tpoint - ВРЕМЯ ОБРАБОТКИ
Tpoint =
2.7125
------------------------
В результате при Tmin получили:
ВВЕДИТЕ ПАРАМЕТРЫ c1,c2,L1,L2 В КОМАНДНОЙ СТРОКЕ
ВВЕДИТЕ c1: 17
ВВЕДИТЕ c2: 25
ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): 19
ВВЕДИТЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: 27
------------------------------------------------------------------------------
Tmin - МИНИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ
Tmin =
1.8207
Tmax - МАКСИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ
Tmax =
42
-----------------------------------------------------------------------------------------
ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: 1.8207
----------------------------------------------------------------------------------------
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ
t - СЛУЧАЙНО НАЗНАЧАЕМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ В СОД
t =
1.8207
ALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ В 1-Й ПОДСИСТЕМЕ
ALPHA1 =
1.0000
ALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ ВО 2-Й ПОДСИСТЕМЕ
ALPHA2 =
1.0000
kvadrat - КООРДИНАТЫ ДОПУСТИМОГО РАБОЧЕГО ПРЯМОУГОЛЬНИКА НА ПЛОСКОСТИ:
kvadrat =
1 1
19 27
---------------------------------------------------------------------------------------------------------
ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ОБРАБОТКИ ДАННЫХ ПО ЗАДАННЫМ КОЭФФИЦИЕНТАМ ПОТЕРЬ
ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: 1
ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: 1
Tpoint - ВРЕМЯ ОБРАБОТКИ
Tpoint =
1.8207
----------------------------------------------------------------------
В результате при Tmax получили:
ВВЕДИТЕ ПАРАМЕТРЫ c1,c2,L1,L2 В КОМАНДНОЙ СТРОКЕ
ВВЕДИТЕ c1: 17
ВВЕДИТЕ c2: 25
ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): 19
ВВЕДИТЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: 27
------------------------------------------------------------------------------
Tmin - МИНИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ
Tmin =
1.8207
Tmax - МАКСИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ
Tmax =
42
-----------------------------------------------------------------------------------------
ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: 42
----------------------------------------------------------------------------------------
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ
ALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ В 1-Й ПОДСИСТЕМЕ
ALPHA1 =
19
ALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ ВО 2-Й ПОДСИСТЕМЕ
ALPHA2 =
27
kvadrat - КООРДИНАТЫ ДОПУСТИМОГО РАБОЧЕГО ПРЯМОУГОЛЬНИКА НА ПЛОСКОСТИ:
kvadrat =
1 1
19 27
---------------------------------------------------------------------------------------------------------
ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ОБРАБОТКИ ДАННЫХ ПО ЗАДАННЫМ КОЭФФИЦИЕНТАМ ПОТЕРЬ
ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: 19
ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: 27
Tpoint - ВРЕМЯ ОБРАБОТКИ
Tpoint =
42
----------------------------------------------------------------------
clear,clc,close all
helpdlg('c1 - НОРМАТИВНАЯ ТРУДОЕМКОСТЬ РАБОТЫ С МАССИВАМИ В МИНУТАХ 1-Й ПОДСИСТЕМЫ; c2 - НОР-МАТИВНАЯ ТРУДОЕМКОСТЬ РАБОТЫ С МАССИВАМИ В МИНУТАХ 2-Й ПОДСИСТЕМЫ; L1 - ЧИСЛО ЭВМ В 1-Й ПОД-СИСТЕМЕ; L2 - ЧИСЛО ЭВМ ВО 2-Й ПОДСИСТЕМЕ.','СПРАВКА: НАЗНАЧЕНИЕ ПАРАМЕТРОВ')
disp('ВВЕДИТЕ ПАРАМЕТРЫ c1,c2,L1,L2 В КОМАНДНОЙ СТРОКЕ')
%- АНАЛИЗ ВВОДА НОРМАТИВНОЙ ТРУДОЕМКОСТИ РАБОТЫ 1-Й ПОДСИСТЕМЫ В МИНУТАХ (ЧАСАХ) -
c1=input('ВВЕДИТЕ c1: ');
if c1==0|c1<0|c1==inf|isreal(c1)==0
errordlg('ОШИБКА ВВОДА: c1 ДОЛЖНО БЫТЬ КОНЕЧНЫМ ДЕЙСТВИТЕЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ! ПО-ВТОРИТЕ ВВОД')
c1=input('ВВЕДИТЕ c1: ');
if c1==0|c1<0|c1==inf|isreal(c1)==0
errordlg(' ПОВТОРНАЯ ОШИБКА ВВОДА! ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%-------------------------------------------------------------------------------------------------------------------
% -АНАЛИЗ ВВОДА НОРМАТИВНОЙ ТРУДОЕМКОСТИ РАБОТЫ 2-Й ПОДСИСТЕМЫ В МИНУТАХ (ЧАСАХ) -
c2=input('ВВЕДИТЕ c2: ');
if c2==0|c2<0|c2==inf|isreal(c2)==0
errordlg('ОШИБКА ВВОДА: c2 ДОЛЖНО БЫТЬ КОНЕЧНЫМ ДЕЙСТВИТЕЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ! ПО-ВТОРИТЕ ВВОД')
c2=input('ВВЕДИТЕ c2: ');
if c2==0|c2<0|c2==inf|isreal(c2)==0
errordlg(' ПОВТОРНАЯ ОШИБКА ВВОДА! ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%---------------------------------------------------------------------------------------------------------------
% ---------------------- АНАЛИЗ ВВОДА ЧИСЛА ЭВМ 1-Й ПОДСИСТЕМЫ ------------------
L1=input('ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): ');
if L1==0|L1<0|isreal(L1)==0|L1-fix(L1)~=0
errordlg('ОШИБКА ВВОДА ЧИСЛО ЭВМ ДОЛЖНО БЫТЬ НАТУРАЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
L1=input('ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): ');
if L1==0|L1<0|L1-fix(L1)~=0
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ЧИСЛА ЭВМ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%------------------------------------------------------------------------------------------------------------
%----------------------- АНАЛИЗ ВВОДА ЧИСЛА ЭВМ 2-Й ПОДСИСТЕМЫ ---------------
L2=input('ВВЕДИТЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: ');
if L2==0|L2<0|L2-fix(L2)~=0
errordlg('ОШИБКА ВВОДА ЧИСЛО ЭВМ ДОЛЖНО БЫТЬ НАТУРАЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
L1=input('ВВЕДИТЕ ДОПУСТИМОЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: ');
if L2==0|L2<0|L2-fix(L2)~=0
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ЧИСЛА ЭВМ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%----------------------------------------------------------------------------------------------------
% -------------------------- АНАЛИЗ ИНТЕРВАЛА ВРЕМЕНИ --------------------------
disp('------------------------------------------------------------------------------')
disp('Tmin - МИНИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ')
Tmin=c1/L1+c2/L2
disp('Tmax - МАКСИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ')
Tmax=c1+c2
disp('-----------------------------------------------------------------------------------------')
T=input('ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: ');
if isempty(T)==1|T<=0|isreal(T)==0|T==inf|T<Tmin|T>Tmax
errordlg('ОШИБКА ВВОДА ВРЕМЕНИ ОБРАБОТКИ. ВРЕМЯ ДОЛЖНО БЫТЬ ДЕЙСТВИТЕЛЬНЫМ КОНЕЧНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
T0=input('ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: ');
if isempty(T)==1|T<=0|isreal(T)==0|T==inf|T<Tmin|T>Tmax
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ВРЕМЕНИ ОБРАБОТКИ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
disp('----------------------------------------------------------------------------------------')
t=rand(1);
if T==Tmax
disp('ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ')
t=T;
else
disp('ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ')
disp('t - СЛУЧАЙНО НАЗНАЧАЕМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ В СОД')
t=t*(T-Tmin)+Tmin
end
X1=(t-c2/L2)*L1/c1;
if X1>L1
X1=L1;
end
X2=(t-c2)*L1/c1;
if X2<1
X2=1;
end
ALPHA1=rand(1);
disp('ALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ В 1-Й ПОДСИСТЕМЕ')
ALPHA1=ALPHA1*(X1-X2)+X2
disp('ALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ ВО 2-Й ПОДСИСТЕМЕ')
ALPHA2=(t-c1*ALPHA1/L1)*L2/c2
ALPHA1t=[0 t*L1/c1];
ALPHA2t=[t*L2/c2 0];
disp('kvadrat - КООРДИНАТЫ ДОПУСТИМОГО РАБОЧЕГО ПРЯМОУГОЛЬНИКА НА ПЛОСКОСТИ: ')
kvadrat=[1 1;L1 L2]
kvadratx=[1 L1 L1 1 1];
kvadraty=[L2 L2 1 1 L2];
Tminx=[Tmin*L1/c1 0];
Tminy=[0 Tmin*L2/c2];
Tmaxx=[Tmax*L1/c1 0];
Tmaxy=[0 Tmax*L2/c2];
Tx=[0 T*L1/c1];
Ty=[T*L2/c2 0];
figure(1),
plot(kvadratx,kvadraty,'linewidth',2), grid on, hold on,
plot(Tx,Ty,'k:','linew',2),
plot(Tminx,Tminy,'m:','linew',2),
plot(Tmaxx,Tmaxy,'g:','linew',2)
plot(ALPHA1t,ALPHA2t,'r','linew',2),
plot(ALPHA1,ALPHA2,'r*','linew',2),
xlabel('\bf\alpha_1'),ylabel('\bf\alpha_2'),
legend('Valid square','T extreme','Tmin','Tmax','T','Random point'),
set(gcf,'color','w'),hold off
%--------------------- ВЫБОР ТОЧКИ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ ------------------------
disp('---------------------------------------------------------------------------------------------------------')
disp('ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ОБРАБОТКИ ДАННЫХ ПО ЗАДАННЫМ КОЭФФИЦИЕНТАМ ПОТЕРЬ')
TALPHA1=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: ');
if (TALPHA1<1)|(TALPHA1>L1)|TALPHA1==inf|isreal(TALPHA1)==0
errordlg('ОШИБКА ВЫБОРА: КООРДИНАТА (TALPHA1) НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПОВТОРИТЕ ВВОД ДЕЙСТВИТЕЛЬНОГО ЧИСЛА')
TALPHA1=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: ');
if TALPHA1<1|TALPHA1>L1|TALPHA1==inf|isreal(TALPHA1)==0
errordlg('ПОВТОРНАЯ ОШИБКА ВЫБОРА: КООРДИНАТА TALPHA1 НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
close all
break
end
end
TALPHA2=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: ');
if (TALPHA2<1)|(TALPHA2>L2)|TALPHA2==inf|isreal(TALPHA2)==0
errordlg('ОШИБКА ВЫБОРА: КООРДИНАТА (TALPHA2) НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПОВТОРИТЕ ВВОД ДЕЙСТВИТЕЛЬНОГО ЧИСЛА')
TALPHA2=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: ');
if TALPHA2<1|TALPHA2>L2|TALPHA2==inf|isreal(TALPHA2)==0
errordlg('ПОВТОРНАЯ ОШИБКА ВЫБОРА: КООРДИНАТА TALPHA2 НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
close all
break
end
end
%--------------------------------------------------------------------------------------------
disp(' Tpoint - ВРЕМЯ ОБРАБОТКИ')
Tpoint=c1*TALPHA1/L1+c2*TALPHA2/L2
xpoint=[0 Tpoint*L1/c1];
ypoint=[Tpoint*L2/c2 0];
disp('---------------------------------------------------------------------------------------')
%----------------- ПРОВЕРКА РЕЗУЛЬТАТОВ -------------------------
if TALPHA1*c1/L1+TALPHA2*c2/L2==Tpoint
helpdlg('РЕЗУЛЬТАТЫ ВЕРНЫ')
else errordlg('РЕЗУЛЬТАТЫ НЕВЕРНЫ')
end
%----------------------------- ГРАФИЧЕСКИЕ ПОСТРОЕНИЯ --------------------------------------
figure(2),
plot(kvadratx,kvadraty,'linew',2), grid on, hold on,
plot(Tx,Ty,'k:','linew',2),
plot(Tminx,Tminy,'m:','linew',2),
plot(Tmaxx,Tmaxy,'g:','linew',2)
plot(ALPHA1t,ALPHA2t,'r','linew',2),
plot(ALPHA1,ALPHA2,'r*','linew',2),
plot(xpoint,ypoint,'k','linew',2),
plot(TALPHA1,TALPHA2,'k*','linew',2),
xlabel('\bf\alpha_1'),ylabel('\bf\alpha_2'),
legend('Volid square','T extreme','Tmin','Tmax','T','Random point','Tpoint','Input point'),
set(gcf,'color','w'),hold off
clear,clc,close all
helpdlg('c1 - НОРМАТИВНАЯ ТРУДОЕМКОСТЬ РАБОТЫ С МАССИВАМИ В МИНУТАХ 1-Й ПОДСИСТЕМЫ; c2 - НОР-МАТИВНАЯ ТРУДОЕМКОСТЬ РАБОТЫ С МАССИВАМИ В МИНУТАХ 2-Й ПОДСИСТЕМЫ; L1 - ЧИСЛО ЭВМ В 1-Й ПОД-СИСТЕМЕ; L2 - ЧИСЛО ЭВМ ВО 2-Й ПОДСИСТЕМЕ.','СПРАВКА: НАЗНАЧЕНИЕ ПАРАМЕТРОВ')
disp('ВВЕДИТЕ ПАРАМЕТРЫ c1,c2,L1,L2 В КОМАНДНОЙ СТРОКЕ')
%- АНАЛИЗ ВВОДА НОРМАТИВНОЙ ТРУДОЕМКОСТИ РАБОТЫ 1-Й ПОДСИСТЕМЫ В МИНУТАХ (ЧАСАХ) -
c1=input('ВВЕДИТЕ c1: ');
if c1==0|c1<0|c1==inf|isreal(c1)==0
errordlg('ОШИБКА ВВОДА: c1 ДОЛЖНО БЫТЬ КОНЕЧНЫМ ДЕЙСТВИТЕЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ! ПО-ВТОРИТЕ ВВОД')
c1=input('ВВЕДИТЕ c1: ');
if c1==0|c1<0|c1==inf|isreal(c1)==0
errordlg(' ПОВТОРНАЯ ОШИБКА ВВОДА! ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%-------------------------------------------------------------------------------------------------------------------
% -АНАЛИЗ ВВОДА НОРМАТИВНОЙ ТРУДОЕМКОСТИ РАБОТЫ 2-Й ПОДСИСТЕМЫ В МИНУТАХ (ЧАСАХ) -
c2=input('ВВЕДИТЕ c2: ');
if c2==0|c2<0|c2==inf|isreal(c2)==0
errordlg('ОШИБКА ВВОДА: c2 ДОЛЖНО БЫТЬ КОНЕЧНЫМ ДЕЙСТВИТЕЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ! ПО-ВТОРИТЕ ВВОД')
c2=input('ВВЕДИТЕ c2: ');
if c2==0|c2<0|c2==inf|isreal(c2)==0
errordlg(' ПОВТОРНАЯ ОШИБКА ВВОДА! ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%---------------------------------------------------------------------------------------------------------------
% ---------------------- АНАЛИЗ ВВОДА ЧИСЛА ЭВМ 1-Й ПОДСИСТЕМЫ ------------------
L1=input('ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): ');
if L1==0|L1<0|isreal(L1)==0|L1-fix(L1)~=0
errordlg('ОШИБКА ВВОДА ЧИСЛО ЭВМ ДОЛЖНО БЫТЬ НАТУРАЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
L1=input('ВВЕДИТЕ ЧИСЛО ЭВМ 1-Й ПОДСИСТЕМЫ (L1): ');
if L1==0|L1<0|L1-fix(L1)~=0
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ЧИСЛА ЭВМ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%------------------------------------------------------------------------------------------------------------
%----------------------- АНАЛИЗ ВВОДА ЧИСЛА ЭВМ 2-Й ПОДСИСТЕМЫ ---------------
L2=input('ВВЕДИТЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: ');
if L2==0|L2<0|L2-fix(L2)~=0
errordlg('ОШИБКА ВВОДА ЧИСЛО ЭВМ ДОЛЖНО БЫТЬ НАТУРАЛЬНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
L1=input('ВВЕДИТЕ ДОПУСТИМОЕ ЧИСЛО ЭВМ 2-Й ПОДСИСТЕМЫ L2: ');
if L2==0|L2<0|L2-fix(L2)~=0
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ЧИСЛА ЭВМ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
%----------------------------------------------------------------------------------------------------
% -------------------------- АНАЛИЗ ИНТЕРВАЛА ВРЕМЕНИ --------------------------
disp('------------------------------------------------------------------------------')
disp('Tmin - МИНИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ')
Tmin=c1/L1+c2/L2
disp('Tmax - МАКСИМАЛЬНО ДОПУСТИМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ')
Tmax=c1+c2
disp('-----------------------------------------------------------------------------------------')
T=input('ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: ');
if isempty(T)==1|T<=0|isreal(T)==0|T==inf|T<Tmin|T>Tmax
errordlg('ОШИБКА ВВОДА ВРЕМЕНИ ОБРАБОТКИ. ВРЕМЯ ДОЛЖНО БЫТЬ ДЕЙСТВИТЕЛЬНЫМ КОНЕЧНЫМ ЧИСЛОМ БОЛЬШЕ НУЛЯ. ПОВТОРИТЕ ВВОД')
T0=input('ВВЕДИТЕ ВРЕМЯ ОБРАБОТКИ T_предел: ');
if isempty(T)==1|T<=0|isreal(T)==0|T==inf|T<Tmin|T>Tmax
errordlg('ПОВТОРНАЯ ОШИБКА ВВОДА ВРЕМЕНИ ОБРАБОТКИ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
break
end
end
disp('----------------------------------------------------------------------------------------')
t=rand(1);
if T==Tmax
disp('ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ')
t=T;
else
disp('ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ ДЛЯ ЗАДАННОГО ВРЕМЕНИ ОБРАБОТКИ')
disp('t - СЛУЧАЙНО НАЗНАЧАЕМОЕ ВРЕМЯ ОБРАБОТКИ ДАННЫХ В СОД')
t=t*(T-Tmin)+Tmin
end
X1=(t-c2/L2)*L1/c1;
if X1>L1
X1=L1;
end
X2=(t-c2)*L1/c1;
if X2<1
X2=1;
end
ALPHA1=rand(1);
disp('ALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ В 1-Й ПОДСИСТЕМЕ')
ALPHA1=ALPHA1*(X1-X2)+X2
disp('ALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ ВО 2-Й ПОДСИСТЕМЕ')
ALPHA2=(t-c1*ALPHA1/L1)*L2/c2
ALPHA1t=[0 t*L1/c1];
ALPHA2t=[t*L2/c2 0];
disp('kvadrat - КООРДИНАТЫ ДОПУСТИМОГО РАБОЧЕГО ПРЯМОУГОЛЬНИКА НА ПЛОСКОСТИ: ')
kvadrat=[1 1;L1 L2]
kvadratx=[1 L1 L1 1 1];
kvadraty=[L2 L2 1 1 L2];
Tminx=[Tmin*L1/c1 0];
Tminy=[0 Tmin*L2/c2];
Tmaxx=[Tmax*L1/c1 0];
Tmaxy=[0 Tmax*L2/c2];
Tx=[0 T*L1/c1];
Ty=[T*L2/c2 0];
figure(1),
plot(kvadratx,kvadraty,'linewidth',2), grid on, hold on,
plot(Tx,Ty,'k:','linew',2),
plot(Tminx,Tminy,'m:','linew',2),
plot(Tmaxx,Tmaxy,'g:','linew',2)
plot(ALPHA1t,ALPHA2t,'r','linew',2),
plot(ALPHA1,ALPHA2,'r*','linew',2),
xlabel('\bf\alpha_1'),ylabel('\bf\alpha_2'),
legend('Valid square','T extreme','Tmin','Tmax','T','Random point'),
set(gcf,'color','w'),hold off
%--------------------- ВЫБОР ТОЧКИ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ ------------------------
disp('---------------------------------------------------------------------------------------------------------')
disp('ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ОБРАБОТКИ ДАННЫХ ПО ЗАДАННЫМ КОЭФФИЦИЕНТАМ ПОТЕРЬ')
% TALPHA1=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: ');
% if (TALPHA1<1)|(TALPHA1>L1)|TALPHA1==inf|isreal(TALPHA1)==0
% errordlg('ОШИБКА ВЫБОРА: КООРДИНАТА (TALPHA1) НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПОВТОРИТЕ ВВОД ДЕЙСТВИТЕЛЬНОГО ЧИСЛА')
% TALPHA1=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA1 - КОЭФФИЦИЕНТ ПОТЕРЬ 1-Й ПОДСИСТЕМЫ) ПО ОСИ АБСЦИСС: ');
% if TALPHA1<1|TALPHA1>L1|TALPHA1==inf|isreal(TALPHA1)==0
% errordlg('ПОВТОРНАЯ ОШИБКА ВЫБОРА: КООРДИНАТА TALPHA1 НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
% close all
% break
% end
% end
% TALPHA2=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: ');
% if (TALPHA2<1)|(TALPHA2>L2)|TALPHA2==inf|isreal(TALPHA2)==0
% errordlg('ОШИБКА ВЫБОРА: КООРДИНАТА (TALPHA2) НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПОВТОРИТЕ ВВОД ДЕЙСТВИТЕЛЬНОГО ЧИСЛА')
% TALPHA2=input('ВВЕДИТЕ КООРДИНАТУ (TALPHA2 - КОЭФФИЦИЕНТ ПОТЕРЬ 2-Й ПОДСИСТЕМЫ) ПО ОСИ ОРДИНАТ: ');
% if TALPHA2<1|TALPHA2>L2|TALPHA2==inf|isreal(TALPHA2)==0
% errordlg('ПОВТОРНАЯ ОШИБКА ВЫБОРА: КООРДИНАТА TALPHA2 НЕ В РАБОЧЕМ ПРЯМОУГОЛЬНИКЕ. ПРОГРАММА ПРЕРЫВАЕТСЯ')
% close all
% break
% end
% end
%--------------------------------------------------------------------------------------------
x = t/c1; y = t/c2;
%так как прямая проходит через точки (0;y) и (x;0) ее уравнение
%высчитывается по формуле Y = kX + B. B = y; k = -B/x = -y/x;
%Y = -y/x*X + y;% УРАВНЕНИЕ ПРЯМОЙ
% подставим в уравнение прямой Y = 1
% найдем Х
X1 = -x*(1 - y)/y;%пересечение с осью Y = 1
X2 = -x*(1/L2 - y)/y;%пересечение с осью Y = 1/L20
Y1 = - (x/y);
if (X1 >= 1/L1 && X1 <= 1)
XL = X1;
else
XL = 1/L1;
end
if (X2 >= 1/L2 && X2 <= 1)
XR = X2;
else
XR = 1;
end
if (Y1 > 1/L2 && Y1 < 1)
YMIN = Y1;
else
YMIN = 1/L2;
end
a = rand(1)*(XR-XL)+XL;
TALPHA1 = a * L1
b = (T - c1 * a) / c2;
TALPHA2 = b * L2
disp(' Tpoint - ВРЕМЯ ОБРАБОТКИ')
Tpoint=c1*TALPHA1/L1+c2*TALPHA2/L2
xpoint=[0 Tpoint*L1/c1];
ypoint=[Tpoint*L2/c2 0];
disp('---------------------------------------------------------------------------------------')
%----------------------------- ГРАФИЧЕСКИЕ ПОСТРОЕНИЯ --------------------------------------
figure(2),
plot(kvadratx,kvadraty,'linew',2), grid on, hold on,
plot(Tx,Ty,'k:','linew',2),
plot(Tminx,Tminy,'m:','linew',2),
plot(Tmaxx,Tmaxy,'g:','linew',2)
plot(ALPHA1t,ALPHA2t,'r','linew',2),
plot(ALPHA1,ALPHA2,'r*','linew',2),
plot(xpoint,ypoint,'k','linew',2),
plot(TALPHA1,TALPHA2,'k*','linew',2),
xlabel('\bf\alpha_1'),ylabel('\bf\alpha_2'),
legend('Volid square','T extreme','Tmin','Tmax','T','Random point','Tpoint','Input point'),
set(gcf,'color','w'),hold off