| 查看: 1223 | 回复: 10 | |||
tmw1987铁杆木虫 (正式写手)
|
[交流]
【求助】函数求零点！！！苦恼挣扎求解中…… 已有4人参与
|
|
大侠： 小弟我遇到如下图所示的周期性函数，需要求解前四十个零点值对应的x值。我用fzero函数求解，但是只能一个个求解，感觉如此强大的matlab肯定有更加高明的方法可以求解。 函数和图形都已知，如何一次性求解多个零点值！ 请赐教！！！  |
回复此楼» 猜你喜欢
生命口会评
已经有6人回复
-
27届辽宁大学应届毕业生申博
已经有3人回复
-
关于如何从代码看上不上会
已经有23人回复
-
大龄残疾硕士的一点执念
已经有25人回复
-
今年E04面上
已经有22人回复
-
chemdraw
已经有6人回复
-
祈祷青基必中
已经有10人回复
-
b口会评
已经有6人回复
-
信息学部会评时间
已经有3人回复
-
在职考研
已经有7人回复
tmw1987
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4627.6
- 散金: 1526
- 红花: 9
- 帖子: 982
- 在线: 540.4小时
- 虫号: 815020
- 注册: 2009-07-24
- 性别: GG
- 专业: 功能与智能高分子
2楼2010-06-16 18:12:37
3楼2010-06-16 18:49:46
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
4楼2010-06-16 19:30:06
tmw1987
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4627.6
- 散金: 1526
- 红花: 9
- 帖子: 982
- 在线: 540.4小时
- 虫号: 815020
- 注册: 2009-07-24
- 性别: GG
- 专业: 功能与智能高分子
5楼2010-06-16 19:37:09
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
6楼2010-06-16 20:14:45
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
用全是1模拟 图像太那个...... 7楼2010-06-16 21:44:39
tmw1987
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4627.6
- 散金: 1526
- 红花: 9
- 帖子: 982
- 在线: 540.4小时
- 虫号: 815020
- 注册: 2009-07-24
- 性别: GG
- 专业: 功能与智能高分子
8楼2010-06-16 22:13:41
tmw1987
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4627.6
- 散金: 1526
- 红花: 9
- 帖子: 982
- 在线: 540.4小时
- 虫号: 815020
- 注册: 2009-07-24
- 性别: GG
- 专业: 功能与智能高分子
9楼2010-06-16 22:34:23
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
★ ★
robert2020(金币+2):辛苦了~~ 2010-06-17 13:43:41
robert2020(金币+2):辛苦了~~ 2010-06-17 13:43:41
|
一时也没什么好方法 还是原始那种 在一定的精度里寻找 >> myfun=inline('(1-0.2040)*sin((0.022-0.02)*x/sqrt(1.6e-7)-x*0.02/sqrt(1e-7))-(0.2040+0.02)*sin((0.022-0.02)*x/sqrt(1.6e-7)-x*0.02/sqrt(1e-7)) ','x') myfun = Inline function: myfun(x) = (1-0.2040)*sin((0.022-0.02)*x/sqrt(1.6e-7)-x*0.02/sqrt(1e-7))-(0.2040+0.02)*sin((0.022-0.02)*x/sqrt(1.6e-7)-x*0.02/sqrt(1e-7)) >> [x,Fval,exitflag]=fsolve(myfun,[0:0.05:1]') Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. x = 0 0.053937 0.10787 0.16181 0.21575 0.26969 0.37756 -0.16181 0.32362 0.4315 0.48543 0.53937 0.59331 0.64724 0.70118 0.75512 0.80906 0.86299 0.91693 0.80906 1.1327 Fval = 0 -7.005e-017 1.401e-016 -2.1015e-016 2.802e-016 1.5072e-009 -2.5225e-015 2.1015e-016 4.4846e-015 4.6247e-015 -6.3045e-016 7.005e-016 1.2616e-015 8.406e-016 -2.9428e-015 9.807e-016 9.8141e-016 1.1208e-015 -3.223e-015 5.0457e-015 5.6111e-016 exitflag = 1 |
10楼2010-06-16 23:05:01