MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.

来源:学生作业帮助网 编辑:作业帮 时间:2024/10/04 18:12:41
MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.
xTn@$_I\Q` ;hlS7vcZRJ%@xjˆ!}  J‹;=s;W.cw:j[%;k~|zF&7 ^U[sK3?L@f,<:~|Faz߭e.S )PCĔyCLXjU7&ri'GUF0Sn 4Э Ɇ@ ,G8 1B# hcSRd]ZQG@R:T YhQ{ZtB@Ґ&YۦuEu+**PڮNt`aTO;~}E#53 J_d 6݈X4.@#<ɼ,.Sf<|5;C\р{5 gSfdl ȍTť4yںN:!(CnSR*p6 K/)M_Pwu=\YhN< v䚤 F2GV^t7Z.;q.z(;1_Ic1&-d.+3ʂ>" U7#:/KEj5EiS;P 6z;Y~+0@/7Ww

MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.
MATLAB求一道定积分方程

已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.

MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.
clc;clear
format long
f=@(t) t./(exp(t)-1);
t=@(x) 1+4./x.*(1./x.*quad(f,0,x)-1)-0.65;
x1=1;
x2=10;
tol=1e-6;
while abs(x1-x2)>tol
    if feval(t,x1)*feval(t,(x1+x2)/2)>0
        x1=(x1+x2)/2;
    elseif feval(t,x1)*feval(t,(x1+x2)/2)<0
        x2=(x1+x2)/2;
    else
        break
    end
end
 x=(x1+x2)/2;
disp('精确到小数点后6位的解为')
disp(x)
t=feval(t,x)+0.65;
disp('将x带入式子结果为')
disp(t)
结果为
精确到小数点后6位的解为
   9.437618821859360

将x带入式子结果为
   0.649999994395076