EDA365电子论坛网

标题: 二重积分求值 [打印本页]

---

作者: uaidenp    时间: 2020-7-1 15:09
标题: 二重积分求值
计算二重积分的时候，结果就是一个式子，不是积分值，这是为什么呢？求教

代码如下：
, I' u' v: D4 V0 M
! F% X& p+ y. gclear all
close all
5 i2 ~+ X0 C6 g: Y- F% Uclc
$ ?) }% O5 K3 p! f0 ~* M, E3 xsyms x y
omega=1;
, {- L% ?- _! h$ Rmu=0.001;
- q8 r& u9 R& g! ?( }2 O6 ~alpha=0.13;
beta=0.0031;
3 @8 B- \' g! Ju=0.001;v=-0.02;
tau=2.2;tau1=1.5;D=0.11;
tau01=1.42;tau02=1.5;tau03=1.56;tau04=1.9;
omega2=1-u*cos(omega*tau)-v*omega*sin(omega*tau01);
# C5 Y6 |/ B$ v1 t  L: ueta=1+(1/mu)*(-(u/omega)*sin(omega*tau)+v*cos(omega*tau01));
/ u2 W" A5 L7 L# `5 ~0 L0 R4 FK=D/(omega2.*(1+omega2.*(tau1.^2)));
- b! c7 }8 m% `# p: f  n& pf1=(1./(omega2.*K)).*exp(-(mu./(1536.*K)).*(15.*beta.*(x.^2+y.^2./omega2).^4-32.*alpha.*(x.^2+y.^2./omega2).^3+96.*(x.^2+y.^2./omega2).^2-768.*eta.*(x.^2+y.^2./omega2)));
$ D& D" S" g! k( q6 o7 G4 @! ?1 K$ mF1=int(int(f1,x,0,10),y,-10,10)
9 ]9 `0 H. Y: r0 s7 e2 A0 @% P


- R& s1 ]6 Z+ e/ p* n结果：
! A6 y; {. [* n+ gF1=
int(int((1029484437153*exp((76368780973521777*(x^2 + (2251799813685248*y^2)/2297649917499219)^3)/922337203685477580800 - (17623564840043487*(x^2 + (2251799813685248*y^2)/2297649917499219)^2)/9223372036854775808 - (546330510041348097*(x^2 + (2251799813685248*y^2)/2297649917499219)^4)/590295810358705651712000 - (3488543873076974433998624506509*x^2)/81129638414606681695789005144064 - (1162847957692324811332874835503*y^2)/27593854166072889506733609713664))/34359738368, x, 0, 10), y, -10, 10)

/ _# U: r5 [/ r( C0 c2 G. S
: a8 q$ ]0 {* P7 N

---

作者: llbnmo    时间: 2020-7-1 15:28

积分结果是个关于y的函数，如何转化成具体数值呢

---

作者: lupkpu9    时间: 2020-7-2 10:16
double(F1)
6 X# Y5 g, o2 ^. q/ F第一句加上这个

---

作者: uaidenp    时间: 2020-7-2 10:20

lupkpu9 发表于 2020-7-2 10:16
. x5 f% c" A( f9 edouble(F1)
6 T1 d4 S# y" g; J" T: I4 z6 B第一句加上这个

. @! w. R* D1 r) y) t" _8 j3 k可以了，谢谢哦
: `: Q5 A. m: Y0 m

---

作者: 勇往直前11    时间: 2020-7-3 09:53
来学习一下

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2