| 查看: 2721 | 回复: 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人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB求解代数方程组问题，请帮忙
已经有9人回复
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
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
投票:
20