24小时热门版块排行榜    

查看: 1420  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

pdeliuyong

银虫 (小有名气)

[交流] 【求助】求方程的正根

>> syms a b c t;
>> eq1=a*(cos(2*pi*(t+0.0001)))^2+b*(cos(2*pi*(t+0.0001)+c))^2-54.1852^2;
>> eq2=a*(cos(2*pi*(t+0.0002)))^2+b*(cos(2*pi*(t+0.0002)+c))^2-53.95478^2;
>> eq3=a*(cos(2*pi*(t+0.0003)))^2+b*(cos(2*pi*(t+0.0004)+c))^2-53.72272^2;
>> eq4=a*(cos(2*pi*(t+0.0004)))^2+b*(cos(2*pi*(t+0.0004)+c))^2-53.48904^2;
>> S=solve(eq1,eq2,eq3,eq4,a,b,c,t);
>> S.a,S.b,S.c,S.t
我想求上面4个方程的正根,上面的程序要如何添加

[ Last edited by woshilsh on 2009-5-22 at 13:03 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pdeliuyong

银虫 (小有名气)

恢复2楼

按照你的要求改了程序
>> syms a b c t;
eq1=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0001)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0001)+(sqrt(c))^2))^2-54.1852^2;
eq2=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0002)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0002)+(sqrt(c))^2))^2-53.95478^2;
eq3=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0003)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0003)+(sqrt(c))^2))^2-53.72272^2;
eq4=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0004)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0004)+(sqrt(c))^2))^2-53.48904^2;
S=solve(eq1,eq2,eq3,eq4,a,b,c,t);
S.a,S.b,S.c,S.t
Warning: Explicit solution could not be found.
> In solve at 170
??? Error using ==> subsref
Reference to non-existent field 'a'.

Error in ==> sym.subsref at 16
  y = builtin('subsref',struct(x),a);
4楼2009-03-22 10:37:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

fspdlh

金虫 (正式写手)

★ ★
kuhailangyu(金币+2,VIP+0):谢谢参与! 3-22 20:48
把要搞成正根的未知数设为另一个数的平方啊
不过你这四个方程四个未知数刚好啊,如果那样有解吗?
2楼2009-03-22 10:09:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pdeliuyong

银虫 (小有名气)

方程正根

按照你的作了,不行
syms a b c t;
eq1=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0001)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0001)+(sqrt(c))^2))^2-54.1852^2;
eq2=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0002)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0002)+(sqrt(c))^2))^2-54.1852^2;
eq3=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0003)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0003)+(sqrt(c))^2))^2-54.1852^2;
eq4=(sqrt(a))^2*(cos(2*pi*((sqrt(t))^2+0.0004)))^2+(sqrt(b))^2*(cos(2*pi*((sqrt(t))^2+0.0004)+(sqrt(c))^2))^2-54.1852^2;
S=solve(eq1,eq2,eq3,eq4,a,b,c,t);
S.a,S.b,S.c,S.t

ans =

2936.0358990400000038789585232735


ans =

2936.0358990400000038789585232735


ans =

413.11943394705781085783760490125


ans =

-170.00030647677566299912270803713
3楼2009-03-22 10:28:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

说明你这个方程四个值都为正的解不存在

[ Last edited by fspdlh on 2009-3-22 at 10:50 ]
5楼2009-03-22 10:49:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见