24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

1275322766

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 614.5
红花: 1
帖子: 17
在线: 10小时
虫号: 3450917
注册: 2014-09-30
专业: 传动机械学

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

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

1.png

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有159人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

1楼 2016-03-20 10:16:53

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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])