24小时热门版块排行榜    

查看: 1720  |  回复: 5

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
回复此楼

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

有用

» 猜你喜欢

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

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

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的回帖

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的回帖

cjg0527

银虫 (初入文坛)

内容已删除
5楼2013-03-20 13:39:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cjg0527

银虫 (初入文坛)

内容已删除
6楼2013-03-20 13:40:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cjg0527 的主题更新
信息提示
请填处理意见