全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
11537 9
2007-04-04
编程方面的知识较弱,请大侠们指教,问题如下:
一个非线性方程组如下:
y1 = u*s*w1+(1-u)/2;
y2 = u*(1-s)*w2+(1-u)/2;
g1 = (s*w1^(1-d)+(1-s)*(w2*t)^(1-d))^(1/(1-d));
g2 = (s*(w1*t)^(1-d)+(1-s)*w2^(1-d))^(1/(1-d));
w1 = (y1*g1^(d-1)+y2*g2^(d-1)*t^(1-d))^(1/d);
w2 = (y1*g1^(d-1)*t^(1-d)+y2*g2^(d-1))^(1/d);
o1 = w1*g1^(-u);
o2 = w2*g2^(-u);
在给定参数d,u,t和s时可利用1stopt软件求解,我的做法如下:
constant d=5, u=0.4, t=1.7, s=0.2;
Parameters y1, y2, g1, g2, w1, w2, o1, o2;
Function y1 = u*s*w1+(1-u)/2;
y2 = u*(1-s)*w2+(1-u)/2;
g1 = (s*w1^(1-d)+(1-s)*(w2*t)^(1-d))^(1/(1-d));
g2 = (s*(w1*t)^(1-d)+(1-s)*w2^(1-d))^(1/(1-d));
w1 = (y1*g1^(d-1)+y2*g2^(d-1)*t^(1-d))^(1/d);
w2 = (y1*g1^(d-1)*t^(1-d)+y2*g2^(d-1))^(1/d);
o1 = w1*g1^(-u);
o2 = w2*g2^(-u);
现在想考察方程组的解(如o1或o2)与参数s之间的关系,并想将其绘成曲线图的形式,则需要对s在[0,1]区间上选取若干值并在这些取值下分别对方程组求解,手动计算方式很繁琐,请高手指教,1stopt软件中是否有相应的循环命令可以实现对方程组在不同参数下的求解呢?具体应该如何操作?谢谢。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2007-4-7 19:38:00

在新版中可以实现,S从0到1,增幅0.1:

constant d=5, u=0.4, t=1.7;
LoopConstant s=[0:0.1:1];
Parameters y1, y2, g1, g2, w1, w2, o1, o2;
Function y1 = u*s*w1+(1-u)/2;
y2 = u*(1-s)*w2+(1-u)/2;
g1 = (s*w1^(1-d)+(1-s)*(w2*t)^(1-d))^(1/(1-d));
g2 = (s*(w1*t)^(1-d)+(1-s)*w2^(1-d))^(1/(1-d));
w1 = (y1*g1^(d-1)+y2*g2^(d-1)*t^(1-d))^(1/d);
w2 = (y1*g1^(d-1)*t^(1-d)+y2*g2^(d-1))^(1/d);
o1 = w1*g1^(-u);
o2 = w2*g2^(-u);

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-4-8 21:45:00
感谢指教。。不过我试了一下好像不行,1stopt把s也当成参数解出来了。不知道是不是我的软件版本的问题?我用的是1.5版,请问一下您说的新版本是哪一版呢?谢谢
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-4-28 17:52:00

新版本好像是2.0

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-4-28 22:02:00

111806.pdf
大小:(72.99 KB)

 马上下载


其实用lingo和matlab综合起来做,很简单

[此贴子已经被作者于2007-4-28 22:29:30编辑过]

附件列表

111796.pdf

大小:18.73 KB

 马上下载

111797.pdf

大小:18.73 KB

 马上下载

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-4-28 22:50:00

太厉害啦,能不能把程序贴出来啊,学习一下

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群