|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我的代码：
/ o/ ~8 |2 H9 a! Gclear all
5 o8 O8 J$ E) x/ Lsyms xc Rw xt yt b h C; % 定义符号: s, x; w) \4 f4 x) A6 f
eqn = (xt-xc)^2+(Rw+4*h*xt^2/b^2-h)^2-Rw^2 == 0; % 表达式+ ] U: |+ L3 K; ~/ ]7 W* O
col = collect(eqn) % 合并同类项. p/ t7 B/ p# [* y
sol = solve(col,xt); % 求解以xt为自变量得形式+ h7 k) {+ |' [
s = vpa(sol,3) % 限制显示精度
3 r5 O# j) { J! r- ~. Ls = roots(sol) % 显示解
! Q. C+ H! H& O. ~+ X( H9 | D0 G7 x4 B) }
我的想法：" f" K( L3 v; X: L
将多项式(xt-xc)^2+(Rw+4*h*xt^2/b^2-h)^2-Rw^2 == 0,将xt提到前面表示为xt=...的形式
+ X5 z) k0 Q! j; J. W I结果显示太长，请问什么原因，还是本来解就这么多，有没有简化的办法，谢谢) _2 H( a# |9 }
|
|
/1 
发表于 2021-3-24 17:44
收藏
淘帖
支持！
反对！