MATLAB中matlab解方程的函数函数是什么原理

人人网 - 抱歉
哦,抱歉,好像看不到了
现在你可以:
看看其它好友写了什么
北京千橡网景科技发展有限公司:
文网文[号··京公网安备号·甲测资字
文化部监督电子邮箱:wlwh@··
文明办网文明上网举报电话: 举报邮箱:&&&&&&&&&&&&【图文】第7章
MATLAB解方程与函数极值_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
MATLAB解方程与函数极值
||文档简介
总评分4.2|
浏览量4670
&&M​A​T​L​A​B​解​方​程​与​函​数​极​值
大小:119.00KB
登录百度文库,专享文档复制特权,财富值每天免费拿!
你可能喜欢查看: 39854|回复: 29|关注: 0
使用MATLAB solve和fsolve函数解方程组
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
关注者: 1
b=0.0105;c=22.15;T=220.14;F=222.2299;
f1='F*(1-exp(-y/b))*x=1';
f2='x*(c-b/(x^2))=2*(sqrt(((b^2)/(x^2)-c*b)*exp(2*(z-T)/b))-b/x)';
f3='5000+c*b-1/2*(b^2)/(x^2)-b/(x^2)*(z-y)-F^2*b*(y-(3*b)/2+2*b*exp(-y/b)-b/2*exp(-2*y/b))=0';
[x,y,z]=solve(f1,f2,f3)
其运行结果为:
[ empty sym ]
请教大侠,为什么会这样?该怎么改?
<h1 style="color:# 麦片财富积分
关注者: 4
试试这个fsolve解得的
设一个函数answer
function M=answer(x)
b=0.0105;c=22.15;T=220.14;F=222.2299;
M=[F*(1-exp(-x(2)/b))*x(1)-1;x(1)*(c-b/(x(1)^2))-2*(sqrt(((b^2)/(x(1)^2)-c*b)*exp(2*(x(3)-T)/b))-b/x(1));5000+c*b-1/2*(b^2)/(x(1)^2)-b/(x(1)^2)*(x(3)-x(2))-F^2*b*(x(2)-(3*b)/2+2*b*exp(-x(2)/b)-b/2*exp(-2*x(2)/b))];
end
复制代码在matlab命令行输入
x0=[1;2;3];
options=optimset('Display','iter');
[x,fval] = fsolve(@answer,x0,options)
复制代码解得
x =
& & 0.0353
& & 9.7667
& & 3.0010
& & 6.8528
& & 1.0798
& & 0.7003
复制代码
<h1 style="color:# 麦片财富积分
关注者: 1
hasen-chen 发表于
试试这个fsolve解得的
设一个函数answer在matlab命令行输入解得
嗯,确实算出来了!非常感谢!
额,能不能再麻烦您一次告诉我问题出在哪里以及为什么要这样做?
<h1 style="color:# 麦片财富积分
关注者: 4
游愚 发表于
嗯,确实算出来了!非常感谢!
额,能不能再麻烦您一次告诉我问题出在哪里以及为什么要这样做? ...
一般用fsolve解比较好写,用solve解时,上面的等号要用==
<h1 style="color:# 麦片财富积分
关注者: 1
hasen-chen 发表于
一般用fsolve解比较好写,用solve解时,上面的等号要用==
我刚才试了一下,报错了诶:
??? Error using ==& solve&getEqns at 173
' F*(1-exp(-y/b))*x==1 ' is not a valid expression or equation.
Error in ==& solve at 84
[eqns,vars] = getEqns(varargin{:});
<h1 style="color:# 麦片财富积分
关注者: 4
游愚 发表于
我刚才试了一下,报错了诶:
??? Error using ==& solve&getEqns at 173
' F*(1-exp(-y/b))*x==1 ' is no ...
你试试
clc
b=0.0105;c=22.15;T=220.14;F=222.2299;
f1=F*(1-exp(-y/b))*x-1;
f2=x*(c-b/(x^2))-2*(sqrt(((b^2)/(x^2)-c*b)*exp(2*(z-T)/b))-b/x);
f3=5000+c*b-1/2*(b^2)/(x^2)-b/(x^2)*(z-y)-F^2*b*(y-(3*b)/2+2*b*exp(-y/b)-b/2*exp(-2*y/b));
[x,y,z]=solve(f1==0,f2==0,f3==0)
复制代码
<h1 style="color:# 麦片财富积分
关注者: 4
我这边用个程序还没有错,但是运行时间很长,现在还没有运行出来,对于solve我也不太知道是怎么回事
<h1 style="color:# 麦片财富积分
关注者: 1
hasen-chen 发表于
我这边用个程序还没有错,但是运行时间很长,现在还没有运行出来,对于solve我也不太知道是怎么回事 ...
呵呵,好吧,谢谢您的指导!
另外,那个fval是指在该X值得时候所取得最大值吧?而且经过实际检验所得的x值也是不符合实际的诶,难道是我的赋值有问题?
<h1 style="color:# 麦片财富积分
关注者: 6
这道题用Matlab的数值解(fsolve)很难得到正解,初值不好猜。用1stOpt试试,不用初值,比Matlab强很多:
Constant b=0.0105,c=22.15,T=220.14,F=222.2299;
Function F*(1-exp(-y/b))*x=1;
& && && &x*(c-b/(x^2))=2*(sqrt(((b^2)/(x^2)-c*b)*exp(2*(z-T)/b))-b/x);
& && && &5000+c*b-1/2*(b^2)/(x^2)-b/(x^2)*(z-y)-F^2*b*(y-(3*b)/2+2*b*exp(-y/b)-b/2*exp(-2*y/b))=0;
复制代码结果:
y: 0.72479
<h1 style="color:# 麦片财富积分
关注者: 1
youyouyou 发表于
这道题用Matlab的数值解(fsolve)很难得到正解,初值不好猜。用1stOpt试试,不用初值,比Matlab强很多:结 ...
额,不会用1stOpt的孩纸伤不起啊,求指点·········
站长推荐 /2
Powered by求matlab高手解方程函数就是题目中给出的那个函数,你可以设一个初始的ψ值对应一个X3值,然后再令ψ+2对应另一个X3值。。。。就可得到一系列方程。怎么求解呢,未知数这么多,急。
clcclear all%8个未知数,建立8个方程syms x1 x2 y1 r z3 b a a0;[x1,x2,y1,r,z3,b,a,a0]=solve((1*(pi/8)-x2)*(r*cos(pi+0.0349)*sin(b)+r*sin(pi+0.0349)*cos(a)*cos(b)+z3*cos(b)*sin(a)+x1)-(r*sin(pi...
为您推荐:
扫描下载二维码【图文】第7章
MATLAB解方程与函数极值_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
MATLAB解方程与函数极值
上传于||文档简介
&&对&#8203;你&#8203;们&#8203;的&#8203;学&#8203;习&#8203;很&#8203;有&#8203;用&#8203;的&#8203;,&#8203;主&#8203;要&#8203;是&#8203;p&#8203;p&#8203;t&#8203;格&#8203;式&#8203;的&#8203;。
大小:135.50KB
登录百度文库,专享文档复制特权,财富值每天免费拿!
你可能喜欢

我要回帖

更多关于 matlab解方程函数 的文章

 

随机推荐