| 查看: 434 | 回复: 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 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有89人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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