|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 chybeyond 于 2014-6-9 15:46 编辑
我想求一个总矩阵,每个元素又是子矩阵,请大虾帮我看看,主程序如下,三个小调用程序已经上传。总是出现这个错误,谢谢!
??? Conversion to double from sym is not possible
- clc;
- syms xi xo x y i j ii jj iii jjj M N r h l
- M=2;
- N=2;
- mu=0.3;
- L=2;
- R=1;
- Ri=0.2;
- H=0.1;
- t=2*mu/(1-2*mu); %lamta/G
- %%%%无量纲值计算公式
- r=Ri/R;
- h=H/R;
- l=R/L;
- xi=l*sqrt((1-h/2)^2-y^2);
- xo=l*sqrt((1+h/2)^2-y^2);
- %%%%调用各个函数
- T1=T_1(n);
- T2=T_2(n);
- f10=f_10(x,y,i,j,ii,jj,iii,jjj,L,R,Ri,H,M,N,xi,xo,r,h,l);
- %%%%质量子矩阵
- MM=zeros(3);
- m11=zeros((M+1)*(N+1));
- m22=m11;
- m33=m11;
- %%%%求解质量子矩阵内的各个元素
- for iii=1:(M+1)*(N+1)
- for jjj=1:(M+1)*(N+1)
- jj=mod(jjj,N+1)-1; %求余数
- if jj==-1
- jj=N;
- end
- j=mod(iii,N+1)-1;
- if j==-1
- j=N;
- end
- ii=mod(jjj-jj,N)-1;
- if ii==-1
- ii=M;
- end
- i=mod(iii-j,N)-1;
- if i==-1
- i=M;
- end
- T1=subs(T1);
- T2=subs(T2);
- f10=subs(f10);
- m11(iii,jjj)=T1*f10;
- m22(iii,jjj)=T1*f10;
- m33(iii,jjj)=T2*f10;
- end
- end
- %%形成总质量矩阵
- MM=[m11 0 0;0 m22 0;0 0 m33]
复制代码
|
|