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

二重积分求值

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
计算二重积分的时候,结果就是一个式子,不是积分值,这是为什么呢?求教
; w4 n. I* B. G- E/ y! r, n
6 C; k! T. i6 Z! }代码如下:
! c% D% t% i1 l1 K( h& D
/ Z5 j) U1 Z% aclear all0 t, X% k0 a! L: ?
close all  A- {  Q1 g' t9 v% _' b
clc
, `" d' d& l7 m' usyms x y
  F5 i4 N" N5 F9 `- nomega=1;
$ l+ i0 S% D) _1 N) t" l; nmu=0.001;
1 X, J/ i1 r2 Ialpha=0.13;
! x7 X/ S$ M( P- Tbeta=0.0031;( }- F7 O" D; C& o: Y  y1 N" ]" k
u=0.001;v=-0.02; & F; t; R* ]1 M* c( r4 r! j
tau=2.2;tau1=1.5;D=0.11;
5 e' X3 I: M# Q$ ^8 X' G2 o' }tau01=1.42;tau02=1.5;tau03=1.56;tau04=1.9;7 I# b+ G1 P5 C/ P) X; t
omega2=1-u*cos(omega*tau)-v*omega*sin(omega*tau01);) e& Y9 v, U* u7 j7 u" _: b& z  Q
eta=1+(1/mu)*(-(u/omega)*sin(omega*tau)+v*cos(omega*tau01));! Y$ h; P) l, Z0 I/ d
K=D/(omega2.*(1+omega2.*(tau1.^2)));
5 d+ W* r. y( af1=(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)));
# n: @3 ]6 B: N) c# A) W3 `. SF1=int(int(f1,x,0,10),y,-10,10)+ F* P! V1 r% v% t) W& r. D3 s2 u

. c5 _3 D  F4 c5 Y0 g7 x- [
7 T4 X! `2 j' w1 B+ \7 P  ^& }
! V, F3 S, f; P3 z& a6 R结果:& J# n7 ]) O( r
F1=& `& z+ X+ V" B
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)# {  q/ @* ^! Z* c* [- j  |

: `6 I! _% a. A2 V7 p4 y) l6 }3 R0 ], N* s* ]% _8 `

该用户从未签到

2#
发表于 2020-7-1 15:28 | 只看该作者
# d* N, H& \* j7 M4 k1 F7 E  g6 `
积分结果是个关于y的函数,如何转化成具体数值呢

该用户从未签到

3#
发表于 2020-7-2 10:16 | 只看该作者
double(F1)
6 r! U1 W  J3 ]第一句加上这个

点评

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

该用户从未签到

4#
 楼主| 发表于 2020-7-2 10:20 | 只看该作者
lupkpu9 发表于 2020-7-2 10:16
& X2 a. h- ~" k% jdouble(F1). ^& D" P1 w1 e9 @2 e5 q
第一句加上这个

9 }& P' ]* i! e可以了,谢谢哦
# A- Y5 Y- B; m1 c* S

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-5 08:43 , Processed in 0.140625 second(s), 27 queries , Gzip On.

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

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

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