24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1485  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Kevin1987s

新虫 (初入文坛)

[求助] matlab求解非线性方程组,求助!

θ1-4未知,用已知表示未知,怎么编程呢?
这是我编的程序,方程中用l代替l1,m代替l2
>> syms \alpha x y z R r l m \theta_1 \theta_2 \theta_3 \theta_4
eq1=[2^1/2*r*sin(\pi/4-\alpha)*cos\alpha+x-R-l*sin\theta_1]^2+[2^1/2*r*sin(\pi/4-\alpha)*sin\alpha+y]^2+(z+l*cos\theta_1)^2-m^2;
eq2=[-2^1/2*r*cos(\pi/4-\alpha)*sin\alpha+x]^2+[2^1/2*r*cos(\pi/4-\alpha)*cos\alpha+y-R-l*sin\theta_2]^2+(z+l*cos\theta_2)^2-m^2;
eq3=[-2^1/2*r*sin(\pi/4-\alpha)*cos\alpha+x+R+l*sin\theta_3]^2+[-2^1/2*r*sin(\pi/4-\alpha)*sin\alpha+y]^2+(z+l*cos\theta_3)^2-m^2;
eq4=[2^1/2*r*cos(\pi/4-\alpha)*sin\alpha+x]^2+[-2^1/2*r*cos(\pi/4-\alpha)*cos\alpha+y+R+l*sin\theta_4]^2+(z+l*cos\theta_4)^2-m^2;
d=solve(eq1,eq2,eq3,eq4)
d.\theta_1 \theta_2 \theta_3 \theta_4
??? Error using ==> syms at 61
Not a valid variable name.
应该怎么改啊?十分感谢!

方程组
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

Kevin1987s

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by csgt0 at 2012-05-21 13:29:15:
变量名不对吧,别用\应该可以

可能是方程太复杂了,一直Busy。。。。有优化的算法么?谢谢您!
3楼2012-05-21 14:15:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

感谢参与,应助指数 +1
变量名不对吧,别用\应该可以
showmethemoney
2楼2012-05-21 13:29:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

我也不太懂,刚用你的2,4式试了下,能解出x,y,但是也是无穷串
showmethemoney
4楼2012-05-21 14:20:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Kevin1987s

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by csgt0 at 2012-05-21 14:20:32:
我也不太懂,刚用你的2,4式试了下,能解出x,y,但是也是无穷串

限制有效数字呢?在倒数第二行加上vpa(theta_1,4),vpa(theta_2,4),vpa(theta_3,4),vpa(theta_4,4)
5楼2012-05-21 14:28:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见