| 查看: 2481 | 回复: 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 |
回复此楼» 猜你喜欢
专硕085601求调剂
已经有6人回复
-
282求调剂 不挑专业 求收留
已经有6人回复
-
材料专硕 调剂
已经有7人回复
-
285求调剂
已经有6人回复
-
280求调剂
已经有17人回复
-
工科 267求调剂
已经有6人回复
-
调剂
已经有5人回复
-
材料调剂
已经有13人回复
-
315求调剂
已经有3人回复
-
270求调剂
已经有4人回复
physicsapple
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 760.8
- 散金: 1
- 帖子: 61
- 在线: 39.5小时
- 虫号: 1713096
- 注册: 2012-03-24
- 性别: GG
- 专业: 力学
2楼2016-03-12 21:13:54
