求助：编写程序求解积分方程中的参数k的值

方程中已知rw/c=0.38;tITZ=40;Dc=0.488;入=1.08;x范围0-tITZ；
怎么编程求参数k的值

[Dyuman Joshi]

Please show what you have tried yet.

[一凡 李]

[Dyuman Joshi]

I meant, please show the code you have written to try to solve the problem.

[一凡 李]

clear
clc
a0=0.5;
r0=1.08;
w_c=0.38;
faic=1-(1+1.31*a0)/(1+3.2*w_c);
t_ITZ=40;
D_cem=100;
D_c=((1.0482*10^-5*D_cem^2+3.246*10^-4*D_cem+0.0146)*w_c-1.79*10^-7*D_cem^2+5.0429*10^-5*D_cem+1.00564)/(1+3.15*w_c);
y=125/(1+3.15*w_c)-(125-t_ITZ)*D_c;
fzero(@(k) integral(@(x) D_c*(x/t_ITZ).^(1-1.08*(x/t_ITZ).^k),0,t_ITZ)-y,0)

Answer 1

[Torsten]

gamma_wc = 0.38;
t_ITZ = 40;
D_c = 0.488;
lambda = 1.08;
fun = @(k) 125/(1+3.15*gamma_wc)-(125-t_ITZ)*D_c-integral(@(x)D_c*(x/t_ITZ).^(1-lambda*(x/t_ITZ).^k),0,t_ITZ);
k = fsolve(fun,1)

Equation solved.

fsolve completed because the vector of function values is near zero
as measured by the value of the function tolerance, and
the problem appears regular as measured by the gradient.

k = 0.3083

[Dyuman Joshi]

I tried fzero() first, and boy, that was a ride. Then I moved on to fsolve() and vpasolve(), both of which worked.

Though, I think you should have provided some hints/guidelines first instead of posting the full solution.

[一凡 李]

thank U so much