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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 特征方程求解
51/1返回列表
查看: 1218  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭