gkyingkyin
2008-11-08, 11:14
function xielouliang
p=input('请输入燃气管道压强(相对压强，单位：兆帕)=');
Cd=1.0;
p0=1.01325e+005;
M=0.016;
R=8.314;
T=293.15;
r=1.307;
pc=1.837*1.01325e+005;
p1=p*1.0e+006+p0;
radius=[0.0005;0.001;0.0015;0.002;0.0025;0.003;0.0035;0.004;0.0045;0.005;0.0075;0.01];
for i=1:12
A(i)=pi*radius(i)*radius(i);
i=i+1;
end;
for i=1:12
if p1<pc;
Q(i)=Cd*p1*A(i)*[2*r*M*[(p0/p1)^(2/r)-(p0/p1)^((r+1)/r)]/(r-1)/R/T]^(0.5);
else if p1>=pc;
Q(i)=Cd*p1*A(i)*(r*M/R/T*(2/(r+1))^((r+1)/(r-1)))^(1/2);
end;
end;
end;
vpa(Q,5)%管道天然气泄露量估算
运行后用
mcc -B sgl xielouliang.m转成c语言的程序
运行程序在一开始输入变量时候显示Warning:not enough space string is truncated
还有运行结果一闪窗口就没了，怎么让它保留一下啊，我不懂c语言，请大家指点，谢谢:ft: ~~
p=input('请输入燃气管道压强(相对压强，单位：兆帕)=');
Cd=1.0;
p0=1.01325e+005;
M=0.016;
R=8.314;
T=293.15;
r=1.307;
pc=1.837*1.01325e+005;
p1=p*1.0e+006+p0;
radius=[0.0005;0.001;0.0015;0.002;0.0025;0.003;0.0035;0.004;0.0045;0.005;0.0075;0.01];
for i=1:12
A(i)=pi*radius(i)*radius(i);
i=i+1;
end;
for i=1:12
if p1<pc;
Q(i)=Cd*p1*A(i)*[2*r*M*[(p0/p1)^(2/r)-(p0/p1)^((r+1)/r)]/(r-1)/R/T]^(0.5);
else if p1>=pc;
Q(i)=Cd*p1*A(i)*(r*M/R/T*(2/(r+1))^((r+1)/(r-1)))^(1/2);
end;
end;
end;
vpa(Q,5)%管道天然气泄露量估算
运行后用
mcc -B sgl xielouliang.m转成c语言的程序
运行程序在一开始输入变量时候显示Warning:not enough space string is truncated
还有运行结果一闪窗口就没了，怎么让它保留一下啊，我不懂c语言，请大家指点，谢谢:ft: ~~