matlab 7 重积分

cx330_123

已于 2024-04-18 19:45:26 修改

阅读量383

点赞数 5

文章标签： matlab 开发语言

于 2024-04-17 16:30:07 首次发布

本文链接：https://blog.csdn.net/cx330_123/article/details/137814534

版权

1. 求解积分 $\int_{-1}^{1}\int_{-1}^{1}e^{x^2+y^2}dxdy$

%求解顺序可变，此处为下到上的求和，改变x和y的for位置，可变为左到右
%数值解
h1=0.01;h2=0.01;s=0;
for y=-1:h1:1;
    for x=-1:h2:1;
        s=s+exp(x^2+y^2)*h2*h1;
    end
end
s


function z = myfun(x,y)
z=exp(x.^2+y.^2);
end
v=integral2(@myfun,-1,1,-1,1)

2. 曲线积分

第一曲线积分：求曲线质量

$M=\int_{L}^{}\rho (x,y)ds\\ds=\sqrt{(dx)^2+(dy^2)}$

s.t. 求解积分 $\int_{L}^{}\frac{Z^2}{x^2+y^2}ds$ ，其中L是螺旋线

syms t
x=cos(t);%t在[0:2*pi]间
y=sin(t);
z=t;
ds=sqrt(diff(x)^2+diff(y)^2+diff(z)^2);
simplify(ds)  %简化函数，最后ds为sqrt(2)*dt
%法一数值解
syms t%剩下求解t^2×ds的积分
f=sqrt(2)*t^2;
v1=int(f,0,2*pi)%结果类型为sym
%法二数值解
syms t
ff=@(t)sqrt(2).*t.^2;
v2=quad(ff,0,2*pi)%结果类型为num，编程序要函数句柄
%法三精确解
v3=integral(ff,0,2*pi)

第二曲线积分：物理意义来源于物体再给定力场中沿着曲线运动，计算所做的总功 (和方向有关） $W= \int_{L}^{}\vec{F}(x,y)d\vec{s}\\\vec{F}=P(x,y)\vec{i}+Q(x,y)\vec{j}\\d\vec{s}=dx\vec{i}+dy\vec{j}\\W=\int_{L}^{}[P(x,y)dx+Q(x,y)dy]$

3. 曲面积分：

第一曲面积分:对给定密度的函数（面密度）计算曲面质量

$M=\int \int_{s}^{}\rho (x,y,z)dS$

dS可以用切平面dA近似，而 $dA=\frac{1}{cos\gamma }dxdy$ ,cosr是切平面法向量n与z轴的方向余弦值

$M\\=\int \int_{s}^{}\rho (x,y,z)dS\\=\int \int_{D}^{}\rho (x,y,f(x,y))\sqrt{1+\frac{\partial z}{\partial x}+\frac{\partial z}{\partial y}}dxdy$

假设 $x=x(t),y=y(t),z=z(t)$ 在曲面s上，则 $F(x(t),y(t),z(t))=c$

两边求导得 ${F}'_{x}(t)\times{x}'(t)+{F}'_{y}(t)\times{y}'(t)+{F}'_{z}(t)\times{z}'(t)=0$

向量形式可表示为 $({F}'_{x},{F}'_{y},{F}'_{z})\cdot ({x}'(t),{y}'(t),{z}'(t))=0$ ，所以该F的偏导是切面法向量，进行单位化后，可得到余弦值。（ $\vec{n}=\sqrt{1+(\frac{\partial F}{\partial x})^2+(\frac{\partial F}{\partial y})^2}$ ）