学习Matlab强大的符号计算（解方程）

学习Matlab强大的符号计算（解方程）

分类： 学习 其它语言 2011-09-19 10:34 3135人阅读 评论(0) 收藏 举报

matlab

[plain] view plain copy print ?

1. syms x1 x2; %定义符号变量  
2. syms y1 y2;  
3. clc  
4. [x1, x2] = solve('x1 + x2 = y1',...  
5.     'x1^2 + x2^2 = y2',...  
6.     'x1', 'x2') %求x1，x2用y1, y2表示的表达  
7. %求解y1 = 1，y2 = 5 时的函数值  
8. y1 = 1;  
9. y2 = 5;  
10. x1 = eval(x1)   % eval %把符号变量转换为数值型,  
11. x2 = vpa(eval(x2), 3) % vpa  %精度,这里用3位小数精度表示  
12. %{  
13. %运行结果为  
14. x1 =  
15.  1/2*y1-1/2*(-y1^2+2*y2)^(1/2)  
16.  1/2*y1+1/2*(-y1^2+2*y2)^(1/2)  
17. x2 =  
18.  1/2*y1+1/2*(-y1^2+2*y2)^(1/2)  
19.  1/2*y1-1/2*(-y1^2+2*y2)^(1/2)  
20. x1 =  
21.     -1  
22.      2  
23. x2 =  
24.   2.  
25.  -1.  
26. %}  

syms x1 x2; %定义符号变量
syms y1 y2;
clc
[x1, x2] = solve('x1 + x2 = y1',...
    'x1^2 + x2^2 = y2',...
    'x1', 'x2') %求x1，x2用y1, y2表示的表达
%求解y1 = 1，y2 = 5 时的函数值
y1 = 1;
y2 = 5;
x1 = eval(x1)   % eval %把符号变量转换为数值型,
x2 = vpa(eval(x2), 3) % vpa  %精度,这里用3位小数精度表示
%{
%运行结果为
x1 =
 1/2*y1-1/2*(-y1^2+2*y2)^(1/2)
 1/2*y1+1/2*(-y1^2+2*y2)^(1/2)
x2 =
 1/2*y1+1/2*(-y1^2+2*y2)^(1/2)
 1/2*y1-1/2*(-y1^2+2*y2)^(1/2)
x1 =
    -1
     2
x2 =
  2.
 -1.
%}

以前把Matlab当做强大的计算器，才知道它不止计算器那么简单。