救助：请大家帮忙看看程序哪里出错
这是根据单自由度振动系统写的，但是程序运行老出错。请大家帮忙看看m1=3;k1=3;c1=0.6;y0='
M=diag();
K=;
C=;
=size(M);
A=;
T=50;dt=0.01;
for i=1:round(T/dt)+1
t(i)=dt*(i-1);
y(:,i)=expm(A*t(i))*y0;
end
subplot(2,1,1),plot(t,y(1,:))
title('波形图')
subplot(2,1,2),plot(y(1,:),y(2,:))% 相图
title('相图')
=eig(A)
for i=1:4
f=p(i,:)
end
回复 楼主 412 的帖子
报错不是告诉LZ了吗?!Index exceeds matrix dimensions.
p为2*2矩阵, for i=1:4, f=p(i,:); end有误!
页:
[1]