找回密码
 注册
关于网站域名变更的通知
查看: 643|回复: 3
打印 上一主题 下一主题

使用roots()显示求解后的多项式,显示输出已经截断

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-3-24 17:44 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
我的代码:  H4 e& ~: D/ n9 q9 G
clear all
- _( [2 E" z" @( i" A4 w: csyms xc Rw xt yt b h C; % 定义符号6 }' C8 ]1 v4 {$ G* u0 b$ p% }: g; {
eqn = (xt-xc)^2+(Rw+4*h*xt^2/b^2-h)^2-Rw^2 == 0; % 表达式
! q# c( v/ f7 }2 v! Q4 {col = collect(eqn)            % 合并同类项
9 U5 ]0 H' d9 ]+ u0 I+ Psol = solve(col,xt);          % 求解以xt为自变量得形式& _' J! c8 H0 g; @# P
s = vpa(sol,3)                % 限制显示精度
6 p4 i' a9 c+ f9 l1 F6 n6 ss = roots(sol)                % 显示解6 c7 o! Z  |3 B7 r$ _
& h$ m& i" F9 d8 L5 r
我的想法:
) Q9 s3 W* t1 K3 B6 g将多项式(xt-xc)^2+(Rw+4*h*xt^2/b^2-h)^2-Rw^2 == 0,将xt提到前面表示为xt=...的形式) ?4 @! ^7 ]$ |; p
结果显示太长,请问什么原因,还是本来解就这么多,有没有简化的办法,谢谢# M4 q) O$ |' D, Z$ ~

该用户从未签到

2#
发表于 2021-3-24 18:10 | 只看该作者
你的写法值得商榷。首先,高次多项式构造的方程,求解析解很困难(五次以上没有求根公式,你这个是四次的也很麻烦);然后,你这个方程里面很多变量都没有给定具体数值,用vpa也顶多在系数、幂次上的有确切数值的地方有些效果,对于未赋值的那些变量是不会有效果的;最后,你写的这个roots(sol)根本就不对,这等于是"将sol(实际为一元四次方程的四个根)依次作为某一元三次多项式的系数,并求解该一元三次方程的根",这个复杂系数的方程的根,除非某些特例,只会更加复杂

该用户从未签到

3#
发表于 2021-3-24 19:02 | 只看该作者
来学习一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-23 22:49 , Processed in 0.140625 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表