1个回答
展开全部
题主给出的积分方程,可以使用双循环语句和数值积分函数来求解T(x,t)。求解思路如下:
首先,根据x,t的范围值,将x,t值细分,即使用linspace函数,把x和t等分50(这个可以根据你的要求而改变)
x=linspace(0,6,50);t=linspace(0,5400,50);
其二,使用for语句,循环二次,即
for i=1:length(x1)
x=x1(i);
for j=1:length(t1)
t=t1(i);
。。。。。。
end
T(ⅰ,j)=T1
end
其三,自定义被积函数,即
func=@(fai)fai_s./(t-fai).^(3/2).*exp(-x^2./exp(a2.*(t-fai)));
其四,使用integral数值积分函数求T(x,t),即
T(x,t)=x/(8.91e-4*sqrt(pi))*integral(func,0,t);
其五,使用surf绘制其【x,t,T(x,t)】三维曲面图,即surf(x, t, T);
按上述思路编程,运行可以得到如下结果。
追问
。。那个能截一下屏吗
追答
如你需要代码,请私信
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询