| 查看: 2699 | 回复: 12 | |||
[求助]
matlab求三元4次符号代数方程组，出现warning且答案不正确，该怎么办？ 已有1人参与
|
是为了求一个截断椭圆的曲线方程，椭圆长半轴、短半轴以及中心点的高度分别是 未知量 a, b 和 m. 截断椭圆.jpg 椭圆方程因此为：  椭圆方程.jpg 另 x=0, 可得到 Rco； 对曲线方程求导后另 x=0， 可得到 tan(theta);  方程组.jpg 三个未知数 a、b、m，三个方程①②③，用matlab 求解符号代数方程的函数 “solve（‘equ1’，‘equ2’，…‘var1’，‘var2’…） ”，却得不到结果，这是为什么？请指教，谢谢 （程序中 theta用 c 表示， Rco 用 r 表示，Req 用 p 表示）  求指教.jpg [ Last edited by yihu_55 on 2014-2-24 at 16:27 ] |
回复此楼» 猜你喜欢
面上意见出来了
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB求解代数方程组问题，请帮忙
已经有9人回复
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2014-02-24 16:28:55
3楼2014-02-24 16:40:49
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
yihu_554楼2014-02-24 16:45:58
5楼2014-02-24 16:58:20
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2014-02-24 17:06:54
7楼2014-02-24 17:16:52
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yihu_55: 金币+15, ★★★很有帮助 2014-02-24 22:13:42
yihu_55: 金币+15, ★★★很有帮助 2014-02-24 22:13:42
|
我试了一下 可以求解出来 你自己对照一下代码 楼主给币币 哈哈 代码： clear; clc; syms a b m c r p; f=tan(c)-r*(b^2)/(a^2)/2; g=r-a*sqrt((b^2-m^2))/b; k=a^2*(3*m*b^2+2*b^3-m^3)/(b^2)-4^p^3; [d,e,s]=solve(f,g,k) 求解结果： d = (b*(-(4*a^4*tan(c)^2 - a^2*b^4)/b^4)^(1/2))/a -(b*(-(4*a^4*tan(c)^2 - a^2*b^4)/b^4)^(1/2))/a e = log((b*(2*a - (-(4*a^4*tan(c)^2 - a^2*b^4)/b^4)^(1/2))*(a + (-(4*a^4*tan(c)^2 - a^2*b^4)/b^4)^(1/2))^2)/a)/(6*log(2)) log((b*(a - (-(4*a^4*tan(c)^2 - a^2*b^4)/b^4)^(1/2))^2*(2*a + (-(4*a^4*tan(c)^2 - a^2*b^4)/b^4)^(1/2)))/a)/(6*log(2)) s = (2*a^2*tan(c))/b^2 (2*a^2*tan(c))/b^2 |
8楼2014-02-24 21:03:02
9楼2014-02-24 22:15:48
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2014-02-24 22:30:20
30