24小时热门版块排行榜

返回列表

oio14644

金虫 (正式写手)

应助: 227 (大学生)
金币: 1008.7
散金: 200
红花: 20
帖子: 674
在线: 59.6小时
虫号: 2581217
注册: 2013-08-05
性别: GG
专业: 环境工程

[求助] matlab 求解一个方程已有2人参与

请问这个方程怎么解？

我已经用matlab求过解析解了求不出来

Explicit solution could not be found.

求教数值的应该怎么做？

X 值域是 0.1到1 求Y

谢谢

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有300人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

MATLAB求解已经有2人回复
matlab 求解已经有1人回复

1楼 2014-07-27 11:54:37

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +1

只能求数值解吧！

还有把你的问题用文本格式给出。你给个图片到很方便，想帮的人还得一个个输入。。。?

赞一下

回复此楼

2楼2014-07-27 14:27:54

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

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 3430
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
oio14644: 金币+50, ★有帮助 2014-07-28 01:30:41

Parameters x[0.1, 1], y;
Function x^2*10^(-2.0588*(x+3*y)^0.5/(1+1.9758*(x+3*y)^0.5))*10^(0.5147*(x+3*y)^0.5/(1+4.1393*(x+3*y)^0.5))*10^(0.5147*(x+3*y)^0.5/(1+4.1393*(x+3*y)^0.5))=3.5/y/10^10;
--------------------------------------------------------
不加x限制，结果有无穷多组；

目标函数值: 0
x: 0.000115458457023581
y: 0.0528939761001652

目标函数值: 0
x: 4.35789128253324E-5
y: 0.660158226930854

目标函数值: 0
x: -8.05412292226326E-5
y: 0.135748411306841

目标函数值: 0
x: 0.000177999439686659
y: 0.0177058509332421

目标函数值: 0
x: 0.000133593794530008
y: 0.0363409793326377

目标函数值: 0
x: 9.50514568389083E-5
y: 0.0878619210580436

---------------------------------------------
加了以后，近似无解。
目标函数值: 0.00160946334974056
x: 0.10000000000019
y: 1.74522020002736E21

赞一下

回复此楼

3楼2014-07-27 14:53:53

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

月只蓝

主管区长 (职业作家)

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

CODE:

function sol_nonl_eq

clear all;clc

global x

format long

y0=0.1;

xspan=0.1:0.1:1;

n=length(xspan);

for i=1:n

    x=xspan(i);

    y(i)=fsolve(@fun,y0);

    y_fval(i)=fun(y(i));

end

[xspan' y' y_fval']

plot(xspan,y,'b-o')

function f=fun(y)

global x

f=y-3.5e-10/(x^2)/(10^(-2.0588*(x+3*y)^0.5/(1+1.9758*(x+3*y)^0.5)))/(10^(-0.5147*(x+3*y)^0.5/(1+4.1393*(x+3*y)^0.5)))^2;

计算结果：
x y 目标函数值
0.100000000000000 0.000000147559815 0.000000025731971
0.200000000000000 0.000000043016195 0.000000003906666
0.300000000000000 0.000000021303785 0.000000001168052
0.400000000000000 0.000000013069400 0.000000000517437
0.500000000000000 0.000000008867610 0.000000000181978
0.600000000000000 0.000000006591207 0.000000000170957
0.700000000000000 0.000000005106442 0.000000000139441
0.800000000000000 0.000000004021441 0.000000000048214
0.900000000000000 0.000000003285544 0.000000000024982
1.000000000000000 0.000000002678200 -0.000000000052131