24小时热门版块排行榜    

查看: 291  |  回复: 3

晓龙追梦

金虫 (小有名气)

[求助] 方程组求解

S=solve('1/(w2*w1)^2-A=0','2/w2^2*k1/w1+2*k2/w2/w1^2-B=0','1/w2^2+4*k2/w2*k1/w1+1/w1^2-C=0','2*k2/w2+2*k1/w1-D=0','k1','k2','w1','w2')
disp('S.k1'),disp(S.k1),disp('S.k2'),disp(S.k2),disp('S.w1'),disp(S.w1),disp('S.w2'),disp(S.w2)
计算结果说没有解析解 怎么办呢?  我的本意是通过ABCD表示四个未知数k1k2w1w2
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jswoo

铁杆木虫 (著名写手)

阿拉多梦

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai(金币+1): 谢谢应助! 2012-03-02 22:07:18
这是符号计算,
大概这个模式:
举例:
f1='1/(w2*w1)^2-A';
solve(f1,'A')
ans =
1/w2^2/w1^2

你试试
2楼2012-03-02 14:23:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dk1013

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai(金币+1): 谢谢参与! 2012-03-02 22:07:30
符号计算还是用mathematica吧?
如果还是无解析解,看来只能数值了。
到得还来别无事 庐山烟雨浙江潮
3楼2012-03-02 17:28:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dk1013

木虫 (正式写手)

【答案】应助回帖

★ ★
xiegangmai(金币+2): 谢谢应助! 2012-03-02 22:07:43
Mathematica 应该有解的
Solve[{1/(w2*w1)^2 - A == 0, 2/w2^2*k1/w1 + 2*k2/w2/w1^2 - B == 0,
  1/w2^2 + 4*k2/w2*k1/w1 + 1/w1^2 - C == 0,
  2*k2/w2 + 2*k1/w1 - D == 0}, {k1, k2, w1, w2}]
只不过这个解太长了,就不在这里贴了
希望帮到你
到得还来别无事 庐山烟雨浙江潮
4楼2012-03-02 17:34:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 晓龙追梦 的主题更新
信息提示
请填处理意见