siyanger 发表于 2006-6-19 09:50

很奇怪的符号矩阵赋值(???)怎么进行动态赋值哪?

用以下方法给符矩阵赋值,有时能给出结果,有时却出现以下的出错信息,不知道这什么哪?<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编辑过]
页: [1]
查看完整版本: 很奇怪的符号矩阵赋值(???)怎么进行动态赋值哪?