怎么用matlab解含有字母系数的方程组的解，举个简单例子_百度知道
怎么用matlab解含有字母系数的方程组的解，举个简单例子
li1127217ye
li1127217ye
采纳数：18
获赞数：299
举个简单例子，解方程组x+A*y=10x-B*y=1其中x,y为变量，A,B为字母系数。只要在Matlab中输入syms x,y,A,B[x y]=solve('x+A*y=10','x-B*y=1','x','y')即可求出解x = (A + 10*B)/(A + B)y = 9/(A + B)对于函数solve的具体用法，可以通过输入help solve来学习。希望我的回答能够解决您的疑问，谢谢。
chenjiechenl
chenjiechenl
采纳数：19
获赞数：136
擅长：暂未定制
例：解方程组a*x+2*y=4,4*x+b*y=7MATLAB代码：z=solve('a*x+2*y=4','4*x+b*y=7')即可。
为你推荐：
其他类似问题
您可能关注的内容
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。matlab&解方程组
最近有多人问如何用matlab解方程组的问题，其实在matlab中解方程组还是很方便的，例如，对于代数方程组Ax=b(A为系数矩阵，非奇异)的求解，MATLAB中有两种方法：
(1)x=inv(A)*b — 采用求逆运算解方程组；
& (2)x=A\B — 采用左除运算解方程组
PS：使用左除的运算效率要比求逆矩阵的效率高很多~
2x1+3x2=13
&&A=[1,2;2,3];b=[8;13];
&&x=inv(A)*b
即二元一次方程组的解x1和x2分别是2和3。
对于同学问到的用matlab解多次的方程组，有符号解法，方法是：先解出符号解,然后用vpa(F,n)求出n位有效数字的数值解.具体步骤如下：
第一步:定义变量syms x y z ...;
第二步:求解[x,y,z,...]=solve('eqn1','eqn2',...,'eqnN','var1','var2',...'varN');
第三步:求出n位有效数字的数值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);...。
如：解二（多）元二（高）次方程组：
x^2+3*y+1=0
y^2+4*x+1=0
解法如下：
&&[x,y]=solve('x^2+3*y+1=0','y^2+4*x+1=0');
&&x=vpa(x,4);
&&y=vpa(y,4);
1.635+3.029*i
1.635-3.029*i
1.834-3.301*i
1.834+3.301*i
二元二次方程组，共4个实数根；
还有的同学问，如何用matlab解高次方程组（非符号方程组）？举个例子好吗？
解答如下：
基本方法是：solve(s1,s2,…,sn,v1,v2,…,vn)，即求表达式s1,s2,…,sn组成的方程组，求解变量分别v1,v2,…,vn。
具体例子如下：
x^2 + x*y + y = 3
x^2 - 4*x + 3 = 0
&& [x,y] = solve('x^2 + x*y + y =
3','x^2 - 4*x + 3 = 0')
运行结果为
即x等于1和3；y等于1和-1.5
&&[x,y] = solve('x^2 + x*y + y =
3','x^2 - 4*x + 3= 0','x','y')
结果一样，二元二方程都是4个实根。
通过这三个例子可以看出，用matlab解各类方程组都是可以的，方法也有多种，只是用到解方程组的函数，注意正确书写参数就可以了，非常方便。
from:http://bbs.seu.edu.cn/pc/pccon.php?id=950&nid=14498&tid=0
2、变参数非线性方程组的求解
对于求解非线性方程组一般用fsolve命令就可以了，但是对于方程组中某一系数是变化的，该怎么求呢？
%定义方程组如下，其中k为变量
function F = myfun(x,k)
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))];
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)
x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组
x1(i)=x(1);
x2(i)=x(2);
plot(k,x1,'-b',k,x2,'-r');
xlabel('k')
legend('x1','x2')
cite from:http://forum.simwe.com/archiver/tid-836299.html
3、非线性方程数值求解
matlab里solve如何使用,是否有别的函数可以代替它.
matlab里我解y=9/17*exp(-1/2*t)*17^(1/2)*sin(1/2*17^(1/2)*t)=0这样的方程为什么只得到0这一个解,如何可以的到1/2*17^(1/2)*t=n*(pi)这样一族解??
在matlab里面solve命令主要是用来求解代数方程（即多项式）的解,但是也不是说其它方程一个也不能解，不过求解非代数方程的能力相当有限，通常只能给出很特殊的实数解。(该问题给出的方程就是典型的超越方程，非代数方程)
从计算机的编程实现角度讲，如今的任何算法都无法准确的给出任意非代数方程的所有解，但是我们有很多成熟的算法来实现求解在某点附近的解。matlab也不例外，它也只能给出任意非代数方程在某点附近的解，函数有两个：fzero和fsolve,具体用法请用help或doc命令查询吧。如果还是不行，你还可以将问题转化为非线性最优化问题，求解非线性最优化问题的最优解，可以用的命令有：fminbnd,
fminsearch, fmincon等等。
*非线性方程数值求解
*单变量非线性方程求解
在MATLAB中提供了一个fzero函数，可以用来求单变量非线性方程的根。该函数的调用格式为：
z=fzero('fname',x0,tol,trace)
其中fname是待求根的函数文件名，x0为搜索的起点。一个函数可能有多个根，但fzero函数只给出离x0最近的那个根。tol控制结果的相对精度，缺省时取tol=eps，trace�指定迭代信息是否在运算中显示，为1时显示，为0时不显示，缺省时取trace=0。
例 求f(x)=x-10x+2=0在x0=0.5附近的根。
步骤如下：
(1) 建立函数文件funx.m。
function fx=funx(x)
fx=x-10.^x+2;
(2) 调用fzero函数求根。
z=fzero('funx',0.5)
**非线性方程组的求解
对于非线性方程组F(X)=0，用fsolve函数求其数值解。fsolve函数的调用格式为：
X=fsolve('fun',X0,option)
其中X为返回的解，fun是用于定义需求解的非线性方程组的函数文件名，X0是求根过程的初值，option为最优化工具箱的选项设定。最优化工具箱提供了20多个选项，用户可以使用optimset命令将它们显示出来。如果想改变其中某个选项，则可以调用optimset()函数来完成。例如，Display选项决定函数调用时中间结果的显示方式，其中‘off’为不显示，‘iter’表示每步都显示，‘final’只显示最终结果。
optimset(‘Display’,‘off’)将设定Display选项为‘off’。
例 求下列非线性方程组在(0.5,0.5) 附近的数值解。
(1) 建立函数文件myfun.m。
function q=myfun(p)
q(1)=x-0.6*sin(x)-0.3*cos(y);
q(2)=y-0.6*cos(x)+0.3*sin(y);
(2) 在给定的初值x0=0.5,y0=0.5下，调用fsolve函数求方程的根。
x=fsolve('myfun',[0.5,0.5]',optimset('Display','off'))
将求得的解代回原方程，可以检验结果是否正确，命令如下：
q=myfun(x)
1.0e-009 *
可见得到了较高精度的结果。
from:http://blog.sina.com.cn/s/blog_56ef652d0100ebew.html
4、fsolve函数解方程
[X,FVAL,EXITFLAG,OUTPUT,JACOB]=FSOLVE(FUN,X0,...) returns the
Jacobian of FUN at X.
FUN can be specified using @:
x = fsolve(@myfun,[2 3 4],optimset('Display','iter'))
where myfun is a MATLAB function such as:
function F = myfun(x)
F = sin(x);
FUN can also be an anonymous function:
x = fsolve(@(x) sin(3*x),[1 4],optimset('Display','off'))
If FUN is parameterized, you can use anonymous functions to capture
problem-dependent parameters. Suppose you want to solve the system
nonlinear equations given in the function myfun, which is
parameterized
by its second argument c. Here myfun is an M-file function such
function F = myfun(x,c)
F = [ 2*x(1) - x(2) - exp(c*x(1))
-x(1) + 2*x(2) - exp(c*x(2))];
To solve the system of equations for a specific value of c, first
assign the
value to c. Then create a one-argument anonymous function that
that value of c and calls myfun with two arguments. Finally, pass
this anonymous
function to FSOLVE:
c = -1; % define parameter first
x = fsolve(@(x) myfun(x,c),[-5;-5])
cite from:
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。当前位置：
&用Matlab求解带有符号变量的方程组，望大神指导
用Matlab求解带有符号变量的方程组，望大神指导
作者 Clinkz
代码如下：
syms x y p1 p2
M1=x^2+(2+p1)*x+1+y
M2=y^2+(1+p2)*y+1+x
[x1,y1]=solve(M1,M2,'x','y')
放在matlab下运行，一直显示处于busy状态，等很久也解不出来，有没有别的函数或者方法解决这个问题。
跪求大神指导
MATLAB可运行，但不存在解析解。
附图1给出部分计算结果。
如果已知p1 p2，可以求数值解。
引用回帖:: Originally posted by 月只蓝 at
MATLAB可运行，但不存在解析解。
附图1给出部分计算结果。
如果已知p1 p2，可以求数值解。
附图1.png 关键就在这，因为我以后要对x和y求关于p1,p2的偏导，所以我想求出x、y关于p1,p2的方程
引用回帖:: Originally posted by 月只蓝 at
MATLAB可运行，但不存在解析解。
附图1给出部分计算结果。
如果已知p1 p2，可以求数值解。
附图1.png 你的程序怎么运行的，我的怎么已知处于busy状态，是不是matlab版本问题
引用回帖:: Originally posted by Clinkz at
你的程序怎么运行的，我的怎么已知处于busy状态，是不是matlab版本问题... R2010a
引用回帖:: Originally posted by 月只蓝 at
R2010a... 运行了多长时间，为什么我这运行一晚上都没有结果？
你求解的式子不是等式，怎么会有结果
24小时热帖
下载小木虫APP
与700万科研达人随时交流24小时热门版块排行榜&&&&
【悬赏金币】回答本帖问题，作者sami512将赠送您 10 个金币
(小有名气)
在线: 59.8小时
虫号: 3962207
注册: 专业: 功能与智能高分子
求助，Matlab怎么求解带参数的方程组已有2人参与
求助各位大神，matlab小白一个，论文需要求解R与a,b,c的函数关系，a,b,c为参数，是否要先求取x,y,z,t的值，网上找了许多程序，都行不通，跪求大神们的程序代码
2*(a+b)*x-(0.5*a+b)*y-(0.5*a+b)*z-a*t=V
-(0.5*a+b) *x+(a+b+c)*y-c*t=0
-(0.5*a+b)*x+(a+b+c)*y-c*t=0
-a*x-c*y-c*z+2*(a+b+c)*t=0
& 猜你喜欢
已经有14人回复
已经有6人回复
已经有24人回复
已经有8人回复
已经有86人回复
已经有6人回复
已经有13人回复
已经有6人回复
已经有72人回复
已经有3人回复
(小有名气)
jjdg: 屏蔽内容
本帖内容被屏蔽
(小有名气)
在线: 40.1小时
虫号: 6172469
注册: 专业: 化学生物学与生物有机化学
★ jjdg: 金币+1, 感谢参与
solve 函数，可以用符号计算
(小有名气)
在线: 59.8小时
虫号: 3962207
注册: 专业: 功能与智能高分子
引用回帖:: Originally posted by Ltd_123 at
solve 函数，可以用符号计算 solve解出来的值都是0
(小有名气)
在线: 25.7小时
虫号: 8588011
注册: 性别: GG专业: 计算机软件
【答案】应助回帖
★ 感谢参与，应助指数 +1jjdg: 金币+1, 感谢参与
x y z t需要给定初值！不确定的话，先模拟出大概位置，之后循环求解。
(小有名气)
在线: 25.7小时
虫号: 8588011
注册: 性别: GG专业: 计算机软件
【答案】应助回帖
引用回帖:: Originally posted by 神级呆呆贼 at
x y z t需要给定初值！不确定的话，先模拟出大概位置，之后循环求解。 我再补充一下，如果想求解析解，可以尝试了解一下mathematica软件
(正式写手)
在线: 145.9小时
虫号: 1476844
注册: 性别: GG专业: 结构工程
【答案】应助回帖
感谢参与，应助指数 +1
第三、四个方程怎么一样？是否打错了？如果不一样，这不就是四元一次方程组么，初中知识吧，直接求解方程啊，不需要用MATLAB算吧，将x、y、z、t用a、b、c、R表示，然后将x的值代入第一个方程，不就知道R、V、a、b、c的关系啦。
(小有名气)
在线: 59.8小时
虫号: 3962207
注册: 专业: 功能与智能高分子
引用回帖:: Originally posted by 神级呆呆贼 at
我再补充一下，如果想求解析解，可以尝试了解一下mathematica软件... 已解决，谢谢大神
(小有名气)
在线: 59.8小时
虫号: 3962207
注册: 专业: 功能与智能高分子
引用回帖:: Originally posted by Ltd_123 at
solve 函数，可以用符号计算 已解决，谢谢大神
(小有名气)
在线: 59.8小时
虫号: 3962207
注册: 专业: 功能与智能高分子
引用回帖:: Originally posted by lwjjtt at
第三、四个方程怎么一样？是否打错了？如果不一样，这不就是四元一次方程组么，初中知识吧，直接求解方程啊，不需要用MATLAB算吧，将x、y、z、t用a、b、c、R表示，然后将x的值代入第一个方程，不就知道R、V、a、b、 ... 第3.4个有一个写错了现在解决了，谢谢你
相关版块跳转
第一性原理
我要订阅楼主
的主题更新
小木虫,学术科研互动社区,为中国学术科研免费提供动力
违规贴举报删除请发送邮件至：
广告投放与宣传请联系 李想 QQ：
QQ：&&邮箱：
Copyright &
MuChong.com, All Rights Reserved. 小木虫 版权所有 上传我的文档
 上传文档
 下载
 收藏
粉丝量：25
该文档贡献者很忙，什么也没留下。
 下载此文档
用matlab求解非线性方程组的几种方法之程序
下载积分：1200
内容提示：用matlab求解非线性方程组的几种方法之程序
文档格式：PDF|
浏览次数：676|
上传日期： 10:54:07|
文档星级：
全文阅读已结束，如果下载本文需要使用
 1200 积分
下载此文档
该用户还上传了这些文档
用matlab求解非线性方程组的几种方法之程序
关注微信公众号}