3楼: Originally posted by hppdyx at 2014-03-15 10:29:56
function pottery
%%
syms B2 A3 B3 k1 k2 k3 a b z y1 y2 y3 y4 me2 me3;
y1=int((B2*sin(k2*z))^2,z,0,a)+int((A3*exp(k3*z)+B3*exp(-k3*z))^2,z,a,b)-1;
y2=B2*sin(k2*a)-A3*exp(k3*a)-B3*exp(k3*a);
y3=B2 ...

function pottery
%%
syms B2 A3 B3 k1 k2 k3 a b z y1 y2 y3 y4 me2 me3;
y1=int((B2*sin(k2*z))^2,z,0,a)+int((A3*exp(k3*z)+B3*exp(-k3*z))^2,z,a,b)-1;
y2=B2*sin(k2*a)-A3*exp(k3*a)-B3*exp(k3*a);
y3=B2*k2*cos(k2*a)/me2-(A3*k3*exp(k3*a)-B3*k3*exp(-k3*a))/me3;
% y4=A3*exp(k3*b)+B3*exp(-k3*b);
% y5=k3*sin(k2*a)/cos(k2*a)*me2+k2*me3*(exp(k3*b-k3*a)-exp(-k3*b+k3*a))/(exp(k3*b-k3*a)+exp(-k3*b+k3*a));
S = solve(y1==0, y2==0, y3==0, B2, B3, A3);
S = [S.B2; S.B3; S.A3];
display(S);
% S.B2; S.B3; S.A3;