matlab求解方程

1.命令行求解

syms a b c y x

[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[x,y])

2.在编辑器中求解

clear
syms A1 A2 B1 B2 r1 r2 t41 t42
k=0.01
n=1;
for d=0:300
for q=1:10
eq1=A1+A2+B1+B2-sqrt(2)*(1+r1);
eq2=A1-A2+B1-B2-sqrt(2)*(1+r2);
eq3=1+r1-t41;
eq4=1+r2-t42;
eq5=A1*exp(i*(q-k)*d)+A2*exp(i*(q+k)*d)+B1*exp(-i*(q+k)*d)+B2*exp(-i*(q-k)*d);
eq6=A1*exp(i*(q-k)*d)-A2*exp(i*(q+k)*d)+B1*exp(-i*(q+k)*d)-B2*exp(-i*(q-k)*d);
eq7=t41*(k-q)+r1*(q+k)+(A1*(k-q)-A2*(q+k)+B1*(q+k)+B2*(q-k))/sqrt(2)+(k-q)*t41;
eq8=t42*(-k-q)+r2*(q-k)+(A1*(k-q)+A2*(q+k)+B1*(q+k)-B2*(q-k))/sqrt(2)-(k+q)*t42;
[A01(n),A02(n),B01(n),B02(n),r01(n),r02(n),t041(n),t042(n)]=solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,A1,A2,B1,B2,r1,r2,t41,t42);
n=n+1;
end
end
disp(' A01 A02 B01 B02 r01 r02 t041 t042 ')
X=[double(A01);double(A02);double(B01);double(B02);double(r01);double(r02);double(t041);double(t042)]'