板翅式换热器（锯齿形）matlab的

 试设计某实验平台所需要使用的一台气水换热器。

设计校核计算通过matlab实现，请求斧正！

clear,clc
%申明所有压力单位统一为Pa
%% 预备计算
%%预备计算
%%假设换热器效率n=0.8，任务书中要求的是大于等于0.7，x假设为为冷测水的进口流量
%初始值设定
t1_1=90;t2_1=25;n=0.8
x=890;
t1_11 = t1_1 - n*(t1_1 - t2_1);%%计算出口气温
%%由查表得出的空气的物性参数
cp1 = 1003+0.02*43+4*0.00001*43^2;
%%由查表得出的水的物性参数
cp2 = 4184.4-0.6964*21+1.06964*21+1.036*10^(-2)*21^2;
qm1=(1000+20*35)/3600;%36为学号
qm2=x/3600;
W1=qm1*cp1;
W2=qm2*cp2;
C=W1/W2;%计算C
%由冷测进口流量和C*的关系可判断，C*在x=1600以内是永远大于0.5的
t2_11=t2_1+W1/W2*(t1_1-t1_11);
%Wmax侧（冷侧）取算数平均温度
tm2=(t2_1+t2_11)/2
%两流体的平均温差为
deta_t1m=((t1_1-tm2)-(t1_11-tm2))/log((t1_1-tm2)/(t1_11-tm2));
%热侧平均温度为
tm1=tm2+deta_t1m;
%根据平均温度求物性参数
mu1=1.50619*10^(-6)*(tm1+273)^1.5;%空气的动力黏度miu
lam1=2.456*10^(-4)*(tm1+273)^0.823;%空气导热系数λ1
Pr1=mu1*cp1/lam1;%计算热侧的普朗特数Pr1
mu2=10^(230.298/(tm2+126.203)-4.5668);%水的动力黏度
lam2=0.598+1.373*10^(-3)*tm2-5.333*10^(-6)*(tm2^2);%水的导热系数
Pr2=mu2*cp2/lam2%计算冷侧的普朗特数Pr2
mixsum=[x;mu1;lam1;Pr1;mu2;lam2;Pr2]';
%% 产品结构规划及其计算
%参数说明     冷端          热端
%间断长度     ls1            ls2
%翅片间距     Pf1            Pf2
%板间距       s1             s2 
%翅片厚度     Deltaf1        Deltaf2
%翅片高度     h1             h2
%翅片层数     m1             m2
%翅片厚度           Deltap        
%侧板厚度           Deltas       
%封条宽度     bs1            bs2
%热侧进口宽度     l2
%冷侧出口         l1
%热侧出口         l3

    ls1=240       ;     ls2=60     ;   %间断长度 
     Pf1=50     ;      Pf2=40       %翅片间距
     s1= 150       ;     s2= 50    ;    %板间距 
     Deltaf1=5.5    ;    Deltaf2=3.35   ; % 翅片厚度