| 查看: 2462 | 回复: 1 | ||
physicsapple铜虫 (小有名气)
|
[求助]
疲劳sn曲线的三参数法方程求解，二分法
|
|
各位高手，我是matlab小白，拼凑了一段matlab程序来求一个方程的根。是用二分法求H（S0）=0的根，但是一直运行不出来。 代码如下： syms S0;%定义变量S0 n=5; S1=[100 80 60 40 30 ];%应力 y=log10([1534 53187 153854 754156]); x=log10(S1-S0); A=x./(S1-S0); B=y./(S1-S0); C=1./(S1-S0); Lx0=sum(A)-(1/n)*sum(x)*sum(C); Ly0=sum(B)-(1/n)*sum(y)*sum(C); Lxx=sum(x.*x)-(1/n)*sum(x)*sum(x); Lyy=sum(y.*y)-(1/n)*sum(y)*sum(y);%L相关式子都是S0含有S0的 H=Lx0/Lxx-Ly0/Lyy; f=inline('');%用二分法求H=0的根，这一步我不太懂，不知道括号里面放什么 a=0; b=100; dlt=0.5*1.0e-5; k=1; while abs(b-a)>dlt c=(a+b)/2; if f(c)==0 break; elseif f(c)*f(b)<0 a=c; else b=c; end fprintf('k=%d,S0=%.5f\n',k,c); k=k+1; end |
回复此楼» 猜你喜欢
求推荐博导
已经有4人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
为什么nbs上溴 没有产物点出现呢
已经有5人回复
-
存款400万可以在学校里躺平吗
已经有34人回复
-
求助一下有机合成大神
已经有4人回复
-
最失望的一年
已经有16人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
疑惑？
已经有5人回复
physicsapple
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 760.8
- 散金: 1
- 帖子: 61
- 在线: 39.5小时
- 虫号: 1713096
- 注册: 2012-03-24
- 性别: GG
- 专业: 力学
2楼2016-03-12 21:13:54