24小时热门版块排行榜

返回列表

牛欢

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 26
帖子: 15
在线: 7.9小时
虫号: 5157169
注册: 2016-10-24
专业: 宏观经济学

[求助] 非线性方程求解已有1人参与

我对matlab没什么认识，各位大神能不能帮我写出具体的代码，关于方程(1-h)^0.461*(0.47880*h+0.32256)=0.28385的解，画一下图。谢谢

@月只蓝发自小木虫Android客户端

回复此楼

» 猜你喜欢

什么时候开奖？已经有13人回复
散金币祈福已经有94人回复
物理学I论文润色/翻译怎么收费? 已经有187人回复
青基已中已经有1人回复
散金币祈福已经有186人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有2人回复

1楼 2017-05-06 10:31:46

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

月只蓝

主管区长 (职业作家)

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

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

format long

[x fval]=fsolve(@(h) (1-h).^0.461.*(0.47880*h+0.32256)-0.28385,    [-0.5 0.5])

h=linspace(-1,1,500);

y=(1-h).^0.461.*(0.47880*h+0.32256)-0.28385;

figure(1)

plot(h,y,'b',[min(h) max(h)],[0 0],'r')

计算结果：
解：x =

-0.108293105549786 0.875531533029126

残差fval =

1.0e-06 *

-0.000000000055511 -0.326191971322665

说明一下，(1-h)^0.461，若h>1则没有实数解。

附图.png

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2017-05-06 12:14:47