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

登录注册

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

24小时热门版块排行榜

返回列表

kuntenghunya

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 98.8
散金: 15
帖子: 25
在线: 47.8小时
虫号: 1818843
注册: 2012-05-15
性别: MM
专业: 传热传质学

[求助] 特征方程求解

求高手帮助，解特征方程

特征方程.jpg

回复此楼

» 收录本帖的淘帖专辑推荐

北大科研分享

» 猜你喜欢

26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有1人回复
如何从铁电相到顺电相。已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有53人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有26人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有24人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有1人回复
PVA溶解已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2013-07-23 11:01:50

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

月只蓝

主管区长 (职业作家)

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

【答案】应助回帖

★
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-07-23 11:45:01

x=solve('x/(tan(x))=1-1.5')

x =

                              [    0.    ]
                              [          ]
                              [-1.836597203]
                              [          ]
                              [1.836597203 ]

赞一下

回复此楼

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

2楼2013-07-23 11:07:08

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

kuntenghunya

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 98.8
散金: 15
帖子: 25
在线: 47.8小时
虫号: 1818843
注册: 2012-05-15
性别: MM
专业: 传热传质学

引用回帖:

2楼: Originally posted by 月只蓝 at 2013-07-23 11:07:08
x=solve('x/(tan(x))=1-1.5')

x =





                     ...

多谢。楼主是用matlab求解的吗？我刚试了一下，得到的结果和上面的不一样

赞一下

回复此楼

3楼2013-07-23 11:20:29

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

月只蓝

主管区长 (职业作家)

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

【答案】应助回帖

引用回帖:

3楼: Originally posted by kuntenghunya at 2013-07-23 11:20:29
多谢。楼主是用matlab求解的吗？我刚试了一下，得到的结果和上面的不一样...

是MATLAB，你的结果不一样，那你的结果是多少呢？

Maple的结果和MATLAB是一致的：

附图2.jpg

赞一下

回复此楼

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

4楼2013-07-23 13:43:39

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

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
kuntenghunya: 金币+5, ★★★很有帮助 2013-07-23 18:10:38
csgt0: 金币+2, 多谢指导 2013-07-24 16:59:04

在excel文件中，修改 K值，比如 1，就得到第一个根 1.836597203，修改成 10，就得到第10个根 29.8618724，等等。当然，根是对称的。 -29.8618724 也是根。前十个根是（除了平凡解 x=0）

1  1.836597203
2  4.815842318
3  7.917052685
4  11.04082982
5  14.17243207
6  17.30764054
7  20.44480347
8  23.58314331
9  26.72224637
10 29.8618724

赞一下

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : Newton05.xls

2013-07-23 16:17:42, 16 K

5楼2013-07-23 16:18:30

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

月只蓝

主管区长 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
kuntenghunya: 金币+5, ★★★很有帮助 2013-07-23 18:10:15
csgt0: 金币+2, 多谢指导 2013-07-24 16:59:18

5楼的虫友说得很对，应当考虑多根的情况。附图表示了在[-1000,1000]，μ存在多解。
根据实际中对μ值的要求，给定μ值的取值区间，可按下程序解得μ(以[-100,100]为例)：

function sol_nol
clear all;clc
format short
n=100;  %μ的取值区间，现在是[-100,100]，可根据需要调节n的大小
x0=-n:1:n;
for i=1:length(x0)
x_sol(i)=fsolve(@fun,x0(i),optimset('Display','Off'));
end
[x0' x_sol']

function y=fun(x)
y=x/(tan(x))+0.5;

结果：
  μ初值 μ最终解
-100.0000  -98.9652
  -99.0000  -98.9652
  -98.0000  -98.9652
  -97.0000  -95.8238
  -96.0000  -95.8238
  -95.0000  -95.8238
  -94.0000  -92.6824
  -93.0000  -92.6824
  -92.0000  -92.6824
  -91.0000  -89.5410
  -90.0000  -89.5410
  -89.0000  -89.5410
  -88.0000  -89.5410
  -87.0000  -86.3996
  -86.0000  -86.3996
  -85.0000  -86.3996
  -84.0000  -83.2582
  -83.0000  -83.2582
  -82.0000  -83.2582
  -81.0000  -80.1169
  -80.0000  -80.1169
  -79.0000  -80.1169
  -78.0000  -76.9755
  -77.0000  -76.9755
  -76.0000  -76.9755
  -75.0000  -73.8342
  -74.0000  -73.8342
  -73.0000  -73.8342
  -72.0000  -70.6929
  -71.0000  -70.6929
  -70.0000  -70.6929
  -69.0000  -67.5516
  -68.0000  -67.5516
  -67.0000  -67.5516
  -66.0000  -67.5516
  -65.0000  -64.4104
  -64.0000  -64.4104
  -63.0000  -64.4104
  -62.0000  -61.2692
  -61.0000  -61.2692
  -60.0000  -61.2692
  -59.0000  -58.1281
  -58.0000  -58.1281
  -57.0000  -58.1281
  -56.0000  -54.9870
  -55.0000  -54.9870
  -54.0000  -54.9870
  -53.0000  -51.8459
  -52.0000  -51.8459
  -51.0000  -51.8459
  -50.0000  -48.7050
  -49.0000  -48.7050
  -48.0000  -48.7050
  -47.0000  -45.5641
  -46.0000  -45.5641
  -45.0000  -45.5641
  -44.0000  -45.5641
  -43.0000  -42.4233
  -42.0000  -42.4233
  -41.0000  -42.4233
  -40.0000  -39.2826
  -39.0000  -39.2826
  -38.0000  -39.2826
  -37.0000  -36.1421
  -36.0000  -36.1421
  -35.0000  -36.1421
  -34.0000  -33.0019
  -33.0000  -33.0019
  -32.0000  -33.0019
  -31.0000  -29.8619
  -30.0000  -29.8619
  -29.0000  -29.8619
  -28.0000  -26.7222
  -27.0000  -26.7222
  -26.0000  -26.7222
  -25.0000  -23.5831
  -24.0000  -23.5831
  -23.0000  -23.5831
  -22.0000  -23.5831
  -21.0000  -20.4448
  -20.0000  -20.4448
  -19.0000  -20.4448
  -18.0000  -17.3076
  -17.0000  -17.3076
  -16.0000  -17.3076
  -15.0000  -14.1724
  -14.0000  -14.1724
  -13.0000  -14.1724
  -12.0000  -11.0408
  -11.0000  -11.0408
  -10.0000  -11.0408
-9.0000 -7.9171
-8.0000 -7.9171
-7.0000 -7.9171
-6.0000 -4.8158
-5.0000 -4.8158
-4.0000 -4.8158
-3.0000 -1.8366
-2.0000 -1.8366
-1.0000 -1.8366
      0       0
1.0000 1.8366
2.0000 1.8366
3.0000 1.8366
4.0000 4.8158
5.0000 4.8158
6.0000 4.8158
7.0000 7.9171
8.0000 7.9171
9.0000 7.9171
10.0000 11.0408
11.0000 11.0408
12.0000 11.0408
13.0000 14.1724
14.0000 14.1724
15.0000 14.1724
16.0000 17.3076
17.0000 17.3076
18.0000 17.3076
19.0000 20.4448
20.0000 20.4448
21.0000 20.4448
22.0000 23.5831
23.0000 23.5831
24.0000 23.5831
25.0000 23.5831
26.0000 26.7222
27.0000 26.7222
28.0000 26.7222
29.0000 29.8619
30.0000 29.8619
31.0000 29.8619
32.0000 33.0019
33.0000 33.0019
34.0000 33.0019
35.0000 36.1421
36.0000 36.1421
37.0000 36.1421
38.0000 39.2826
39.0000 39.2826
40.0000 39.2826
41.0000 42.4233
42.0000 42.4233
43.0000 42.4233
44.0000 45.5641
45.0000 45.5641
46.0000 45.5641
47.0000 45.5641
48.0000 48.7050
49.0000 48.7050
50.0000 48.7050
51.0000 51.8459
52.0000 51.8459
53.0000 51.8459
54.0000 54.9870
55.0000 54.9870
56.0000 54.9870
57.0000 58.1281
58.0000 58.1281
59.0000 58.1281
60.0000 61.2692
61.0000 61.2692
62.0000 61.2692
63.0000 64.4104
64.0000 64.4104
65.0000 64.4104
66.0000 67.5516
67.0000 67.5516
68.0000 67.5516
69.0000 67.5516
70.0000 70.6929
71.0000 70.6929
72.0000 70.6929
73.0000 73.8342
74.0000 73.8342
75.0000 73.8342
76.0000 76.9755
77.0000 76.9755
78.0000 76.9755
79.0000 80.1169
80.0000 80.1169
81.0000 80.1169
82.0000 83.2582
83.0000 83.2582
84.0000 83.2582
85.0000 86.3996
86.0000 86.3996
87.0000 86.3996
88.0000 89.5410
89.0000 89.5410
90.0000 89.5410
91.0000 89.5410
92.0000 92.6824
93.0000 92.6824
94.0000 92.6824
95.0000 95.8238
96.0000 95.8238
97.0000 95.8238
98.0000 98.9652
99.0000 98.9652
  100.0000 98.9652

附图1.jpg

赞一下

回复此楼

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

6楼2013-07-23 17:02:53

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

kuntenghunya

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 98.8
散金: 15
帖子: 25
在线: 47.8小时
虫号: 1818843
注册: 2012-05-15
性别: MM
专业: 传热传质学

引用回帖:

6楼: Originally posted by 月只蓝 at 2013-07-23 17:02:53
5楼的虫友说得很对，应当考虑多根的情况。附图表示了在，μ存在多解。
根据实际中对μ值的要求，给定μ值的取值区间，可按下程序解得μ(以为例)：

function sol_nol
clear all;clc
format short
n=100; %μ ...

多谢多谢

回复此楼

7楼2013-07-23 18:10:23

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

kuntenghunya

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 98.8
散金: 15
帖子: 25
在线: 47.8小时
虫号: 1818843
注册: 2012-05-15
性别: MM
专业: 传热传质学

引用回帖:

5楼: Originally posted by pippi6 at 2013-07-23 16:18:30
在excel文件中，修改 K值，比如 1，就得到第一个根 1.836597203，修改成 10，就得到第10个根 29.8618724，等等。当然，根是对称的。 -29.8618724 也是根。前十个根是（除了平凡解 x=0）

1 1.836597203
2 4.8 ...

多谢

回复此楼

8楼2013-07-23 18:10:53

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

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

【答案】应助回帖

★
csgt0: 金币+1, 谢谢 2013-07-24 16:59:37

6楼的matlab程序很全面啊。不过对于k-th 根，迭代初值可选择为 1.8+k*pi。
顺便补充一点，我是用Newton法，使用excel做的。如果按原来的形式写 f= x/tan(x)-(1-Bi), 除了第一个根外，其他的根不容易收敛。写成
f= 1/tan(x)-(1-Bi)/x，
df/dx=-1/sin(x)^2+(1-Bi)/x^2,
x= x-f/(df/dx)
则要好得多。

赞一下

回复此楼

9楼2013-07-23 18:52:51

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

kuntenghunya

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 98.8
散金: 15
帖子: 25
在线: 47.8小时
虫号: 1818843
注册: 2012-05-15
性别: MM
专业: 传热传质学

引用回帖:

9楼: Originally posted by pippi6 at 2013-07-23 18:52:51
6楼的matlab程序很全面啊。不过对于k-th 根，迭代初值可选择为 1.8+k*pi。
顺便补充一点，我是用Newton法，使用excel做的。如果按原来的形式写 f= x/tan(x)-(1-Bi), 除了第一个根外，其他的根不容易收敛。写成
f= ...

多谢给予帮助，超越方程的第一个根满足我的需要了

赞一下

回复此楼

10楼2013-07-23 19:33:42

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

相关版块跳转我要订阅楼主 kuntenghunya 的主题更新

返回列表