clear,format compact
c=3*10.^9;
er=2.54;
w1=input('w1=');
zh=input('zh=');
zl=input('zl=');
l1=input('l1=');
c2=input('c2=');
l3=input('l3=');
c4=input('c4=');
l5=input('l5=');
c6=input('c6=');
cf=input('cf=');
vh=c;
vl=c/sqrt(er);
yl=1/zl;
yh=1/zh;
g1='w1*l1=zh*sin(w1*s1/vh)+zl*s2*w1/(2*vl)';
g2='w1*c2=yl*s2*w1/vl+2*w1*cf+yh*s1*w1/(2*vh)+yh*s3*w1/(2*vh)';
g3='w1*l3=zh*sin(w1*s3/vh)+zl*s2*w1(2*vl+zl*s4*w1/2*vl)';
g4='w1*c4=yl*s4*w1/vl+2*w1*cf+yh*s3*w1/(2*vh)+yh*s5*w1/(2*vh)';
g5='w1*l5=zh*sin(w1*s5/vh)+zl*s2*w1/(2*vl)+zl*s4*w1/(2*vl)+zl*s6*w1/(2*vl)';
g6='w1*c6=yl*s6*w1/vl+2*w1*cf+yh*s5*w1/(2*vh)+yh*s5*w1/(2*vh)';
[x1,x2,x3,x4,x5,x6]=solve(g1,g2,g3,g4,g5,g6);
x1v=vpa(x1,6);
x2v=vpa(x2,6);
x3v=vpa(x3,6);
x4v=vpa(x4,6);
x5v=vpa(x5,6);
x6v=vpa(x6,6);

这些是上面参数，w1=1.971*10.^9;zh=150;zl=10;l1=4.826*10.^-9;c2=2.329*10.^-12;
l3=8.616*10.^-9;c4=2.610*10.^-12;
l5=8.902*10.^-9;c6=2.610*10.^-12;cf=0.465*10.^-12;
求出上面方程里的未知数，S1,S2,S3,S4,S5,S6