应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » MATLAB求解非线性方程组
51/1返回列表
查看: 1718  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cjg0527

银虫 (初入文坛)

[求助] MATLAB求解非线性方程组

我想求解如图所示的二元一次非线性方程组，t1与t0是未知数，其他均为常数。用如下的方法求解，系统提示不能求解，我以前从来没有接触过MATLAB，不知道用什么方法才能求解，所有的工作的都卡着这个地方了，非常着急，也请大家多多帮助。
Rw0=12e-3
Vc=10e-6/60
Rg=3e-3
a=6e-3
ae=10e-6/60
nw=300*3.14*2/60
f1='((Rw0+Vc*t1)*cos(2*3.14*nw*t1)-a-Vc*t1+(a+Vc*t0)*cos(2*3.14*nw*(t1-t0)))^2+((Rw0+Vc*t1)*sin(2*3.14*nw*t1)-(a+Vc*t0)*sin(2*3.14*nw*(t1-t0)))^2-Rg^2=0'
f2='Rw0+Vc*t0-ae-(Rw0+Vc*t1)*cos(2*3.14*nw*(t1-t0))=0'
[t1,t0]=solve(f1,f2,'t1,t0')

1.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

有用

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-03-19 13:17:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shizijun

木虫 (著名写手)

自定义头衔

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2013-03-20 20:54:26
改写了一点，可以算。不过貌似跟你给出的f2的表达式不太一致

clear all;clc;
syms t1 t0
Rw0=12e-3;
Vc=10e-6/60;
Rg=3e-3;
a=6e-3;
ae=10e-6/60;
nw=300*3.14*2/60;
f1=((Rw0+Vc*t1)*cos(2*3.14*nw*t1)-a-Vc*t1+(a+Vc*t0)*cos(2*3.14*nw*(t1-t0)))^2+((Rw0+Vc*t1)*sin(2*3.14*nw*t1)-(a+Vc*t0)*sin(2*3.14*nw*(t1-t0)))^2-Rg^2;
f2=Rw0+Vc*t0-ae-(Rw0+Vc*t1)*cos(2*3.14*nw*(t1-t0));
[t0,t1]=solve(f1,f2,'t1,t0')
赞一下 回复此楼
个人签名
4楼2013-03-19 18:55:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

xxyyxz

银虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2013-03-20 20:54:10
第二个方程f2何来+Vc*t0？Rw（t0）不等于Rw0,而且没有发现gm和Rw（t1）的表达式，请补充
赞一下 回复此楼
2楼2013-03-19 18:17:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mechzhen

铜虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2013-03-20 20:54:19
那是因为f1='((Rw0+Vc*t1)*cos(2*3.14*nw*t1)-a-Vc*t1+(a+Vc*t0)*cos(2*3.14*nw*(t1-t0)))^2+((Rw0+Vc*t1)*sin(2*3.14*nw*t1)-(a+Vc*t0)*sin(2*3.14*nw*(t1-t0)))^2-Rg^2=0' ，这样写得话系统是不认Rw0等参数的吧？！
赞一下 回复此楼
奔跑-追赶-超越-在路上！
3楼2013-03-19 18:41:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭