24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Q莹儿

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 192.2
散金: 20
红花: 1
帖子: 65
在线: 22.1小时
虫号: 4508623
注册: 2016-03-15
专业: 金属材料的微观结构

[求助] MATLAB求解复杂方程组，求助！！！

MATLAB新手。想用MATLAB求解复杂二元多次方程组。方程式如下所述：

d/x=-1.131*[log(a/x)]^3+13.635*[log(a/x)^2]-30.594[[log(a/x)]]+29.267；
b/ac=(-1.40557y^3+0.77526y^2+0.15830y-0.06831)[ log(a/x)]^3+(17.93006y^3-9.22091y^2-2.37733n+0.86295[[log(a/x)] ^2+(-79.99715y^3+40.55620y^2+9.00157y-2.54543[[log(a/x)]+(122.65069y^3-63.88418y^2-9.58936y+6.20045)
其中，a,b,c,d为参量，他们的数值在Excel表格中，x,y为未知数，我想利用MATLAB，求解出上述方程的解。在a,b,c,d参量变化的情况下，求解出x,y的值。
我看到网上说可以使用m文件自定义函数，也可以使用solve语句，但是我用solve语句写出来如下，
syms x y a b c d
a= xlsread('20190718.xlsx','sheet2','B2:B6');
b= xlsread('20190718.xlsx','sheet2','C2:C6');
c= xlsread('20190718.xlsx','sheet2','A2:A6');
d= xlsread('20190718.xlsx','sheet2','D2

6');
[x,y]=solve('-1.131*[log(a/x)]^3+13.635*[log(a/x)^2]-30.594[[log(a/x)]]+29.267==d/x,b/ac==(-1.40557y^3+0.77526y^2+0.15830y-0.06831)[ log(a/x)]^3+(17.93006y^3-9.22091y^2-2.37733n+0.86295[[log(a/x)] ^2+(-79.99715y^3+40.55620y^2+9.00157y-2.54543[[log(a/x)]+(122.65069y^3-63.88418y^2-9.58936y+6.20045)',x,y)
x=simple(x),y=simple(y)

系统提示'"对无效变量名或定义数字的字符串的支持将在将来的版本中删除。创建符号表达式，首先创建符号变量，然后对其使用操作。""
我不知道该怎么办了，求助，非常感谢！

回复此楼

» 猜你喜欢

求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有243人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复

1楼2019-07-22 13:09:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

你的a b c d每个都是包含5个数的数组你是想对于每组a b c d都解一次方程组一共是解5次对吧？
提供几点思路：
1 不要用solve函数这个函数不适合求复杂方程的数值解改用fsolve函数
2 由于a b c d有5组如果想一次求出所有的解需要写个循环比如用for i=1:5
上述两点我以前写过类似的代码回在别人的帖子里具体翻一下我回复记录
3 MATLAB的fsolve函数需要提供初值如果求解困难可以试试1stOpt软件来解

发自小木虫Android客户端