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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 解带参数方程，参数在变
71/1返回列表
查看: 2285  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 21:13 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 21:07 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:56 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:34 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:29 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 5/250 2026-11-10 20:01 by Z7qfCsdIiD1T
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 5/250 2026-11-10 19:56 by Z7qfCsdIiD1T
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭