版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

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

3209的烟屁

金虫 (小有名气)

应助: 3 (幼儿园)
金币: 789.3
红花: 1
帖子: 131
在线: 57.4小时
虫号: 1003622
注册: 2010-04-23
专业: 钢铁冶金

[求助] 解带参数方程，参数在变

function f=fafun(x)
syms M01 M02 M03 M04 M05 M06 T N0 R
M01=30;
M02=30;
M03=5;
M04=5;
M05=0;
M06=30;
M11=M01/60;
M12=M02/56;
M13=M03/102;
M14=M04/40.3;
M15=M05/62;
M16=M06/78;
M21=M11/(M11+M12+M13+M14+M15+M16);
M22=M12/(M11+M12+M13+M14+M15+M16);
M23=M13/(M11+M12+M13+M14+M15+M16);
M24=M14/(M11+M12+M13+M14+M15+M16);
M25=M15/(M11+M12+M13+M14+M15+M16);
M26=M16/(M11+M12+M13+M14+M15+M16);
N31=0.5*M21/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N32=0.6875*M22/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N33=0.3542*M23/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N34=0.4583*M24/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N35=0.6736*M25/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N36=0.7444*M26/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
T=1623.15;
N11=(243.2+0.031*T)/1000;
N12=(791-0.0935*T)/1000;
N13=(1024-0.177*T)/1000;
N14=(1770-0.636*T)/1000;
N15=(438-0.116*T)/1000;
N16=(1604.6-0.72*T)/1000;
N0=6.02*10^23;
R=8.314;
N21=R*T/( N0^(1/3)*(27.516*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N22=R*T/( N0^(1/3)*(20.7*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N23=R*T/(N0^(1/3)*(28.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N24=R*T/(N0^(1/3)*(16.1*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N25=R*T/(N0^(1/3)*(33.0*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N26=R*T/(N0^(1/3)*(31.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
f=[N31*exp((x(7)-N11)/N21)-x(1);
N32*exp((x(7)-N12)/N22)-x(2);
N33*exp((x(7)-N13)/N23)-x(3);
N34*exp((x(7)-N14)/N24)-x(4);
N35*exp((x(7)-N15)/N25)-x(5);
N36*exp((x(7)-N16)/N26)-x(6);
N31*exp((x(7)-N11)/N21)+N32*exp((x(7)-N12)/N22)+N33*exp((x(7)-N13)/N23)+N34*exp((x(7)-N14)/N24)+N35*exp((x(7)-N15)/N25)+N36*exp((x(7)-N16)/N26)-1];
M文件如上所示，我想让M05+M06=30，所以就是M05=30,29，28....不想一个一个去改，有没有类似for M05=30 if M05>1 M05=M05-1这样的办法直接能输出30个结果？还有可能变三个量，所以想请教怎么办？

回复此楼

» 猜你喜欢

瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有130人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
物理学I论文润色/翻译怎么收费? 已经有124人回复
PVA溶解已经有2人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复
什么时候开奖？已经有6人回复
散金币祈福已经有94人回复

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

微分方程组参数拟合的问题（多参数动力系统参数估计）已经有11人回复
如何用得到的吸附数据进行Langmuir-Freundlich方程拟合，而得出相关参数的具体值已经有5人回复
请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。已经有17人回复
求高版本1stopt，解决一个问题，微分方程组参数拟合已经有3人回复
动力学方程参数估计方法已经有14人回复
用matlab解带参数的一元三次方程--求解已经有3人回复
拜求，这种带参数的方程怎么解比较好？已经有4人回复
origin8.0拟合方程中参数怎么求已经有5人回复
一阶微分方程求解，并优化参数值已经有14人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
matlab拟合方程参数时初值的选择已经有15人回复
matlab解微分方程已经有10人回复
matlab-常微分方程参数估计已经有12人回复
求回归方程参数估计的最优化算法matlab代码已经有9人回复
求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊已经有5人回复
请教Matlab高手解多元非线性方程组中的参数已经有6人回复
【求助】matlab 求解微分方程中的未知参数已经有20人回复
求非线性方程中的参数问题已经有13人回复
【求助】如何从Murnaghan方程的结果中得到优化的晶格参数？已经有11人回复
【求助】不仅控制方程里含有未知参数，边界条件里也有已经有5人回复

1楼 2013-04-27 10:44:04

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

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
计算强帖: 2
应助: 367 (硕士)
贵宾: 0.457
金币: 7812
散金: 8897
红花: 20
沙发: 7
帖子: 2810
在线: 515.9小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

function f3209
global M05 M06
options=optimset('Display','off');
x0=[0.4;0.2;0.007;0.005;0.16;0.14;0.4];
for i=1:31
M05=i-1;
M06=31-i;
[x(:,i),fval(:,i),exitflag(i),output(:,i),Jacobian(i,:,

]=fsolve(@fafun,x0,options);
end
x
fval
exitflag
output
end

function f=fafun(x)
global M05 M06
M01=30;
M02=30;
M03=5;
M04=5;
M11=M01/60;
M12=M02/56;
M13=M03/102;
M14=M04/40.3;
M15=M05/62;
M16=M06/78;
M21=M11/(M11+M12+M13+M14+M15+M16);
M22=M12/(M11+M12+M13+M14+M15+M16);
M23=M13/(M11+M12+M13+M14+M15+M16);
M24=M14/(M11+M12+M13+M14+M15+M16);
M25=M15/(M11+M12+M13+M14+M15+M16);
M26=M16/(M11+M12+M13+M14+M15+M16);
N31=0.5*M21/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N32=0.6875*M22/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N33=0.3542*M23/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N34=0.4583*M24/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N35=0.6736*M25/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N36=0.7444*M26/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
T=1623.15;
N11=(243.2+0.031*T)/1000;
N12=(791-0.0935*T)/1000;
N13=(1024-0.177*T)/1000;
N14=(1770-0.636*T)/1000;
N15=(438-0.116*T)/1000;
N16=(1604.6-0.72*T)/1000;
N0=6.02*10^23;
R=8.314;
N21=R*T/( N0^(1/3)*(27.516*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N22=R*T/( N0^(1/3)*(20.7*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N23=R*T/(N0^(1/3)*(28.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N24=R*T/(N0^(1/3)*(16.1*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N25=R*T/(N0^(1/3)*(33.0*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N26=R*T/(N0^(1/3)*(31.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
f=[N31*exp((x(7)-N11)/N21)-x(1);
N32*exp((x(7)-N12)/N22)-x(2);
N33*exp((x(7)-N13)/N23)-x(3);
N34*exp((x(7)-N14)/N24)-x(4);
N35*exp((x(7)-N15)/N25)-x(5);
N36*exp((x(7)-N16)/N26)-x(6);
N31*exp((x(7)-N11)/N21)+N32*exp((x(7)-N12)/N22)+N33*exp((x(7)-N13)/N23)+N34*exp((x(7)-N14)/N24)+N35*exp((x(7)-N15)/N25)+N36*exp((x(7)-N16)/N26)-1];
end

赞一下

回复此楼

showmethemoney

7楼2013-05-02 15:05:46

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

查看全部 7 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
计算强帖: 2
应助: 367 (硕士)
贵宾: 0.457
金币: 7812
散金: 8897
红花: 20
沙发: 7
帖子: 2810
在线: 515.9小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
3209的烟屁: 金币+10, ★★★很有帮助 2013-04-28 11:28:23

1、用矩阵M06=0:30,M01=30-M06;后面都加上点
2、用循环，用数组存储结果

赞一下

回复此楼

showmethemoney

2楼2013-04-27 17:15:11

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

3209的烟屁

金虫 (小有名气)

应助: 3 (幼儿园)
金币: 789.3
红花: 1
帖子: 131
在线: 57.4小时
虫号: 1003622
注册: 2010-04-23
专业: 钢铁冶金

引用回帖:

2楼: Originally posted by csgt0 at 2013-04-27 17:15:11
1、用矩阵M06=0:30,M01=30-M06;后面都加上点
2、用循环，用数组存储结果

大哥不要生气，我想再问一下具体应该怎么写，自己改了几下也不行，希望把改的地方写一下好吗，跪谢！

赞一下

回复此楼

3楼2013-04-28 11:30:48

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

3209的烟屁

金虫 (小有名气)

应助: 3 (幼儿园)
金币: 789.3
红花: 1
帖子: 131
在线: 57.4小时
虫号: 1003622
注册: 2010-04-23
专业: 钢铁冶金

引用回帖:

2楼: Originally posted by csgt0 at 2013-04-27 17:15:11
1、用矩阵M06=0:30,M01=30-M06;后面都加上点
2、用循环，用数组存储结果

M01=30;
M02=30;
M03=5;
M04=5;
for M05=0:30
if M05<=30
M05=M05+1;
end
M06=30-M05;
这是我改的，运行不了啊，问题是这些代码是M文件中的，解这个方程还需要在命令窗口中输入
x0=[0.4;0.2;0.007;0.005;0.16;0.14;0.4];
options=optimset('Display','off');
[x,fval,exitflag,output,Jacobian]=fsolve(@fafun,x0,options)
我需要改变M01-M05得到很多数据点，这样一个一个来太慢了

赞一下

回复此楼

4楼2013-04-28 16:13:46

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

查看全部 7 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 估计是周四 +3	archvillain 2026-08-18	3/150	2026-08-21 01:48 by jnhyjjm
[基金申请] 基金啊基金 +3	longfie172 2026-08-20	3/150	2026-08-21 00:04 by wuchongjun
[基金申请] 科研孤儿太难了 +14	我4大白菜 2026-08-20	15/750	2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧！ +8	大萍1987 2026-08-20	10/500	2026-08-20 21:48 by 大萍1987
[基金申请] 时间戳又变了 +12	wuchongjun 2026-08-20	17/850	2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护，明天很可能放榜！ +10	zju2000 2026-08-16	11/550	2026-08-20 20:02 by cl479861084
[基金申请] 只有每年这种时候来逛逛小木虫 +21	yaoyewhu2008 2026-08-20	22/1100	2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3	snowwithsea 2026-08-19	12/600	2026-08-20 15:56 by lfy8008
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5	瞬息宇宙 2026-08-15	5/250	2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4	yuleib84 2026-08-19	5/250	2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14	布布和一二 2026-08-19	15/750	2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10	wangze12014 2026-08-14	12/600	2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4	奶牛小黑 2026-08-15	10/500	2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6	Tide man 2026-08-14	7/350	2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6	CrisMessi 2026-08-18	6/300	2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜？ +5	Shxjjxjkx 2026-08-18	5/250	2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4	wwm09 2026-08-17	6/300	2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14	archvillain 2026-08-15	26/1300	2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4	且听虎啸 2026-08-17	5/250	2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7	Tide man 2026-08-14	8/400	2026-08-15 01:11 by arzu_hma