24小时热门版块排行榜

>论坛更新日志 (2057)
>导师招生 (368)
>硕博家园 (172)
>虫友互识 (126)
>招聘信息布告栏 (92)
>考博 (59)
>博后之家 (23)
>基金申请 (7)
>生物科学 (6)
>考研 (6)
>教师之家 (5)
>晶体 (4)
>论文投稿 (4)
>电化学 (3)
>找工作 (3)
>公派出国 (3)

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

abo

木虫 (正式写手)

应助: 5 (幼儿园)
金币: 4383.6
散金: 429
红花: 2
帖子: 600
在线: 338.5小时
虫号: 65387
注册: 2005-04-24
性别: GG
专业: 自动化

[交流] ［求助］１０１个今币求解

'(a*cos(2*pi*t))^2+(b*cos(2*pi*t+n))^2=(54.18525)^2',
'(a*cos(2*pi*(t+0.013)))^2+(b*cos(2*pi*(t+0.013)+n))^2=(53.01687)^2',
'(a*cos(2*pi*(t+0.007)))^2+(b*cos(2*pi*(t+0.007)+n))^2=(51.56303)^2',
'(a*cos(2*pi*(t+0.016)))^2+(b*cos(2*pi*(t+0.016)+n))^2=(50.81597)^2')

以上方程组如何解决啊
下面是用ＭＡＴＬＡＢ编写的程序

syms('a','b','n','t')
[a,b,n,t]=solve('(a*cos(2*pi*t))^2+(b*cos(2*pi*t+n))^2=(54.18525)^2','(a*cos(2*pi*(t+0.013)))^2+(b*cos(2*pi*(t+0.013)+n))^2=(53.01687)^2','(a*cos(2*pi*(t+0.007)))^2+(b*cos(2*pi*(t+0.007)+n))^2=(51.56303)^2', '(a*cos(2*pi*(t+0.016)))^2+(b*cos(2*pi*(t+0.016)+n))^2=(50.81597)^2')

结果：Warning: Explicit solution could not be found.
不要很精确的解，如何求解啊
谢谢

回复此楼

» 猜你喜欢

课题组招2027级博士上海工程技术大学激光智能制造方向已经有4人回复
现代”学阀”该如何界定已经有8人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有3人回复
实名举报：华南理工大学物理与光电学院副院长李志远婚内两次出轨女学生（博士）已经有16人回复
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告已经有11人回复
上海工程技术大学激光智能制造课题组招收博士研究生已经有10人回复
申博发邮件已经有11人回复
各位大神，目前国内有哪些比较好用的逆合成软件？已经有12人回复
我的奶奶已经有3人回复
国社科系统bug了，是不是要放榜了？已经有9人回复

1楼 2009-03-24 09:12:36

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

fspdlh

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

用http://muchong.com/bbs/viewthread.php?tid=1311004&fpage=1上的粒子群优化算法和遗传算法搞了一下，也不好搞，再顶一下这个帖子，期待高手，这方程不会是无解吧

function y=fun(x)
a=x(1);
b=x(2);
t=x(3);
n=x(4);
y=0;
y=y+( (a*cos(2*pi*t))^2          + (b*cos(2*pi*t+n))^2       - (54.18525)^2       )^2;
y=y+( (a*cos(2*pi*(t+0.013)))^2 + (b*cos(2*pi*(t+0.013)+n))^2 - (53.01687)^2       )^2;
y=y+( (a*cos(2*pi*(t+0.007)))^2 + (b*cos(2*pi*(t+0.007)+n))^2 - (51.56303)^2       )^2;
y=y+( (a*cos(2*pi*(t+0.016)))^2 + (b*cos(2*pi*(t+0.016)+n))^2 - (50.81597)^2'    )^2;

>> [x y]=YSPSO(@fun,30,2.8,1.3,10000,4)

x =

  1.0e+002 *

1.56387432102746
  -0.61124242527700
0.15221726954985
0.10463512319258

y =

3.330657386605306e+004

>> [x y]=YSPSO(@fun,40,2.05,2.05,10000,4)

x =

  1.0e+002 *

1.34814607821837
  -1.00091039452162
  -0.22801667138287
  -0.15051297762441

y =

3.330657386602951e+004

>> [x y]=AsyLnCPSO(@fun,10,2.5,0.5,0.5,2.5,0.7,100,4)

x =

-58.94847244528419
  -4.51791509340913
  16.06831925267567
-43.81674017186624

y =

3.808986350835356e+004

>> [x y]=SecVibratPSO(@fun,40,0.78,1,1,10000,4)

x =

  1.0e+002 *

  -1.89352933572822
  -0.80524286901973
0.03722412069556
  -0.21252264587985

y =

3.431387208971817e+004

>> [x y]=BreedPSO(@fun,40,2,2,0.7,0.9,0.2,10000,4)

x =

  1.0e+002 *

1.84970839650448
  -0.51723006683137
  -3.62265092903016
  -5.92246460256721

y =

3.365643828984751e+004