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

二重积分求值

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
计算二重积分的时候,结果就是一个式子,不是积分值,这是为什么呢?求教+ T; p4 H" z0 C; a9 [

1 W) D' ]( g7 k" `6 U代码如下:! V  f/ N4 t& }4 b9 @. U4 Y. i/ v
& c# v# U* K  O8 T2 U- ?0 _
clear all
& q) \. Y+ b- s) u7 ~( \2 \3 Aclose all  _0 i0 J1 i/ B. w! D- J
clc
5 |! T7 a, e! k8 R4 wsyms x y
" ^; G3 n3 r" `; A" [  e0 ]1 Womega=1;
2 Y" [' h! ?3 ?0 P) c$ ~- Emu=0.001;
: I" J  |$ B; {) _0 q# K6 h( Walpha=0.13;
- B6 D; Y1 V$ j5 [  ~: N8 Abeta=0.0031;4 _3 c/ F5 |/ f: G$ |  b- R/ r
u=0.001;v=-0.02;
2 q9 G1 u2 c; L4 D9 P* P- b+ T, \tau=2.2;tau1=1.5;D=0.11;- e" N6 U- f' E* H
tau01=1.42;tau02=1.5;tau03=1.56;tau04=1.9;
$ e- J# n- i! g) d6 o- e' L' Uomega2=1-u*cos(omega*tau)-v*omega*sin(omega*tau01);8 T0 M: u# I/ l7 @
eta=1+(1/mu)*(-(u/omega)*sin(omega*tau)+v*cos(omega*tau01));
8 S* G; n) m; i* PK=D/(omega2.*(1+omega2.*(tau1.^2)));  @8 z  q. t1 B- g1 r
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)));8 _, o; y3 e4 h! z% S. ^
F1=int(int(f1,x,0,10),y,-10,10)$ i5 i$ ?" t* V6 J2 K
% _% N1 I- l5 w

! F6 \0 A/ Z/ ?* U$ m
  K1 O- d$ {$ n$ U2 M/ L结果:
+ `$ b* E* s- ]; p, w% l7 CF1=2 V% p% ~5 S+ B) E3 M5 _
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)! @" ?$ j6 U0 g  E2 R

7 m8 C" }, T. N; W( b6 F
  s8 N8 o/ G* w% _2 F! H: [* L

该用户从未签到

2#
发表于 2020-7-1 15:28 | 只看该作者
0 F  g% A# n- X( ^8 a. z
积分结果是个关于y的函数,如何转化成具体数值呢

该用户从未签到

3#
发表于 2020-7-2 10:16 | 只看该作者
double(F1)
  h" N* k5 W3 A  \* @- S( R第一句加上这个

点评

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

该用户从未签到

4#
 楼主| 发表于 2020-7-2 10:20 | 只看该作者
lupkpu9 发表于 2020-7-2 10:16
; {5 l# u) x  P% }7 fdouble(F1)& h- [% Q' V- l
第一句加上这个

0 I' ~) ]1 q5 e2 G& r可以了,谢谢哦
9 ?8 y2 b9 g0 Z) Y+ l( o+ P( T+ ^) c* e

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-10-8 01:44 , Processed in 0.125000 second(s), 24 queries , Gzip On.

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

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

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