Question

如何用matlab的for 循环语句求一个 含参数的 二重积分的 数值积分 急急急??!!!

注意我用过int来做过两次符号积分,但是这个积分用符号积分积不出显式来
Warning: Explicit integral could not be found.
Warning: Explicit integral could not be found.
被积函数如下:
(a^2+x^2-2*a*x*cos(b-y)))*exp(i*2*b)*exp(-a^2)*(a/2)^2 *a
其中a b 为参数,x ,y 为积分变量,
在线等,急急急急,求大神啊!!!!

能够解决的大神,小弟再追加分 !!!

Answer 1

clear all;clc;
syms a b x y
f=(a^2+x^2-2*a*x*cos(b-y))*exp(j*2*b)*exp(-a^2)*(a/2)^2 *a;
jf=int(int(f,y),x)

jf =

1/4*exp(2*i*b)*exp(-a^2)*a^3*(a^2*y*x+1/3*x^3*y-a*x^2*sin(-b+y))

追问

哦,不好意思,被积函数应该是这样,我刚刚仔细检查了下,被积函数应该是这样:
exp(a^2+x^2-2*a*x*cos(b-y))*exp(i*2*y)*exp(-x^2/0.13^2)*(a/2) *a
积分区域为x: [0,0.009] y:[0,2*pi]
非常谢谢了.

追答

改了就没有解析解了

追问

那请问一下,怎么用数值积分来计算上面的积分式呢?

追答

你好，数值积分的话，要给出a和b的值的
还要定义一个函数，算了，用匿名函数吧：
clear all;clc;
a=2;b=1;
syms x y
f=@(x,y)exp(a^2+x.^2-2*a*x.*cos(b-y)).*exp(j*2*y).*exp(-x.^2/0.13^2)*(a/2)*a;
jf=dblquad(f,0,0.009,0,2*pi)

jf =
-1.3838e-004 +3.0236e-004i

追问

这样子,其实我也用过dblaquad积过,我的目的是积出含a ,b 这两个参数的数值表达式,然后画出Z=f(a,b) 的图形,不知道有没有什么方法呢?
其实说了这么多,我就是想请问, 如何进行含参数的 数值积分 ,而且还是个二重的....

追答

知道你的意思，但是你的函数不能得出解析解
符号积分积不出来，只能数值积分