qy469110098
2011-03-15, 16:43
x=linspace(1,10)
A=[0 -1 -2
1 0 -1
2 1 0]
c1=1/360*pi
D=5
d=8
R1=(1-erf(-sqrt(2)*pi/c1))^(-1).*exp(j*D*sin(x).*A)./(1+c1^2/2*(D*cos(x).*A).^2)
Y=d.*R1
E=eig(Y)
e1=E(1)
e2=E(2)
e3=E(3)
y=1/2*(e1^2*(1-sqrt(e1/(1+e1)))/(e1-e2)/(e1-e3)+e2^2*(1-sqrt(e2/(1+e2)))/(e2-e1)/(e2-e3)+e3^2*(1-sqrt(e3/(1+e3)))/(e3-e2)/(e3-e1))
plot(x,y)
运行出错
Matrix dimensions must agree
求大侠解释
A=[0 -1 -2
1 0 -1
2 1 0]
c1=1/360*pi
D=5
d=8
R1=(1-erf(-sqrt(2)*pi/c1))^(-1).*exp(j*D*sin(x).*A)./(1+c1^2/2*(D*cos(x).*A).^2)
Y=d.*R1
E=eig(Y)
e1=E(1)
e2=E(2)
e3=E(3)
y=1/2*(e1^2*(1-sqrt(e1/(1+e1)))/(e1-e2)/(e1-e3)+e2^2*(1-sqrt(e2/(1+e2)))/(e2-e1)/(e2-e3)+e3^2*(1-sqrt(e3/(1+e3)))/(e3-e2)/(e3-e1))
plot(x,y)
运行出错
Matrix dimensions must agree
求大侠解释