matlab中怎么解方程中参数怎么引用前面求出结果.function []=jiachunxll(P,d)%d为泄露孔直径毫米,P为液体压力,T0为天然气温度℃p=792;%甲醇密度A=(d/2000)^2*pi;P0=10130;%环境压力C0=0.5;g=9.8;hl=2;%液体距泄漏
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/12 11:49:18
![matlab中怎么解方程中参数怎么引用前面求出结果.function []=jiachunxll(P,d)%d为泄露孔直径毫米,P为液体压力,T0为天然气温度℃p=792;%甲醇密度A=(d/2000)^2*pi;P0=10130;%环境压力C0=0.5;g=9.8;hl=2;%液体距泄漏](/uploads/image/z/10211967-63-7.jpg?t=matlab%E4%B8%AD%E6%80%8E%E4%B9%88%E8%A7%A3%E6%96%B9%E7%A8%8B%E4%B8%AD%E5%8F%82%E6%95%B0%E6%80%8E%E4%B9%88%E5%BC%95%E7%94%A8%E5%89%8D%E9%9D%A2%E6%B1%82%E5%87%BA%E7%BB%93%E6%9E%9C.function+%5B%5D%3Djiachunxll%28P%2Cd%29%25d%E4%B8%BA%E6%B3%84%E9%9C%B2%E5%AD%94%E7%9B%B4%E5%BE%84%E6%AF%AB%E7%B1%B3%2CP%E4%B8%BA%E6%B6%B2%E4%BD%93%E5%8E%8B%E5%8A%9B%2CT0%E4%B8%BA%E5%A4%A9%E7%84%B6%E6%B0%94%E6%B8%A9%E5%BA%A6%E2%84%83p%3D792%3B%25%E7%94%B2%E9%86%87%E5%AF%86%E5%BA%A6A%3D%28d%2F2000%29%5E2%2Api%3BP0%3D10130%3B%25%E7%8E%AF%E5%A2%83%E5%8E%8B%E5%8A%9BC0%3D0.5%3Bg%3D9.8%3Bhl%3D2%3B%25%E6%B6%B2%E4%BD%93%E8%B7%9D%E6%B3%84%E6%BC%8F)
matlab中怎么解方程中参数怎么引用前面求出结果.function []=jiachunxll(P,d)%d为泄露孔直径毫米,P为液体压力,T0为天然气温度℃p=792;%甲醇密度A=(d/2000)^2*pi;P0=10130;%环境压力C0=0.5;g=9.8;hl=2;%液体距泄漏
matlab中怎么解方程中参数怎么引用前面求出结果.
function []=jiachunxll(P,d)%d为泄露孔直径毫米,P为液体压力,T0为天然气温度℃
p=792;%甲醇密度
A=(d/2000)^2*pi;
P0=10130;%环境压力
C0=0.5;
g=9.8;
hl=2;%液体距泄漏口高度
A0=0;%防火堤内液池面积
u=1.5;%风速
Dz=pi*10*p;
T0=293.15;%环境温度
for t=1:50
Qm=p*A*C0*(2*((P-P0)/p+g*hl))^0.5;%泄漏量
Hv=1169690;%蒸发热j/kg————
Cp=2.534;%热熔kj/(kg*K)
Tb=337.75;%泄漏液体的沸点K
Tt=20+273;%泄漏液体的储存温度K
a=1.29*10^(-7);%水泥地面热扩散系数
K=1.1;%水泥地面导热系数
b=5.285*10^(-3);%大气稳定程度e时系数
n=0.3;
P1=12.8*1000;%液体表面蒸汽压Pa
M=32;
R=8.314;
if Tt >= Tb
Fv=Cp*(Tt-Tb)/Hv;%闪蒸比例
else
Fv=0
end
Q1=Qm*Fv;%闪蒸蒸发率kg/s
D=5*Fv*Qm%带到空气中的液体量kg/s
Q2=K*A0*(T0-Tb)/(Hv*(pi*a*1)^0.5)
r=(A0/pi)^0.5;%液池半径
Q3=b*P1*M/(R*T0)*u^((2-n)*(2+n))*r^((4+n)/(2+n))%————
Ds=Qm-D-Q2-Q3;
Dh=D+Q2+Q3;%损失的液体质量
s0=(Dz-Qm*t)/10
if s0>=pi/2
y=solve('s0-pi/2=y*(1-y^2)^0.5+asin(y)','y')
y=1+y
else
y=solve('pi/2-s0=y*(1-y^2)^0.5+asin(y)','y')
y=1-y
end
hl=y
A0=Ds/(0.015*p)
if A0 >= 112
A0=112
end
end
程序中s0上面计算出的结果下面解方程一直引用不了
matlab中怎么解方程中参数怎么引用前面求出结果.function []=jiachunxll(P,d)%d为泄露孔直径毫米,P为液体压力,T0为天然气温度℃p=792;%甲醇密度A=(d/2000)^2*pi;P0=10130;%环境压力C0=0.5;g=9.8;hl=2;%液体距泄漏
将y=solve('s0-pi/2=y*(1-y^2)^0.5+asin(y)','y')代码中,去掉丹引号,即
y=solve(s0-pi/2-(y*(1-y^2)^0.5+asin(y)),'y')
同理,y=solve('pi/2-s0=y*(1-y^2)^0.5+asin(y)','y')改为
y=solve(pi/2-s0-(y*(1-y^2)^0.5+asin(y)),'y')