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

Question

0 个投票

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

4 个评论
显示 2更早的评论隐藏 2更早的评论

Dyuman Joshi 2024-1-24

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

一凡李 2024-1-24

在 MATLAB Online 中打开

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)

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question

Answer 1

Torsten 2024-1-24

移动：Torsten 2024-1-24

在 MATLAB Online 中打开

1 个投票

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

2 个评论
显示无隐藏无

Dyuman Joshi 2024-1-24

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.

一凡李 2024-1-24

thank U so much

请先登录，再进行评论。

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

4 个评论
显示 2更早的评论隐藏 2更早的评论

回答（1 个）

2 个评论
显示无隐藏无

类别

产品

版本

标签

Community Treasure Hunt

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

4 个评论 显示 2更早的评论 隐藏 2更早的评论

回答（1 个）

2 个评论 显示 无 隐藏 无

类别

产品

版本

标签

另请参阅

Community Treasure Hunt

4 个评论
显示 2更早的评论隐藏 2更早的评论

2 个评论
显示无隐藏无