全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
11268 16
2006-09-12

各位好,请不吝赐教。

%题干(求5个未知数的非线性方程组)

% 已知p0=(bv0+((roe1-r)./(1+r)).*bv0+((roe2-r)./(1+r).^2).*bv1+((roe2-r)./(r.*(1+r).^2)).*bv1;

并假定P0BVjROEj为已知5×1的列向量

p0=[5;4;7;8;9],bv0=[3;2;4;5;7], bv1=[2.5;3;3.5;5.2;6],roe1=[0.1;0.15;0.08;0.07;0.12], roe1=[0.12;0.10;0.10;0.05;0.15],

% 如何求解r(r5×1的列向量)

%编写mycost2000try.m函数

function F = mycost2000try(r);

F=[5;4;7;8;9]-([3;2;4;5;7]+(([0.1;0.15;0.08;0.07;0.12]-r)./(1+r).^2).*[3;2;4;5;7]+(( [0.12;0.10;0.10;0.05;0.15]-r)./(r.*(1+r).^2)).* [2.5;3;3.5;5.2;6];

%在命令行输入如下命令:

one=ones(5);

r0=0.1*one(:,1); % 初始点

options = optimset('Display','off'); % Turn off Display

[r,fval] = fsolve(@mycost2000try,r0, options)

问题1:为什么无法求解?

问题2:本题只是我的一个简化的例子,实际例子p0bv0等应该是(1000×1)的列向量,相应地求解的r也为1000×1的列向量,这时应如何编写和调用m函数,为什么在mycost2000try .m先定义如下函数:

F=(p0-(bv0+((roe1-r)./(1+r)).*bv0+((roe2-r)./(1+r).^2).*bv1+((roe2-r)./(r.*(1+r).^2)).*bv1);

然后在命令行中定义p0bv0等是不可识别的?

请各位不吝赐教。

二维码

扫码加我 拉你入群

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

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

全部回复
2006-9-13 14:48:00

可以以论坛现金作为奖励

二维码

扫码加我 拉你入群

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

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

2006-9-14 05:45:00

纯粹个人建议

能不用fsolve就不用fsolve,更何况你要求这么高的维数。你可以试试把方程变成多项式形式,求多项式的解比较方便快速。

二维码

扫码加我 拉你入群

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

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

2006-9-15 23:39:00

用1stOpt试了一下,不知对否?

Constant p0(1:5)=[5,4,7,8,9];
Constant bv0(1:5)=[3,2,4,5,7];
Constant bv1(1:5)=[2.5,3,3.5,5.2,6];
Constant roe1(1:5)=[0.1,0.15,0.08,0.07,0.12];
Constant roe2(1:5)=[0.12,0.10,0.10,0.05,0.15];
Function For(i=1:5)(p0=(bv0+((roe1-r)/(1+r))*bv0+((roe2-r)/(1+r)^2)*bv1+((roe2-r)/(r*(1+r)^2))*bv1));

结果:

r1: 0.0662098312457287
r2: 0.0606946081232086
r3: 0.0534039089862346
r4: 0.0320727048208649
r5: 0.110382787884127

二维码

扫码加我 拉你入群

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

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

2006-9-16 10:05:00

dingd,你这个程序如何运行?是m文件脚本吗?1stOpt是什么意思?不好意思,由于学的主要是SASSTATAmatlab刚学,还请赐教。不知是否已在银行开户,我将以100论坛现金略表心意。

二维码

扫码加我 拉你入群

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

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

2006-9-17 10:16:00
1stOpt是个国产优化、拟合、数值分析计算软件,网址:http://www.7d-soft.cn。1stOpt的优化拟合、方程求解等能力远强于Matlab、SAS、SPSS、Orgin等,建议试试!
二维码

扫码加我 拉你入群

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

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

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

分享

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