| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 431 人关注过本帖
标题:下面程序本来要得到24*24的矩阵的,结果是3*3的矩阵,请问问题出在哪里?
只看楼主 加入收藏
Juneus
Rank: 1
等 级:新手上路
帖 子:7
专家分:4
注 册:2013-1-9
结帖率:100%
收藏
已结贴  问题点数:5 回复次数:3 
下面程序本来要得到24*24的矩阵的,结果是3*3的矩阵,请问问题出在哪里?
下面程序本来要得到24*24的矩阵的,结果是3*3的矩阵,请问问题出在哪里?

E=210E9;NU=0.3;
a=35/3;b=45/4;c=10;
A1=NU/(1-NU);A2=(1-2*NU)/(2*(1-NU));
A3=E*(1-NU)/((1+NU)-(1-2*NU));
x=[1.667,1.667,1.667,1.667,-10,-10,10,-10];
y=[-30,-18.75,-18.75,-30,-30,-18.75,-18.75,-30];
z=[20,20,10,10,20,20,10,10];

% Ke=zeros(24,24);
for i=1;8;
    for j=1;8;
        ei=x(i)/a;
        ej=x(j)/a;
        fi=y(i)/b;
        fj=y(j)/b;
        gi=z(i)/c;
        gj=z(j)/c;
        G11=b*b*c*c*ei*ej*(3+fi*fj)*(3+gi+gj)+A2*a*a*(3+ei*ej)*...
            (3*b*b*gi*gj+fi*fj*(3*c*c+(b*b+c*c)*gi*gj));
        G12=3*a*b*c*c*(A1*ei*fj+A2*fi*ej)*(3+gi*gj);
        G13=3*a*b*b*c*(A1*ei*gj+A2*gi*ej)*(3+fi*fj);
        G21=3*a*b*c*c*(A1*ej*fi+A2*fj*ei)*(3+gi*gj);
        G22=a*a*c*c*fi*fj*(3+ei*ej)*(3+gi*gj)+A2*b*b*(3+fi*fj)*...
            (3*a*a*gi*gj+ei*ej*(3*c*c+(a*a+c*c)*gi*gj));
        G23=3*a*a*b*c*(A1*fi*gj+A2*gi*fj)*(3+ei*ej);
        G31=3*a*b*b*c*(A1*gi*ej+A2*ei*gj)*(3+fi*fj);
        G32=3*a*a*b*c*(A1*gi*fj+A2*fi*gj)*(3+ei*ej);
        G33=a*a*b*b*gi*gj*(3+ei*ej)*(3+fi*fj)+A2*c*c*(3+gi*gj)*...
            (3*a*a*fi*fj+ei*ej*(3*b*b+(a*a+b*b)*fi*fj));
        Ke(3*i-2,3*j-2)=A3/(72*a*b*c)*G11;
        Ke(3*i-2,3*j-1)=A3/(72*a*b*c)*G12;
        Ke(3*i-2,3*j)=A3/(72*a*b*c)*G13;
        Ke(3*i-1,3*j-2)=A3/(72*a*b*c)*G21;
        Ke(3*i-1,3*j-1)=A3/(72*a*b*c)*G22;
        Ke(3*i-1,3*j)=A3/(72*a*b*c)*G23;
        Ke(3*i,3*j-2)=A3/(72*a*b*c)*G31;
        Ke(3*i,3*j-1)=A3/(72*a*b*c)*G32;
        Ke(3*i,3*j)=A3/(72*a*b*c)*G33;
    end
end
2013-01-16 20:20
zklhp
Rank: 20Rank: 20Rank: 20Rank: 20Rank: 20
来 自:china
等 级:贵宾
威 望:254
帖 子:11485
专家分:33241
注 册:2007-7-10
收藏
得分:5 
% Ke=zeros(24,24);

是不是这里不注释就好了呢

不会这个 瞎说的。。。
2013-01-16 21:40
Juneus
Rank: 1
等 级:新手上路
帖 子:7
专家分:4
注 册:2013-1-9
收藏
得分:0 
嗯,已经知道哪里错了,与注释无关,是for语句中的冒号打成分号了
2013-01-16 21:59
Juneus
Rank: 1
等 级:新手上路
帖 子:7
专家分:4
注 册:2013-1-9
收藏
得分:0 
额,注释那行确实得把%去掉
2013-01-16 22:03
快速回复:下面程序本来要得到24*24的矩阵的,结果是3*3的矩阵,请问问题出在哪里 ...
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.017527 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved