matlab中用for 循环+条件语句编程 Illegal use of reserved keyword "else". 5

functionF=Eigenvalue_Equation(delt,M,Tg,amax)%求地震力F(前三阶振型下)%Misamassmatrix质量矩阵%delt柔度... function F= Eigenvalue_Equation(delt,M,Tg,amax)
%求地震力F(前三阶振型下)
% M is a mass matrix 质量矩阵
% delt 柔度矩阵
% V is a eigenvector特征向量矩阵
% D is a eigenvalue特征值对角阵
%Tg 特征周期
%amax 水平地震影响系数最大值
ik=inv(delt);
[V,D]=eig(ik,M);
d=diag(sqrt(D));
for i=1:3
[d1(i),j]=min(d);
xgd(:,i)=V(:,j);
d(j)=max(d)+1;
end
w=d1; %所求自振频率
x=xgd; %所求结构主振型
for i=1:3
T(i)=0.7*2*pi/w(i);
end
for j=1:3
Mj=(x(:,j))'*M* x(:,j);
fi(:,j)=x(:,j)/sqrt(Mj); %振型的正交归一化
n=length(x(:,j));
I=ones(n,1);
zhcan(j)=(fi(:,j))'*M*I/(( fi(:,j))'*M* fi(:,j));
end
for i=1:3
if T(i)<=0.1
a(i)=0.45*amax+5.5*amax*T(i);
elseif(T(i)>0.1)&(T(i)<=Tg)
a(i)=amax;
elseif(T(i)>Tg)&(T(i)<=5*Tg)
l(i)=(Tg/T(i));
a(i)=l(i).^0.9*amax;
elseif(T(i)>5*Tg)&(T(i)<=6)
a(i)=[0.2.^0.9-0.02*(T(i)-5*Tg)*amax;
else a(i)==NaN; 【%此处Illegal use of reserved keyword "else".】另end老是回不到最前面而且不变蓝!!!
end
end
for j=1:3
for i=1:n
F(i,j)=a(j)*zhcan(j)*fi(i,j)*G(i);
end
end
展开
 我来答
沙冰巧001
2011-03-15 · TA获得超过824个赞
知道小有建树答主
回答量:827
采纳率:0%
帮助的人:742万
展开全部
速度大于从动盘角速度循环继续
if t<0.2
Mm=Mm(t+h);
a1=a1_1; %a1_1要有参数的,改成a1_1(t,w1)之类的
w1=w1+a1*h;
a2=a2_1; %同上
w2=w2+a2*h;
else Mm=beta*TN;
a1=a1_2; % 同上
w1=w1+a1*h;
a2=a2_2; %同上
w2=w2+a2*h;
end

另外,看看你这些函数是不是都另存到单独的m文件里了
hangeneral
2011-03-26 · TA获得超过146个赞
知道小有建树答主
回答量:164
采纳率:0%
帮助的人:119万
展开全部
else后面不能是判断语句,换行是很好的习惯。
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
Carrot_hy
2011-03-14 · TA获得超过1836个赞
知道小有建树答主
回答量:374
采纳率:0%
帮助的人:232万
展开全部
哥哥,matlab和c不一样,换一行试试嘛
这样,
else
a(i)=NaN;
或者

加一个逗号
else, a(i) = NaN;

还有赋值,不要加两个=号
本回答被网友采纳
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 1条折叠回答
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式