24小时热门版块排行榜    

查看: 1250  |  回复: 7

1275322766

新虫 (初入文坛)

[求助] matlab解方程 已有1人参与

此方程怎么在matlab编程求解,我用solve就得一个负值,好像是错的!求大神帮忙写个程序!本人不善用matlab!

matlab解方程
1.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1275322766

新虫 (初入文坛)

方程中D为150,d为45
2楼2016-03-20 10:18:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuansen

木虫 (著名写手)

明天有时间帮你看看,没道理解不出来啊

发自小木虫Android客户端
3楼2016-03-21 00:52:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuansen

木虫 (著名写手)

晕,还没时间看了,又要拖到明天了!

发自小木虫Android客户端
4楼2016-03-24 00:12:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
1275322766: 金币+10, ★★★很有帮助 2016-03-25 09:34:54
CODE:
function solve_n_eqs
clear all;clc
d=45;
D=150;
n=5000;   
x=linspace(-0.0005,0.001,n);
y=0.1248*x.^4.*( 1./(pi*d-0.2*x)-1./(pi*D-0.2*x)   )+0.2496*x.^3.*log(abs( (pi*D-0.2*x)./(pi*d-0.2*x)  ))+1.248*pi*x.^2*(D-d)-0.2*x.^3*(D^3-d^3)    ;
figure(1)
plot(x,y,'b',[min(x) max(x)],[0 0],'r'),axis([min(x) max(x) min(y) max(y)]);

format long

[x fval]=fsolve(@(x) 0.1248*x.^4.*( 1./(pi*d-0.2*x)-1./(pi*D-0.2*x)   )+0.2496*x.^3.*log(abs( (pi*D-0.2*x)./(pi*d-0.2*x)  ))+1.248*pi*x.^2*(D-d)-0.2*x.^3*(D^3-d^3),[1e-5 1e-3])

x =

  1.0e-003 *

   0.001238815125239   0.632354921061680


残差fval =

  1.0e-005 *

   0.000063053264693  -0.145570151313058
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
5楼2016-03-24 13:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1275322766

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by 月只蓝 at 2016-03-24 13:09:52
function solve_n_eqs
clear all;clc
d=45;
D=150;
n=5000;   
x=linspace(-0.0005,0.001,n);
y=0.1248*x.^4.*( 1./(pi*d-0.2*x)-1./(pi*D-0.2*x)   )+0.2496*x.^3.*log(abs( (pi*D-0.2*x)./(pi*d-0.2*x)   ...

结果不是这样的,应该在该在120~150之间吧,这和我想要的结果相差甚远!不过还是非常感谢你的!
6楼2016-03-25 09:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
6楼: Originally posted by 1275322766 at 2016-03-25 09:35:35
结果不是这样的,应该在该在120~150之间吧,这和我想要的结果相差甚远!不过还是非常感谢你的!...

抱歉,5楼的方程有个地方写错了。现检查后,重新的计算结果:
x =

  1.0e+002 *

  -2.328129109456768   1.979382267670229


fval =

  1.0e-008 *

  -0.372529029846191   0.372529029846191
x有两个根,分别为-233和198左右,也不在120~150之间。
matlab解方程-1
附图1.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
7楼2016-03-25 12:28:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuansen

木虫 (著名写手)

专家有图有码,不得不服!!

发自小木虫Android客户端
8楼2016-04-03 11:21:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1275322766 的主题更新
信息提示
请填处理意见