| 查看: 1931 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
解带参数方程，参数在变
|
|||
|
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个结果？还有可能变三个量，所以想请教怎么办？ |
回复此楼» 猜你喜欢
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有208人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有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人回复
|
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 |
7楼2013-05-02 15:05:46
2楼2013-04-27 17:15:11
3楼2013-04-28 11:30:48
4楼2013-04-28 16:13:46
20