高人可以帮我高试一下吗 为什么第四行有错啊？？？？
clearxs=2*0.287;
fid=fopen('c\\:临时文件\\:EL centro','r');
dzhbo=fscanf(fid,'%f');
ag=dzhbo*0.01*xs;
dt=0.02;
ndzh=400;
cn=3;
m0=*1e+3;
k0=*1e+5;
l=diag(ones(cn));
m=diag(m0);
=matrixju(k0,cn);
=eig(ik,m);
d=diag(sqrt(d));
for i=1:cn;
=min(d);
xgd(:,i)=x(:,j);
d(j)=max(d)+1;
end
w=d1;
x=xgd;
a1=2*w(1)*w(2)*(0.05*w(2)-0.07*w(1))/(w(2)^2-w(1)^2);
a2=2*(0.07*w(2)-0.05*w(1))/(w(2)^2-w(1)^2);
for j=1:cn
x(:,j)=x(:,j)/x(cn,j);
znb0(j)=(a1+a2*w(j)^2)/2/w(j);
zhcan(j)=(x(:,j))'*m*l/((x(:,j))'*m*x(:,j));
=zxzj(znb0(j),w(j),ag);
end
for i=1:cn;
disp1=0;
accel1=0;
for j=1:cn
disp0=zhcan(j)*dltacceler(j,:)*x(i,j);
accel0=zhcan(j)*dltacceler(j,:)*x(i,j);
disp1=disp1+disp0;
accel1=accel1+accel0;
end
disp(i,:)=disp1;
accel(i,:)=accel1;
end
t=0:dt:ndzh*dt;
subplot(2,2,1)
plot(t,disp(3,:)*1e+3,'b')
subplot(2,2,2)
plot(t,accel(3,:),'g')
subplot(2,2,3)
plot(t,dzhbo,'r')
回复 楼主 pzstar 的帖子
首先检查下第3行的fid是否为-1, 再检查档案格式是否对! "："是否需要？
页:
[1]