[分享]符号表达

suffer

<P>j=input('请输入j的值:');<BR>k_1(1)=5;k_1(2)=5;k_1(3)=8;l=100;w=200;<BR>A=[];C=[];Q=[];W=[];<BR>kl=0;<BR>%这里A和C代表的是公式(4-79)中的Anj2和Cmj2，Q和W代表的是Anj1和Cmj1<BR>for tp=1:j<BR>    for tt=1:k_1(tp)<BR>           A=[A;'A' num2str(tp) num2str(tt) num2str(2) ];<BR>           C=[C;'C' num2str(tp) num2str(tt) num2str(2) ];<BR>           Q=[Q;'Q' num2str(tp) num2str(tt) num2str(1) ];<BR>           W=[W;'W' num2str(tp) num2str(tt) num2str(1) ];<BR>    end<BR>   kl=kl+k_1(tp);<BR>end<BR>for fd=1:kl<BR>    dk{fd}=strcat(A(fd,1),A(fd,2),A(fd,3),A(fd,4));<BR>    dc{fd}=strcat(C(fd,1),C(fd,2),C(fd,3),C(fd,4));<BR>    db{fd}=strcat(Q(fd,1),Q(fd,2),Q(fd,3),Q(fd,4));<BR>    de{fd}=strcat(W(fd,1),W(fd,2),W(fd,3),W(fd,4));<BR>end<BR>tk=1;<BR>sk=0;<BR>for i=1:j<BR>   sk=sk+k_1(i);<BR>   fe{i}=dk(:,[tk:sk]);<BR>   fs{i}=dc(:,[tk:sk]);<BR>   fa{i}=db(:,[tk:sk]);<BR>   fb{i}=de(:,[tk:sk]);<BR>   tk=tk+k_1(i);<BR>end<BR>%这样可以实现选取所需的变量名例如A112等，使用四个大括号<BR>%fe,fs矩阵里面包含的是Anj2和Cmj2,fa和fb里面包含的是Anj1和Cmj1<BR>for i=1:j<BR>      for n=1:k_1(i)<BR>          aa=num2str(4./1);<BR>          ffe=sym(strcat(aa,'*',fe{i}{n}));<BR>          ffs=sym(strcat(aa,'*',fs{i}{n}));<BR>          ffa=sym(strcat(aa,'*',fa{i}{n}));<BR>          ffb=sym(strcat(aa,'*',fb{i}{n}));<BR>            Bj1(:,n,i)=[0;0;0;ffe;ffs;0];<BR>            Bj2(:,n,i)=[0;0;0;ffa;ffb;0];<BR>    end<BR>end</P>