一个MATLAB二重积分的问题要用MATLAB画一个表达式中含二重积分的曲线,但是不会写程序.看到您之前回答过一个类似的问题,很强大,希望能得到您的帮助.表达式如下,其中参数设定为a=0.05,b=0.1,t=1
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/05 14:50:34
![一个MATLAB二重积分的问题要用MATLAB画一个表达式中含二重积分的曲线,但是不会写程序.看到您之前回答过一个类似的问题,很强大,希望能得到您的帮助.表达式如下,其中参数设定为a=0.05,b=0.1,t=1](/uploads/image/z/5905628-44-8.jpg?t=%E4%B8%80%E4%B8%AAMATLAB%E4%BA%8C%E9%87%8D%E7%A7%AF%E5%88%86%E7%9A%84%E9%97%AE%E9%A2%98%E8%A6%81%E7%94%A8MATLAB%E7%94%BB%E4%B8%80%E4%B8%AA%E8%A1%A8%E8%BE%BE%E5%BC%8F%E4%B8%AD%E5%90%AB%E4%BA%8C%E9%87%8D%E7%A7%AF%E5%88%86%E7%9A%84%E6%9B%B2%E7%BA%BF%2C%E4%BD%86%E6%98%AF%E4%B8%8D%E4%BC%9A%E5%86%99%E7%A8%8B%E5%BA%8F.%E7%9C%8B%E5%88%B0%E6%82%A8%E4%B9%8B%E5%89%8D%E5%9B%9E%E7%AD%94%E8%BF%87%E4%B8%80%E4%B8%AA%E7%B1%BB%E4%BC%BC%E7%9A%84%E9%97%AE%E9%A2%98%2C%E5%BE%88%E5%BC%BA%E5%A4%A7%2C%E5%B8%8C%E6%9C%9B%E8%83%BD%E5%BE%97%E5%88%B0%E6%82%A8%E7%9A%84%E5%B8%AE%E5%8A%A9.%E8%A1%A8%E8%BE%BE%E5%BC%8F%E5%A6%82%E4%B8%8B%2C%E5%85%B6%E4%B8%AD%E5%8F%82%E6%95%B0%E8%AE%BE%E5%AE%9A%E4%B8%BAa%3D0.05%2Cb%3D0.1%2Ct%3D1)
一个MATLAB二重积分的问题要用MATLAB画一个表达式中含二重积分的曲线,但是不会写程序.看到您之前回答过一个类似的问题,很强大,希望能得到您的帮助.表达式如下,其中参数设定为a=0.05,b=0.1,t=1
一个MATLAB二重积分的问题
要用MATLAB画一个表达式中含二重积分的曲线,但是不会写程序.看到您之前回答过一个类似的问题,很强大,希望能得到您的帮助.
表达式如下,其中参数设定为a=0.05,b=0.1,t=1,x的定义域是(0.,正无穷).
一个MATLAB二重积分的问题要用MATLAB画一个表达式中含二重积分的曲线,但是不会写程序.看到您之前回答过一个类似的问题,很强大,希望能得到您的帮助.表达式如下,其中参数设定为a=0.05,b=0.1,t=1
最简单的方法是用数值积分,但是上下限一改差别非常大,说明你这个函数的振荡的非常厉害,不适合,具体的上下限你可以看看你的函数的收敛区间再定:
clc
clear
a=0.05;b=0.1;t=1;
xx=linspace(0,50);Q=zeros(size(xx));
for i=1:length(xx)
x=xx(i);
Q(i) = quad2d(@(y,v) sin(4*pi*v/b^2/t).*sinh(v).*exp(-2.*v.^2/b^2/t).*y.^(2*a/b^2).*exp( -2/(b^2*x).*(y.^2+2.*y.*cosh(v)+1) ),0,300,0,300,'AbsTol',1e-8);
end
Q=8*Q.*(pi*b^3*x.^2*sqrt(2*pi*t)).^(-1).*exp(4*pi^2-(a*t)^2/(2*b^2*t));
plot(xx,Q,'k.')
再就是用符号积分,不过这个会非常慢
syms y v
ap=0.05;bp=0.1;tp=1;
a=sym(0.05);b=sym(0.1);t=sym(1);
xx=linspace(0,50);Q=zeros(size(xx));
for i=1:length(xx)
x=sym(xx(i));
Q(i)=str2num(int(int(sin(4*pi*v/b^2/t)*sinh(v)*exp(-2.*v^2/b^2/t).*y^(2*a/b^2)*exp( -2/(b^2*x)*(y^2+2.*y*cosh(v)+1) ),v,0,Inf),y,0,Inf));
end
Q=8*Q.*(pi*bp^3*x.^2*sqrt(2*pi*t)).^(-1).*exp(4*pi^2-(ap*tp)^2/(2*bp^2*tp));
plot(xx,Q,'k.')
再有就是用quadlgk来积分,这个积分可以用无穷作为上下限,但是只能做一维,第二维感觉还是要用数值,使用你可以参考这个函数的帮助