24小时热门版块排行榜    

查看: 1906  |  回复: 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

matlab求二元三次方程
无标题.jpg
回复此楼
坚持就是胜利
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
yangboa007: 金币+2, 有帮助, 感谢,我试过了,求出来的是负数。 2020-08-30 14:55:17
尽量避免符号求解,选择数值求解,用 fsolve 函数试试

发自小木虫Android客户端
数值计算
2楼2020-08-29 11:15:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangboa007 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 7/350 2026-03-17 20:20 by 花125533
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见