| 查看: 1404 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】求方程的正根
|
|||
|
>> 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 ] |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
6楼2009-03-23 20:43:46
2楼2009-03-22 10:09:14
方程正根
|
按照你的作了，不行 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
恢复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
