| 查看: 550 | 回复: 1 | |||
[求助]
Zarray=[Zarray,z(Z)];这个啥意思，matlab的程序，这个语句在程序中的作用
|
|
function main20130210() %P=[1541000 4875180 4871930 ] P=[4000000 300 1541000 3000000 2000000 1541000] Zarray=[]; L=length(P); for ii=1:L z=[]; T=360; Tc=408.1; Pc=3.65*10^6; %单位为pa w=0.176; % 偏心因子 Pr=P(ii)/Pc ; %对比压力,这里得出的Pr是数组； Tr=T/Tc ; %对比温度； z(1)=1; m=0.480+1.574*w-0.176*w^2 ; % 具体的数字； F=(1+m*(1-Tr^0.5))^2/Tr ; % 具体的数字； for j=1:100 h(j)=0.08664*Pr/(z(j)*Tr); zz(j)=1/(1-h(j))-4.9340*F*h(j)/(1+h(j)); t(j)=abs(zz(j)-z(j)); if t(j)<0.0001 z(j)=zz(j); % Z(j)为对应压力下计算得出的压缩因子的值 break %终止循环； else %与if对应的转移语句； z(j+1)=zz(j); end end % z Z=numel(z) ; %中间过程中计算得出的压缩因子的数目。 yz=z(Z); %最后得到的压缩因子的值。取值最后一个数。 Zarray=[Zarray,z(Z)]; end Zarray |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有300人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab语句问题
已经有5人回复
- matlab中如何通过for语句得到已有矩阵的部分数据
已经有3人回复
- 跪求matlab的for循环语句的编写
已经有20人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 计算强帖: 3
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2013-02-18 00:43:27
20