| 查看: 1913 | 回复: 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年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有274人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有29人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有53人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有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人回复
3楼2013-04-28 11:30:48
2楼2013-04-27 17:15:11
4楼2013-04-28 16:13:46
tracuer2011
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1745.2
- 红花: 1
- 帖子: 223
- 在线: 93.2小时
- 虫号: 1260856
- 注册: 2011-04-10
- 性别: GG
- 专业: 光学
★
csgt0: 金币+1, 谢谢 2013-05-02 09:10:58
csgt0: 金币+1, 谢谢 2013-05-02 09:10:58
|
M01=30; M02=30; M03=5; M04=5; x0=[0.4;0.2;0.007;0.005;0.16;0.14;0.4]; options=optimset('Display','off'); for M05=0:30 M06=30-M05; options=optimset('Display','off'); [x,fval,exitflag,output,Jacobian]=fsolve(@fafun,x0,options) Fval(M05+1)=fval; %Fval 即为所求数组 end %%类似于这样行不行 |
5楼2013-04-30 00:09:22
