Текст программы с комментариями

Функция VvodIshodnDannih описанная ниже, обеспечивает считывание значений исходных данных из соответствующих «окон» ввода в интерфейсе, показанном выше на рис.3.2.5.

Для окон ввода (элементов типа Edit Text) предусмотрены тэги – обозначения внутри текста программы, заданные при создании визуального интерфейса. Как видно из текста функции VvodIshodnDannih, T0cold считывается, как описано в строке T0cold=str2double(get(handles.edit1,'String'));, из окна интерфейса (элемента типа Edit Text) с тэгом edit1, тэг был указан при создании визуального интерфейса в «Инспекторе свойств» (Property Inspector). Аналогично, легко определить из текста функции VvodIshodnDannih тэги остальных окон ввода (элементов типа Edit Text). Это edit1, …, edit10. Как видно, нумерация идет подряд (на рис.3.2.5. – сверху вниз по окнам ввода).

function VvodIshodnDannih(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

T0cold=str2double(get(handles.edit1,'String'));

T1cold=str2double(get(handles.edit2,'String'));

Vcold=str2double(get(handles.edit3,'String'));

Ccold=str2double(get(handles.edit4,'String'));

Ft=str2double(get(handles.edit5,'String'));

Kt=str2double(get(handles.edit6,'String'));

T0hot=str2double(get(handles.edit7,'String'));

L=str2double(get(handles.edit8,'String'));

Vhot=str2double(get(handles.edit9,'String'));

Chot=str2double(get(handles.edit10,'String'));

end

Функция VivodArgumentovNaEkran, описанная ниже, обеспечивает вывод на экран, то есть в соответствующие окна вывода (элементы типа Edit Text), которые легко определить из текста функции. Как видно из строк S=sprintf('%g',T0cold); и затем set(handles.edit1,'String',S);, T0cold выводится в окно вывода (элемент типа Edit Text) с тэгом edit1. Аналогично легко определить, что переменные выводятся в окна с тэгами edit1, …, edit10, то есть те же самые, которые служат для ввода исходной информации, и в соответствии с пояснениями (названиями вводимых параметров), написанными на панели интерфейса (рис.3.2.5) рядом с окнами ввода исходных значений.

function VivodArgumentovNaEkran(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

%vivod argumentov na ekran   

S=sprintf('%g',T0cold);

set(handles.edit1,'String',S);

S=sprintf('%g',T1cold);

set(handles.edit2,'String',S);

S=sprintf('%g',T0hot);

set(handles.edit7,'String',S);

S=sprintf('%g',Vcold);

set(handles.edit3,'String',S);

S=sprintf('%g',Ccold);

set(handles.edit4,'String',S);

S=sprintf('%g',Ft);

set(handles.edit5,'String',S);

S=sprintf('%g',Kt);

set(handles.edit6,'String',S);

S=sprintf('%g',L);

set(handles.edit8,'String',S);

S=sprintf('%g',Vhot);

set(handles.edit9,'String',S);

S=sprintf('%g',Chot);

set(handles.edit10,'String',S);

end

Функция ReshitZadachu имеет центральное (наиглавнейшее) значение в ходе решения поставленной задачи (об отыскании результатов исходя из известных начальных данных). Функция ReshitZadachu получает в качестве аргумента функции вектор VectorArg который содержит все основные исходные данные, последовательно расположенные как элементы вектора. Функция ReshitZadachu возвращает (ее возвращаемым значением является) вектор результатов VectorRes, который содержит вычисляемые результаты, которые в дальнейшем другими функциями будут выводиться в качестве результатов в окнах интерфейса (рис.3.2.5) и сохраняться в файл (ниже будет описано).

Кроме того, функция ReshitZadachu строит график в окне графика на рис.3.2.5, подписывает его, подписывает горизонтальную и вертикальную оси. 

function [VectorRes,Lmas,Tmas,n]=ReshitZadachu(VectorArg)

T0cold=VectorArg(1); T1cold=VectorArg(2); Vcold=VectorArg(3); Ccold=VectorArg(4); Ft=VectorArg(5); Kt=VectorArg(6); T0hot=VectorArg(7); L=VectorArg(8); Vhot=VectorArg(9); Chot=VectorArg(10);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function LocInt=funLocInt(T0hot)

   LocInt=Kt*(T0hot-T1cold);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Функция difurThot содержит дифференциальное уравнение (43) в конечно-разностной форме:

function dThot=difurThot(L,T0hot)

   dThot=-(Ft*funLocInt(T0hot))/(Vhot*Chot*L);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Функция uravnenijeT1hot решает с помощью функции MATLAB ode45 методом Рунге–Кутта 4 или 5 порядка уравнение в функции difurThot, описанной выше при граничных условиях [0.1, L ] по длине L змеевика, и получает массивы данных по L и температуре T, после чего выходная температура жидкости в змеевике T1hot приравнивается последнему элементу массива по Т.

function T1hot=uravnenijeT1hot(L,T0hot)       

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

  end   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function equ=equiation(T0cold,T1cold,Vcold,Ccold, Ft,Kt,T0hot,L,Vhot,Chot)

   equ=Vcold*Ccold*(T0cold-T1cold)-Vhot*Chot*(uravnenijeT1hot(L,T0hot)-T0hot);       

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podgotovka k deleniju popolam (poisk T1cold)

%granici poiska

Tleft=T0hot;

Tright=T0cold;

%shag poiska

stepT=(Tright-Tleft)/49;

x(1)=Tleft;

T1cold=x(1);

y(1)=equiation(T0cold,T1cold,Vcold,Ccold, Ft,Kt,T0hot,L,Vhot,Chot);

%poisk intervala smeni znakov

   for i=2:50

       x(i)=x(i-1)+stepT;

       T1cold=x(i);

y(i)=equiation(T0cold,T1cold,Vcold,Ccold, Ft,Kt,T0hot,L,Vhot,Chot);

       if y(i)*y(i-1)<0

           Tleft=x(i-1);

           Tright=x(i);

           break;

       end

   end

%delenije popolam -- naiti T1cold

   count=0;

   while abs(Tleft-Tright)>0.0001

       count=count+1;

       if count>100

           break;

       end

       Tmid=(Tright+Tleft)/2;

       T1cold=Tleft;

yleft=equiation(T0cold,T1cold,Vcold,Ccold, Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tmid;

ymid=equiation(T0cold,T1cold,Vcold,Ccold, Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tright;

yright=equiation(T0cold,T1cold,Vcold,Ccold, Ft,Kt,T0hot,L,Vhot,Chot);

       if yleft==0

           T1cold=Tleft;

           break;

       elseif ymid==0

           T1cold=Tmid;

           break;

       elseif yright==0

           T1cold=Tright;

           break;

       end

       if yleft*ymid<0

           Tright=Tmid;

       elseif yright*ymid<0

           Tleft=Tmid;

       end

       T1cold=Tmid;

   end

Ниже приведена программа решения рассматриваемой задачи. Последовательно вызываются ранее описанные функции. Затем строится график.

Tout_cold=T1cold;

T1hot=uravnenijeT1hot(L,T0hot);

Tout_hot=T1hot;

LocInt=-Kt*(Tout_hot-Tout_cold);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

   plot(Lmas,Tmas,'-k');

   grid;

   xlabel(' Length(m) ');

   ylabel(' T.out.#2(C) (Zmeevik) ');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

VectorRes=[Tout_cold,Tout_hot, LocInt];

end

Функция ReshitZadachuG аналогична описанной выше функции ReshitZadachu, однако в отличие от нее не строит график. Функция ReshitZadachuG предназначена для многократного вызова при расчете точек, например, для трехмерного графика (рис.3.2.5).

function [VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg)

T0cold=VectorArg(1); T1cold=VectorArg(2); Vcold=VectorArg(3); Ccold=VectorArg(4); Ft=VectorArg(5);Kt=VectorArg(6); T0hot=VectorArg(7); L=VectorArg(8); Vhot=VectorArg(9); Chot=VectorArg(10);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function LocInt=funLocInt(T0hot)

   LocInt=Kt*(T0hot-T1cold);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function dThot=difurThot(L,T0hot)

   dThot=-(Ft*funLocInt(T0hot))/(Vhot*Chot*L);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function T1hot=uravnenijeT1hot(L,T0hot)       

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

  end   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function equ=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot)

   equ=Vcold*Ccold*(T0cold-T1cold)-Vhot*Chot*(uravnenijeT1hot(L,T0hot)-T0hot);       

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podgotovka k deleniju popolam (poisk T1cold)

%granici poiska

Tleft=T0hot;

Tright=T0cold;

%shag poiska

stepT=(Tright-Tleft)/49;

x(1)=Tleft;

T1cold=x(1);

y(1)=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

%poisk intervala smeni znakov

   for i=2:50

       x(i)=x(i-1)+stepT;

       T1cold=x(i);

y(i)=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       if y(i)*y(i-1)<0

           Tleft=x(i-1);

           Tright=x(i);

           break;

       end

   end

%delenije popolam -- naiti T1cold

   count=0;

   while abs(Tleft-Tright)>0.0001

       count=count+1;

       if count>100

           break;

       end

       Tmid=(Tright+Tleft)/2;

       T1cold=Tleft;

yleft=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tmid;

ymid=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tright;

yright=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       if yleft==0

           T1cold=Tleft;

           break;

       elseif ymid==0

           T1cold=Tmid;

           break;

       elseif yright==0

           T1cold=Tright;

           break;

       end

       if yleft*ymid<0

           Tright=Tmid;

       elseif yright*ymid<0

             Tleft=Tmid;

       end

       T1cold=Tmid;

   end

Tout_cold=T1cold;

T1hot=uravnenijeT1hot(L,T0hot);

Tout_hot=T1hot;

LocInt=-Kt*(Tout_hot-Tout_cold);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

VectorRes=[Tout_cold,Tout_hot, LocInt];

end

Функция VivodResultNaEkran выводит результаты вычислений в окна вывода на панели «Результаты вычислений» (рис.3.2.5) с тэгами edit15, …, edit17, то есть пронумерованными подряд. Вывод происходит в окна в соответствии с подписями рядом с ними (рис.3.2.5) с названиями выводимых расчетных параметров (результатов вычислений). 

function VivodResultNaEkran(hObject, eventdata, handles) 

global Tout_cold Tout_hot LocInt;  

%vivod resultatov na ekran   

S=sprintf('%g',Tout_cold);

set(handles.edit15,'String',S);

S=sprintf('%g',Tout_hot);

set(handles.edit16,'String',S);

S=sprintf('%g',LocInt);

set(handles.edit17,'String',S);

end

Функция PostrGrafic строит график. Как видно из теста функции global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;, она использует имя (порядковый номер подписанный у соответствующего окна в интерфейсе) NameArg аргумента (отображаемого по горизонтальной оси), использует имя (порядковый номер подписанный у соответствующего окна в интерфейсе) NameVar зависимой переменной (отображаемой по вертикальной оси), шкалу по горизонтальной оси от значения LeftNameVar до значения RightNameVar с шагом StepGr (шаг отображается в разметке горизонтальной оси и в шаге между точками графика).

Кроме того, как видно из текста plot(k1Array,v1Array,'o-');, функция PostrGrafic строит график. Затем, как видно из текста TextVectorArg=[' T.in.#1(C) ','          ','Stream#1(kg/s) ',' C.#1(G/kg*K) ',' Square(m2) ',' Koef.Transmit ',' T.in.#2(C) ',' Length(m) ','Stream#2(kg/s) ',' C.#2(G/kg*K) ']; TextVectorRes=[' T.out.#1(C) ',' T.out.#2(C) ', ' Local.Intens'];, функция подготавливает подписи к горизонтальной оси и к вертикальной оси и размещает их в соответствующих двух массивах. В каждом из массивов подписи из одинакового количества символов, включая пробелы. Последняя подпись в каждом массиве на один символ длиннее.

Затем, как видно из текста xLlabel=' ';

for i=((NameArg-1)*15+1):(NameArg*15+1)

   xLlabel=strcat(xLlabel,TextVectorArg(i));

end

yLlabel=' ';

for i=((NameVar-1)*13+1):(NameVar*13+1)

   yLlabel=strcat(yLlabel,TextVectorRes(i));

end

xlabel(xLlabel);

ylabel(yLlabel); функция подготавливает подпись к горизонтальной оси, учитывая, что длина подписи принимается пятнадцать символов (последняя подпись в массиве делается на один знак длиннее). Функция выбирает подпись - элемент массива с номером NameArg.

Также функция PostrGrafic подготавливает подпись к вертикальной оси, учитывая, что длина подписи принимается тринадцать символов (последняя подпись в массиве делается на один знак длиннее). Функция выбирает подпись - элемент массива с номером NameVar.

Затем функция PostrGrafic подписывает оси подготовленными для этого подписями.

Как видно из текста Array_3D=0; zLlabel=0; NumberSteps_3D=0; LeftNameArg_3D=0; StepGr_3D=0; RightNameArg_3D=0; функция PostrGrafic присваивает значения равные нулю переменным, связанным с построением трехмерного графика. Это обеспечивает согласованную работу с некоторыми другими функциями (будут описаны ниже).

function PostrGrafic(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

global Lmas Tmas n;

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

if NameArg>1

   NameArg=NameArg+1;

end

VectorArg(NameArg)=LeftNameVar;

k1Array(1)=LeftNameVar;

[VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg);

v1Array(1)=VectorRes(NameVar);

for i=2:NumberSteps+1

   k1Array(i)=k1Array(i-1)+StepGr;

   VectorArg(NameArg)=k1Array(i);

   [VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg);

   v1Array(i)=VectorRes(NameVar);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%stroim grafic

plot(k1Array,v1Array,'o-');

set(gca,'XGrid','on');

set(gca,'YGrid','on'); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podpisi k osyam

TextVectorArg=[' T.in.#1(C) ','          ','Stream#1(kg/s) ',' C.#1(G/kg*K) ',' Square(m2) ',' Koef.Transmit ',' T.in.#2(C) ',' Length(m) ','Stream#2(kg/s) ',' C.#2(G/kg*K) '];

TextVectorRes=[' T.out.#1(C) ',' T.out.#2(C) ', ' Local.Intens'];

xLlabel=' ';

for i=((NameArg-1)*15+1):(NameArg*15+1)

   xLlabel=strcat(xLlabel,TextVectorArg(i));

end

yLlabel=' ';

for i=((NameVar-1)*13+1):(NameVar*13+1)

   yLlabel=strcat(yLlabel,TextVectorRes(i));

end

xlabel(xLlabel);

ylabel(yLlabel);

Array_3D=0;

zLlabel=0;

NumberSteps_3D=0;

LeftNameArg_3D=0;

StepGr_3D=0;

RightNameArg_3D=0;

end

Функция PostrGrafic_3D строит трехмерный график. Функция использует для третьей оси (второй оси в горизонтальной плоскости) значение NameArg_3D как имя (номер указанный в подписях к элементам интерфейса) переменной по третьей оси (второй аргумент при построении графика зависимости от двух аргументов), шкалу по третьей оси от LeftNameArg_3D до RightNameArg_3D с шагом StepGr_3D. В конце работы функция присваивает v1Array=0; для совместимости с некоторыми другими функциями (будут описаны ниже)

function PostrGrafic_3D(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

global Lmas Tmas n;

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

if NameArg>1

   NameArg=NameArg+1;

end

if NameArg_3D>1

   NameArg_3D=NameArg_3D+1;

end

for i=1:NumberSteps+1

   if i==1

       k1Array(i)=LeftNameVar; 

   else

       k1Array(i)=k1Array(i-1)+StepGr;           

   end %if

   VectorArg(NameArg)=k1Array(i);

   for j=1:NumberSteps_3D+1

       if j==1

           z1Array(j)=LeftNameArg_3D; 

       else

           z1Array(j)=z1Array(j-1)+StepGr_3D; 

       end %if

       VectorArg(NameArg_3D)=z1Array(j);

       [VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg);

       Array_3D(j,i)=VectorRes(NameVar);

   end %for j

end %for i

[x,y]=meshgrid(LeftNameVar:StepGr:RightNameVar,LeftNameArg_3D:StepGr_3D:RightNameArg_3D);

%x %число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

%y %число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

%Array_3D %число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

mesh(x,y,Array_3D);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podpisi k osyam

TextVectorArg=[' T.in.#1(C) ','          ','Stream#1(kg/s) ',' C.#1(G/kg*K) ',' Square(m2) ',' Koef.Transmit ',' T.in.#2(C) ',' Length(m) ','Stream#2(kg/s) ',' C.#2(G/kg*K) '];

TextVectorRes=[' T.out.#1(C) ',' T.out.#2(C) ', ' Local.Intens'];

xLlabel=' ';

for i=((NameArg-1)*15+1):(NameArg*15+1)

   xLlabel=strcat(xLlabel,TextVectorArg(i));

end

yLlabel=' ';

for i=((NameArg_3D-1)*15+1):(NameArg_3D*15+1)

   yLlabel=strcat(yLlabel,TextVectorArg(i));

end

zLlabel=' ';

for i=((NameVar-1)*13+1):(NameVar*13+1)

   zLlabel=strcat(zLlabel,TextVectorRes(i));

end

xlabel(xLlabel);

ylabel(yLlabel);

zlabel(zLlabel);

v1Array=0;

end

Функция pushbutton13_Callback вызывается при нажатии кнопки интерфейса с надписью «ПОСТРОИТЬ ГРАФИК». Необходимо присвоить этой кнопке тэг pushbutton13 при создании интерфейса. Функция считывает значения переменных из окон интерфейса. NameVar из окна «номер зависимой переменной» с тэгом edit31, LeftNameVar из окна «минимальное значение аргумента» с тэгом edit29, RightNameVar из окна «максимальное значение аргумента» с тэгом edit30, NameArg из окна «номер аргумента» с тэгом edit28, StepGr из окна «шаг аргумента» с тэгом edit32, а также с панели для построения трехмерного графика «Объемный график» считывает NameArg_3D из окна «номер аргумента» с тэгом edit36, LeftNameArg_3D из окна «минимальный» с тэгом edit35, RightNameArg_3D из окна «максимальный» с тэгом edit33, StepGr_3D из окна «шаг» с тэгом edit34. Необходимо присвоить этим окнам соответствующие тэги.

Если NameArg_3D равен нулю, то вызывается функция PostrGrafic для построения плоского графика. По умолчанию в окне панели для построения трехмерного графика «Объемный график» NameArg_3D «номер аргумента» (рис.3.1.1) стоит ноль. Его надо указать при создании интерфейса как значение ‘String’.

Если NameArg_3D не равен нулю, то вызывается функция PostrGrafic_3D для построения трехмерного графика.

function pushbutton13_Callback(hObject, eventdata, handles)

global PG;

if PG==1       

else

   clear global;

end %if

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

NameVar=str2double(get(handles.edit31,'String'));

LeftNameVar=str2double(get(handles.edit29,'String'));

RightNameVar=str2double(get(handles.edit30,'String'));

NameArg=str2double(get(handles.edit28,'String'));

StepGr=str2double(get(handles.edit32,'String'));

VvodIshodnDannih(hObject, eventdata, handles);

NumberSteps=round((RightNameVar-LeftNameVar)/StepGr);

NameArg_3D=str2double(get(handles.edit36,'String'));

if NameArg_3D==0

   PostrGrafic(hObject, eventdata, handles);

else

LeftNameArg_3D=str2double(get(handles.edit35,'String'));

RightNameArg_3D=str2double(get(handles.edit33,'String'));

   StepGr_3D=str2double(get(handles.edit34,'String'));

   NumberSteps_3D=round((RightNameArg_3D-LeftNameArg_3D)/StepGr_3D);

   PostrGrafic_3D(hObject, eventdata, handles);

end %if

end

Функция pushbutton3_Callback вызывается кнопкой «ВЫЧИСЛИТЬ». Необходимо этой кнопке присвоить тэг pushbutton3. Она обеспечивает вызов функции VvodIshodnDannih то есть ввод исходных данных из окон интерфейса, вызов функции ReshitZadachu то есть вычисление, и вызов функции VivodResultNaEkran то есть вывод результатов в окна интерфейса.

function pushbutton3_Callback(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

VvodIshodnDannih(hObject, eventdata, handles);

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

[VectorRes,Lmas,Tmas,n]=ReshitZadachu(VectorArg);

Tout_cold=VectorRes(1); Tout_hot=VectorRes(2); LocInt=VectorRes(3);

VivodResultNaEkran(hObject, eventdata, handles);

end

Функция pushbutton2_Callback вызывается кнопкой «ВЫЧИСЛИТЬ И СОХРАНИТЬ», необходимо присвоить этой кнопке тэг pushbutton2, она обеспечивает вычисление результатов и сохранение в файл.

function pushbutton2_Callback(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

%Vichislit resultati

pushbutton3_Callback(hObject, eventdata, handles);

%otkrit fail dlya zapisi

[f,p]=uiputfile('C:\MATLAB701\work\S292Res.txt','Окно выбора');

KudaZapisat=strcat(p,f);

f=fopen(KudaZapisat,'wt');

%gotovim dannije dlja avtomaticheskogo schitivanija

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

VectorRes=[Tout_cold,Tout_hot, LocInt];

fprintf(f,'%d\n',n);

for i=1:n

   fprintf(f,'%g\t',Tmas(i));           

end

for i=1:n

   fprintf(f,'%g\t',Lmas(i));    

end   

fprintf(f,'\n');

for i=1:3

   fprintf(f,'%g\t',VectorRes(i));

end

for i=1:10

   fprintf(f,'%g\t',VectorArg(i));

end       

%sohranit v fail resultati vichislenij

fprintf(f,'\nS292_TEPLOOBM_SMESHENIYE_VITESNENIYE_02\n');

fprintf(f,'RESULTATI VICHISLENIY\n');

TextString=strcat('Tout#1(C)\t','Tout2(C)\t','LocInt\n');

fprintf(f,TextString);

fprintf(f,'%g\t\t',VectorRes(1));

fprintf(f,'%g\t\t',VectorRes(2));

fprintf(f,'%g\n',VectorRes(3));

%sohranit v fail ishodnije argumenti

fprintf(f,'ARGUMENTI VICHISLENIY\n');

TextString=strcat(' T.in.#1(C) ','Predpolag Tout#1','Stream#1(kg/s) ',' C.#1(G/kg*K) ',' Square(m2) ',' Koef.Transmit ',' T.in.#2(C) ',' Length(m) ','Stream#2(kg/s) ',' C.#2(G/kg*K) ');

fprintf(f,TextString);

for i=1:10

   fprintf(f,'%g\t\t\t',VectorArg(i));

end

fclose(f);

end

Функция pushbutton14_Callback вызывается кнопкой «ОТКРЫТЬ ИЗ ФАЙЛА», необходимо этой кнопке присвоить тэг pushbutton14, функция обеспечивает открытие из файла исходных данных и результатов и вывод их в окна интерфейса.

function pushbutton14_Callback(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

%Otkrit is faila

[f,p]=uigetfile('C:\MATLAB701\work\S292Res*.txt','Окно выбора');

OtkudaChitat=strcat(p,f);

f=fopen(OtkudaChitat,'rt');

%chitat is faila

%chitat matricu T(n,m)

n=fscanf(f,'%d',1);

for i=1:n

   Tmas(i)=fscanf(f,'%g',1);       

end

%chitat stolbec L(n)

for i=1:n

   Lmas(i)=fscanf(f,'%g',1);

end   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for i=1:3

   VectorRes(i)=fscanf(f,'%g',1);

end

for i=1:10

   VectorArg(i)=fscanf(f,'%g',1);

end

fclose(f);

%vivod na ekran

T0cold=VectorArg(1); T1cold=VectorArg(2); T0hot=VectorArg(7); Vcold=VectorArg(3); Ccold=VectorArg(4); Ft=VectorArg(5); Kt=VectorArg(6); L=VectorArg(8); Vhot=VectorArg(9); Chot=VectorArg(10);

Tout_cold=VectorRes(1); Tout_hot=VectorRes(2); LocInt=VectorRes(3);

VivodResultNaEkran(hObject, eventdata, handles);

VivodArgumentovNaEkran(hObject, eventdata, handles);  

%postroit grafic

plot(Lmas,Tmas,'-k');

grid;

xlabel('L (m)');

ylabel('T (C)');  

end

Функция pushbutton8_Callback вызывается при нажатии на кнопку «ПОСТРОИТЬ ГРАФИК И СОХРАНИТЬ», необходимо этой кнопке присвоить тэг pushbutton8, функция обеспечивает построение графика и его сохранение в файл.

function pushbutton8_Callback(hObject, eventdata, handles)

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

global PG;

%postroit grafic

PG=1;

pushbutton13_Callback(hObject, eventdata, handles);

PG=0;

%Sohranit grafic v fail

[f,p]=uiputfile('C:\MATLAB701\work\S292Graf.txt','Окно выбора');

KudaZapisat=strcat(p,f);

f=fopen(KudaZapisat,'wt');

%%%%%%%%%%%%%%%%%%%

if Array_3D==0

  %gotovim dannije dlja avtomaticheskogo chitivanija

  fprintf(f,'%d\t',NumberSteps);

  for i=1:(NumberSteps+1)

      fprintf(f,'%g\t',k1Array(i));

  end

  fprintf(f,'\n');

  for i=1:(NumberSteps+1)

      fprintf(f,'%g\t',v1Array(i));

  end       

  %podpisi k osyam grafica

  fprintf(f,'%s\t%s\n',xLlabel,yLlabel);

else

  %gotovim dannije dlja avtomaticheskogo chitivanija

  NS=0;

  fprintf(f,'%d\n',NS);

fprintf(f,'%d\t%d\t%d\t%d\t%d\t%d\n',LeftNameVar,StepGr,RightNameVar,LeftNameArg_3D,StepGr_3D,RightNameArg_3D);

  fprintf(f,'\n');

  for i=1:(NumberSteps_3D+1)

      for j=1:(NumberSteps+1)           

         fprintf(f,'%g\t',Array_3D(i,j));

      end %for j

     fprintf(f,'\n');

  end

  fprintf(f,'\n');

  %podpisi k osyam grafica

  fprintf(f,'%s\t%s\t%s\n',xLlabel,yLlabel,zLlabel); 

end %if

%zapisivajem argumenti

VvodIshodnDannih(hObject, eventdata, handles);

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

for i=1:10

  fprintf(f,'%g\t',VectorArg(i));

end 

%podpisi k osyam grafica

fclose(f);

end

Функция pushbutton15_Callback вызывается кнопкой «ОТКРЫТЬ ГРАФИК ИЗ ФАЙЛА», необходимо присвоить этой кнопке тэг pushbutton15, функция обеспечивает открытие графика из файла (то есть построение графика) и выводит в окна интерфейса значения исходных данных.

function pushbutton15_Callback(hObject, eventdata, handles)

clear global;

global k1Array v1Array xLlabel yLlabel;

global LeftNameVar NumberSteps StepGr NameVar NameArg;

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

%otkrit grafic is faila

[f,p]=uigetfile('C:\MATLAB701\work\S292Graf*.txt','Окно выбора');

OtkudaChitat=strcat(p,f);

f=fopen(OtkudaChitat,'rt');

%chitat is faila

NumberSteps=fscanf(f,'%d',1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if NumberSteps==0

  LeftNameVar=fscanf(f,'%g',1);

  StepGr=fscanf(f,'%g',1);

  RightNameVar=fscanf(f,'%g',1);

  LeftNameArg_3D=fscanf(f,'%g',1);

  StepGr_3D=fscanf(f,'%g',1);

  RightNameArg_3D=fscanf(f,'%g',1);

[x,y]=meshgrid(LeftNameVar:StepGr:RightNameVar,LeftNameArg_3D:StepGr_3D:RightNameArg_3D);

%число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

   NumberSteps_3D=round((RightNameArg_3D-LeftNameArg_3D)/StepGr_3D);

   NumberSteps=round((RightNameVar-LeftNameVar)/StepGr);

   for i=1:(NumberSteps_3D+1)

       for j=1:(NumberSteps+1)           

          Array_3D(i,j)=fscanf(f,'%g',1);

       end %for j

   end % for i

   %podpisi k osyam grafica

   xLlabel=fscanf(f,'%s',1);

   yLlabel=fscanf(f,'%s',1);

   zLlabel=fscanf(f,'%s',1);

   %stroim grafic

   mesh(x,y,Array_3D);

   xlabel(xLlabel);

   ylabel(yLlabel);

   zlabel(zLlabel);

else

   for i=1:NumberSteps+1

       k1Array(i)=fscanf(f,'%g',1);

   end

   for i=1:NumberSteps+1

       v1Array(i)=fscanf(f,'%g',1);

   end

   %podpisi k osyam grafica

   xLlabel=fscanf(f,'%s',1);

   yLlabel=fscanf(f,'%s',1);

   %stroim grafic

   plot(k1Array,v1Array,'o-');

   set(gca,'XGrid','on');

   set(gca,'YGrid','on'); 

   xlabel(xLlabel);

   ylabel(yLlabel);

end %if

for i=1:10

  VectorArg(i)=fscanf(f,'%g',1);

end

fclose(f);

%vivod na ekran

T0cold=VectorArg(1); T1cold=VectorArg(2); T0hot=VectorArg(7); Vcold=VectorArg(3); Ccold=VectorArg(4); Ft=VectorArg(5); Kt=VectorArg(6); L=VectorArg(8); Vhot=VectorArg(9); Chot=VectorArg(10);

 %vivod na ekran

VivodArgumentovNaEkran(hObject, eventdata, handles); 

end

Для начинающих освоение МАТЛАБа, можно рекомендовать такую же по содержанию программу, но без визуального интерфейса.

Ниже изложено ее содержание.

Для наиболее начальной стадии изучения МАТЛАБа можно рекомендовать ее выполнить сначала не всю, а только следующие ее части:

ИСХОДНЫЕ ДАННЫЕ

ВЫЧИСЛИТЬ

ВЫВОД РЕЗУЛЬТАТОВ

Ниже будет пояснено, какие несколько функций для этого надо подготовить.

ИСХОДНЫЕ ДАННЫЕ:

Функция DATA, необходимая для задания исходных данных, так как графического интерфейса с окнами ввода данных не предусмотрено. Данные помещаются в файл DATA, соответствующие переменные объявляются глобальными. В дальнейшем именно к файлу DATA будут обращаться различные функции, нуждающиеся в исходных данных, вместо прежнего, описанного выше, обращения к функции ввода данных из окон графического интерфейса.

function DATA

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

T0cold=500; % T nachaln №1 (C)

T1cold=100; % predpolagaem T konechn №1 (C)

Vcold=1; % rashod №1 (kg/s)

Ccold=10000; % teploemkost №1 (J/(kg*K))

Ft=10; % poverhn teploobm (m2)

Kt=100; % koef teploperedachi

T0hot=100;% T nachaln №2 (C)

L=100; % dlina zmeevika (m)

Vhot=1;% rashod №2 (kg/s)

Chot=10000;% teploemkost №2 (J/(kg*K))

end

ВЫЧИСЛИТЬ:

function [VectorRes,Lmas,Tmas,n]=ReshitZadachu(VectorArg)

T0cold=VectorArg(1); T1cold=VectorArg(2); Vcold=VectorArg(3); Ccold=VectorArg(4); Ft=VectorArg(5); Kt=VectorArg(6); T0hot=VectorArg(7); L=VectorArg(8); Vhot=VectorArg(9); Chot=VectorArg(10);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function LocInt=funLocInt(T0hot)

   LocInt=Kt*(T0hot-T1cold);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function dThot=difurThot(L,T0hot)

   dThot=-(Ft*funLocInt(T0hot))/(Vhot*Chot*L);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function T1hot=uravnenijeT1hot(L,T0hot)       

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

  end   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function equ=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot)

   equ=Vcold*Ccold*(T0cold-T1cold)-Vhot*Chot*(uravnenijeT1hot(L,T0hot)-T0hot);       

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podgotovka k deleniju popolam (poisk T1cold)

%granici poiska

Tleft=T0hot;

Tright=T0cold;

%shag poiska

stepT=(Tright-Tleft)/49;

x(1)=Tleft;

T1cold=x(1);

y(1)=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

%poisk intervala smeni znakov

   for i=2:50

       x(i)=x(i-1)+stepT;

       T1cold=x(i);

y(i)=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       if y(i)*y(i-1)<0

           Tleft=x(i-1);

           Tright=x(i);

           break;

       end

   end

%delenije popolam -- naiti T1cold

   count=0;

   while abs(Tleft-Tright)>0.0001

       count=count+1;

       if count>100

           break;

       end

       Tmid=(Tright+Tleft)/2;

       T1cold=Tleft;

yleft=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tmid;

ymid=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tright;

yright=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       if yleft==0

           T1cold=Tleft;

           break;

       elseif ymid==0

           T1cold=Tmid;

           break;

       elseif yright==0

           T1cold=Tright;

           break;

       end

       if yleft*ymid<0

           Tright=Tmid;

       elseif yright*ymid<0

           Tleft=Tmid;

       end

       T1cold=Tmid;

   end

Tout_cold=T1cold;

T1hot=uravnenijeT1hot(L,T0hot);

Tout_hot=T1hot;

LocInt=-Kt*(Tout_hot-Tout_cold);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

   plot(Lmas,Tmas,'-k');

   grid;

   xlabel(' Length(m) ');

   ylabel(' T.out.#2(C) (Zmeevik) ');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

VectorRes=[Tout_cold,Tout_hot, LocInt];

end

Функция GLAV_Vichislit_pushbutton3_Callback не связана ни с какими элементами (кнопками) графического интерфейса, так как графического интерфейса нет, но выполняет те же самые задачи, что и раньше (как было описано в предыдущем варианте программы, где был графический интерфейс). Однако, вместо обращения к функции ввода данных из окон графического интерфейса, теперь функция GLAV_Vichislit_pushbutton3_Callback обращается к функции DATA, чтобы получить исходные данные. Объявленные в функции GLAV_Vichislit_pushbutton3_Callback глобальные переменные в результате обращения к функции DATA становятся содержащими те же значения, как одинаковые с ними по названию глобальные переменные в функции DATA.

function GLAV_Vichislit_pushbutton3_Callback

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

DATA;

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

[VectorRes,Lmas,Tmas,n]=ReshitZadachu(VectorArg);

Tout_cold=VectorRes(1); Tout_hot=VectorRes(2); LocInt=VectorRes(3);

VivodResultNaEkran;

end

ВЫВОД РЕЗУЛЬТАТОВ:

Функция VivodResultNaEkran выводит результаты в Command Window (так как графического интерфейса не предусмотрено).

function VivodResultNaEkran

global Tout_cold Tout_hot LocInt;

disp('РЕЗУЛЬТАТЫ ВЫЧИСЛЕНИЙ');

disp(['Температура выходная потока №1 (Tout_cold)=' num2str(Tout_cold,'%10.2f') ' С']);

disp(['Температура выходная потока №2 (Tout_hot)=' num2str(Tout_hot,'%10.2f') ' С']);

disp(['Локальная интенсивность теплопередачи (LocInt)=' num2str(LocInt,'%10.2f') ' ']);

end

Функция VivodArgumentovNaEkran выводит значения исходных данных в Command Window (так как графического интерфейса не предусмотрено).

function VivodArgumentovNaEkran

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

%vivod argumentov na ekran

disp('ИСХОДНЫЕ ДАННЫЕ');

disp(['Температура начальная потока №1 (T0cold)=' num2str(T0cold,'%10.2f') ' С']);

disp(['Предполагаемая температура конечная потока №1 (T1cold)=' num2str(T1cold,'%10.2f') ' С']);

disp(['Расход потока №1 (Vcold)=' num2str(Vcold,'%10.2f') ' кг/с']);

disp(['Теплоемкость потока №1 (Ccold)=' num2str(Ccold,'%10.2f') ' Дж/(кг*град)']);

disp(['Поверхность теплообмена (Ft)=' num2str(Ft,'%10.2f') ' м2']);

disp(['Коэф теплопередачи (Kt)=' num2str(Kt,'%10.2f') ' ']);

disp(['Температура начальная потока №2 (T0hot)=' num2str(T0hot,'%10.2f') ' С']);

disp(['Длина змеевика (L)=' num2str(L,'%10.2f') ' м']);

disp(['Расход потока №2 (Vhot)=' num2str(Vhot,'%10.2f') ' кг/с']);

disp(['Теплоемкость потока №2 (Chot)=' num2str(Chot,'%10.2f') ' Дж/(кг*град)']);

end

Для более подробного изучения, следует подготовить нижеследующие функции. Они позволят построить графики – плоские и трехмерные, а также сохранить результаты вычислений и графики в файл и открыть их из файла. Все они не связаны с графическим интерфейсом и для получения исходных данных обращаются к функции DATA, а для вывода значений аргументов и вычисленных результатов – к вышеописанным функциям VivodArgumentovNaEkran и VivodResultNaEkran, выводящим в Command Window.

function [VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg)

T0cold=VectorArg(1); T1cold=VectorArg(2); Vcold=VectorArg(3); Ccold=VectorArg(4); Ft=VectorArg(5); Kt=VectorArg(6); T0hot=VectorArg(7); L=VectorArg(8); Vhot=VectorArg(9); Chot=VectorArg(10);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function LocInt=funLocInt(T0hot)

   LocInt=Kt*(T0hot-T1cold);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function dThot=difurThot(L,T0hot)

   dThot=-(Ft*funLocInt(T0hot))/(Vhot*Chot*L);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function T1hot=uravnenijeT1hot(L,T0hot)       

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

  end   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function equ=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot)

   equ=Vcold*Ccold*(T0cold-T1cold)-Vhot*Chot*(uravnenijeT1hot(L,T0hot)-T0hot);       

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podgotovka k deleniju popolam (poisk T1cold)

%granici poiska

Tleft=T0hot;

Tright=T0cold;

%shag poiska

stepT=(Tright-Tleft)/49;

x(1)=Tleft;

T1cold=x(1);

y(1)=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

%poisk intervala smeni znakov

   for i=2:50

       x(i)=x(i-1)+stepT;

       T1cold=x(i);

y(i)=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       if y(i)*y(i-1)<0

           Tleft=x(i-1);

           Tright=x(i);

           break;

       end

   end

%delenije popolam -- naiti T1cold

   count=0;

   while abs(Tleft-Tright)>0.0001

       count=count+1;

       if count>100

           break;

       end

       Tmid=(Tright+Tleft)/2;

       T1cold=Tleft;

yleft=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tmid;

ymid=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       T1cold=Tright;

yright=equiation(T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot);

       if yleft==0

           T1cold=Tleft;

           break;

       elseif ymid==0

           T1cold=Tmid;

           break;

       elseif yright==0

           T1cold=Tright;

           break;

       end

       if yleft*ymid<0

           Tright=Tmid;

       elseif yright*ymid<0

           Tleft=Tmid;

       end

       T1cold=Tmid;

   end

Tout_cold=T1cold;

T1hot=uravnenijeT1hot(L,T0hot);

Tout_hot=T1hot;

LocInt=-Kt*(Tout_hot-Tout_cold);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   [Lmas,Tmas]=ode45(@difurThot,[0.1 L], T0hot);

   n=length(Tmas);

   T1hot=Tmas(n);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

VectorRes=[Tout_cold,Tout_hot, LocInt];

end

 

function PostrGrafic_3D

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

global Lmas Tmas n;

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

if NameArg>1

   NameArg=NameArg+1;

end

if NameArg_3D>1

   NameArg_3D=NameArg_3D+1;

end

for i=1:NumberSteps+1

   if i==1

       k1Array(i)=LeftNameVar; 

   else

       k1Array(i)=k1Array(i-1)+StepGr;           

   end %if

   VectorArg(NameArg)=k1Array(i);

   for j=1:NumberSteps_3D+1

       if j==1

           z1Array(j)=LeftNameArg_3D; 

       else

           z1Array(j)=z1Array(j-1)+StepGr_3D;               

       end %if

       VectorArg(NameArg_3D)=z1Array(j);

       [VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg);

       Array_3D(j,i)=VectorRes(NameVar);

   end %for j

end %for i

[x,y]=meshgrid(LeftNameVar:StepGr:RightNameVar,LeftNameArg_3D:StepGr_3D:RightNameArg_3D);

%x %число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

%y %число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

%Array_3D %число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

mesh(x,y,Array_3D);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podpisi k osyam

TextVectorArg=[' T.in.#1(C) ','          ','Stream#1(kg/s) ',' C.#1(G/kg*K) ',' Square(m2) ',' Koef.Transmit ',' T.in.#2(C) ',' Length(m) ','Stream#2(kg/s) ',' C.#2(G/kg*K) '];

TextVectorRes=[' T.out.#1(C) ',' T.out.#2(C) ', ' Local.Intens'];

xLlabel=' ';

for i=((NameArg-1)*15+1):(NameArg*15+1)

   xLlabel=strcat(xLlabel,TextVectorArg(i));

end

yLlabel=' ';

for i=((NameArg_3D-1)*15+1):(NameArg_3D*15+1)

   yLlabel=strcat(yLlabel,TextVectorArg(i));

end

zLlabel=' ';

for i=((NameVar-1)*13+1):(NameVar*13+1)

   zLlabel=strcat(zLlabel,TextVectorRes(i));

end

xlabel(xLlabel);

ylabel(yLlabel);

zlabel(zLlabel);

v1Array=0;

end

 

function PostrGrafic

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

global Lmas Tmas n;

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

if NameArg>1

   NameArg=NameArg+1;

end

VectorArg(NameArg)=LeftNameVar;

 k1Array(1)=LeftNameVar;

[VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg);

v1Array(1)=VectorRes(NameVar);

for i=2:NumberSteps+1

   k1Array(i)=k1Array(i-1)+StepGr;

   VectorArg(NameArg)=k1Array(i);

   [VectorRes,Lmas,Tmas,n]=ReshitZadachuG(VectorArg);

   v1Array(i)=VectorRes(NameVar);

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%stroim grafic

plot(k1Array,v1Array,'o-');

set(gca,'XGrid','on');

set(gca,'YGrid','on'); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%podpisi k osyam

TextVectorArg=[' T.in.#1(C) ','          ','Stream#1(kg/s) ',' C.#1(G/kg*K) ',' Square(m2) ',' Koef.Transmit ',' T.in.#2(C) ',' Length(m) ','Stream#2(kg/s) ',' C.#2(G/kg*K) '];

TextVectorRes=[' T.out.#1(C) ',' T.out.#2(C) ', ' Local.Intens'];

xLlabel=' ';

for i=((NameArg-1)*15+1):(NameArg*15+1)

   xLlabel=strcat(xLlabel,TextVectorArg(i));

end

yLlabel=' ';

for i=((NameVar-1)*13+1):(NameVar*13+1)

   yLlabel=strcat(yLlabel,TextVectorRes(i));

end

xlabel(xLlabel);

ylabel(yLlabel);

Array_3D=0;

zLlabel=0;

NumberSteps_3D=0;

LeftNameArg_3D=0;

StepGr_3D=0;

RightNameArg_3D=0;

end

 

function GLAV_Vichislit_Sohranit_pushbutton2_Callback

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

%Vichislit resultati

GLAV_Vichislit_pushbutton3_Callback;

%otkrit fail dlya zapisi

[f,p]=uiputfile('C:\MATLAB701\work\S292Res.txt','Окно выбора');

KudaZapisat=strcat(p,f);

f=fopen(KudaZapisat,'wt');

%gotovim dannije dlja avtomaticheskogo schitivanija

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

VectorRes=[Tout_cold,Tout_hot, LocInt];

fprintf(f,'%d\n',n);

for i=1:n

   fprintf(f,'%g\t',Tmas(i));           

end

for i=1:n

   fprintf(f,'%g\t',Lmas(i));    

end   

fprintf(f,'\n');

for i=1:3

   fprintf(f,'%g\t',VectorRes(i));

end

for i=1:10

   fprintf(f,'%g\t',VectorArg(i));

end       

%sohranit v fail resultati vichislenij

fprintf(f,'\nS292_TEPLOOBM_SMESHENIYE_VITESNENIYE_02\n');

fprintf(f,'RESULTATI VICHISLENIY\n');

TextString=strcat('Tout#1(C)\t','Tout2(C)\t','LocInt\n');

fprintf(f,TextString);

fprintf(f,'%g\t\t',VectorRes(1));

fprintf(f,'%g\t\t',VectorRes(2));

fprintf(f,'%g\n',VectorRes(3));

%sohranit v fail ishodnije argumenti

fprintf(f,'ARGUMENTI VICHISLENIY\n');

TextString=strcat(' T.in.#1(C) ','Predpolag Tout#1','Stream#1(kg/s) ',' C.#1(G/kg*K) ',' Square(m2) ',' Koef.Transmit ',' T.in.#2(C) ',' Length(m) ','Stream#2(kg/s) ',' C.#2(G/kg*K) ');

fprintf(f,TextString);

for i=1:10

   fprintf(f,'%g\t\t\t',VectorArg(i));

end

fclose(f);

end

 

function GLAV_PostrGrafic_Sohranit_pushbutton8_Callback

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

  global PG;

%postroit grafic

PG=1;

GLAV_PostrGrafic_pushbutton13_Callback;

PG=0;

%Sohranit grafic v fail

[f,p]=uiputfile('C:\MATLAB701\work\S292Graf.txt','Окно выбора');

KudaZapisat=strcat(p,f);

f=fopen(KudaZapisat,'wt');

%%%%%%%%%%%%%%%%%%%

if Array_3D==0

  %gotovim dannije dlja avtomaticheskogo chitivanija

  fprintf(f,'%d\t',NumberSteps);

  for i=1:(NumberSteps+1)

      fprintf(f,'%g\t',k1Array(i));

  end

  fprintf(f,'\n');

  for i=1:(NumberSteps+1)

      fprintf(f,'%g\t',v1Array(i));

  end       

  %podpisi k osyam grafica

  fprintf(f,'%s\t%s\n',xLlabel,yLlabel);

else

  %gotovim dannije dlja avtomaticheskogo chitivanija

  NS=0;

  fprintf(f,'%d\n',NS);

fprintf(f,'%d\t%d\t%d\t%d\t%d\t%d\n',LeftNameVar,StepGr,RightNameVar,LeftNameArg_3D,StepGr_3D,RightNameArg_3D);

  fprintf(f,'\n');

  for i=1:(NumberSteps_3D+1)

      for j=1:(NumberSteps+1)           

         fprintf(f,'%g\t',Array_3D(i,j));

      end %for j

      fprintf(f,'\n');

  end

  fprintf(f,'\n');

  %podpisi k osyam grafica

  fprintf(f,'%s\t%s\t%s\n',xLlabel,yLlabel,zLlabel);       

end %if

%zapisivajem argumenti

DATA;

VectorArg=[T0cold,T1cold,Vcold,Ccold,Ft,Kt,T0hot,L,Vhot,Chot];

for i=1:10

  fprintf(f,'%g\t',VectorArg(i));

end 

fclose(f);

end

 

function GLAV_PostrGrafic_pushbutton13_Callback

global PG;

if PG==1       

else

   clear global;

end %if

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

global NameVar LeftNameVar RightNameVar NameArg StepGr LeftNameArg_3D StepGr_3D RightNameArg_3D NameArg_3D;

global k1Array v1Array Array_3D xLlabel yLlabel zLlabel NumberSteps NumberSteps_3D;

DATA_Graf;

DATA;

NumberSteps=round((RightNameVar-LeftNameVar)/StepGr);

if NameArg_3D==0       

   PostrGrafic;

else

   NumberSteps_3D=round((RightNameArg_3D-LeftNameArg_3D)/StepGr_3D);

   PostrGrafic_3D;

end %if

end

 

function GLAV_Otkrit_Grafic_pushbutton15_Callback

clear global;

global k1Array v1Array xLlabel yLlabel;

global LeftNameVar NumberSteps StepGr NameVar NameArg;

global T0cold T1cold Vcold Ccold Ft Kt T0hot L Vhot Chot;

global Tout_cold Tout_hot LocInt;

global Lmas Tmas n;

%otkrit grafic is faila

[f,p]=uigetfile('C:\MATLAB701\work\S292Graf*.txt','Окно выбора');

OtkudaChitat=strcat(p,f);

f=fopen(OtkudaChitat,'rt');

%chitat is faila

NumberSteps=fscanf(f,'%d',1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if NumberSteps==0

  LeftNameVar=fscanf(f,'%g',1);

  StepGr=fscanf(f,'%g',1);

  RightNameVar=fscanf(f,'%g',1);

  LeftNameArg_3D=fscanf(f,'%g',1);

  StepGr_3D=fscanf(f,'%g',1);

  RightNameArg_3D=fscanf(f,'%g',1);

[x,y]=meshgrid(LeftNameVar:StepGr:RightNameVar,LeftNameArg_3D:StepGr_3D:RightNameArg_3D);

%число столбцов=NumberSteps+1 число строк=NumberSteps_3D+1

   NumberSteps_3D=round((RightNameArg_3D-LeftNameArg_3D)/StepGr_3D);

   NumberSteps=round((RightNameVar-LeftNameVar)/StepGr);

   for i=1:(NumberSteps_3D+1)

       for j=1:(NumberSteps+1)           

          Array_3D(i,j)=fscanf(f,'%g',1);

       end %for j

   end % for i