|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
qdzb
k=[1,1,1,1,1]*4.2e+7;
m0=[1,1,1,1,1]*2.0e+5;
g1=1.5e+7;
g2=2.01e+7;
nn=2;
cn=5;
as=3e-2;
t=1.3e-2;
O=zeros(cn);
I=eye(cn);
k0=z(k cn);
m=diag(mo);
[x,d]=eig(k0,m);
d=sqrt(d);
w=sort(diag(d));
a=2*w(1)*w(2)*(0.05*w(2)-0.07*w(1))/(w(2)^2-w(1)^2);
b=2*(0.07*w(2)-0.05*w(1))/(w(2)^2-w(1)^2);
co=a*m+b*k0;
nd=[1,1,1,1,1];
kd1=nn*g1*as/t;
cd1=nn*g2*as/w(1)/t;
kd=matrixju(kd1*nd,cn);
cd=matrixju(cd1*nd,cn);
A=[O I;-inv(m)*k0-inv(m)*c0];
B=[O;I];
R=inv(m)*[kd cd];
II=diag(I);
Ca=zeros(cn,2*cn);
Ca(cn,2*cn)=1;
Cw=zeros(cn,2*cn);
Cw(cn,cn)=1;
D=diag(O);
[kcju]=matrixju(korc,cn)
kcju=zeros(cn);
for i=1:cn-1
kcju(i,i)=korc(i)+korc(i+1);
kcju(i,i+1)=-korc(i+1);
kcju(i+1,i)=-korc(i+1);
end
kcju(cn,cn)=korc(cn);
qdzb怎么定义呀,怎样进行dzb模块的定义?请各位帮忙。我这个输进去,运行就是出错。
[ 本帖最后由 eight 于 2007-12-10 15:53 编辑 ] |
|