maple 求maple求解方程组组解析解

这类方程应该算是不定方程,在给定条件下只有一组解,下面的Maple代码怎么运行不了呢?或者能用其他相关软件算出来也行,比如Matlab,R语言等solve({a^2+b^2+169+sqrt(c-13)-24*a-10*b = 0},{a, b, c})
assuming a&0, b&0, c&0;
solve([1/(cos(a)^2)+1/(sin(a)^2*sin(b)^2*cos(b)^2) = 9,
a&0, a&Pi/2, b&0, b&Pi/2], [a,b] );
MATLAB 的 Symbolic Math Toolbox 可以通过 solve 用解析方法解代数方程,不过如果条件是不定的,可能不容易快速得到结果。也可以用 vpasolve (我用的R2013a) 来用数值方法解,但没用过不知道不定条件下如何。其实 MATLAB 最强大的就是数值计算,这类简单问题可以用数值逼近方法来解决。当然只是一种方法,也许不能得到解析结果,也许运气不好遇到局部极值啥的,权作参考。对于这两题,都可以将变量限定在某个范围内,设定一个精度,用数值模拟求得最小结果来解。第二题相对简单。clear;
range = [0 : 0.001 : pi / 2];
for a = range
for b = range
index_a = find(a == range);
index_b = find(b == range);
c(index_a, index_b) = 1 / (cos(a) ^ 2) + ...
1 / (sin(a) ^ 2 * sin(b) ^ 2 * cos(b) ^ 2) - 9;
solution = abs(c);
[y, x] = find(solution == min(min(solution)));
solution_a = range(y);
solution_b = range(x);
solution_nv = c(y, x);
fprintf('a = %8.5f\nb = %8.5f\nnearest equation value = %8.5f\n', ...
radtodeg(solution_a), radtodeg(solution_b), solution_nv);
得到结果:a = 54.71747
b = 44.97719
nearest equation value =
第二题涉及到三角形三边关系,某一边的条件,以及变量范围的问题,姑且假设三边均不超过200,在粗糙的精度下运气好也能得到正确答案:clear;
range_a = [1 : 1 : 200];
range_b = [1 : 1 : 200];
range_c = [13 : 1 : 200];
for a = range_a
for b = range_b
for c = range_c
index_a = find(a == range_a);
index_b = find(b == range_b);
index_c = find(c == range_c);
if index_a + index_b & index_c || ...
index_a + index_c & index_b || ...
index_b + index_c & index_a
d(index_a, index_b, index_c) = a ^ 2 + b ^ 2 + 169 + ...
sqrt(c - 13) - 24 * a - 10 * b;
d(index_a, index_b, index_c) = 10000;
solution = abs(d);
[z, y, x] = find(solution == min(min(min(solution))));
solution_a = range_a(z);
solution_b = range_b(y);
solution_c = range_c(x);
solution_nv = d(z, y, x);
fprintf(['a = %8.5f\nb = %8.5f\nc = %8.5f\n'...
'nearest equation value = %8.5f\n'],...
solution_a, solution_b, solution_c, solution_nv);
得到结果:a = 12
nearest equation value =
原来恰巧就是直角三角形。
我不知道 maple 怎么算,mathematica 是可以算的。&br&第一个问题:&br&&div class=&highlight&&&pre&&code class=&language-text&&Solve[{a^2 + b^2 + 169 + Sqrt[c - 13] - 24 a - 10 b == 0, a &= 0,
b &= 0, c &= 0}, {a, b, c}]
&/code&&/pre&&/div&&br&能得出正确答案&br&第二个问题:&br&&div class=&highlight&&&pre&&code class=&language-text&&Solve[{1/(a b (1 - b)) + 1/(1 - a) == 9, a &= 0, b &= 0, a &= 1,
b &= 1}, {a, b}]
&/code&&/pre&&/div&&br&也能得出正确答案(你应该能看出这个形式和原问题是什么关系吧?)&br&&br&既然你会手算,答案我就不贴了。另外……虽然题主说不是作业题,但是这种凑得很好的题目很难让人相信不是作业题。
我不知道 maple 怎么算,mathematica 是可以算的。第一个问题:Solve[{a^2 + b^2 + 169 + Sqrt[c - 13] - 24 a - 10 b == 0, a &= 0,
b &= 0, c &= 0}, {a, b, c}]能得出正确答案第二个问题:Solve[{1/(a b (1 - b)) + 1/(1 - a) == 9, a &= 0, b &= 0, a…
数学问题果断使用Mathematica。第一个用FindInstance,第二个用Solve。&br&&img src=&/c874c990fdb2bfaacbcb9f12_b.jpg& data-rawwidth=&861& data-rawheight=&358& class=&origin_image zh-lightbox-thumb& width=&861& data-original=&/c874c990fdb2bfaacbcb9f12_r.jpg&&
数学问题果断使用Mathematica。第一个用FindInstance,第二个用Solve。
已有帐号?
无法登录?
社交帐号登录
写代码的 / 跑步爱好者 / 大脑骇客 / 登山爱好者Maple解方程组后,用eval验证解,不知错在哪里,请指教f := {(x-y)^2-3*(x-y) = 10, x^2+2*x*y+y^2 = 9}:g := solve(f):g {x = -5/2, y = -1/2}, {x = 1, y = -4}, {x = 1/2, y = 5/2}, {x = 4, y = -1}然后验证:eval(f, g) Error, invalid input: eval received ({x = -5/2, y = -1/2}, {x = 1, y = -4}, {x = 1/2, y = 5/2}, {x = 4, y = -1}), which is not valid for its 2nd argument, eqns但是用subs验证却是对的subs(g, f) {9 = 9, 10 = 10}解了几个其他方程组也是这样的过程,都没有错,不知道为什么会出现这样的结果,高手指教呀!
张起灵1053
不对应 可以这样> eval(f,g[1]); 一个一个测试
为您推荐:
扫描下载二维码maple解复杂方程组_百度知道
maple解复杂方程组
com/zhidao/pic/item/eb086e061d95f725.1e-2*[sin(c+d+e)-sin(e)] = 0}.1e-2*(sin(e+c)-sin(e))/sin(2*e+c+d)-0://d,α=20.jpg" target="_blank" title="点击查看大图" class="ikqb_img_alink"><img class="ikqb_img" src="/zhidao/wh%3D600%2C800/sign=c50af43eb1003af34defd466051aea64/eb086e061d95f725.2e-2*sin(d)/sin(e+d)+0.728e-1*(1/a-1/b)请问一下这个方程组怎么求解啊.2e-2*sin(e)/sin(e+d)]/sin(2*e+c+d)+0.hiphotos, e = 20.1e-2*[sin(e+c)-sin(e)]*[sin(e+d)-1]/sin(2*e+c+d)+sin(e+c)*[sin(e+d)-1]*[1-0.baidu,麻烦大神给出算式代码好么://d
提问者采纳
应该要有6个方程才能求解.jpg" />你这个问题涉及到6个未知量.jpg" esrc="/zhidao/pic/item/bf096b63fc43ecf81a4c510fa25f.jpg" target="_blank" title="点击查看大图" class="ikqb_img_alink"><img class="ikqb_img" src="http<a href="http://g://g,你这才4个.baidu.hiphotos./zhidao/wh%3D450%2C600/sign=fb9d168aee98b9/bf096b63fc43ecf81a4c510fa25f./zhidao/wh%3D600%2C800/sign=53f13e32e2dde711edfe223/bf096b63fc43ecf81a4c510fa25f://g
没有,有些要赋值,R1=0.001,R2=0.002,阿尔法=20
提问者评价
来自团队:
为您推荐:
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁(Albertjun)
第三方登录:}

我要回帖

更多关于 maple求解方程组 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信