| 查看: 2016 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者yangboa007将赠送您 98 个金币 | ||
yangboa007铁杆木虫 (职业作家)
|
[求助]
matlab求二元三次方程 已有1人参与
|
|
|
有一组二元三次方程，想要求得a和ε的值，用matlab运算，发现只有复数解，由于a和ε本身可以允许有误差，方程中6770.132和3941.435也可以有正负500的误差，所以我想用matlab求a和ε的约数，a的值在0~1之间、ε的值在1000~3000之间，只要能代入方程偏差不大即可。相当于这组数无限趋近，能在我的接受范围。还请大家赐教，谢谢！ 以下是二元三次方程（绿色图片为二元三次方程组）： matlab程序代码如下，用b代替ε： [a,b]=solve('6770.132=(4.17*a^2*(1-b)^2/b^3*1.4*10^-5*0.25)+(0.29*a*(1-b)/b^3*8.37*0.25^2)','3941.435=(4.17*a^2*(1-b)^2/b^3*1.4*10^-5*0.196)+(0.29*a*(1-b)/b^3*8.703*0.196^2)','a','b') 得出的结果如下： 警告: support of character vectors that are not valid variable names or define a number will be removed in a future release. to create symbolic expressions, first create symbolic variables and then use operations on them. > in sym>convertexpression (line 1559) in sym>convertchar (line 1464) in sym>tomupad (line 1216) in sym (line 179) in solve>geteqns (line 405) in solve (line 225) 警告: support of character vectors that are not valid variable names or define a number will be removed in a future release. to create symbolic expressions, first create symbolic variables and then use operations on them. > in sym>convertexpression (line 1559) in sym>convertchar (line 1464) in sym>tomupad (line 1216) in sym (line 179) in solve>geteqns (line 405) in solve (line 225) 警告: do not specify equations and variables as character vectors. instead, create symbolic variables with syms. > in solve>geteqns (line 445) in solve (line 225) a = -2153.8654367915894618926415071084 - 3123.0895646907707350229180631831 - 1192.7113132876986000825947006744i - 3123.0895646907707350229180631831 + 1192.7113132876986000825947006744i b = - 0.36131483557811412999971700246919 - 1.1551442818685524875109192327699e-40i 0.1806574177890570649998585012346 + 0.31290782637484434808955229418758i 0.1806574177890570649998585012346 - 0.31290782637484434808955229418758i  无标题.jpg |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有59人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有7人回复
-
散金币祈福
已经有94人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
yangboa007: 金币+2, ★有帮助, 感谢，我试过了，求出来的是负数。 2020-08-30 14:55:17
感谢参与，应助指数 +1
yangboa007: 金币+2, ★有帮助, 感谢，我试过了，求出来的是负数。 2020-08-30 14:55:17
|
尽量避免符号求解，选择数值求解，用 fsolve 函数试试 发自小木虫Android客户端 |
2楼2020-08-29 11:15:09
40