24小时热门版块排行榜    

查看: 492  |  回复: 7

wy515077587

金虫 (小有名气)

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

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

» 猜你喜欢

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

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

wy515077587

金虫 (小有名气)

自己顶一下
2楼2014-08-23 13:44:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子


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

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的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by wy515077587 at 2014-08-24 10:05:37
这是有限元法计算程序的一句,别人编的程序,学习中......

嗯,不要全信别人的。自己也可以试试。或者百度某个函数,或者用Matlab的help查询函数的用法。都挺好的
7楼2014-08-24 10:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子


小木虫: 金币+0.5, 给个红包,谢谢回帖
这样看来就是赋值了,但是不要跟Matlab自带函数重名啊,容易引起误会。可以用换个名字啊
8楼2014-08-24 11:36:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wy515077587 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见