| 查看: 1335 | 回复: 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 ] |
回复此楼» 猜你喜欢
生物学303求调剂，一志愿华中农微生物，六级已过，有科研有文章，党员
已经有4人回复
-
318求调剂
已经有5人回复
-
0703调剂
已经有11人回复
-
化学357分，考研调剂
已经有10人回复
-
0854求调剂
已经有4人回复
-
272求调剂
已经有3人回复
-
调剂求助
已经有13人回复
-
332求调剂
已经有8人回复
-
341求调剂
已经有3人回复
-
一志愿电子科技大学085600材料与化工 329分求调剂
已经有10人回复
方程正根
|
按照你的作了，不行 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楼2009-03-22 10:09:14
恢复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
5楼2009-03-22 10:49:04
