24小时热门版块排行榜

>论坛更新日志 (4518)
>导师招生 (1020)
>休闲灌水 (807)
>虫友互识 (397)
>文献求助 (364)
>博后之家 (274)
>招聘信息布告栏 (199)
>论文投稿 (159)
>考博 (128)
>硕博家园 (117)
>找工作 (108)
>第一性原理 (93)
>生物科学 (38)
>绿色求助(高悬赏) (37)
>基金申请 (36)
>考研 (30)

返回列表

voice2001

铜虫 (著名写手)

应助: 7 (幼儿园)
金币: 8067
散金: 12
帖子: 2233
在线: 302.7小时
虫号: 288175
注册: 2006-10-21
专业: 海洋物理学

[求助] 【求助】matlab2007的solve求解三角函数，计算结果错误已有1人参与

matlab程序求解三角函数程序如下：

syms x
solve('tan(x-pi/4)-2*x/3=0')
return

matlab2007运行结果如下：
-.13791712807073613045661524523992-1.5323481937969730123523907483188*i

matlab2014运行结果如下：
1.6044193659902407984267577001949

请问
1.这是什么原因？
2.应该如何修改程序使结果一致？

回复此楼

1楼 2017-03-02 13:44:39

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68174.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3292.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
voice2001: 金币+5, ★★★★★最佳答案 2017-03-03 14:06:30

>> x=[-.13791712807073613045661524523992-1.5323481937969730123523907483188*i 1.6044193659902407984267577001949];
>> y=tan(x-pi/4)-2*x/3

y =

1.0e-016 *

0.2776 0

两个解都能使得原方程残差几乎为0，何来“计算结果错误”的说法？只是其中一个解是复数而已。

原方程与直线y=0存在多个交点，说明原方程存在多个解。
比如想求得x=1附近的解：
format long
x0=1;
[x fval]=fsolve(@(x) tan(x-pi/4)-2*x/3,x0)

计算结果：
x =

1.604419385283544

fval =

2.850402580634182e-008

如果想求得某一个范围之内，所有的实数解，可参考：http://muchong.com/bbs/viewthread.php?tid=9917857

untitled.png

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2017-03-02 16:21:46

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

voice2001

铜虫 (著名写手)

应助: 7 (幼儿园)
金币: 8067
散金: 12
帖子: 2233
在线: 302.7小时
虫号: 288175
注册: 2006-10-21
专业: 海洋物理学

引用回帖:

2楼: Originally posted by 月只蓝 at 2017-03-02 16:21:46
>> x=;
>> y=tan(x-pi/4)-2*x/3

y =

1.0e-016 *

0.2776 0

两个解都能使得原方程残差几乎为0，何来“计算结果错误”的说法？只是其中一个解是复数而已。

原方程与直线y=0 ...

【追问】在matlab2007上只得到一个复数解，在matlab2014上只得到一个实数解，而不是同时得到多个解呢？

赞一下

回复此楼

3楼2017-03-02 16:44:01

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