Question

如何在MATLAB中求解这个复杂的超越方程，要有具体代码

此超越方程为变幅杆的频率方程，已知R1=25, R2=12.5, R3=10, L1=30, L2=20, K=24.8 ,求解L3=？本人刚开始学MATLAB，还望各位好心人多多帮忙！！！





Answer 1

该超越方程可以用matlab的vpasolve（）函数命令求得。求解过程如下：

syms L3

R1=25;R2=12.5;R3=10;L1=30;L2=20;k=24.8;

eq1=(R2^2+R3^2*(tan(k*L2))^2)*(1/tan(k*L1)+(R1-R2)/(k*L1*R2));

eq2=(R2^2*(tan(k*L2))^2+R3^2)*tan(k*(L2+L3));

eq3=(R3^2-R2^2)*(1-tan(k*(L2+L3))/tan(k*L1)-(R1-R2)*tan(k*(L2+L3))/(k*L1*R2))*tan(k*L2);

L3=vpasolve(eq1-eq2+eq3)

执行结果为

L3=-0.03476≈0

追问

很感谢您的热心回答，不过我刚才按照您的代码运行了一下，出错了，麻烦您给看一下，谢谢！

追答

你把代码复制编辑器里，然后保存为一个文件，再运行。

追问

按照您说的，试了一下，果然可以，但是我想问您一下为啥直接在命令执行窗口就会出错，还有就是运算结果不对，应该是0.036，还得麻烦您给解答一下，长度忘记单位了，是毫米。

追答

你有没有注意到多了这些符号﹥﹥

追问

程序运行没问题，但是运算结果不对，咋回事，正确值应该是36e-3米，还希望您再帮忙看一下，谢谢

追答