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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 特征方程求解
51/1返回列表
查看: 1293  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kuntenghunya

铁虫 (初入文坛)

[求助] 特征方程求解

求高手帮助，解特征方程
特征方程求解
特征方程.jpg
回复此楼

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

北大科研分享

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-07-23 11:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
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
附图1.jpg
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2013-07-23 17:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★
感谢参与，应助指数 +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

铁虫 (初入文坛)
引用回帖:
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的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

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

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

Maple的结果和MATLAB是一致的：
特征方程求解-2
附图2.jpg
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2013-07-23 13:43:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +4 小蜡新笔 2026-03-26 4/200 2026-03-26 12:33 by ms629
[考研] 299求调剂 +4 15188958825 2026-03-25 4/200 2026-03-25 22:56 by 418490947
[考研] 333求调剂 +6 87639 2026-03-21 11/550 2026-03-25 16:17 by 87639
[考研] 材料与化工考研调剂 +9 孅華 2026-03-22 9/450 2026-03-25 13:09 by cmz0325
[考研] 359求调剂 +3 王了个楠 2026-03-25 3/150 2026-03-25 12:50 by Dyhoer
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 一志愿211 初试270分 求调剂 +5 谷雨上岸 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 一志愿深大，0703化学，总分302，求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
信息提示
关闭
请填处理意见
关闭