tgr.m

来自「已知过热蒸汽压力和焓」· M 代码 · 共 45 行

%超临界下过热蒸汽的焓：28~23 340~610

P1=[23 24 25 26 27 28 23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28 ...
     23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28 ...
      23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28 ...
       23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28 ...
        23 24 25 26 27 28  23 24 25 26 27 28  23 24 25 26 27 28 ]';
    P=1.0e+6*P1;
T=[340 340 340 340 340 340 350 350 350 350 350 350 360  360  360  360  360  360 370  370  370  370  370  370 380  380  380  380  380  380 ...
     390  390  390  390  390  390 400  400  400  400  400  400 410  410  410  410  410  410 420  420  420  420  420  420 430  430  430  430  430  430 ...
     440  440  440  440  440  440 450  450  450  450  450  450 460  460  460  460  460  460  470   470   470   470   470   470 480  480  480  480  480  480 ...
     490  490  490  490  490  490 500  500  500  500  500  500  510   510   510   510   510   510 520  520  520  520  520  520  530  530  530  530  530  530 ...
     540  540  540  540  540  540  550  550  550  550  550  550  560  560  560  560  560  560  570  570  570  570  570  570  580  580 580 580 580 580 ...
     590   590  590  590  590  590 600  600 600 600 600 600 610  610 610 610 610 610]';
I=[1564.7 1562.0 1559.5 1557.2 1555.0 1552.9...
    1634.1 1630.1 1626.4 1623.0 1619.8 1616.8...
    1714.1 1707.3 1701.3 1696.0 1691.1 1686.6...
    1821.0 1804.8 1792.3 1782.1 1773.4 1765.8...
    2373.2 2028.0 1939.5 1904.6 1882.2 1865.5...
    2584.2 2505.3 2400.3 2250.6 2103.2 2028.4 ...
    2693.6 2642.0 2583.2 2514.9 2433.9 2339.3 ...
    2774.7 2734.7 2691.3 2643.8 2591.4 2533.4...
    2841.4 2808.2 2773.0 2735.5 2695.4 2652.5...
    2899.2 2870.6 2840.6 2809.2 2776.2 2741.4...
    2951.0 2925.7 2899.4 2872.2 2843.8 2814.3...
    2998.2 2975.5 2952.1 2927.2 2903.0 2877.2 ...
    3042.1 3021.5 3000.3 2978.5 2956.1 2933.2 ...
    3083.4 3064.4 3045.0 3025.2 3004.9 2984.1 ...
    3122.5 3104.9 3087.0 3068.8 3050.2 3031.2...
    3159.9 3143.5 3126.9 3110.0 3092.8 3075.4...
    3195.8 3180.5 3165.0 3149.3 3133.3 3117.1...
    3230.5 3216.2 3201.6 3186.9 3172.0 3156.8...
    3264.2 3250.7 3237.0 3223.2 3209.2 3195.0 ... 
    3297.0 3284.2 3271.3 3258.2 3245.0 3231.7...
    3329.0 3316.9 3304.7 3292.3 3279.9 3267.3...
    3360.4 3348.9 3337.3 3325.6 3313.8 3301.8 ...
    3391.2 3380.2 3369.2 3358.0 3346.8 3335.5...
    3421.5 3411.0 3400.5 3389.8 3379.1 3368.4...
    3451.3 3441.3 3431.2 3421.1 3410.8 3400.6...
    3480.7 3471.1 3461.5 3451.8 3442.0 3432.2...
    3509.7 3500.5 3491.2 3481.9 3472.6 3463.2...
    3538.3 3529.4 3520.7 3511.7 3502.7 3493.7]';
x=[ones(size(P)) ,log10(P) ,log10(T)];
   a=x\log10(I)