matlab替换变量名,MATLAB只是简单地把表达式里的变量名替换成数值，而不给出结果...

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

%1.1.2.1快群

%堆芯慢化能力

%kesi_macr=zeros(5,1);

%for i=1:5

%kesi_macr=kesi(1,i)*N(1,i)*micro_s1(1,i)

%end

kesi_macr=kesi_H*N_H*micro_s1_H+kesi_O*N_O*micro_s1_O+kesi_Zr*N_Zr*micro_s1_Zr+kesi_5*N_5*micro_s1_5+kesi_8*N_8*micro_s1_8;

%快群宏观输运截面

macr_tr1=N_5*micro_s1_5*0.9972+N_8*micro_s1_8*0.9972+N_H*micro_s1_H*0.3386+N_O*micro_s1_O*0.9583+N_Zr*micro_s1_Zr*0.9927

%快群宏观扩散系数

D1=1/(3*macr_tr1)

%1.1.2.1热群

L_sigama_th=N_5*micro_a_5+N_8*micro_a_8+N_O*micro_a_O+N_H*micro_a_H+N_Zr*micro_a_Zr

L_sigama_kTM=L_sigama_th*sqrt(293/TM)       %式(2-80)

Tn=TM*(1+1.4*L_sigama_kTM/kesi_macr)        %式(2-79)

ga_5=1.03888*(Tn-273)^(-0.0184)

ga_8=0.9855*(Tn-273)^(0.0043)

gf_5=1.053*(Tn-273)^(-0.0226)

micro_a2_5=micro_a_5/1.128*sqrt(293/Tn)*ga_5;

micro_a2_8=micro_a_8/1.128*sqrt(293/Tn)*ga_8;

micro_a2_Zr=micro_a_Zr/1.128*sqrt(293/Tn);

micro_a2_H=micro_a_H/1.128*sqrt(293/Tn);

micro_a2_O=micro_a_O/1.128*sqrt(293/Tn);

micro_f2_5=micro_f_5/1.128*sqrt(293/Tn)*gf_5;

%堆芯的热中子宏观吸收截面

macro_a2=N_5*micro_a2_5+N_8*micro_a2_8+N_H*micro_a2_H+N_O*micro_a2_O+N_Zr*micro_a2_Zr

%堆芯的热中子宏观裂变截面

macro_f2=N_5*micro_f2_5

%L_H2O

L_H2O=2.719/den_H20_311*((t+273)/295)^(1/2)*(1.0614/(1+0.0614*((t+273)/295)^(1/2)))^(1/2)

%水的热群扩散系数

D2_H2O=L_H2O^2*(N_H*micro_a2_H+N_H2O*V_H2O/V*micro_a2_O)

macr_tr2_H2O=1/(3*D2_H2O)

问题是：为什么从

L_sigama_th=N_5*micro_a_5+N_8*micro_a_8+N_O*micro_a_O+N_H*micro_a_H+N_Zr*micro_a_Zr

L_sigama_kTM=L_sigama_th*sqrt(293/TM)       %式(2-80)

Tn=TM*(1+1.4*L_sigama_kTM/kesi_macr)        %式(2-79)

开始计算的结果如下呢：

L_sigama_kTM =

.30589717806208658168495092291557e-5*17115595^(1/2)

Tn =

11683/20+.45936265372953731364559710690310e-1*17115595^(1/2)

ga_5 =

6493/6250/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(23/1250)

ga_8 =

1971/2000*(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(43/10000)

gf_5 =

1053/1000/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(113/5000)

macro_a2 =

.13711453320619852340828168976573e-1*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(23/1250)+.17678485256047348144979250378501e-2*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)*(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(43/10000)+.84913479087977694805938554307674e-3*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)

macro_f2 =

.11909787189519441142160261181845e-1*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(113/5000)

L_H2O =

0.0054

D2_H2O =

69386399805210663738033738022873015464822941615415/3138550867693340381917894711603833208051177722232017256448*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)

macr_tr2_H2O =

3138550867693340381917894711603833208051177722232017256448/60990645428780173425731655722105380593579365679949785*293^(1/2)*(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)