TA的每日心情 | 奋斗
2022-1-21 15:15 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0， 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc). 求出b的表达式。$ w3 d! X* E) h4 G4 T
代码如下：
: r, [ p" i1 _7 x, `, Z, D% Define the given equations1 u2 u- U8 Q0 m* J
syms pc qc N alpha po b pp qq q_ p_, F9 ?+ q5 C$ [' ^5 V0 r
eq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0;
/ L, ^* p- w5 K; a6 peq2 = q_ == qc + b * (qq - qc);7 S, {! S# h. w8 j) [ i- ?$ p
eq3 = p_ == pc + b * (pp - pc);8 r, ], j4 h+ C. |4 l _
- X0 r2 u; e- c& I8 h+ J' p6 v) e
% Solve for b; K0 x6 {( v: L8 j3 x6 w
eq4 = subs(eq1, [q_, p_], [eq2, eq3]);
8 ^# Y# A& ~7 ~1 ^! vdisp(eq4);
6 l' e/ W* w6 K( i0 ~. ]9 r/ K+ ub = solve(eq4, b);1 D8 G9 I2 `+ M8 T0 {
H/ U, r5 _* G运行会报错，求各位大佬帮我看看是哪里出错了，感谢！！% e6 T5 z4 M4 a* T2 K- W
|
|
/1 
发表于 2023-2-16 13:45
收藏
淘帖
支持！
反对！