【原创】变参数非线性方程组的求解！

studyboy

对于求解非线性方程组一般用fsolve命令就可以了，但是对于方程组中某一系数是变化的，该怎么求呢？

%定义方程组如下，其中k为变量
function F = myfun(x,k)
H=0.32; Pc0=0.23; W=0.18;
F=[Pc0+H*(1+1.5*(x(1)/W-1)-0.5*(x(1)/W-1)^3)-x(2);  x(1)-k*sqrt(x(2))];

%求解过程
H=0.32; Pc0=0.23; W=0.18;
x0 = [2*W; Pc0+2*H];           % 取初值
options = optimset('Display','off');
k=0:0.01:1;             % 变量取值范围[0 1]
for i=1:1:length(k)
kk=k(i);                        
x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组
x1(i)=x(1); x2(i)=x(2);
end
plot(k,x1,'-b',k,x2,'-r'); xlabel('k'); legend('x1','x2')

很简单，大家见笑了！

[ 本帖最后由 ChaChing 于 2009-9-9 22:09 编辑 ]

andy_3656

不用谦虚哦，大家到这里来就是相互学习的啊！

无水1324

不要谦虚了
大哥太感谢了
我也作了一个，但是总是出问题，好好研究一下你的，虽然别人说太简单了
但是对于我这种初学，不深的人来说，
简直太好了。。。

x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组

这个怎么解释，有点问题哦

studyboy

不好意思，前几天不在。
其实这一问题的关键求解就在于：
x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组。
这里把kk作为第二个参数代入非线性方程组，然后求得某一kk值时的解，释放kk。然后继续循环kk。这样就可以得到整个kk取值区间上的解了。
具体内容可以参考Matlab help中的fslove中 的变参数部分。

jaguar

x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组。
这句话好像matlab 6.5 里通不过呀，
Line: 15 Column: 13
"identifier" expected, "(" found.
能否请解释一下

qqzjyywx

??? Error: File: D:\MATLAB6p5\work\Untitled2222.m Line: 16 Column: 13
"identifier" expected, "(" found.
同样的问题怎么回事 阿？

studyboy

我是在7的环境下面运行的，具体的6.5的情况不太清楚。不好意思。

welcoemjimy

不错 我是刚学的  半天纔農明白

[ 本帖最后由 无水1324 于 2007-7-23 12:06 编辑 ]

wugang326

如果KK值事先不能确定,都要靠前一次的计算结果迭代产生,应该怎么办啊

suffer

原帖由 jaguar 于 2006-8-7 21:33 发表
x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组。
这句话好像matlab 6.5 里通不过呀，
Line: 15 Column: 13
"identifier" expected, "(" found.
能否请解释一下

6.5的调用格式有点不一样，只用修改一下这一句就行了

1. x = fsolve(@myfun, x0, options,kk);%求解非线性方程组

复制代码

suffer

原帖由 wugang326 于 2006-10-19 18:07 发表
如果KK值事先不能确定,都要靠前一次的计算结果迭代产生,应该怎么办啊

这里再把完整的程序帖一下，实际上只是改动了一个句子：

1. H=0.32; Pc0=0.23; W=0.18;
   
2. x0 = [2*W; Pc0+2*H];           % 取初值
   
3. options = optimset('Display','off');
   
4. k=0:0.01:1;             % 变量取值范围[0 1]
   
5. for i=1:1:length(k)
   
6. kk=k(i);                        
   
7. x = fsolve(@myfun, x0, options,kk);%求解非线性方程组
   
8. x1(i)=x(1); x2(i)=x(2);
   
9. end
   
10. plot(k,x1,'-b',k,x2,'-r'); xlabel('k'); legend('x1','x2')
    

复制代码

1. function F = myfun(x,k)
   
2. H=0.32; Pc0=0.23; W=0.18;
   
3. F=[Pc0+H*(1+1.5*(x(1)/W-1)-0.5*(x(1)/W-1)^3)-x(2);
   
4.     x(1)-k*sqrt(x(2))];

复制代码

[ 本帖最后由 ChaChing 于 2009-2-26 13:58 编辑 ]

michael1003

将k设为全局变量，应该也可以解决这个问题吧。

chbo76

好，不错。
谢谢。:lol

花如月

好文章，俺给顶起来

huanhuan023

为什么同样的程序到我机子上就不能运行呢？
是我装的软件的问题吗？疑惑中……我的也是7.0版本……