找回密码
 注册
关于网站域名变更的通知
查看: 499|回复: 4
打印 上一主题 下一主题

二重积分求值

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-7-1 15:09 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
计算二重积分的时候,结果就是一个式子,不是积分值,这是为什么呢?求教/ o% m. h) G* f; m0 F9 d: M
# z: y9 W1 P( `4 |7 E0 ?  k" [$ z! q& h
代码如下:
" o0 C! s& T3 ~- y* g& e" K
# x1 D) ?: J2 @clear all
! q, p- z7 Q3 z8 ]1 ^close all
/ t- r1 ~! J; I6 Rclc
2 e3 \7 @  p$ M' |7 Fsyms x y6 ~& s8 W% e0 A' O$ e# [
omega=1;
4 }; d" t* o6 }% Amu=0.001;
' o5 M% _0 G1 e% \9 q1 G/ Z2 v" U2 W) Jalpha=0.13;  H5 l. _7 |3 r: [
beta=0.0031;" m7 v4 X! ^$ d7 p2 ?* T
u=0.001;v=-0.02;
2 Y3 p% s. ^4 Z8 W: K( dtau=2.2;tau1=1.5;D=0.11;5 u) g6 }9 |) |8 W6 G7 G% m! ~
tau01=1.42;tau02=1.5;tau03=1.56;tau04=1.9;, u3 _" a2 i* t, i
omega2=1-u*cos(omega*tau)-v*omega*sin(omega*tau01);0 T: L1 a) m5 W5 Q& @2 `+ h
eta=1+(1/mu)*(-(u/omega)*sin(omega*tau)+v*cos(omega*tau01));
6 e, W9 g. s/ NK=D/(omega2.*(1+omega2.*(tau1.^2)));4 B' j% ~0 W, b4 M/ l) k  z1 X
f1=(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)));
! J( j) h& A  \+ n9 R5 GF1=int(int(f1,x,0,10),y,-10,10)+ {. G4 a* ?5 p7 S$ \" M

/ d5 p1 N9 u  B6 X$ O$ _
& ^& a  K7 ~  g6 I2 w4 s
; c( x. y; f  u# q/ b7 ]结果:
, |1 i) K" |3 a3 B4 ?3 gF1=+ K( Z7 h6 l" x* _' k" n
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)
8 z( o  }1 E8 n- m
8 N/ y6 A) [& s% n3 a
- W( O* [6 Y* S  Z  L& I7 [

该用户从未签到

2#
发表于 2020-7-1 15:28 | 只看该作者

9 C- i* U! n, @$ s积分结果是个关于y的函数,如何转化成具体数值呢

该用户从未签到

3#
发表于 2020-7-2 10:16 | 只看该作者
double(F1)4 Q: D2 P( @' X( N( h0 U1 W
第一句加上这个

点评

可以了,谢谢哦  详情 回复 发表于 2020-7-2 10:20

该用户从未签到

4#
 楼主| 发表于 2020-7-2 10:20 | 只看该作者
lupkpu9 发表于 2020-7-2 10:16
( \) d; e, s, B3 g% mdouble(F1)4 i( v: I1 @' y, `2 o
第一句加上这个
& n! u* r) \# \' A8 h
可以了,谢谢哦/ n5 u" s) }) Z. y! ~

该用户从未签到

5#
发表于 2020-7-3 09:53 | 只看该作者
来学习一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-6-26 16:48 , Processed in 0.078125 second(s), 27 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表