24小时热门版块排行榜    

查看: 493  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见