Labview and Matlab for PT

管路中之損失除管路本身外，其出口、入口及管徑擴張之損失均與動能有關，即V²/(2g)之係數。其程式內容如下

% pipefittings.m
% Local losses at pipe fittings
%
while 1
    k=menu('請輸入接合件之型式','驟然放大','驟然逸出','擴 大 管');
    Q=input('請輸入管路流量(m^3/s)[10]: ','s');
    if isempty(Q), Q='10';end
    D1=input('請輸入入口管徑(m)[1]: ','s');
    if isempty(D1), D1='1';end
    d1=str2num(D1);
    ke=str2num(Q)./(pi*d1.^2/4)/(2*9.81);
    switch k
        case 1
            HL=ke;
        case 2
            HL=ke/2;
        case 3
           D2=input('請輸入出口管徑(m)[2]: ','s');
           if isempty(D2), D2='2';end
           HL=(1-(d1./str2num(D2)).^2).*ke;
    end
    disp(['此接頭損失為 ' num2str(HL) ' (m).'])
    A=input('是否繼續執行(Y/N)[Y]:','s');
    if isempty(A), A='Y';end
    if upper(A)~='Y', break;end
    disp('請由選單選擇轉換單位。。。')
end

管路配件損頭

經過認養後

% pipefittings.m
% Local losses at pipe fittings
%
while 1
    k=menu('請輸入接合件之型式','驟然放大','驟然逸出','擴 大 管');
    switch k
        case 1
            prompt = {'請輸入管路流量(m^3/s)[10]','請輸入入口管徑(m)[1]:'};
            dlg_title = '您所選擇的是驟然放大';
            num_lines = 1;
            def = {'10','1'};
            answer = inputdlg(prompt,dlg_title,num_lines,def);
            Q=answer{1,1};d1=answer{2,1};
            ke=Q./(pi*d1.^2/4)/(2*9.81);
            HL=ke;
        case 2
            prompt = {'請輸入管路流量(m^3/s)[10]','請輸入入口管徑(m)[1]:'};
            dlg_title = '您所選擇的是驟然逸出';
            num_lines = 1;
            def = {'10','1'};
            answer = inputdlg(prompt,dlg_title,num_lines,def);
            Q=answer{1,1};d1=answer{2,1};
            ke=Q./(pi*d1.^2/4)/(2*9.81);
            HL=ke/2;
        case 3
            prompt = {'請輸入管路流量(m^3/s)[10]','請輸入入口管徑(m)[1]:','請輸入出口管徑(m)[2]:'};
            dlg_title = '您所選擇的是擴大管';
            num_lines = 1;
            def = {'10','1','2'};
            answer = inputdlg(prompt,dlg_title,num_lines,def);
            Q=answer{1,1};d1=answer{2,1};d2=answer{3,1};
            ke=Q./(pi*d1.^2/4)/(2*9.81);
            HL=(1-(d1./d2).^2).*ke;
    end
    figure('Position',[400 400 400 300],'Name','結果')
    uicontrol('style','text','position',[160 200 80 80],'string','此接頭損失為')
    uicontrol('style','text','position',[100 90 200 100],'string',HL)
    result=questdlg('是否繼續執行？','程式執行詢問','Yes','No','No');
    if strcmp(result,'No');
        break
    end
end

或是點此網頁