声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 682|回复: 0

[编程技巧] 帮我看下用matlab求解非线性方程时哪里出了问题

[复制链接]
发表于 2009-1-8 20:35 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
程序如下:其中s1,s2,s3,s4,s5,s6,s7,s8是方程;k2 y1 y2 m1 m2 x1 x2为相应的7个变量,但是求解时,有如下显示,没求出解,需要高手指导,谢谢!!!
syms k2 y1 y2 m1 m2 x1 x2;
>> s1=k2+8.5847142026029796983461767280000;
>> s3=y1-(-2.-2.040*k2)/(-.2996e-1-.200e-3*k2);
>> s4=y2-(1+.20e-1*k2)/(-.2996e-1-.200e-3*k2);
s5=m1-(-.6e-1*(-2+.10e-1*y1)*(1+.10e-1*y2)+1.960+.1920e-1*y2-.200e-3*y1+.40e-3*x1*(1+.10e-1*y2))/(.6e-1*(1+.10e-1*y2)^2+.40e-1+.400e-3*y2);
s6=m2-((-.6e-1-.120e-2*k2)*(1+.10e-1*y2)-.20e-1-.400e-3*k2+.80e-3*x2*(1+.10e-1*y2))/(.6e-1*(1+.10e-1*y2)^2+.40e-1+.400e-3*y2);
s7=x1-2.5*(-1020000.+6.*k2*m1*y2+3.*m2*y2*y1+800.*m2*m1*y2+3.*m2*y2^2*m1-20400.*k2+40000.*m1-1040000.*m2+300.*y1+40800.*x2-10200.*m2*y2+300.*m1*y2+800.*k2*m1+6.*k2*y1+50000.*m2*m1+400.*m2*y1-400.*x2*m1-4.*x2*y1)/(5.*m2*y2+1498.+10.*k2+500.*m2);
s8=x2-.625*(30000.+1200.*k2+80000.*m2+600.*m2*y2+12.*k2^2+1600.*k2*m2+12.*k2*m2*y2+50000.*m2^2+800.*m2^2*y2+3.*m2^2*y2^2)/(1499.+10.*k2+500.*m2+5.*m2*y2);
>> solve(s1,s3,s4,s5,s6,s7,s8,' k2 y1 y2 m1 m2 x1 x2')

ans =
    k2: [2x1 sym]
    m1: [2x1 sym]
    m2: [2x1 sym]
    x1: [2x1 sym]
    x2: [2x1 sym]
    y1: [2x1 sym]
    y2: [2x1 sym]
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-12 01:11 , Processed in 0.062425 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表