emwave 发表于 2009-9-5 12:59:53

符号方程组的求解，谢谢！

本帖最后由 ljelly 于 2009-9-8 09:00 编辑

符号方程组的matlab程序如下所示，
我在求解时，发现程序不能得出结果。
注释部分为该方程组的一个正确解，
请高手帮忙看看，如何解决这个问题，感激不尽！

clear;clc
syms S2 T2 S3 T3;
R = 0.1;
k = 4;
theta_m = 1.3351769;
S1 = 2;
T1 = 1.163;
% S2 = 2;
% S3 = 3.9055;
% T2 = 0.82026;
% T3 = 0.958;
E1 = 2*T2^2*T3^2*(2*(S1+S2)-T1^2) - T1^4*T3^2*(S2+S3) + 4*S3*T2^4 - 5/k/cos(theta_m);
E2 = (S2+S3)*( 4*(S1+S2)*(2*S3*T2^2-T3^2*(T1^2-(S1+S2))) - S3*T1^4*(S2+S3) - 4*S3*T1^2*T2^2 ) - 20/k/cos(theta_m)^3;
E3 = S3*(S1+S2)*(S2+S3)^2*((S1+S2)-T1^2) - 4/k/cos(theta_m)^5;
E4 = T1^2*T2^2*T3^2 - 1/(2*k*sin(theta_m)*sqrt(10^(R/10)-1));
=solve(E1,E2,E3,E4)

emwave 发表于 2009-9-6 16:17:31

请高手现身。。。

debateshang 发表于 2009-9-6 17:57:50

不是高手，但你可以试试fsolve
刚才 用fsolve试了试结果为
1.9892    3.9232    0.8204    0.9619
代码就不写了，help fsolve就ok了

emwave 发表于 2009-9-7 21:11:19

我试过了，可以
感谢debateshang友的帮助！
页: [1]
查看完整版本: 符号方程组的求解，谢谢！