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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab解方程
81/1返回列表
查看: 1316  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298-一志愿中国农业大学-求调剂 +11 手机用户 2026-03-17 12/600 2026-03-23 23:51 by 热情沙漠
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 070300化学求调剂 +8 苑豆豆 2026-03-20 8/400 2026-03-23 20:57 by baobaoye
[考研] 一志愿中国石油大学（华东） 本科齐鲁工业大学 +4 石能伟 2026-03-17 4/200 2026-03-23 17:51 by 17862566385
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 08工学调剂 +7 用户573181 2026-03-20 11/550 2026-03-23 15:47 by 我爱学习学习使�
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 070300，一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 305分求调剂（食品工程） +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 初始318分求调剂（有工作经验） +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
信息提示
关闭
请填处理意见
关闭