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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 解带参数方程，参数在变
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 1943  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

3209的烟屁

金虫 (小有名气)

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

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个结果？还有可能变三个量，所以想请教怎么办？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-27 10:44:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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的烟屁

金虫 (小有名气)
引用回帖:
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的烟屁

金虫 (小有名气)
引用回帖:
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的回帖

tracuer2011

木虫 (小有名气)
★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3209的烟屁

金虫 (小有名气)
引用回帖:
5楼: Originally posted by tracuer2011 at 2013-04-30 00:09:22
M01=30;
M02=30;
M03=5;
M04=5;
x0=;
options=optimset('Display','off');
for M05=0:30
     M06=30-M05;
     options=optimset('Display','off');
     =fsolve(@fafun,x0,options)
     Fval(M05+1)= ...

试了一下，还是不行。M01-m06是定义在M文件中的，这样在命令窗口给行吗，是不是这个循环要放在M文件的代码中啊？
赞一下 回复此楼
6楼2013-04-30 17:08:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图
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的回帖
相关版块跳转 我要订阅楼主 3209的烟屁 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-05 9/450 2026-04-05 22:24 by dongzh2009
[考研] 338求调剂 +3 我想上岸ii 2026-04-05 3/150 2026-04-05 19:59 by nepu_uu
[考研] 考研调剂生寻找导师 +3 顾瞻考研啊 2026-04-05 3/150 2026-04-05 18:18 by 啵啵啵0119
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 0703化学 +12 妮妮ninicgb 2026-04-04 13/650 2026-04-05 10:46 by 啊俊！
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 调剂求助 +10 想换手机不想解� 2026-04-02 13/650 2026-04-05 09:41 by sam3303
[考研] 341求调剂 +3 学无止境，冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 085600，321分求调剂 +10 大馋小子 2026-04-04 11/550 2026-04-05 08:25 by 544594351
[考研] 311分 22408 求调剂 +3 bing_bot 2026-04-03 3/150 2026-04-05 00:43 by chongya
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 282求调剂 不挑专业 求收留 +7 Yam. 2026-03-30 8/400 2026-04-03 14:12 by zhangdingwa
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 材料化工340求调剂 +5 jhx777 2026-03-30 5/250 2026-04-02 12:45 by smileboy2006
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
信息提示
关闭
请填处理意见
关闭