24小时热门版块排行榜

返回列表

physicsapple

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 760.8
散金: 1
帖子: 61
在线: 39.5小时
虫号: 1713096
注册: 2012-03-24
性别: GG
专业: 力学

[求助] 疲劳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

回复此楼