Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2010-09-17, 15:43   #1
xiaoxiong5227
初级会员
 
注册日期: 2010-09-17
年龄: 37
帖子: 1
声望力: 0
xiaoxiong5227 正向着好的方向发展
默认 solve解非线性方程组

代码如下:
clear;clc
syms x1 x2 x3 x4 x5 x6 t2
y1=-221/14+6/35*x1+9/70*x3+3*x5*(1-t2)^2*t2;
y2=-69/14+6/35*x2+9/70*x4+3*x6*(1-t2)^2*t2;
y3=-233/14+9/70*x1+6/35*x3+3*x5*(1-t2)*t2^2;
y4=-7/2+9/70*x2+6/35*x4+3*x6*(1-t2)*t2^2;
y5=30*(1-t2)^3+3*(1-t2)^2*t2*x1+3*(1-t2)*t2^2*x3+90*t2^3-6233770349685199/140737488355328;
y6=-10*(1-t2)^3+3*(1-t2)^2*t2*x2+3*(1-t2)*t2^2*x4-30*t2^3-2480330921882157/281474976710656;
d=solve(y1,y2,y3,y4,y5,y6);
x1=vpa(d.x1,6)
x2=vpa(d.x2,6)
x3=vpa(d.x3,6)
x4=vpa(d.x4,6)
x5=vpa(d.x5,6)
x6=vpa(d.x6,6)
f=30*(1-t2)^3+3*(1-t2)^2*t2*x1+3*(1-t2)*t2^2*x3+90*t2^3-6233770349685199/140737488355328;
f=vpa(f,6)
t2=0:0.0001:1
plot(t2,subs(f))



为什么画出来的图形不是0 啊,是不是非线性方程组不能这样解啊?
xiaoxiong5227 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 10:09


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.