matlab求一个二重积分,最终给出一个数值就行,

来源：学生作业帮编辑：百度作业网作业帮分类：数学作业时间：2024/07/31 19:52:42

matlab求一个二重积分,最终给出一个数值就行,
syms h z ;
gam=8;
u=0.33;
a=50;
b=50;
z1=z-h;
z2=z+h;
R1=(a^2+b^2+z1^2)^0.5;
R2=(a^2+b^2+z2^2)^0.5;
r1=(a^2+z1^2)^0.5;
r2=(a^2+z2^2)^0.5;
r3=(b^2+z1^2)^0.5;
r4=(b^2+z2^2)^0.5; %定义积分符号
m1=(gam/pi)*atan((a*b)/(z1*R1)); %第一部分
m2=(gam/pi)*atan((a*b)/(z2*R2)); %第二部分
m3=(gam*a*b*z1*(R1^2+z1^2))/(pi*(1-u)*2*r1^2*r3^2*R1); %第三部分
m4=(gam*a*b*(h+(3-4*u)*z)*(R2^2+z2^2))/(pi*(1-u)*2*r2^2*r4^2*R2); %第四部分
m5=(a*b*h*z*z2*(2*R2^2+r4^2))/(R2^3*r4^4); %第五部分
m6=(a*b*h*z*z2*(2*R2^2+r2^2))/(R2^3*r2^4); %第六部分
f=@(z,h)m1+m2+m3+m4+m5+m6; %被积分函数
F=dblquad(@(z,h)f,0,10),h),10,50); %求积分,先对h积分,后对z积分

clc;
clear all;
syms h z;
gam=8;
u=0.33;
a=50;
b=50;
z1=z-h;
z2=z+h;
R1=(a^2+b^2+z1^2)^0.5;
R2=(a^2+b^2+z2^2)^0.5;
r1=(a^2+z1^2)^0.5;
r2=(a^2+z2^2)^0.5;
r3=(b^2+z1^2)^0.5;
r4=(b^2+z2^2)^0.5; %定义积分符号
m1=(gam/pi)*atan((a*b)/(z1*R1)); %第一部分
m2=(gam/pi)*atan((a*b)/(z2*R2)); %第二部分
m3=(gam*a*b*z1*(R1^2+z1^2))/(pi*(1-u)*2*r1^2*r3^2*R1); %第三部分
m4=(gam*a*b*(h+(3-4*u)*z)*(R2^2+z2^2))/(pi*(1-u)*2*r2^2*r4^2*R2); %第四部分
m5=(a*b*h*z*z2*(2*R2^2+r4^2))/(R2^3*r4^4); %第五部分
m6=(a*b*h*z*z2*(2*R2^2+r2^2))/(R2^3*r2^4); %第六部分
f=m1+m2+m3+m4+m5+m6; %被积分函数的符号变量
ff=eval(['@(h,z)',vectorize(char(f)),';']);%转为函数
F=dblquad(ff,0,10,10,50) %求积分,先对h积分,后对z积分
运行之后
F =
2.7928e+003
也就是2792.8
再问：谢谢大神啊，这个积分能不能得到解释解呢？就是对h积分一次以后，能不能得到不含h却含有z的被积表达式？积分两次后是否能得到解释式呢
再答：由于你的式子非常复杂，你说的对h积分是无法得到不含h却含有z的被积表达式，这个我已经试过了。积分两次也不能，得到的结果都是同时含有h和z的表达式。实际上，我们上面得到的h和z都积分的结果是数值解。

求一个matlab数值积分,计算以下等式在三角形区域内A（0,0）,B（1,0）,C（1,1）的积分值(二重积分) matlab求二重积分 MATLAB如何求二重积分 MATLAB问题,给出一个数,求所有组合求EXCEL中,当结果等于2个数值中的任何一个,则最终结果等于2个数值中的另外一个,函数公式怎么写? 如何用matlab做二重积分的数值计算 matlab求数值积分已知一个正定且对称矩阵的特征值,用matlab怎么求此矩阵,给出一系列矩阵也行在matlab里面,对于一个方程组,里面有用符号表示的参数,参数值在前面已经给出来了,怎么得到数值解? 马克思的智商是多少请给出具体的事实和充足的理由，以及智商的具体数值（实在不行就给出一个范围）。 matlab中solve求的解怎么就一个呢 MATLAB中定义一个未知矩阵,在到入一维数值时,默认为行矩阵,怎么设置使得matlab输入数值后默认为列矩阵