Question

MATLAB程序运行提示“Error using + Too many output arguments.”是怎么回事？

clc;clear all;
En=31.25*10^9;Ea=12.5*10^9;%岩石弹性模量、剪切模量
EV=3.5*10^9;EM=2.0*10^9;%节理法向、切向刚度
NV=150*10^6;NM=1.29*10^6;%节理法向、切向黏度
w=100*pi;f=50;%波频率
u=En/(2*Ea)-1;%泊松比
Pa=0;Sb=asin(sin(Pa)*sqrt((1-2*u)/(2*(1-u))));%P波入射角、S波反射角
p=2650;Cp=5830;Cs=3830;%岩石密度、P波波速、S波波速
Zp=p*Cp;Zs=p*Cs;
Dt=0.001;%时间间距
EZ=((En+EV)*En*EV+En*(NV^2)*(w^2))/((En+EV)^2+(NV^2)*(w^2))+(((En^2)*NV*w)/((En+EV)^2+(NV^2)*(w^2)))*i;
GZ=(Ea*(EM^2)+(NM^2)*(Ea+EM)*(w^2))/(EM^2+(NM^2)*(w^2))+(((EM^2)*NM*w)/(EM^2+(NM^2)*(w^2)))*i;%复模量
A=[Zp*cos(2*Sb);Zp*sin(2*Sb)*tan(Sb)*atan(Pa)];
B=[Zp*cos(2*Sb),-Zs*sin(2*Sb);-Zp*sin(2*Sb)*tan(Sb)*atan(Pa),-Zs*cos(2*Sb)];
C=[Zp*cos(2*Sb),Zs*sin(2*Sb);Zp*sin(2*Sb)*tan(Sb)*atan(Pa),-Zs*cos(2*Sb)];
D=[EZ*Dt*cos(Pa);GZ*Dt*sin(Pa)];
E=[-EZ*Dt*cos(Pa),EZ*Dt*sin(Sb);-GZ*Dt*sin(Pa),GZ*Dt*cos(Sb)];
F=[-EZ*Dt*cos(Pa),-EZ*Dt*sin(Sb);-GZ*Dt*sin(Pa),GZ*Dt*cos(Sb)];%参数定义
VIp(1)=0.01;VTp(1)=0;VTs(1)=0;%初始赋值
for i=1:1:100000
[VRp(i);VRs(i)]=-(B^(-1))*A*VIp(i)+(B^(-1))*C*[VTp(i);VTs(i)];
[VTp(i+1);VTs(i+1)]=(C^(-1))*D*VIp(i)+(C^(-1))*E*[VRp(i);VRs(i)]+(C^(-1))*(C+F)*[VTp(i);VTs(i)];
t=Dt*i;
Tkc(i)=VTp(i)/VIp(i);
Rkc(i)=VRp(i)/VIp(i);
plot(Tkc(i),t,'*');
hold on;%绘制拟合曲线
plot(Rkc(i),t,'0');
end
运行上述一段程序，出现错误是什么情况，谢谢解答

Answer 1

Multiple left-hand sides must be separated by commas.是这个错误吧。[VRp(i);VRs(i)]=-(B^(-1))*A*VIp(i)+(B^(-1))*C*[VTp(i);VTs(i)];多变量赋值怎么会有分号。。。应该用逗号。

追问

您说得对，改成逗号之后提示为输出参量过多是什么问题，谢谢

追答

你好好分析下你左右赋值的维度是不是一样的，比如一个三维的列向量是不能赋给一个四维的列向量的。

追问

左边和右边都是二行一列的矩阵，可是左边赋值用逗号不是变成一行二列了？

追答

你把右边也改成逗号
[VRp(i),VRs(i)]=-(B^(-1))*A*VIp(i)+(B^(-1))*C*[VTp(i),VTs(i)];

追问

改了矩阵行列就乱了啊，系数ABCDEF也是矩阵

追答

设置个中间变量
Y=-(B^(-1))*A*VIp(i)+(B^(-1))*C*[VTp(i);VTs(i)];

VRp(i)=Y(1)
VRs(i)=Y(2)