matlab+微分方程组的解法_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
文档贡献者
评价文档：
喜欢此文档的还喜欢
matlab+微分方程组的解法
m​a​t​l​a​b​微​分​方​程​组​的​解​法
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
大小：586.50KB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢第七讲Matlab在求解优化问题中的应用(2)[1]_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
文档贡献者
评价文档：
喜欢此文档的还喜欢
第七讲Matlab在求解优化问题中的应用(2)[1]
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
大小：708.50KB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢求Matlab大神！我用fsolve函数求解六元非线性方程组，但是方程组中的系数是通过其它函数求出，然后通过全_百度知道
求Matlab大神！我用fsolve函数求解六元非线性方程组，但是方程组中的系数是通过其它函数求出，然后通过全
但是调用fsolve函数所求的结果却不对！，是什么原因啊！求解~~不胜感激，经验证！！，我所求的系数确实传递到这个函数中了！局变量global变量传递过来的
提问者采纳
初值的问题
那请问一下，初值的设定是不是要尽可能地接近实际值啊？
提问者评价
其他类似问题
fsolve的相关知识
为您推荐：
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁&& 查看话题
matlab fsolve得出的解全部为复数解
我用matlab中的FSOLVE解一个六元非线性方程组，为什么得出的解全部为复数解，可是理论上应该解为实数，本人Matlab新手，代码如下，大神们帮忙看一下吧！
function=fsolveyalinjieliu
d=0.045;D=0.05;k=1.4;ff=0.004;L=100;P1==300;Pa=101325;M=29;R=8.314;
x0=unifrnd(0,1,1,6);
options=optimset('MaxFunEvals',20000,'MaxIter',5000);
=fsolve(@f,x0,options);
function F=f(x)
似乎是5个方程6个未知数，无穷组解。 有些虚部非常小的，应该是误差，其实就是实数解，忽略虚部即可。 : Originally posted by dingd at
似乎是5个方程6个未知数，无穷组解。 本来是六个方程的，由于得出的是复数解，然后我就删掉了一个试试:D : Originally posted by nono2009 at
有些虚部非常小的，应该是误差，其实就是实数解，忽略虚部即可。 这样可以么 但平方应该是小于零的吧 试着改改初始值，fsolve的算法貌似有些问题，很容易陷入局部最优。如果初始值改变，解也跟着改变，就表示算法陷入局部最优了。这时候就可以考虑用其他的函数试试。 : Originally posted by duyoudu at
试着改改初始值，fsolve的算法貌似有些问题，很容易陷入局部最优。如果初始值改变，解也跟着改变，就表示算法陷入局部最优了。这时候就可以考虑用其他的函数试试。 谢谢你的回答，解非线性方程组可以用solve函数么？ : Originally posted by xiangzigood at
本来是六个方程的，由于得出的是复数解，然后我就删掉了一个试试:D... 把完整方程都贴出来看看。 : Originally posted by xiangzigood at
谢谢你的回答，解非线性方程组可以用solve函数么？... 应该是没问题的。fsolve的算法是个通用算法，可以解答线性和非线性的，虽然效率不如专门的求解函数高。但这个局部最优的问题是存在的，可以考虑函数工具箱中提供的全局求解的函数，具体是哪个你可以看看帮助。 : Originally posted by dingd at
把完整方程都贴出来看看。... function=fsolveyalinjieliu
d=0.04;D=0.05;k=1.4;ff=0.004;L=;Pa=101325;M=29;R=8.314;
x0=unifrnd(0,1,1,6);
options=optimset('MaxFunEvals',20000,'MaxIter',5000);
=fsolve(@f,x0,options);
function F=f(x)
这个应该比较完整，六个方程六个未知数 : Originally posted by duyoudu at
应该是没问题的。fsolve的算法是个通用算法，可以解答线性和非线性的，虽然效率不如专门的求解函数高。但这个局部最优的问题是存在的，可以考虑函数工具箱中提供的全局求解的函数，具体是哪个你可以看看帮助。... 全局求解？您能给推荐一个不，刚解除matlab不久。。 : Originally posted by duyoudu at
应该是没问题的。fsolve的算法是个通用算法，可以解答线性和非线性的，虽然效率不如专门的求解函数高。但这个局部最优的问题是存在的，可以考虑函数工具箱中提供的全局求解的函数，具体是哪个你可以看看帮助。... 您能给推荐一个全局求解的函数不？刚解除matlab不久，不是很熟悉 用1stOpt吧，其它软件得到正解的概率几乎为0：
Parameter x(6);
Constant d=0.045,D1=0.05,k=1.4,ff=0.004,L=100,P1==300,Pa=101325,M=29,R=8.314;
Function ((k+1)/2)*ln((x2^2*x3)/(x1^2*x4))-(1/(x1^2)-1/(x2^2))+4*ff*L/D1;
& && &&&(D1^2/D1^2)*x5*sqrt((2*k*M)*((Pa/x5)^(2/k)-(Pa/x5)^((k+1)/k))/(R*x6*(k-1)))-x2*x5*sqrt((k*M)/(R*x6));
& && &&&(D1^2/D1^2)*x5*sqrt((2*k*M)*((Pa/x5)^(2/k)-(Pa/x5)^((k+1)/k))/(R*x6*(k-1)))-x1*P1*sqrt((k*M)/(R*T1));
& && &&&1+((k-1)*x1^2)/2-x3;
& && &&&1+((k-1)*x2^2)/2-x4;
& && &&&x6-(x3*T1)/x4;
x5: 370304
x6: 291. : Originally posted by dingd at
用1stOpt吧，其它软件得到正解的概率几乎为0：
Parameter x(6);
Constant d=0.045,D1=0.05,k=1.4,ff=0.004,L=100,P1==300,Pa=101325,M=29,R=8.314;
Function ((k+1)/2)*ln((x2^2*x3)/(x1^2*x4))-(1/ ... 您好，请问是只得到这一组正实数解么？而且网上只能下到1stopt 1.5，您是用的什么版本呢？穷学生表示买不起哈 只得到这一组实数解。 : Originally posted by xiangzigood at
全局求解？您能给推荐一个不，刚解除matlab不久。。... Global Optimization Toolbox,你可以试试这里面的函数。不过如何定义你自己的问题是一个问题，这个就真的只有靠你自己看帮助了。matlab求解含有字母的非线性方程组syms x1 x2 x3 x4 x5 x6 x20 D;f1 =x1*(67/100*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/1879*x3)-D);f2 =D*(x20-x2)-(11/5+/(x2+7/25)*(1-1/2039*x2)*(1-2/29/50*x2/(x2+))*x1;f3 =(-269/100+4535_百度作业帮
拍照搜题，秒出答案
matlab求解含有字母的非线性方程组syms x1 x2 x3 x4 x5 x6 x20 D;f1 =x1*(67/100*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/1879*x3)-D);f2 =D*(x20-x2)-(11/5+/(x2+7/25)*(1-1/2039*x2)*(1-2/29/50*x2/(x2+))*x1;f3 =(-269/100+4535
matlab求解含有字母的非线性方程组syms x1 x2 x3 x4 x5 x6 x20 D;f1 =x1*(67/100*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/1879*x3)-D);f2 =D*(x20-x2)-(11/5+/(x2+7/25)*(1-1/2039*x2)*(1-2/29/50*x2/(x2+))*x1;f3 =(-269/100+00*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/59/100*x2/(x2+31/2))*x1-D*x3;f4 =*x2/(x2+67/50)+--*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/580/151*x2/(x2+/100*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/1879*x3)*x4;f5 =/5-*x2/(x2+7/25)*(1-1/2039*x2)*(1-2//100*x2/(x2+241/100))*x4/(53/100+3/0*x5+x4)-*(58/25-*x2/(x2+7/25)*(1-1/2039*x2)*(1-2//50*x2/(x2+27/20))*x5/(7/50+x5*(1+500/209*x5))-67/100*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/1879*x3)*x5;f6 =*(58/25-*x2/(x2+7/25)*(1-1/2039*x2)*(1-2//50*x2/(x2+27/20))*x5/(7/50+x5*(1+500/209*x5))-*x6+*x3-67/100*x2/(x2+7/25)*(1-1/2039*x2)*(1-2/1879*x3)*x6;[x1 x2 x3 x4 x5 x6]=solve(f1,f2,f3,f4,f5,f6,'x1','x2','x3','x4','x5','x6');解不出来
我直接拷进去并没有报错,S=solve('d=k*(d/tan(a)-r*n/num)+h','sin(a)=asin(n*sin(acot(k)))','a','k')S = a:[2x1 sym]k:[2x1 sym]aa =atan(k*d*num/(d*num+k*r*n-h*num))kk =(-sin(sin(a))^2+n^2)^(1/2)/sin(sin(a))-(-sin(sin(a))^2+n^2)^(1/2)/sin(sin(a))你是不是没有清除变量,或者是matlab 2008a版的,那个版本有些函数有改动,可以试一下[a,k]=solve('d=k*(d/tan(a)-r*n/num)+h','sin(a)=asin(n*sin(acot(k)))')}