24小时热门版块排行榜    

查看: 2961  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

请叫我呼呼呼

铜虫 (小有名气)

[求助] fsolve解方程组有问题 已有2人参与

function ChemEqui
x0=[1,1,1,1,1,1,1,1,1,1];
x=fsolve(@equi,x0);
function f=equi(x)
f1=(x(1)-x(2))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7))-0.959313*(1-2*x(1)-x(3))^2;
f2=(x(2)-x(3))*(x(2)+x(4)+x(7))-1.003773*(x(1)-x(2))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
f3=(x(3)-x(4))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7))-0.964235*(x(2)-x(3))*(1-2*x(1)-x(3));
f4=(x(2)+x(4)+x(7))*(x(4)-x(5))-0.994909*(x(3)-x(4))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
f5=x(5)-x(6)-x(7)-1.002046*(x(4)-x(5));
f6=(x(6)-2*x(8))*x(6)-1.052621*(x(5)-x(6)-x(7))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
f7=x(7)*(x(2)+x(4)+x(7))-1.043431*(x(5)-x(6)-x(7))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
f8=x(8)^2-0.988734*(x(6)-2*x(8))^2;
f=[f1 f2 f3 f4 f5 f6 f7 f8];

上面的写的,计算出来后有警告Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm instead.   请问一下什么情况
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
jjdg: 金币+2, 感谢参与 2017-07-09 11:32:28
1stOpt求解:
CODE:
Parameter x(8);
Function
(x(1)-x(2))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7))-0.959313*(1-2*x(1)-x(3))^2;
(x(2)-x(3))*(x(2)+x(4)+x(7))-1.003773*(x(1)-x(2))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
(x(3)-x(4))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7))-0.964235*(x(2)-x(3))*(1-2*x(1)-x(3));
(x(2)+x(4)+x(7))*(x(4)-x(5))-0.994909*(x(3)-x(4))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
x(5)-x(6)-x(7)-1.002046*(x(4)-x(5));
(x(6)-2*x(8))*x(6)-1.052621*(x(5)-x(6)-x(7))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
x(7)*(x(2)+x(4)+x(7))-1.043431*(x(5)-x(6)-x(7))*(1-x(1)-x(3)+x(2)+x(4)+x(6)+x(7));
x(8)^2-0.988734*(x(6)-2*x(8))^2;

x1: -5.49608040273917E-17
x2: 1.80094386556405E-22
x3: 1
x4: -5.06683228281383E-22
x5: 7.04504319547718E-23
x6: 2.921739051712E-22
x7: 1.98623664044073E-22
x8: 1.53554561331567E-13

x3=1,其它均为0即为最好解。
4楼2017-07-09 09:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2017-07-08 17:25:55
没有使用默认的算法

使用了LM算法

参考
https://cn.mathworks.com/help/optim/ug/fsolve.html
ส็็็็็็็็็็็็็็็็็็็็
2楼2017-07-08 12:32:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

Matlab专家


jjdg: 金币+1, 感谢参与 2017-07-09 11:32:18
本帖仅楼主可见
3楼2017-07-08 21:45:09
已阅   申请程序强帖   回复此楼   编辑   查看我的主页

匿名

用户注销 (初入文坛)

本帖仅楼主可见
5楼2018-11-28 14:39:57
已阅   申请程序强帖   回复此楼   编辑   查看我的主页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 招08考数学 +3 laoshidan 2026-03-20 7/350 2026-03-22 12:54 by 在风落中
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
信息提示
请填处理意见