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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab解方程
81/1返回列表
查看: 1321  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

1275322766

新虫 (初入文坛)

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

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

matlab解方程
1.png
回复此楼

» 猜你喜欢
1楼 2016-03-20 10:16:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1275322766

新虫 (初入文坛)
方程中D为150，d为45
赞一下 回复此楼
2楼2016-03-20 10:18:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuansen

木虫 (著名写手)
明天有时间帮你看看，没道理解不出来啊

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

yuansen

木虫 (著名写手)
晕，还没时间看了，又要拖到明天了！

发自小木虫Android客户端
赞一下(1人) 回复此楼
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
赞一下(2人) 回复此楼
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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
[考研] 300求调剂，材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 材料学硕，求调剂 6+3 糖葫芦888ll 2026-03-22 6/300 2026-03-24 15:05 by 123456zjm
[考研] 293求调剂 +6 加一一九 2026-03-24 6/300 2026-03-24 14:29 by JourneyLucky
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 求材料，环境专业调剂 +3 18567500178 2026-03-18 3/150 2026-03-23 23:50 by 热情沙漠
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 328求调剂，英语六级551，有科研经历 +7 生物工程调剂 2026-03-17 12/600 2026-03-23 18:18 by YMU施老师
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境） +3 我爱学电池 2026-03-23 3/150 2026-03-23 17:16 by AZMK
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123！
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
信息提示
关闭
请填处理意见
关闭