版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

zhmdream

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2535.9
红花: 1
帖子: 776
在线: 292.2小时
虫号: 591487
注册: 2008-09-02
性别: GG
专业: 凝聚态物性I:结构、力学和

[求助] 求助Matlab解方程问题

大家好，请大家一个Matlab解方程的问题，程序：
a=input('输入a:');
solve('a^2/(a+x)+a/(x-a)=9');
ans

运行结果：
输入a:5

ans =

(1/18*a+1/18+1/18*(a^2-34*a+361)^(1/2))*a
(1/18*a+1/18-1/18*(a^2-34*a+361)^(1/2))*a

我定义了a=5，但计算结果中a怎么是个字母啊?
是不是我没有定义a为数值？还是其他原因？
望大家帮忙，谢谢啦。

回复此楼

人生如梦，岁月无情~~~

1楼 2011-10-09 19:59:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sudo

木虫 (正式写手)

程序强帖: 16
应助: 6 (幼儿园)
金币: 1297.6
散金: 1486
红花: 20
帖子: 588
在线: 641小时
虫号: 1211394
注册: 2011-02-24
性别: GG
专业: 文艺美学

【答案】应助回帖

zhmdream(金币+10): 这样虽然可以，但是对于输入任意的a，还是手动改，比较麻烦 2011-10-10 10:19:02
zhmdream(金币+5): 谢谢 2011-10-10 11:05:26
zhmdream(金币+45): 谢谢啦！问题解决差不多了。我之后再翻翻书，消化消化。看来我还是对很多函数和命令理解太少了。呵呵 2011-10-10 18:31:10

SOLVE Symbolic solution of algebraic equations.

solve是用来解出符号解的啦，所以solve里面字符串里面的a和你定义的a并不是同一个

你可以用subs函数

CODE:

A=solve('a^2/(a+x)+a/(x-a)=9');

subs(A, 5)

这样就把5“代入到符号表达式”中，求出数值解来了

PS：如果solve不能求出一些方程的解析解，那么应该考虑用其它数值方法

赞一下

回复此楼

2楼2011-10-10 08:59:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sudo

木虫 (正式写手)

程序强帖: 16
应助: 6 (幼儿园)
金币: 1297.6
散金: 1486
红花: 20
帖子: 588
在线: 641小时
虫号: 1211394
注册: 2011-02-24
性别: GG
专业: 文艺美学

呼唤变通...

CODE:

a=input('输入a:');

A=solve('a^2/(a+x)+a/(x-a)=9');

subs(A, a)

你是不是想这样？

赞一下

回复此楼

3楼2011-10-10 10:36:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhmdream

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2535.9
红花: 1
帖子: 776
在线: 292.2小时
虫号: 591487
注册: 2008-09-02
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

3楼: Originally posted by sudo at 2011-10-10 10:36:37:

呼唤变通...

CODE:

a=input('输入a:');

A=solve('a^2/(a+x)+a/(x-a)=9');

subs(A, a)

你是不是想这样？

谢谢啦！呵呵，其实我是想编一个这样的程序：

c=input('请输入c:')
a=x+1;
b=2*x+3;
A=solve('a*b+2*x+a-b=c');
subs(A,a);
subs(A,b);
subs(A,c);
ans

a、b包含了x，然后求解x。但是运行还是出现了以前的问题：

请输入c:2
c =2
ans =-1/2*a*b-1/2*a+1/2*b+1

呵呵，我是刚学Matlab，不是很懂，所以想一步一步弄明白。谢谢帮忙。帮我看看问题出在哪儿？剩下50金币全给你了。

赞一下

回复此楼

人生如梦，岁月无情~~~

4楼2011-10-10 11:04:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sudo

木虫 (正式写手)

程序强帖: 16
应助: 6 (幼儿园)
金币: 1297.6
散金: 1486
红花: 20
帖子: 588
在线: 641小时
虫号: 1211394
注册: 2011-02-24
性别: GG
专业: 文艺美学

【答案】应助回帖

★ ★ ★ ★ ★
微尘、梦想(金币+5): 谢谢参与应助~ 2011-10-10 19:48:48

引用回帖:

4楼: Originally posted by zhmdream at 2011-10-10 11:04:43:
谢谢啦！呵呵，其实我是想编一个这样的程序：

c=input('请输入c:')
a=x+1;
b=2*x+3;
A=solve('a*b+2*x+a-b=c');
subs(A,a);
subs(A,b);
subs(A,c);
ans

a、b包含了x，然后求解x。但是运行还是出现了 ...

嗯，慢慢来吧

CODE:

c=input('请输入c:');

syms x;

a=x+1;

b=2*x+3;

f=a*b+2*x+a-b; %构造方程式的左边

%将方程式构造为字符串

equation = [char(f) ' = ' num2str(c,9)]; %num2str把数值转为字符串，第二个参数是有效数字

%计算符号解（经过上面的运算equation里面只有一个变量x）

A=solve(equation, 'x')

%输出数值解

double(A)