| 查看: 1204 | 回复: 10 | |||
tmw1987铁杆木虫 (正式写手)
|
[交流]
【求助】函数求零点！！！苦恼挣扎求解中…… 已有4人参与
|
|
大侠： 小弟我遇到如下图所示的周期性函数，需要求解前四十个零点值对应的x值。我用fzero函数求解，但是只能一个个求解，感觉如此强大的matlab肯定有更加高明的方法可以求解。 函数和图形都已知，如何一次性求解多个零点值！ 请赐教！！！  |
回复此楼» 猜你喜欢
有谁可曾问过你过的还好吗？
已经有21人回复
-
E0414, 我的本子有没有希望？
已经有7人回复
-
一篇论文同时出现在两个期刊，一模一样，这算不算学术不端，请各位老师斧正。
已经有12人回复
-
希望面上有个好结果
已经有7人回复
-
三区计算机方向期刊推荐
已经有5人回复
-
sci论文二审求助
已经有5人回复
-
函评
已经有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
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 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