Question

matlab求解一元多次隐式函数

各位大神，本人在写毕业论文的过程中，需求解一个一元多次的方程，试了很多办法，终无结果，毕业时间越来越紧张，在此拜托各位大神帮帮忙，本人必有重谢：程序代码都写好了，需求解式中的Fai值。可随时联系我，QQ：1013346516
clear,clc

syms Fai

L_2=6.5*10^(-3);
L_3=41*10^(-3);

L_hj=35.826665*10^(-3);
L_gh=7.502767*10^(-3);

x_j1=15.999999*10^(-3);
y_j1=-20.828428*10^(-3);

x_j2=10.037947*10^(-3);
y_j2=24.270621*10^(-3);

x_j3=-26.037948*10^(-3);
y_j3=-3.442193*10^(-3);

x_G1=-26.702767*10^(-3);
y_G1=-27.5*10^(-3);

x_G2=37.167081*10^(-3);
y_G2=-9.375274*10^(-3);

x_G3=-10.464317*10^(-3);
y_G3=36.875275*10^(-3);

ARFA_1=0*2*pi/360;
ARFA_2=120*2*pi/360;
ARFA_3=240*2*pi/360;
L_gh1=7.502767*10^(-3);
L_gh2=7.502767*10^(-3);
L_gh3=7.502767*10^(-3);

d_1=120*10^(-6);%该值可变
d_2=10*10^(-6);%该值可变
d_3=50*10^(-6);%该值可变

R_1=x_j1*cos(Fai)-y_j1*sin(Fai)-(L_gh1-d_1)*cos(ARFA_1)-x_G1;
R_2=x_j2*cos(Fai)-y_j2*sin(Fai)-(L_gh2-d_2)*cos(ARFA_2)-x_G2;
R_3=x_j3*cos(Fai)-y_j3*sin(Fai)-(L_gh3-d_3)*cos(ARFA_3)-x_G3;

S_1=x_j1*sin(Fai)+y_j1*cos(Fai)-(L_gh1-d_1)*sin(ARFA_1)-y_G1;
S_2=x_j2*sin(Fai)+y_j2*cos(Fai)-(L_gh2-d_2)*sin(ARFA_2)-y_G2;
S_3=x_j3*sin(Fai)+y_j3*cos(Fai)-(L_gh3-d_3)*sin(ARFA_3)-y_G3;

C_41=2*(R_2-R_1);
C_42=2*(S_2-S_1);
C_43=(R_1-R_2)*(R_1+R_2)+(S_1-S_2)*(S_1+S_2);

C_51=2*(R_3-R_1);
C_52=2*(S_3-S_1);
C_53=(R_1-R_3)*(R_1+R_3)+(S_1-S_3)*(S_1+S_3);

((C_52*C_43-C_53*C_42)/(C_52*C_41-C_42*C_51)+R_1)^2+((C_53*C_41-C_43*C_51)/(C_52*C_41-C_42*C_51)+S_1)^2-L_hj^2=0;%求解该方程

Answer 1

你的问题，可以直接用solve（）函数命令求得。即

Fai=solve(((C_52*C_43-C_53*C_42)/(C_52*C_41-C_42*C_51)+R_1)^2+((C_53*C_41-C_43*C_51)/(C_52*C_41-C_42*C_51)+S_1)^2-L_hj^2)

另外，可以将sin（）改为sind（）,cos（）改为cosd（）,前者的单位为弧度，后者的单位为角度

运行结果

Fai =-.1308  （°）