实现了参数的传递程序如下:syms k unw=2;ng=1;d=17;x1=4;z1=3;h=1;[k,u]=solve(nw^2/(ng^2)+h^2/ng^2*(k-u)^2-(z1^2/(x1-k)^2)-1,(ng^2*(d^2+u^2)-u^2)*(k-u)^2-(h*u)^2)运行结果:k =5.70377268258286107727264771430512.27253227917761092390431295578145.69657811689502986076014670178622.2737497376940458328574334799876u =5.386**2243936879520194506351032.14628048588996587257629556934916.05261424920096922705765587064792.4158590962999236974110230724868不过现在又有一个新的问题,上面的系数变量中前两个nw,ng其实可以是固定不变的,所以我直接写到方程里,这样变成如下:syms k ud=17;x1=4;z1=3;h=1;[k,u]=solve(4+h^2/(k-u)^2-(z1^2/(x1-k)^2)-1,((d^2+u^2)-u^2)*(k-u)^2-(h*u)^2)k =4.79841291699915956872494193203953.4015034167698993289352485030634.71571359061161051172221932610333.453751018098001328107450582844u =5.09831372431160704177025080279193.61409738031801803699370153450444.45372950224429881662654047465313.2618759615370012543237033282415 20210311
function f=fx(x)f(1)=1+cos(x(1))*cosh(x(1))-x(2);f(2)=2.4*x(1)*(sin(x(1))*cosh(x(1))-cos(x(1))*sinh(x(1)))-x(2);[x,f]=fsolve(@fx,[1,30])x=1.0328 1.8109f=1.0e-010*-0.0955 0.9024>>.