| 查看: 2480 | 回复: 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人回复
-
070300化学求调剂
已经有11人回复
-
抱歉
已经有5人回复
-
环境285分，过六级，求调剂
已经有9人回复
-
求材料调剂，一志愿郑州大学289分
已经有11人回复
-
化学调剂
已经有9人回复
-
303求调剂
已经有4人回复
-
化工调剂303分，过四级
已经有25人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
298分 070300求调剂
已经有11人回复
physicsapple
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 760.8
- 散金: 1
- 帖子: 61
- 在线: 39.5小时
- 虫号: 1713096
- 注册: 2012-03-24
- 性别: GG
- 专业: 力学
2楼2016-03-12 21:13:54
