MATLAB运行出错怎么解决

落在沙滩的鱼

我运行原书上的M文件
结果出错：??? Error using ==> le
Not enough input arguments.

Error in ==> le at 18
  [varargout{1:nargout}] = builtin('le', varargin{:});

Error in ==> shicheng at 5
m=[2.762 2.760 2.300]*le+3;
请问怎么解决？

下面是源程序：
%MAIN PROGRAM
%response analysis of structure-elastic time history analysis method
%structure parameter
h=[4000,3300,3300];   %各层层高 mm
m=[2.762 2.760 2.300]*le+3;
k0=[2.485 1.921 1.522]*le+5;
le=diag(ones(cn));
cn=length(m);         %总层数
%earthquake parameter
ct=1.4;              %Wilson法sit值
dt=0.02;
xs=200/max(abs(dzhbo));  %调整地震输入加速度值
ag=dzhbo*0.01*xs;   %地震波，400步，步长为0.02S，单位m/s2
ndzh=400;
agl=ag(1:ndzh);
ag2=ag(2:ndzh+1);
agtao=ct*(ag2-ag1);
%initial value
chsh=zeros(cn,1);
wyi1=chsh;
sdu1=chsn;
jsdu1=chsh;
wyimt=chsh;
sdumt=chsh;
jsdumt=chsh;
unit=ones(cn,1);
m=diag(m);
[ik]=matrixju(k0,cn);
[x,d]=eig(ik,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);
c0=a*m+b*ik;
for i=1:ndzh
    kxin=ik+(3/(ct*dt))*c0+(6/(ct*ct*dt))*m;
    dpxin=-m*unit*agtao(i)+m*(6/(ct*dt)*sdu1+3*jsdu1)+c0*(3*sdu1+ct*dt/2*jsdu1);
    dxtao=kxin\dpxin;
    dtjsdu=6*dxtao/(ct*(ct^2*dt^2))-6*sdu1/(ct*ct*dt)-(3/ct)*jsdu1;
    jsdu=jsdu1+dtjsdu;
    dtsdu=(dt/2)*(jsdu+jsdu1);
    sdu=sdu1+dtsdu;
    dtwyi=dt*sdu1+(1/3)*dt^2*jsdu1+(dt^2/6)*jsdu;
    wyi=wyi1+dtwyi;
    jsdu=-m\(m*unit*ag2(i)+c0*sdu+ik*wyi);  %调整加速度
    wyi1=wyi;
    sdu1=sdu;
    jsdu1=jsdu;
    wyimt=[wyimt wyi*1000];
    sdumt=[sdumt sdu];
    jsdumt=[jsdumt jsdu];
end
t=0:dt:ndzh*dt;
subplot(2,2,1)
plot(t,wyimt(3,:),'r-')
subplot(2,2,2)
plot(t,jsdumt(3,:),'r-')

%子程序
%This sub-program is for matrix aggregation of system
function[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)

mklh

m=[2.762 2.760 2.300]*le+3;
这个的意思应该是科学计数法, 你改成乘以1000就好了

落在沙滩的鱼

非常感谢，问题解决了

神瑛侍者

>> %MAIN PROGRAM
%response analysis of structure-elastic time history analysis method
%structure parameter
h=[4000,3300,3300]; %各层层高 mm
m=[2.762 2.760 2.300]*l000;
k0=[2.485 1.921 1.522]*le+5;
le=diag(ones(cn));
cn=length(m); %总层数
%earthquake parameter
ct=1.4; %Wilson法sit值
dt=0.02;
xs=200/max(abs(dzhbo)); %调整地震输入加速度值
ag=dzhbo*0.01*xs; %地震波，400步，步长为0.02S，单位m/s2
ndzh=400;
agl=ag(1:ndzh);
ag2=ag(2:ndzh+1);
agtao=ct*(ag2-ag1);
%initial value
chsh=zeros(cn,1);
wyi1=chsh;
sdu1=chsn;
jsdu1=chsh;
wyimt=chsh;
sdumt=chsh;
jsdumt=chsh;
unit=ones(cn,1);
m=diag(m);
[ik]=matrixju(k0,cn);

??? Undefined function or variable 'l000'.
初学者，望高手指点！多谢！

ChaChing

l000 or 1000:@)

神瑛侍者

ChaChing 发表于 2011-1-3 22:56

                               
登录/注册后可看大图

l000 or 1000:@)

嗯嗯，多谢啦，现在看来当时的问题真是太初级了~