关于Mathematica解方程组的问题，输出时没有错误，解不出来的原因_百度知道
关于Mathematica解方程组的问题，输出时没有错误，解不出来的原因
如题，我在一个三元一次方程组时，应该能解出来的，但是老是得不到结果。由于我刚刚接触Mathematica这个软件，不懂，所以我把解得编程给大家，希望大家能够帮帮我，看看问题出在哪里？f = 60 D
arfA = 300 D
费了好半天劲才查出来你的逗号是中文的！请看区别:，,左边的是中文全角的。
采纳率：59%
为您推荐：
其他类似问题
您可能关注的内容
mathematica的相关知识
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。100分，谁会用Mathematica解这个参数方程组啊!写出代码也行，能解决的话再加一百。_百度知道
100分，谁会用Mathematica解这个参数方程组啊!写出代码也行，能解决的话再加一百。
syms a b c x1 y1 z1eq1=sym('cosb*cosc*x+sina*sinb*cosc*y-cosa*sinc*y+cosa*sinb*cosc*z+sina*sinc*z=x1') eq2=sym('cosb*sinc*x+sina*sinb*sinc*y+cosa*cosc*y+cosa*sinb*sinc*z-sina*cosc*z=y1')eq3=sym('-sinb*x+sina*cosb*y+cosa*cosb*z=z1')[...
其实原问题是这样的：以球心建立三维直角坐标系，球面上有一点A（x,y,z）由于某种原因球发生了绕球心的旋转，A点坐标变为（x1,y1,z1）。问题是：怎么旋转坐标系使A点在旋转后的坐标系内的坐标也是（x,y,z）.
我是想这样解决的：设绕xyz轴旋转的角度一次为abc，利用旋转矩阵可得式子（x1,y1,z1，1）=（x,y,z，1）Gx*Gy*Gz ,
进而求出旋转角度abc。。。。这样解没问题吧，可怎么就是解不出啊？？？？？？？？？？
我有更好的答案
不懂Matlab代码，你最好能直接给原方程。不过……也大致猜得到含义：这是关于a，b，c的三个三角方程是吧？那么，如果1 你希望求得的是符号解而非数值解2 你还不能对a，b，c乃至参数x, y, z, x1, y1, z1的定义域做出比较好的限制那么，我可以负责任的告诉你，这个方程不管用什么软件都是不好解的，因为三角函数本身的周期性会使方程通解变得非常复杂，关于这一部分，你可以参考这里：zhidao.baidu.com/question/.html?oldq=1&from=evaluateTo#reply-box-此外，除了限制未知量的定义域之外，另一个可行的方法，是适当换元，规避三角函数。一个可用的换元是万能公式：(*&注意变形前后的a,&b,&c含义不同&*)eqn&=&{Cos[b]*Cos[c]*x&+&Sin[a]*Sin[b]*Cos[c]*y&-&Cos[a]*Sin[c]*y&+&Cos[a]*Sin[b]*Cos[c]*z&+&Sin[a]*Sin[c]*z&==&x1,&Cos[b]*Sin[c]*x&+&Sin[a]*Sin[b]*Sin[c]*y&+&Cos[a]*Cos[c]*y&+&Cos[a]*Sin[b]*Sin[c]*z&-&Sin[a]*Cos[c]*z&==&y1,&-Sin[b]*x&+&Sin[a]*Cos[b]*y&+&Cos[a]*Cos[b]*z&==&z1}&/.&{Cos[a_]&-&&(1&-&a^2)/(1&+&a^2),&Sin[a_]&-&&2&a/(1&+&a^2)}//Simplify&&&&&&(*&{((1&+&a^2)&(-1&+&b^2)&(-1&+&c^2)&x&+&2&(-(1&+&b^2)&c&y&+&b&z&-&b&c^2&z&+&2&a&(b&(y&-&c^2&y)&+&c&z&+&b^2&c&z)&+&a^2&((1&+&b^2)&c&y&-&b&z&+&b&c^2&z)))/((1&+&a^2)&(1&+&b^2)&(1&+&c^2))&==&x1,&(-(1&+&b^2)&((-1&+&a^2)&y&+&2&a&z)&+&(1&+&b^2)&c^2&((-1&+&a^2)&y&+&2&a&z)&+&2&c&(-(1&+&a^2)&(-1&+&b^2)&x&+&2&b&(2&a&y&+&z&-&a^2&z)))/((1&+&a^2)&(1&+&b^2)&(1&+&c^2))&==&y1,&(-2&(1&+&a^2)&b&x&+&2&a&y&+&z&-&a^2&z&+&b^2&(-2&a&y&-&z&+&a^2&z))/((1&+&a^2)&(1&+&b^2))&==&z1}&*)显然，即使化成了多项式方程，这个方程组的阶数还是很高，高阶多项式（尤其是方程组）的求解也是很难的，但是，如果能够对参数的定义域也做出限制，确切地说，是给出参数的具体数值的话，这个方程的求解还是有希望的，比如：(*&这里随机生成一组参数&*){x,&y,&z,&x1,&y1,&z1}&=&RandomReal[{1,&2},&6](*&本次生成的参数是：{1.953,&1.659,&1.198}&*)(*&注意上面我已经给方程组取名为&eqn&了&*)Solve[eqn,&{a,&b,&c}](*&然后我们要换元换回来，注意这里只是简单地使用了反三角函数，要得到全部解，需要使用{Solve[Tan[d/2]&==&a,&d],&Solve[Tan[e/2]&==&b,&e],&Solve[Tan[f/2]&==&c,&f]}/.%&*){2&ArcTan[a],&2&ArcTan[b],&2&ArcTan[c]}&/.&%(*{{1.16575&+&18.1583&I,&-1.5708&-&1.&I,&-2.93477&-&&&&18.2374&I},&{1.16575&-&18.1583&I,&-1.5708&+&&&&1.&I,&-2.93477&+&18.2374&I},&{-1.97584&-&&&&18.1583&I,&-1.5708&-&1.&I,&&&0.206818&+&18.2374&I},&{-1.97584&+&18.1583&I,&-1.5708&+&&&&1.&I,&&&0.206818&-&18.2374&I},&{-0.810092&-&0.774135&I,&-0.538095&-&&&&35.6682&I,&-0.958933&-&0.890268&I},&{-0.810092&+&&&&0.774135&I,&-0.538095&+&35.6682&I,&-0.958933&+&&&&0.890268&I},&{2.3315&+&0.774135&I,&-2.50062&-&35.8224&I,&&&2.18266&+&0.890268&I},&{2.3315&-&0.774135&I,&-2.50062&+&35.8224&I,&&&2.18266&-&0.890268&I},&{2.3315&+&0.774135&I,&0.353705&-&37.366&I,&&&2.18266&-&0.890268&I},&{2.3315&-&0.774135&I,&0.353705&+&37.366&I,&&&2.18266&+&0.890268&I},&{-0.810092&-&0.774135&I,&&&1.5708&-&38.4263&I,&-0.958933&+&0.890268&I},&{-0.810092&+&&&&0.774135&I,&&&1.5708&+&38.4263&I,&-0.958933&-&0.890268&I},&{1.16575&+&18.1583&I,&&&1.5708&-&1.&I,&&&1.01691&+&18.2374&I},&{1.16575&-&18.1583&I,&&&1.5708&+&1.&I,&&&1.01691&-&18.2374&I},&{-1.97584&-&18.1583&I,&&&1.5708&-&1.&I,&-2.12468&-&18.2374&I},&{-1.97584&+&&&&18.1583&I,&1.5708&+&1.&I,&-2.12468&+&18.2374&I}}*)如你所见，依旧是非常复杂。当然了，使用Reduce求解这个方程的通解并非完全不可能。（毕竟在我写这个答案的这段时间程序一直在运行，并没有返回原式。）但是需要多久那就不知道了，你有兴趣可以试试：(*&这是最基本的代码&*)Reduce[{Cos[b]*Cos[c]*x&+&Sin[a]*Sin[b]*Cos[c]*y&-&Cos[a]*Sin[c]*y&+&&&&&Cos[a]*Sin[b]*Cos[c]*z&+&Sin[a]*Sin[c]*z&==&x1,&&Cos[b]*Sin[c]*x&+&Sin[a]*Sin[b]*Sin[c]*y&+&Cos[a]*Cos[c]*y&+&&&&&Cos[a]*Sin[b]*Sin[c]*z&-&Sin[a]*Cos[c]*z&==&y1,&&-Sin[b]*x&+&Sin[a]*Cos[b]*y&+&Cos[a]*Cos[b]*z&==&z1},&{a,&b,&c}]
我是想用x1,y1，z1表示a,b,c.
是常量，希望求得的是符号解而非数值解，，如果说条件的话：（x1,y1，z1）（
x ， y ， z ）都是球表面上的点球心（0,0,0）半径是330.
a,b,c.的范围0~π。
如果x, y, z是常量那希望能把具体值给我。
sorry ,我都不知道什么是常量了，它代表数值（算是已知的），但是并不是某一个具体的数字。
……那我也无能为力了。真要碰碰运气的话就把这个代码执行个十几小时试试：Reduce[{Cos[b]&Cos[c]&x&+&Sin[a]&Sin[b]&Cos[c]&y&-&Cos[a]&Sin[c]&y&+&Cos[a]&Sin[b]&Cos[c]&z&+&Sin[a]&Sin[c]&z&==&x1,&Cos[b]&Sin[c]&x&+&Sin[a]&Sin[b]&Sin[c]&y&+&Cos[a]&Cos[c]&y&+&Cos[a]&Sin[b]&Sin[c]&z&-&Sin[a]&Cos[c]&z&==&y1,&-Sin[b]&x&+&Sin[a]&Cos[b]&y&+&Cos[a]&Cos[b]&z&==&z1,&x^2&+&y^2&+&z^2&==&330^2,&x1^2&+&y1^2&+&z1^2&==&330^2,&0&&=&a&&=&Pi,&0&&=&b&&=&Pi,&0&&=&c&&=&Pi},&{a,&b,&c},&Reals]
采纳率：87%
来自团队：
lz的方程似乎有问题，x, y, z, x1, y1, z1为常量，数值化后也无解（如果有的话请列出一组）这种方程即使存在解析解，输出的式子也会非常巨大，可以占几个屏幕，又有什么意义呢？
如果是要解关于x,y,z的方程组,其中a,b,c,x1,y1,z1都是参数的话代码如下:(最好自行检查一下下面方程有没有抄错...我很容易抄错的...)(我只用过mathematica8和9,不确定以前版本是否可用)s = Solve[{Cos[b] Cos[c] x + Sin[a] Sin[b] Cos[c] y -
Cos[a] Sin[c] y + Cos[a] Sin[b] Cos[c] z + Sin[a] Sin[c] z == x1, Cos[b] Sin[c] x + Sin[a] Sin[b] Sin[c] y + Cos[a] Cos[c] y + Cos[a] Sin[b] Sin[c] z - Sin[a] Cos[c] z == y1, -Sin[b] x + Sin[a] Cos[b] y + Cos[a] Cos[b] z == z1}, {x, y, z}];X = FullSimplify[s[[1, 1, 2]]]Y = FullSimplify[s[[1, 2, 2]]]Z = FullSimplify[s[[1, 3, 2]]] 结果:X是-z1 Sin[b] + Cos[b] (x1 Cos[c] + y1 Sin[c])Y是Cos[b] Cos[c] (Cos[a] Sec[b] (y1 - x1 Tan[c]) + Sin[a] (z1 Sec[c] + Tan[b] (x1 + y1 Tan[c])))Z是Cos[a] (z1 Cos[b] + Sin[c] (y1 Sin[b] + x1 Tan[a]) + Cos[c] (x1 Sin[b] - y1 Tan[a])) 如果不是要解这个可以再描述的清楚一些.
我是想用x1,y1，z1表示a,b,c.
是常量，您看？
为您推荐：
其他类似问题
mathematica的相关知识
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。框出来的这个方程组是怎么解出来的啊？这上面只有解得，我不会算。化学_百度知道
框出来的这个方程组是怎么解出来的啊？这上面只有解得，我不会算。化学
我有更好的答案
14x:16y=7:1614x=7yy=2x14x+16y=9214x+16*2x=9246x=92x=2y=2x=4
为什么算分子质量氧要乘2
就是y为什么等于2X
采纳率：76%
来自团队：
为您推荐：
其他类似问题
方程组的相关知识
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
下载作业帮安装包
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
这个方程组具体怎么解啊?（1/x+1/y)×12=1 1/x×5+(1/x+1/y)×9=1我解了半天也没解出来,把得数也写出来,
作业帮用户
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
（1/x+1/y)×12=1 1/x×5+(1/x+1/y)×9=1 1/x+1/y=1/12.(1)(5+9)/x+9/y=1.(2)(2)-(1)*9:5/x=1/4x=20代入到（1）：1/20+1/y=1/121/y=1/30y=30所以：x=20,y=30
为您推荐：
其他类似问题
扫描下载二维码mathematica 如何解三角函数方程组_百度知道
mathematica 如何解三角函数方程组
如图，让上面两个式子等于0，然后求解t1、t2的值，我用solve求解时，一直显示“正在运行”，怎么才能快速解出啊，谢谢
啊啦，不知道怎么选择满意答案的同学，你又来问了啊，首先给你点建议，对于这种程序问题，你最好是直接给我们代码，你这样贴一张图，你传图传的累，我们把公式抄下来抄的更累。然后嘛，虽然我知道这次你肯定又不会把我选为满意（并且你很可能还在心里以为你已经把我选为满意答案了因为你已经点过赞同按钮了可你却没意识到那个和选为满意是不一样的……），但是你的这个问题有点意思所以我来答一下。首先，Solve应该是解不了这个问题的，这点你可以看下自带帮助里的更多信息，虽然Solve可以解的方程类型非常广泛，但它主攻的还是多项式方程，在其他方程上常有问题。对于这些方程，应该用Reduce。然后，让我们试试Reduce，直接输入你的两个方程，照样算了半天算不出来，这是怎么回事呢？Reduce也应付不了这个方程组吗？为了了解在这个方程组里发生了什么，让我们观察一下它的图形：ContourPlot[{94 (91 + 36 Cos[t1 - t2] - 75 Sin[t1]) - 5400 Sin[t2] ==
0, 94 (91 + 75 Cos[\[Pi]/10 + t1] - 36 Cos[\[Pi]/10 + t1 - t2]) -
5400 Cos[t2] == 0}, {t1, -5 \[Pi], 5 Pi}, {t2, -5 Pi, 5 Pi}]运行之后你将看到，这个方程组的解是周期性的，那么，很可能是它的解太过复杂，所以迟迟得不出结果，那么，让我们限制一下求解的定义域：Reduce[{94 (91 + 36 Cos[t1 - t2] - 75 Sin[t1]) - 5400 Sin[t2] == 0,
94 (91 + 75 Cos[\[Pi]/10 + t1] - 36 Cos[\[Pi]/10 + t1 - t2]) -
5400 Cos[t2] == 0, -Pi & t1 & Pi, -Pi & t2 & Pi}, {t1, t2}]求解后，产生了一个很大的输出，点击“显示全部”，你将看到这个解的长度和复杂度是相当惊人的，因此Reduce迟迟得不出结果也是不难理解的，可以猜测，只要有足够的时间并且你的电脑有足够的内存，Reduce应该可以得到这个方程的全部解。
采纳率：87%
来自团队：
为您推荐：
其他类似问题
您可能关注的内容
三角函数的相关知识
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。}