24小时热门版块排行榜

返回列表

当前主题已经存档。

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

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

建议用数值方法，t、t+0.007、t+0.013、t+0.016差别不大，得出的四个cos值差别也不会太大，其大小顺序和等式右边值不同，这地方可以好好讨论一下，我想可以得出：
１：a,b平方和在50的平方左右
２：8个cos值都在1附近
３：t和0左右，周期是1，n在pi左右，周期是pi，且范围都很窄
我写了这样的一个函数：
function y=fun(a,b,t,n)
y(1,

=(a.*cos(2.*pi.*t)).^2 +(b.*cos(2.*pi.*t+n)).^2 -(54.18525).^2;
y(2,

=(a.*cos(2.*pi.*(t+0.013))).^2+(b.*cos(2.*pi.*(t+0.013)+n)).^2-(53.01687).^2;
y(3,

=(a.*cos(2.*pi.*(t+0.007))).^2+(b.*cos(2.*pi.*(t+0.007)+n)).^2-(51.56303).^2;
y(4,

=(a.*cos(2.*pi.*(t+0.016))).^2+(b.*cos(2.*pi.*(t+0.016)+n)).^2-(50.81597).^2;
可以作二维的或是三维的图，找归律，找零点，我搞了半天也没搞出来，只能把四个函数值控制在了正负200以内，不过这个思路我觉得应该还是可以的

这是我找到的比较好的一组解，周围值都差不多
>> y=fun(47,47,0.0966,3.14)

y =

47.9831
-173.5101
140.5550
-26.9664

[ Last edited by fspdlh on 2009-3-25 at 08:55 ]