| 查看: 2563 | 回复: 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 |
回复此楼» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有14人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组-博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
physicsapple
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 760.8
- 散金: 1
- 帖子: 61
- 在线: 39.5小时
- 虫号: 1713096
- 注册: 2012-03-24
- 性别: GG
- 专业: 力学
2楼2016-03-12 21:13:54