| 查看: 2063 | 回复: 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 |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有293人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31030.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5603
- 在线: 856.8小时
- 虫号: 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
50