24小时热门版块排行榜    

查看: 620  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wy515077587

金虫 (小有名气)

[交流] 一句程序含义 已有2人参与

disp(1:jdx*jdy,1:5)=1;表示什么意思呢?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wy515077587

金虫 (小有名气)

引用回帖:
3楼: Originally posted by hytao2012 at 2014-08-23 16:18:04
你这程序没问题?感觉这句话与本意违背。
这句话的实际效果是创建一个5*(jdx*jdy)的矩阵,且令每一个值都等于1
但是实际要表达什么的确让人费解

这是有限元法计算程序的一句,别人编的程序,学习中...
4楼2014-08-24 10:05:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

hytao2012

铁杆木虫 (正式写手)

木头虫子


小木虫: 金币+0.5, 给个红包,谢谢回帖
你这程序没问题?感觉这句话与本意违背。
这句话的实际效果是创建一个5*(jdx*jdy)的矩阵,且令每一个值都等于1
但是实际要表达什么的确让人费解
3楼2014-08-23 16:18:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbl37379

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这是matlab?为啥还弄个disp,难道disp是个行向量,你看看,是不是准备画图看迭代收敛,创建一个行向量是为了控制某个变量的长度,确实让人费解

[ 发自小木虫客户端 ]
他山之石,可以攻玉
5楼2014-08-24 10:43:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wy515077587

金虫 (小有名气)

引用回帖:
5楼: Originally posted by xbl37379 at 2014-08-24 10:43:54
这是matlab?为啥还弄个disp,难道disp是个行向量,你看看,是不是准备画图看迭代收敛,创建一个行向量是为了控制某个变量的长度,确实让人费解

clear
syms E poisson density t lx ly
jdx=11;               %number of nodes in x direction
jdy=11;               %number of nodes in y direction

k(1:110*5,1:110*5)=0;     %system stiffness matrix
m(1:110*5,1:110*5)=0;     %system mass matrix
D_k(1:110*5,1:110*5)=0;     %刚度矩阵的导数
D_m(1:110*5,1:110*5)=0;     %质量矩阵的导数

%prepare the arrays containing numbers of nodes needed to describe this problem
en(1:100,1:4)=0;       %element node  
for ni=1:jdx-1
    for nj=1:jdy-1
        en(ni+(nj-1)*(jdx-1),1)=ni+(nj-1)*jdx;
        en(ni+(nj-1)*(jdx-1),2)=ni+1+(nj-1)*jdx;
        en(ni+(nj-1)*(jdx-1),4)=ni+nj*jdx;
        en(ni+(nj-1)*(jdx-1),3)=ni+1+nj*jdx;
    end
end
%constraints for boundary conditions
disp(1:jdx*jdy,1:5)=1;      % node displacement
constraints=1:jdx:jdx*jdy; % constraints for boudaries
disp(constraints,=0;
dof=0;                  
这是源程序的一部分,研究中...
6楼2014-08-24 10:50:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
请填处理意见