| 查看: 543 | 回复: 1 | ||||
huier168木虫 (正式写手)
金虫
|
[求助]
matlab编程,求助,如下: 已有1人参与
|

wurongjun
专家顾问 (职业作家)
-

专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
huier168: 金币+30, ★有帮助 2015-09-28 12:52:37
感谢参与,应助指数 +1
huier168: 金币+30, ★有帮助 2015-09-28 12:52:37
|
问题比较麻烦 首先,t>0不能任意取,必须使0.5-cos(t)>0!(对角线的元素必须大于零) 在这个范围内再算二阶子式: A = [ -cos(t)+1/2, 13/10-sin(t), 2/5] [ 13/10-sin(t), ta+1, 7/10] [ 2/5, 7/10, ta-2] >> D2=factor(det(A(1:2,1:2))) D2 = -cos(t)*ta-cos(t)+1/2*ta-119/100+13/5*sin(t)-sin(t)^2 由D2>0可得 >> solve(D2,'ta') ans = -1/50*(100*cos(t)+119-260*sin(t)+100*sin(t)^2)/(2*cos(t)-1) ta大于上式,但是基本上上式递增到无穷! 再看行列式大于零 >> D=factor(det(A)) D = -cos(t)*ta^2+cos(t)*ta+249/100*cos(t)+1/2*ta^2-47/20*ta+2703/1000+13/5*sin(t)*ta-144/25*sin(t)-sin(t)^2*ta+2*sin(t)^2 >> D=collect(D,ta) D = (-cos(t)+1/2)*ta^2+(-47/20-sin(t)^2+13/5*sin(t)+cos(t))*ta+2703/1000+249/100*cos(t)-144/25*sin(t)+2*sin(t)^2 >> D>0的范围大致与D2的情况相似! 也就是说没有最大值! |

2楼2015-09-13 17:23:43













回复此楼