查看: 139|回复: 1|关注: 0

[已答复] matlab两层迭代算不出结果,无报错

[复制链接]

新手

5 麦片

财富积分


050


1

主题

1

帖子

0

最佳答案
发表于 2019-2-11 12:22:35 | 显示全部楼层 |阅读模式
本帖最后由 xcs111 于 2019-2-11 12:35 编辑

H=6;
[T0,P0]=ambient(H);
k=1.4;
kg=1.33;
const=-kg/(kg-1);
EFFb=0.98;
EFFm=0.98;
deltaH=43000;%KJ/kg
ccloss=0.03;
exloss=0.215;
P5=P0/(1-exloss);
ERfpt=2;
ERfptacl=1;
ERmap=1.5;
ERcal=1.1;
RPM=50500;
    BETA=1;
    while(abs(ERfptacl-ERfpt)/ERfpt>0.01)
        BETA=BETA+0.01;
        [M1,PR,EFF]=compr_map(BETA,RPM);
        NDM1=M1;
        T1=T0;
        P1=P0;
        P2=P1*PR;
        T2=T1*(1+(PR^((k-1)/k)-1)/EFF);
        PWRc=Air_H(T2)-Air_H(T1);
        Cp2=Air_Cp(T2);        
        f=0.019;
        while (abs(ERmap-ERcal)/ERmap>0.01)
            f=f+0.0001;
            T3=T2;
            t=1100;
            while (abs(T3-t)/t>0.002)
                t=T3;
                AAA=(1+f)*Gas_Cp(t,f);
                T3=(Cp2*T2+EFFb*f*deltaH)/AAA;
            end
            P3=P2*(1-ccloss);
            PWRggt=PWRc/(EFFm*(1+f));
            T4=T3;
            tt=800;
            pwrggt=PWRggt/2;
            while (abs(PWRggt-pwrggt)/PWRggt>0.01)
                tt=T4;
                pwrggt=Gas_H(T3,f)-Gas_H(tt,f);
                T4=tt-1;
            end
            NDM3=NDM1*(1+f)*((T3/T1)^0.5)*(P1/P3);
            [ERggt,EFFggt]=GGT_map(NDM3);
            ERmap=ERggt;
            ERcal=(1-(1-T4/T3)/EFFggt)^(const);
        end
        P4=P3/ERcal;
        NDM4=NDM3*((T4/T3)^0.5)*(P3/P4);
        [ERfpt,EFFfpt]=FPT_map(NDM4);
        ERfptacl=P4/P5;
    end

无报错,主要就是两层while循环,但是运行一直显示busy,得不到结果(结果以disp任意参数为例都可)。
如果取消内层迭代(即内层while的条件删除,外层每循环一次,内层里的内容只执行一次),只外层迭代可以计算,且计算很快(能够得到不满足内层while条件的结果);如果删掉外层while判断条件,只内层迭代也可执行,(能够得到不满足外层while条件的结果);但是两层嵌套就是算不出来,算了一天一夜都是busy。
几个函数文件都在附件里,应该没有问题。
matlab7.5.0
实在找不出原因了,希望帮忙解决一下,无比感谢。



Air_Cp.m

374 Bytes, 下载次数: 2

Air_H.m

380 Bytes, 下载次数: 2

ambient.m

210 Bytes, 下载次数: 1

compr_map.m

1.48 KB, 下载次数: 1

二维插值

FPT_map.m

301 Bytes, 下载次数: 1

插值

Gas_Cp.m

476 Bytes, 下载次数: 1

Gas_H.m

442 Bytes, 下载次数: 2

GGT_map.m

313 Bytes, 下载次数: 1

插值

入门

97 麦片

财富积分


50500


0

主题

86

帖子

4

最佳答案
  • 关注者: 6
发表于 2019-2-13 19:03:24 | 显示全部楼层
第一个循环的条件一直满足所以结果出不来
203.png
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /2 下一条

快速回复 返回顶部 返回列表