很奇怪的符号矩阵赋值（???）怎么进行动态赋值哪？

siyanger

用以下方法给符矩阵赋值，有时能给出结果，有时却出现以下的出错信息，不知道这什么哪？<br>出错信息<br>??? Conversion to double from sym is not possible.<br><br>syms U1 V1 U2 V2 U3 V3 U4 V4 <br>syms Ksp Csp Krp Crp Ks Kr Ka i12   %  i12=rbp1/rbp2<br>for (i=1:5)<br>    for(j=1:5)<br>        if(i==1) % 1<br>            if(j==1) H(i,j) =1;<br>            else H(i,j) =0;<br>            end<br>        elseif(i==2) % 2<br>            if(j==1) H(i,j) =(U1+U4)*Ksp;<br>            elseif(j==2) H(i,j) =2*(U1+U4)*Csp;<br>            elseif(j==3) H(i,j) =U1*Ksp;<br>            elseif(j==4) H(i,j) =2*U1*Csp; <br>            else H(i,j) = 0;<br>            end<br>       elseif(i==3) % 2<br>            if(j==1) H(i,j) =U2*Ksp;<br>            elseif(j==2) H(i,j) =2*U2*Csp;<br>            elseif(j==3) H(i,j) =(U2+U4)*Ksp;<br>            elseif(j==4) H(i,j) =2*(U2+U4)*Csp; <br>            else H(i,j) = 0;<br>            end  <br>        elseif(i==4) % 2<br>            if(j==1) H(i,j) =U3*Ksp;<br>            elseif(j==2) H(i,j) =2*U3*Csp;<br>            else H(i,j) = 0;<br>            end   <br>        else<br>            H(i,j) = 0;<br>        end<br>      end<br> end<br><br><FONT color=#e6421a>怎么一会可出结果一会又提示出错呢？</FONT>  大家有没有遇到过啊？          <br><FONT color=#e61a6b>??? Conversion to double from sym is not possible.<br></FONT>有时会出结果，不知道为什么？<br><br>

[此贴子已经被作者于2006-6-19 15:41:35编辑过]