24小时热门版块排行榜    

查看: 737  |  回复: 4
当前主题已经存档。

abing2807

银虫 (小有名气)

[交流] 【求助】请问 cos(x)ch(x)=1 如何解?

这个方程有无穷多个解,x1=4.730041,x2=......
我想求出前20阶的解,如何解呢?

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

redcrag

木虫 (著名写手)


wenzhenzhong(金币+1,VIP+0)::),感谢参与 5-6 00:05
自己编程还是用软件?
用软件(如Matlab、Mathmatica等)是比较容易的,定义一个cos(x)ch(x)-1的函数,画图看一下它的零点,然后调用求零的函数,可以设置初始猜测值或者给定求解区间求解。
http://redcrag.spaces.live.com/
2楼2009-05-06 00:00:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

看看这个吧,把方程改一下就可以了
http://muchong.com/bbs/viewthread.php?tid=1251344
function x=solvefun(a,b,tol)

%--------------------------------------------------------------------------
range=a:1:b;
if range(end)     range=[range,b];
end
tag=[abs(diff(sign(fun(range)))) 0];
range=range(find(tag>0))';
range=[range range+1];
%--------------------------------------------------------------------------
n=size(range,1);
tol=tol/10;
x=zeros(n,1);
for i=1:n
    xmin=range(i,1);
    xmax=range(i,2);
    xcur=(xmin+xmax)/2;
    while xmax-xmin>tol
        if sign(fun(xmin))*sign(fun(xcur))>0
            xmin=xcur;
        else
            xmax=xcur;
        end
        xcur=(xmin+xmax)/2;
    end
    x(i)=xcur;
end
%--------------------------------------------------------------------------
function y=fun(x)
y=cos(x).*cosh(x)-1;

>> solvefun(0,200,0.0001)

ans =

    0.0000
    4.7300
    7.8532
   10.9956
   14.1372
   17.2788
   20.4204
   23.5619
   26.7035
   29.8451
   32.9867
   36.1283
   39.2699
   42.4115
   45.5531
   48.6947
   51.8363
   54.9779
   58.1195
   61.2611
   64.4027
   67.5442
   70.6858
   73.8274
   76.9690
   80.1106
   83.2522
   86.3938
   89.5354
   92.6770
   95.8186
   98.9602
  102.1018
  105.2434
  108.3849
  111.5265
  114.6681
  117.8097
  120.9513
  124.0929
  127.2345
  130.3761
  133.5177
  136.6593
  139.8009
  142.9425
  146.0841
  149.2257
  152.3672
  155.5088
  158.6504
  161.7920
  164.9336
  168.0752
  171.2168
  174.3584
  177.5000
  180.6416
  183.7832
  186.9248
  190.0664
  193.2080
  196.3495
  199.4911

>>

[ Last edited by fspdlh on 2009-5-7 at 00:13 ]
3楼2009-05-06 23:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhfzh

木虫 (正式写手)

呵呵

可以把方程化为复数方程.(用COS(X)转换)

[ Last edited by zhfzh on 2009-5-16 at 10:46 ]
4楼2009-05-16 10:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小作神童


小木虫(金币+0.5):给个红包,谢谢回帖交流
可以把方程化为复数方程.(用COS(X)转换)
5楼2009-08-09 13:26:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 abing2807 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 20:23 by Kddn1pB6pAzX
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 20:22 by Kddn1pB6pAzX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 20:17 by Kddn1pB6pAzX
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 20:11 by Kddn1pB6pAzX
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 13/650 2026-08-01 20:02 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jzg3xk1po8 2026-07-31 13/650 2026-08-01 20:02 by Kddn1pB6pAzX
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 3/150 2026-08-01 19:00 by Kddn1pB6pAzX
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 3/150 2026-08-01 18:57 by Kddn1pB6pAzX
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 3/150 2026-08-01 18:51 by Kddn1pB6pAzX
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[博后之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:02 by Kddn1pB6pAzX
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见