急!matlab solve用法

想知道matlab中solve 和 subs的基本用法,现在有几个方程,想通过这些方程导出y关于x的表达式,可以做到么 谢谢了~~

在MATLAB中,solve函数主要是用来求解代数方程(多项式方程)的符号解析解。

用法示例:

solve(eq)

solve(eq, var)

solve(eq1, eq2, ..., eqn)

g = solve(eq1, eq2, ..., eqn, var1, var2, ..., varn)

其中,eq代表一个符号表达式或字符串,var代表一个变量名称

详细的解释:

g=solve(eq)

函数求代数方程的符号解析解。参量eq表示符号表达式或字符串。若eq是一符号表达式或一没有等号的字符串,则函数对方程的默认变量求解方程eq=0,默认变量由命令findsym(eq)确定。若输出参量g为单一变量,则对于有多重解的非线性方程,g为一行向量。

g=solve(eq,var)

扩展资料

solve使用注意事项

1、solve解非代数方程的能力较弱,最好结合其他方式求解非代数方程

2、如果解得是一个方程组,而且采用了形如[a,b]=solve(a+b==1, 2*a-b==4,a,b) 的格式,那么,在MATLAB R2014a中没问题,可以保证输出的a,b就等于相应的解。

3、但是在R2012b等早先版本中不能保证输出的顺序就是你声明变量时的顺序。所以最好采用g=solve(a+b==1, 2*a-b==4,a,b)这种单输出格式,这样输出的是一个结构体,g.a和g.b就是对应的解。

温馨提示:答案为网友推荐,仅供参考
第1个回答  2008-10-25
matlab的运行结果是:
0 -8.0000
-8.5000 - 0.8660i -0.0000 + 0.0000i
-8.5000 + 0.8660i -0.0000 - 0.0000i
上面的结果自身就是不正确的,我分析的原因是:matlab先求出了解析解,接着代入系数值,由于matlab运算精度低,系数又太多,公式极其复杂,从而导致结果失真过于严重

mathematica运行结果显示:mathematica先求出解析解,接着精确代入系数值计算,也出现了数量级为-16的误差,有部分失真,但是这个失真很容易发现,而得到修正

mathematica执行消去y的运算结果显示原方程组只能由三组解

借助mathematica的结果分析,可以得出:matlab的运行结果中,并没有丢失一组解,而是解的失真太为严重达到了0.1,
可以用下面的程序执行数值解:
function s=sss(p)
syms x y
a=p(1);
b=p(2);
c=p(3);
d=p(4);
e=p(5);
f=p(6);
g=p(7);
a1=simplify(c*x^2+c*x*y+(c*g-c*a-c*b+1)*x-a-e);
a2=simplify(d*y^2+d*x*y+(d*g-d*a-d*b+1)*y-b-f);
[x0,y0]=solve(a1,a2);
s=[x0 y0]

这个程序运行速度快且结果准确,
结果为
[ 0, 0]
[ 0, -8]
[ -13, 0]

注:可以看到新的程序输出结果是矩阵形式,而老程序的输出不是,问题很可能和eval的使用有关.
第2个回答  2018-07-19

主要是这个命令上的错误,['num2str(T0)*(0.9)^x=',num2str(Tend)]

最后一句代码应该这样改一下就可以了。

Time1=ceil(double(solve([num2str(T0),'*(0.9)^x=',num2str(Tend)]))) %计算迭代的次数

修改后执行结果

第3个回答  2018-06-28
1、solve 这个命令用来求解符号方程、方程组等。
x=solve('eqn1','eqn2',...,'var1','var2',...)
例子:
解方程组:x^2+y-6=0;y^2+x-6=0
程序设计:[x,y]=solve('x^2+y-6','y^2+x-6','x','y')
2、subs可以把符号表达式里面的符号变量替换为新的变量
F=subs(f,old,new)
例子:
用符号表达式计算x^2+x-1在x=2时的值
程序设计:g=sym('x^2+x-1');
G=subs(g,'x',2)
第4个回答  推荐于2018-02-27
1、solve 这个命令用来求解符号方程、方程组等。
x=solve('eqn1','eqn2',...,'var1','var2',...)
例子:
解方程组:x^2+y-6=0;y^2+x-6=0
程序设计:[x,y]=solve('x^2+y-6','y^2+x-6','x','y')
2、subs可以把符号表达式里面的符号变量替换为新的变量
F=subs(f,old,new)
例子:
用符号表达式计算x^2+x-1在x=2时的值
程序设计:g=sym('x^2+x-1');
G=subs(g,'x',2)本回答被提问者和网友采纳
相似回答